论文标题
蒙特卡洛集合统计的敏感性和研究深度
Sensitivity and depth of investigation from Monte Carlo ensemble statistics
论文作者
论文摘要
对于许多地球物理测量值,例如直流电流或电磁诱导方法,信息随深度逐渐消失。解释从此类测量值估计的模型时,必须考虑这一点。因此,测量灵敏度分析和确定研究深度是地球物理数据处理过程中的标准步骤。在基于确定性梯度的反转中,最常用的灵敏度度量是差异灵敏度,因为这些反转需要计算Jacobian矩阵。相比之下,蒙特卡洛反转方法中可能不容易获得差异敏感性,因为这些方法不一定包括对正时问题的线性化。取而代之的是,先前的合奏用于模拟正向响应的合奏。然后,根据贝叶斯推论更新了先前的合奏。我们建议使用先前的合奏和正向响应集合之间的协方差来构建灵敏度度量。在Monte Carlo方法中,此协方差的估计不需要对正向模型的其他计算。通过先前合奏的方差使这种协方差归一化,人们获得了简化的回归系数。我们使用简单的前向模型研究了这种简化的回归系数和差异灵敏度之间的差异。对于线性向前模型,简化的回归系数等于差异灵敏度,除了采样误差的影响和先前分布的相关结构的影响。在非线性情况下,对简单的非线性正向模型和频域电磁向前模型分析了简化回归系数作为灵敏度度量的行为。 [...]
For many geophysical measurements, such as direct current or electromagnetic induction methods, information fades away with depth. This has to be taken into account when interpreting models estimated from such measurements. For that reason, a measurement sensitivity analysis and determining the depth of investigation are standard steps during geophysical data processing. In deterministic gradient-based inversion, the most used sensitivity measure, the differential sensitivity, is readily available since these inversions require the computation of Jacobian matrices. In contrast, differential sensitivity may not be readily available in Monte Carlo inversion methods, since these methods do not necessarily include a linearization of the forward problem. Instead, a prior ensemble is used to simulate an ensemble of forward responses. Then, the prior ensemble is updated according to Bayesian inference. We propose to use the covariance between the prior ensemble and the forward response ensemble for constructing sensitivity measures. In Monte Carlo approaches, the estimation of this covariance does not require additional computations of the forward model. Normalizing this covariance by the variance of the prior ensemble, one obtains a simplied regression coefficient. We investigate differences between this simplified regression coefficient and differential sensitivity using simple forward models. For linear forward models, the simplied regression coefficient is equal to differential sensitivity, except for the influences of the sampling error and of the correlation structure of the prior distribution. In the non-linear case, the behaviour of the simplified regression coefficient as sensitivity measure is analysed for a simple non-linear forward model and a frequency-domain electromagnetic forward model. [...]