基于垂向一阶导数与解析信号比值的欧拉反演方法

常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。     

分离三分量远震接收函数的多道最大或然性反褶积方法

迄今为止 ,广泛应用的分离接收函数方法主要基于单个台站的观测数据 ,并依赖等效震源假定的有效性。当地壳内存在明显的间断面时 ,等效震源假定不适于宽频带地震波形数据 ,接收函数的估计将在不同程度上遭到破坏。解决这个问题的出路在于避免使用等效震源假定。依据多道最大或然性反褶积原理 ,给出了利用多道观测数据分离接收函数的方法。我们的方法避免了等效震源假定 ,这对于在复杂构造环境条件下改善接收函数的估计具有实际价值。合成地震图的数值检验结果验证了给出的方法。利用该方法 ,我们得到了跨越大别山造山带 ,从江西大冶附近的大箕铺 (30°2 0′N ,1 1 5°0 3′E)到河南兰考附近的崔林 (34°4 0′N ,1 1 4°4 9′E) ,总长度约 5 0 0km的流动地震台阵剖面各台站的三分量接收函数估计。与基于单台数据的接收函数径向分量比较表明 ,两种方法给出的结果可能存在不可忽略的差别     

瞬变电磁法三维拟地震成像信息提取技术

在瞬变电磁法中,为提高对地下电性界面的分辨率,可移植地震勘探中的诸多技术对界面直接成像,以克服反演中的不适定性.首先为使得波形在高低不同的电性基底所产生的虚拟反射波振幅有可比性,在将原有瞬变电磁场逆变换成虚拟波场数据的基础之上,需要求取均衡系数;接着通过对虚拟波场的离散数据求取反褶积,可基本抵消波场变换的展宽作用和大地滤波影响,以达到对电性界面精确定位的效果;然后,对波形做低通滤波处理后运用克希霍夫积分公式,将虚拟地震波场由地面向地下反向逆时偏移成像.把波场分析原理用于对瞬变电磁场的解释,就形成了瞬变电磁偏移方法.从波场角度看,它拓展和丰富了瞬变电磁场的内涵.使用本方法,可更加容易得到从前用常规解释方法处理无法提取的信息.在延拓成像中可更精确的确定界面深度,并为在瞬变电磁法中合理引入拟地震处理解释方法提供重要的理论依据和实验参照.     

地表一致性反褶积在地震勘探中的应用及效果

地震资料的反褶积处理是通过改造地震激发子波,进而消除地震激发子波在传播过程中所受的虚反射、层间多次反射和大地滤波等影响的一种地震勘探资料处理方法.反褶积的方法很多,如:脉冲反褶积、预测反褶积、地表一致性反褶积.它们之间主要区别之一在于对地震子波的假设和估计地震子波的方法.所以在处理过程中应根据不同的区域资料特征采取不同的反褶积方法.本文以河南省某煤预查区地震勘探为实例,着重总结和比较地表一致性反褶积技术在地震资料处理中的应用效果.应用研究表明,适当选择时窗和自相关步长进行自相关分析,地表一致性反褶积能够展宽频谱,压缩地震子波,并能校正地震信号的相位谱,输出零相位子波,可以较大程度地提高地震资料的分辨率,提高勘探能力.     

基于欧拉反褶积方法计算川滇交界重力变化场源特征

本文以川滇地区2012-2014年来的多期流动重力观测网络数据为基础,采用欧拉反褶积方法对引起重力变化的场源深度和空间分布规律进行了反演和解释.本文通过构建理论模型。得到了最优的反演参数,并对实际数据进行了计算和分析.结果表明:选择构造指数(SI)为1时适合对流动重力数据进行反演.实际资料的反演结果显示在2012-2014年间的重力变化场源位置更集中于垂直于昭通断裂的北西向分布,场源深度分布在30±10 km范围,形态与鲁甸地震发震构造位置相似.这可能与2014年8月3日鲁甸M_S 6.5地震前的中上地壳的深部应力场变化相关.依据本文研究结果和方法,可以开展由场至源的定量研究,为流动重力变化信号的反演和解释提供一条新思路.     

