叠前三参数非高斯反演方法研究

针对地球物理反演中广泛采用的"噪声高斯分布假设",本文研究了叠前地震资料中噪声的非高斯分布特征,提出了针对非高斯噪声的地震叠前非高斯反演概念和思想,构造了能同时压制高斯和非高斯噪声的混合范数作为反演目标函数,采用改进的Powell算法进行求解,有效地抑制了叠前地震资料中的高斯和非高斯混合噪声.模型试算和实际地震数据的反演结果验证了方法的正确性和算法的可靠性.     

九江-瑞昌矿集区的3D结构及对区域找矿的启示

本文分析了九瑞地区重、磁场的分布特征.运用最新的方法技术对已有的重磁资料重新进行处理:采用重磁多尺度边缘检测方法,对九瑞矿集区区域重力和航磁数据进行了边缘检测,并根据检测结果重新厘定了断裂系统的展布位置.在整理、分析九瑞地区地质、地层物性资料的基础上,对实测的1∶5万重磁数据进行较细致的准三维反演.同时将重磁三维物性反演应用到岩浆岩空间结构研究中,获得了矿集区地层结构及岩浆岩三维空间形态特征.根据反演所得磁化率强弱,分析了岩体的基性程度,为寻找与火山岩、侵入岩体有关的金属矿产提供了指示信息.最终建立的模型给出了地下地层的分布特征、控矿构造的展布规律、与成矿相关岩体的三维形态以及已知矿点的空间分布特征,为在九瑞矿集区的深部寻找隐伏矿体提供了新的信息.     

2008年新疆乌恰6.8级地震序列震源特征及帕米尔东北缘应力场研究

本文使用新疆区域数字地震台站记录的宽频带长周期数字波形资料,在时间域反演了2008年10月5日新疆乌恰6.8级地震的强余震及其周围先后发生的52次中等强度地震的矩张量解,结合Harvard大学在该区域的地震矩张量结果,研究了帕米尔东北缘的应力场分区特征.研究结果显示,位于印度板块向欧亚板块推挤的前缘及向北凸出的弧型构造的最北缘的卡兹克阿尔特弧形活动褶皱-逆断裂带,以逆冲推覆活动为主,并有部分走滑类型的地震,基本不存在正断层类型的地震;该弧型构造近东西走向的顶部(文中的西区)与其北西走向的东侧(文中的东区)的局部应力场最大主压应力方向不同,分别为NW、NNE方向,显示出在承受印度板块向欧亚板块俯冲作用的同时,东区也更多的受到了塔里木块体顺时针旋转作用的影响.位于帕米尔陆内俯冲和变形作用强烈、碰撞造成深源地震带东段的南区,地震以走滑错动为主, 逆断、正断层都有,显示出相对复杂的应力状态.位于帕米尔高原内部的西区和南区的应力场最大主压应力方向一致,由北向南,由最大主压应力轴接近水平,过渡为最大主张应力轴接近水平,一定程度揭示了板块俯冲的状态.结合南区和西区的地震深度差异及机制解中断层面的倾角,推测在中帕米尔的东部,由北向南的板块俯冲至150~170km深度,俯冲角度为60°左右.     

测井、井间地震与地面地震数据联合约束下的地质统计学随机建模方法研究

利用能够整合测井信息与井间地震信息的地质统计学随机模拟方法,结合传统的地质统计学反演思路,得到了一种能够同时整合测井、井间地震与地面地震三种先验信息的地质统计学反演与储层建模方法.由于井间射线信息、测井信息与地面地震数据在随机反演与建模过程当中都得到了尊重,因此与传统地质统计学反演仅利用了测井与地面地震数据相比,本文的地质统计学反演与建模方法更充分地利用了先验信息,有效提高了反演的精度,降低了随机建模中的多解性.基于理论数据的测试证实了上述观点.     

Frequency‐domain acoustic‐elastic coupled waveform inversion using the Gauss‐Newton conjugate gradient method

