基于截断牛顿法的VTI介质声波多参数全波形反演

不同类别参数间的相互耦合使多参数地震全波形反演的非线性程度显著增加,地震波速度与各向异性参数取值数量级的巨大差异也会使反演问题的性态变差.合理使用Hessian逆算子可以减弱这两类问题对反演的影响,提高多参数反演的精度,而截断牛顿法是一种可以比较准确地估计Hessian逆算子的优化方法.本文采用截断牛顿法在时间域进行了VTI介质的声波双参数同时反演的研究.不同模型的反演试验表明,在VTI介质声波双参数同时反演中,截断牛顿法比有限内存BFGS(Limited-memory Broyden-Fletcher-Goldfarb-Shanno,L-BFGS)法能更准确地估计Hessian逆算子,进而较好地平衡两类不同参数的同时更新,得到了比较精确的反演结果.     

航空重力测量中牛顿中心差分器对确定载体加速度的影响

讨论不同点数的牛顿中心差分器对航空重力测量中确定载体加速度的影响。理论分析表明,当航空重力测量的有效频段位于牛顿中心差分器的通带时,不同点数的牛顿中心差分器所得结果经低通滤波后几乎没有差异。实验表明,当采样间隔为1 s、载体平均速度为155 m/s、所需重力数据的分辨率为10 km时,低通滤波后不同点数的牛顿中心差分器输出结果之间的较差优于2.1 mGal。     

3D simultaneous joint PP‐PS prestack seismic inversion at Schiehallion field,United Kingdom Continental Shelf

A workflow for simultaneous joint PP‐PS prestack inversion of data from the Schiehallion field on the United Kingdom Continental Shelf is presented and discussed. The main challenge, describing reasonable PS to PP data registration before any prestack or joint PP‐PS inversion, was overcome thanks to a two‐stage process addressing the signal envelope, then working directly on the seismic data to estimate appropriate time‐variant time‐shift volumes. We evaluated the benefits of including PS along with PP prestack seismic data in a joint inversion process to improve the estimated elastic property quality and also to enable estimation of density compared with other prestack and post‐stack inversion approaches. While the estimated acoustic impedance exhibited a similar quality independent of the inversion used (PP post‐stack, PP prestack or joint PP‐PS prestack inversion) the shear impedance estimation was noticeably improved by the joint PP‐PS prestack inversion when compared to the PP prestack inversion. Finally, the density estimated from joint PP and PS prestack data demonstrated an overall good quality, even where not well‐controlled. The main outcome of this study was that despite several data‐related limitations, inverting jointly correctly processed PP and PS data sets brought extra value for reservoir delineation as opposed to PP‐only or post‐stack inversion.     

Stability of source mechanisms inverted from P‐wave amplitude microseismic monitoring data acquired at the surface

We study the stability of source mechanisms inverted from data acquired at surface and near‐surface monitoring arrays. The study is focused on P‐wave data acquired on vertical components, as this is the most common type of acquisition. We apply ray modelling on three models: a fully homogeneous isotropic model, a laterally homogeneous isotropic model and a laterally homogeneous anisotropic model to simulate three commonly used models in inversion. We use geometries of real arrays, one consisting in surface receivers and one consisting in ‘buried’ geophones at the near‐surface. Stability was tested for two of the frequently observed source mechanisms: strike‐slip and dip‐slip and was evaluated by comparing the parameters of correct and inverted mechanisms. We assume these double‐couple source mechanisms and use quantitatively the inversion allowing non‐double‐couple components to measure stability of the inversion. To test the robustness we inverted synthetic amplitudes computed for a laterally homogeneous isotropic model and contaminated with noise using a fully homogeneous model in the inversion. Analogously amplitudes computed in a laterally homogeneous anisotropic model were inverted in all three models. We show that a star‐like surface acquisition array provides very stable inversion up to a very high level of noise in data. Furthermore, we reveal that strike‐slip inversion is more stable than dip‐slip inversion for the receiver geometries considered here. We show that noise and an incorrect velocity model may result in narrow bands of source mechanisms in Hudson's plots.     

江西省冻雨垂直温度层结分析及预报

利用1960-2013年江西省89个气象观测站常规观测资料,采用统计分析和插值处理方法,对江西106次冻雨过程垂直温度层结进行分析,并建立冻雨预报思路。结果表明：江西省冻雨发生通常伴随1000-850 hPa剧烈降温和700 hPa增温,925-700 hPa存在明显逆温。融化层是冻雨发生的必要条件,87.7 %样本融化层位于700 hPa,9.4 %样本融化层位于850 hPa;一般以1000 hPa温度<1.0 ℃、925 hPa和850 hPa温度<-2.0 ℃和700 hPa温度≥0.0 ℃作为江西冻雨预报的温度阈值标准,同时也应注意融化层位于850 hPa的情况。地面气温越低,越有利于冻雨形成,冻雨发生时最低气温为0.0 ℃以下,平均气温为1.0 ℃以下。冷暖空气持续在27°-28 °N交汇,且鄱阳湖以南至抚河流域平均气温＜1.0 ℃和最低气温＜0.0 ℃出现比其他地区多,造成该区域冻雨发生频次比其他地区明显偏多。另外,通过大气逆温层结、地面气温与冻雨的对应关系,建立基于模式资料和常规观测资料的自动诊断方法。     