由精密控制人工震源观测到的绵竹5.6级地震前后波速变化

地震危险区断层的应力状态变化及其精确测量一直是地球物理学关注的热点问题之一,主动震源探测技术的发展为解决上述问题带来了可能.汶川大地震后,结合汶川地震断裂带科学钻探,在四川省绵竹市九龙镇--跨越龙门山断裂带前山断裂,利用10 t精密控制人工震源进行了连续监测试验,在试验观测期间的2009年6月30日凌晨2时3分,当地发生了M_s5.6级地震,地震前后观测数据的对比分析表明,该地震发生前后,穿过断裂带的直达S波走时发生了时延为5~9 ms的微弱变化,射线路径上的平均波速相对下降了约0.3%,该时延主要是由地震同震效应引起的断裂带内应力调整和物性变化造成的.所使用的精密控制人工震源性能稳定,适于开展断层监测.这一结果对于发展断裂带主动震源动态监测技术具有重要意义.     

中国东北地区重力场研究——利用改进的三方向小子域滤波划分主构造线及大地构造单元

利用重力异常进行大地构造单元的划分是区域重力资料处理与解释的重要任务.为了精细地划分中国东北地区大地构造单元及其分界线位置,对三方向小子域滤波进行了改进,探讨了改进的三方向小子域滤波的合理性.利用小子域及改进的三方向小子域滤波对研究区布格重力异常进行处理,并与欧拉反褶积反演结果进行对比分析.结果表明,改进的三方向小子域滤波在识别平面断裂的应用中具有更好的效果.综合分析东北地区布格重力异常特征、小子域和改进的三方向小子域滤波结果,划分出延续性较好、规模较大的断裂18条.对研究区重力异常进行了分区,描述了一、二级异常的分区特点,给出了东北地区大地构造单元的基本格架.黑河-齐齐哈尔-白城等重要构造线平面位置的确定突破了传统的重力异常分区模式,指出东北地区NE向构造具有走滑特征,且走滑距离具有由西向东逐渐增大的特点.     

识别断裂构造的DCT法

为了提高欧拉反褶积(Euler法)的反演精度,提出基于DCT的Euler法确定断层断点位置的反演方法,从理论上给出了基于DCT的Euler齐次方程,指出异常导数的计算精度是影响Euler法反演精度的核心问题.台阶模型实验证实:用DCT法计算的重力异常垂向一阶导数的最大误差为0.460×10-9/s2,均方差为0.189×10-9/s2,水平导数最大误差为0.182×10-9/s2,均方差为0.028×10-9/s2,与DFT法相比,其计算的异常导数具有很高的精度.通过对正演模型的计算,确定了基于DCT的Euler法模型反演的位置为台阶断面的中心点位,与理论点位相对应;在反演模型实验中,证实了DCT法具有较高的反演精度,同时在实际应用中,进一步验证了该法的可行性、实用性和准确性.     

最小熵反褶积的微机程序

讨论了最小熵反褶积的基本原理,并给出BASIC微机程序,利用该程序处理了合成地震记录。      

海上地震资料子波提取鬼波压制技术

海洋地震勘探由于受到海面强反射界面的影响,记录数据中存在鬼波,致使数据频谱中产生陷波效应,限制数据的频带宽度,影响数据的分辨率,干扰地层反射的识别.近年来发展的鬼波压制方法,能够较好的压制地震数据中的鬼波,然而现有鬼波压制技术容易产生延续相位,且计算效率较低.为进一步提高鬼波压制效果和工作效率,本文对海上地震资料特性进行分析,研究认为,海上地震资料中的海底反射信号真实记录了海底有效反射、震源鬼波和检波点鬼波的形态,可以利用提取的海底反射子波对数据进行确定性子波反褶积处理,能够有效的压制震源鬼波和检波点鬼波.模型海上实际数据测试结果表明,该方法具有鬼波压制效果好、计算效率高的优势,处理之后的数据频带得到有效拓宽,低频和高频信息均得到加强.