控制点与检查点数量和比例对水深遥感反演精度的影响分析

在水深遥感反演半经验模型中,水深控制点和检查点选取是一项非常重要的工作,前者用于建立实测水深值与遥感影像灰度值之间的数量关系,后者用于评价水深反演精度。但前人在相关研究中并未就控制点和检查点的选取数量和比例给出有实验依据的建议,亦未开展其对水深反演精度的影响分析。针对这个问题,作者使用6020个实测水深点和World View-2、GF-1 WFV、Landsat8 OLI 3种多光谱影像,基于三波段水深反演模型开展了81组实验,分析比较了不同数量和比例的控制点与检查点对反演结果的影响,评价过程利用平均绝对误差(Mean Absolute Error,MAE)、平均相对误差(Mean Relative Error,MRE)和反映实测水深与反演水深相关性的决定系数R2等3种参数。结果表明:(1)当控制点数量达到31个时,水深反演精度即趋于稳定;(2)检查点数量在30个时,其评价指标已可以代表模型反演精度;(3)控制点和检查点的数量比例对反演精度并无影响,但控制点选取过多则会产生反演精度降低的现象,采用人工选取检查点时剔除浅水区的异常点即可有效避免这一问题。     

基于叠前深度偏移的AVO反演在神狐海域试采区的应用

为了更好地研究神狐海域试采区天然气水合物的展布形态和赋存状况,基于高精度层析速度反演后的叠前深度偏移地震资料进行AVO反演技术研究。为在天然气水合物富集区有效地进行AVO研究,AVO处理前运用叠前深度偏移技术,通过高精度层析反演不断迭代速度模型;然后通过偏移来改善资料的质量,结合实测钻井资料进行AVO正演分析、反演可行性分析、属性交会分析;最后将具有高横向分辨能力的地震资料和高垂向分辨能力的测井资料与地质相结合进行叠前同时反演。该方法可反演出纵、横波阻抗等属性,直接进行岩性和流体识别,结果表明,在AVO反演中引入叠前深度偏移技术,为反演提供了较高质量的叠前共反射点道集和精细的层速度模型,改善了AVO技术在试采区的应用效果,提高了反演的可靠性。反演出的各种属性为AVO分析提供了丰富的信息,能够比较准确地反映地下水合物矿体的展布情况,可用于后续水合物储层的分析及综合研究,对未来的水合物勘探开采具有重要的现实意义。     

Combining discrete cosine transform and convolutional neural networks to speed up the Hamiltonian Monte Carlo inversion of pre-stack seismic data

Markov chain Monte Carlo algorithms are commonly employed for accurate uncertainty appraisals in non-linear inverse problems. The downside of these algorithms is the considerable number of samples needed to achieve reliable posterior estimations, especially in high-dimensional model spaces. To overcome this issue, the Hamiltonian Monte Carlo algorithm has recently been introduced to solve geophysical inversions. Different from classical Markov chain Monte Carlo algorithms, this approach exploits the derivative information of the target posterior probability density to guide the sampling of the model space. However, its main downside is the computational cost for the derivative computation (i.e. the computation of the Jacobian matrix around each sampled model). Possible strategies to mitigate this issue are the reduction of the dimensionality of the model space and/or the use of efficient methods to compute the gradient of the target density. Here we focus the attention to the estimation of elastic properties (P-, S-wave velocities and density) from pre-stack data through a non-linear amplitude versus angle inversion in which the Hamiltonian Monte Carlo algorithm is used to sample the posterior probability. To decrease the computational cost of the inversion procedure, we employ the discrete cosine transform to reparametrize the model space, and we train a convolutional neural network to predict the Jacobian matrix around each sampled model. The training data set for the network is also parametrized in the discrete cosine transform space, thus allowing for a reduction of the number of parameters to be optimized during the learning phase. Once trained the network can be used to compute the Jacobian matrix associated with each sampled model in real time. The outcomes of the proposed approach are compared and validated with the predictions of Hamiltonian Monte Carlo inversions in which a quite computationally expensive, but accurate finite-difference scheme is used to compute the Jacobian matrix and with those obtained by replacing the Jacobian with a matrix operator derived from a linear approximation of the Zoeppritz equations. Synthetic and field inversion experiments demonstrate that the proposed approach dramatically reduces the cost of the Hamiltonian Monte Carlo inversion while preserving an accurate and efficient sampling of the posterior probability.     

单轴压缩与劈裂荷载下灰岩声发射b值特性研究

通过开展灰岩巴西劈裂和单轴压缩声发射试验,探讨两种加载方式下岩石破裂声发射b值特性及b值计算影响因素。结果表明,累积计算b值拟合度高,误差较小;步距为5 dB时不同门槛值下b值随时间变化的总体规律基本一致,选择5 dB的步距和40 dB的门槛值计算b值比较合理,能够降低非岩石破裂信号对b值变化的影响。由动态b值的波动程度可以看出,在劈裂荷载下灰岩损伤演化过程可分为3个阶段,即(0～40%)?c(?c为峰值强度)、(40%～90%)?c和(90%～100%)?c,而在单轴加载下可分为：(0～80%)?c和(80%～100%)?c两个阶段,且不同阶段岩石内部损伤的程度不同。在裂纹扩展阶段,相比单轴加载,劈裂荷载下动态b值增加较为稳定,但在整个破坏过程中劈裂荷载下动态b值波动更大,在单轴压缩和劈裂加载方式下b值的大小取决于灰岩破裂面上的结构性质及破裂模式。     

基于微观尺度的页岩损伤破裂数值试验