2010年2月辽宁地区冻雨天气诊断分析

应用常规天气资料,从天气、气候及大气物理量场等方面,对2010年2月下旬发生在辽宁地区历史罕见的冰冻雨雪灾害性天气发生发展的成因、主要特征以及造成严重灾害的主要原因进行诊断分析。结果表明：500 hPa低涡分裂的小槽东移是冻雨天气的环流背景;700 hPa源源不断的水汽输送和低层辐合作用是降水产生的基本条件。地面强冷空气是冻雨天气的触发机制;高层为冷的冰晶层,对流层中低层有温度t>0℃的暖的融化层,近地面1-2 km为气温t<0℃的冷层,同时地面达到 0℃以下,这种复杂的逆温层是冻雨产生的天气条件。高层暖平流是维持中空暖层结构的重要条件之一,对暖层的建立和破坏起较大作用,暖层消失冻雨天气也随之结束。中低空水汽饱和度、地面温度和中低空上升、下沉运动与冻雨的强度有密切关系。     

基于探地雷达的混凝土无损检测反演成像方法

本文针对水利水电工程中的混凝土构件无损探测问题,提出了一种新的反演成像方法.传统方法主要依赖人工经验观察偏移探测数据,精度无法保证.这种方法以Maxwell方程的TM问题为数学模型,在混凝土的表面放置发射源与接收器,根据探地雷达的接收数据,采用同伦优化方法结合收敛速度较快的阻尼高斯牛顿方法作为反演算法,其成像结果更具定量化和可视性.该成像方法克服了传统雷达剖面图只能近似地反映混凝土中埋藏物及缺陷的深度、大小等缺点,不仅能够探测混凝土内部异物的位置,而且能够比较准确地确定异物的大小及属性.通过计算机模拟以及对实际资料的处理,验证了该反演成像结果直观、可靠、具有一定的应用价值.     

解释位场全张量数据的张量局部波数法及其与常规局部波数法的比较

全张量探测技术以其信息量大、精度高、干扰小等优点在地球物理领域中得到广泛应用.本文提出采用张量局部波数法来进行位场全张量数据的解释,首先给出了张量局部波数的定义,然后推导出利用张量局部波数法进行反演的基本公式.本文方法在进行张量数据反演时无需事先知道场源体的类型(构造指数)即可获得场源体的位置信息,且可根据位置参数对场源体的类型进行估计.通过理论模型证明张量局部波数法可以很好地完成位场全张量数据的反演工作,并将其与常规局部波数法进行对比,证明全张量局部波数法的反演结果更加准确,即使在测点分布不合理的情况下,张量局部波数法仍可以获得准确的结果.最后应用张量局部波数法对美国得克萨斯州实测重力数据进行了反演,其反演结果与已有的研究成果相一致.     

塔里木盆地航磁场分析与磁源体结构

本文研究塔里木盆地区域磁异常图的反演及磁源体结构.由于众多异常的叠加和反演固有的多解性,区域磁异常图的准确解释是非常困难的.三维欧拉反褶积是一种确定地质体位置和埋藏深度的自动定量反演方法,比较适用于计算区域磁异常源的埋藏深度.由于大型克拉通沉积盆地地层具有上新下老的规律性,将磁异常源分解为三个深度层次,圈定它们各自的分布区域,便可将它们与形成的地质作用及时代联系起来,为准确解释区域磁异常图提供可靠的依据.本文应用三维欧拉反褶积反演方法,计算出的塔里木盆地深度为2~5 km、5~10 km、10~20 km三个等级的磁异常源,它们与形成的地质作用及时代分别为: 中生代构造运动,海西期玄武岩侵位和太古代结晶基底的变质作用;圈定了它们各自的分布区域.     

海岸效应对近海地区大地电磁测深数据畸变作用研究

在近海地区采集的大地电磁测深数据通常受到海岸效应的影响,使得大地电磁测深数据发生畸变,因而很难利用大地电磁测深资料较为可靠地获得地下深部的电性结构.本文通过正演模拟方法,分析和总结海水深度变化和海底地形变化对近海地区大地电磁测深数据的畸变影响.当测区与海岸线的距离小于目标频率的大地电磁场趋肤深度时,高导海洋的存在会严重影响测区内电磁场的分布.由于海岸效应的影响,大地电磁测深视电阻率曲线和相位曲线均会发生不同程度的畸变,在低频部分,这种畸变作用尤为明显.大地电磁测深一维Occam反演方法和二维非线性共轭梯度反演方法,对近海地区浅部地层具有较好的反演效果.随着海水深度的增加和海底地形的复杂变化,两种反演方法均会出现不同程度的假异常,为地质解释工作造成了影响.近渤海地区的实测大地电磁测深数据在低频部分可能受到海岸效应的影响而导致视电阻率曲线的严重畸变.