We developed a frequency‐domain acoustic‐elastic coupled waveform inversion based on the Gauss‐Newton conjugate gradient method. Despite the use of a high‐performance computer system and a state‐of‐the‐art parallel computation algorithm, it remained computationally prohibitive to calculate the approximate Hessian explicitly for a large‐scale inverse problem. Therefore, we adopted the conjugate gradient least‐squares algorithm, which is frequently used for geophysical inverse problems, to implement the Gauss‐Newton method so that the approximate Hessian is calculated implicitly. Thus, there was no need to store the Hessian matrix. By simultaneously back‐propagating multi‐components consisting of the pressure and displacements, we could efficiently extract information on the subsurface structures. To verify our algorithm, we applied it to synthetic data sets generated from the Marmousi‐2 model and the modified SEG/EAGE salt model. We also extended our algorithm to the ocean‐bottom cable environment and verified it using ocean‐bottom cable data generated from the Marmousi‐2 model. With the assumption of a hard seafloor, we recovered both the P‐wave velocity of complicated subsurface structures as well as the S‐wave velocity. Although the inversion of the S‐wave velocity is not feasible for the high Poisson's ratios used to simulate a soft seafloor, several strategies exist to treat this problem. Our example using multi‐component data showed some promise in mitigating the soft seafloor effect. However, this issue still remains open.     

综合定点前兆异常推测地震危险区域的一种途径

1998年张北6.2级地震前,山西、河北及其附近地区出现的定点形变、流体前兆异常,具有随时间推移从震区外围大范围区域向震区附近迁移、集中的趋势特性.通过构建异常信号传播函数,借助网格搜索—拟牛顿最小二乘法反演汇聚区域即推测可能的地震危险区,所得汇聚区域边缘距实际震中仅25 ～51 km、汇聚时刻与发震时刻仅差数天,且震前中短期阶段汇聚速度加快.     

Yabuki&Matsu’ura方法在汶川M_W7.9地震反演中的应用

2008年5月12日,青藏高原东缘龙门山断裂带发生汶川MW7.9地震,该地震使得北川—映秀断裂、灌县—江油断裂发生了同震破裂。本文主要利用震后通过复测获得的GPS同震形变场,采用Yabuki&Matsu’ura反演计算方法和分段平面断层模型,反演了地震同震滑动分布。结果表明:映秀—北川主破裂带的断层错动,在映秀附近以逆冲滑动为主,而在北川以北,其走滑运动明显大于逆冲,这一结果与震后地质调查结果与通过地震波研究获得的断层破裂特征相一致;反演得到的最大滑动量达到9.3m和9.6m,分别对应于这次地震中地表破坏最为严重的北川和映秀地区;由所获得的滑动分布计算的地震矩为8.07×1020 N.m,对应的震级为MW7.9。研究结果初步显示,Yabuki&Matsu’ura反演方法可适用在内陆地震断层反演计算中。     

扩频微波网桥在新疆测震网络中的应用

以频率为2.4 GHz和5.8 GHz设备为例,对扩频微波网桥在新疆测震网络中的应用范围进行阐述;根据实际工作中无线网桥的安装架设经验,对在不同条件、不同环境下可能出现的问题提出相应解决方案.     

地震转换波对叠前同时反演的影响研究

叠前同时反演是油气探测的一种有效工具.其理论基础是平面P波Zoeppritz方程计算的反射系数的近似,是入射角的函数.叠前同时反演可以利用三项或两项Fatti方程进行反演分析.本文针对实际油田的测井数据,利用反射率法模拟了仅包含P波一次反射记录,包含P波一次反射和P波层间多次波记录以及全波场地震记录,再利用叠前同时反演对合成地震记录进行反演研究.研究结果表明,在大偏移距处P波主要反射受到其它模式波的污染,从而影响了叠前同时反演结果的精度.对于薄互层介质当转换波影响严重时,使用小角度数据的两项AVO反演比使用大角度数据的三项AVO反演更合理可靠.     

基于全区视电阻率的线源FCSEM二维正演模拟

频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.