为正确认识下寒武统页岩微观尺度破裂特征及声发射信号演化规律,基于数字图像处理技术对含石英充填的页岩进行非均匀性表征,采用岩石真实破裂分析系统RFPA^2D-DIP建立数值模型,模拟页岩在不同围压作用下的破坏过程。试验结果表明：在微观尺度结构下,页岩破裂模式大致可归纳为3种形式,分别为倒V形（0、2 MPa）、V形（6 MPa）和倒Z形（10 MPa）。由于石英矿物颗粒分布不规律且具有非均匀性,导致应力分布对页岩微观结构产生显著差异。单轴时,累计声发射（简称AE）表现为“平缓-线性-平缓-稳定”模式;围压为2 MPa时,累计AE表现为“平缓-非线性-激增-平缓-稳定”模式;围压为6、10 MPa时,累计AE表现为“平缓-非线性-稳定”模式。页岩的破裂模式和微观结构破坏规律为页岩气的压裂开采提供重要的理论指导。     

基于同轴差距测量法的南海深水海底沉积物声衰减特性研究

通过分析沉积物声波测量过程能量损失的实质,阐述平行轴差距衰减测量法和垂直轴差距衰减测量法的研究基础——差动式衰减测量方法原理,并推导了其衰减系数公式。根据南海沉积物的柱状样品分装的特点,结合以上两种方法,提出了同轴差距衰减测量法,这种方法具有原理上的合理性并且其可操作性强;运用此方法测量了沉积物在常温和温度控制变化条件下的声波信号,计算了衰减系数,研究了温度对沉积物声波传播能量的影响,得出南海深水海底沉积物具有以下声衰减特性：沉积物含砂量高,声衰减系数大;随着温度的升高,沉积物的声衰减系数变化具有不均匀性,整体呈非线性减小趋势。以上研究将为声学遥测和反演海底沉积物的物理力学特性提供数据和方法支持。     

重磁异常正反演可视化实时方法技术改进

摘 要 给出了2∙5D 多边形截面棱柱体重磁异常的简化表达式及其偏导数;提供了实用的可 视化模型建立、修改、实时正反演中的一些具体方法技术。借助这一操作方便、功能强大的 图形化实时反演系统‚资料解释工作者可以综合利用各种地质资料对计算过程实现引导和控 制‚使其解释思想及技术方案得到充分体现‚从而提高计算速度和质量。     

基于BP神经网络反演分析的隧道塌方机理研究

岩溶区隧道在施工过程中易发生崩塌,针对力学方面的隧道塌方机理分析较多,而针对岩溶软弱破碎带等地层方面的隧道塌方机理研究分析较少。为保证隧道施工的安全性、经济性和可行性,掌握隧道施工中的塌方机理非常有必要。依托贵州某岩溶破碎地层隧道在开挖过程中发生的坍塌现象,结合隧道的监测数据,运用BP神经网络的构建原理,对隧道的地层参数进行反演。将反演土体力学参数输入到FLAC3D有限元软件构建的不同施工方法模型中,对典型断面的崩塌破坏机制和风险进行判断和分析。结果表明：施工方法对隧道开挖的稳定性影响较大,针对围岩等级为Ⅴ级的隧道,采用三台阶七步法和单侧壁导坑法施工较安全,隧道塌方与隧道双向同时开挖没有关系;反演所得的隧道拱顶位移预测值为2.3 cm,地表位移预测值为1.2 cm,与监测数据偏差13%左右,反演结果具有一定的可信度。研究结果对岩溶区软弱破碎地层断面隧道公路建设具有重要指导意义。     

Identifying causes of ground-penetrating radar reflections using time-domain reflectometry and sedimentological analyses

Ground-penetrating radar (GPR) is a geophysical technique widely used to study the shallow subsurface and identify various sediment features that reflect electromagnetic waves. However, little is known about the exact cause of GPR reflections because few studies have coupled wave theory to petrophysical data. In this study, a 100- and 200-MHz GPR survey was conducted on aeolian deposits in a quarry. Time-domain reflectometry (TDR) was used to obtain detailed information on the product of relative permittivity (ɛ_r) and relative magnetic permeability (μ_r), which mainly controls the GPR contrast parameter in the subsurface. Combining TDR data and lacquer peels from the quarry wall allowed the identification of various relationships between sediment characteristics and ɛ_rμ_r. Synthetic radar traces, constructed using the TDR logs and sedimentological data from the lacquer peels, were compared with the actual GPR sections. Numerous peaks in ɛ_rμ_r, which are superimposed on a baseline value of 4 for dry sand, are caused by potential GPR reflectors. These increases in ɛ_rμ_r coincide with the presence of either organic material, having a higher water content and relative permittivity than the surrounding sediment, or iron oxide bands, enhancing relative magnetic permeability and causing water to stagnate on top of them. Sedimentary structures, as reflected in textural change, only result in possible GPR reflections when the volumetric water content exceeds 0·055. The synthetic radar traces provide an improved insight into the behaviour of radar waves and show that GPR results may be ambiguous because of multiples and interference.     

瑞雷波多阶模式频散曲线稀疏正则化反演方法研究

目前瑞雷波多阶模式频散曲线反演中仅考虑数据的拟合,缺乏对模型的约束,不能很好地刻画地层间断面的问题,针对此问题,研究了瑞雷波多阶模式频散曲线稀疏正则化反演方法.正演模拟基于广义反射-透射系数法,数值计算上采用一种快速求根方法,与二等分方法相比,能够在很短的时间内达到最优的收敛效果;反演建模时采用L1范数正则化方法对模型进行稀疏性刻画,使反演结果更加符合地质实际;在反问题的数值实现上,针对稀疏正则化模型提出一种隐式迭代正则化算法,其迭代算子具有非膨胀特性,可以收敛到极小化问题的解.数值实验结果表明,新的反演方案具有计算效率高,模型"逐块"光滑的特性刻画好,对非高斯噪声鲁棒性强的特点.