基于PSO的地震盲反褶积方法

本文基于地层反射系数非高斯的统计特性,在反褶积输出单位方差约束下,将反褶积输出的负熵表示为非多项式函数,作为盲反褶积的目标函数,然后采用粒子群算法优化目标函数寻找最佳反褶积算子,实现地震信号的盲反褶积.数值模拟和实际资料处理结果表明,与传统反褶积方法相比,本文方法同时适应于最小相位子波及混合相位子波的反褶积,能够更好地从地震数据中估计反射系数,有效拓宽地震资料的频谱,得到高分辨率的地震资料.     

欧拉反褶积中构造指数的计算与优化选取

构造指数是欧拉反演中最重要的参数,以往选取构造指数往往根据地质先验信息和人工经验试探选取,这样往往不能快速,准确地确定反演构造指数,从而影响反演效果,甚至反演失败.本文提出利用异常值自动估算构造指数,这样避免了经验选取构造指数产生的不确定因素的干扰,可以快速方便地选出较合理的构造指数,提高反演效率.模型实验显示了该方法良好的反演效果,证明了文中提出的方法的可行性.     

基于欧拉反褶积方法估算华北地区重力变化的等效源参数

利用华北地区2009—2014年绝对重力与相对重力多期重复观测资料, 得到不同时空尺度的华北区域重力场动态变化图像。 采用欧拉反褶积方法, 通过对理论模型试算, 获得最优反演参数, 对引起华北地区重力场变化的场源深度和空间分布规律进行了反演计算和分析。 结果表明: 当构造指数为1时, 适合对流动重力场变化数据进行反演。 实际资料的反演结果在2009—2014年间的场源位置集中于河套断裂带。 本文研究方法可以用于重力场反演和以场求源的定量研究。 本工作为流动重力变化信号的反演和解释提供一条新思路, 也为构建地震重力预报指标体系提供了定量依据。     

非稳态地震稀疏约束反褶积研究（英文）

传统Robinson褶积模型主要受缚于三种不合理的假设,即白噪反射系数、最小相位地震子波与稳态假设,而现代反射系数反演方法（如稀疏约束反褶积等）均在前两个假设上寻求突破的同时却忽视了一个重要事实：实际地震信号具有典型的非稳态特征,这直接冲击着反射系数反演中地震子波不随时间变化的这一基础性假设。本文首先通过实际反射系数测试证实,非稳态效应造成重要信息无法得到有效展现,且对深层影响尤为严重。为校正非稳态影响,本文从描述非稳态方面具有普适性的非稳态褶积模型出发,借助对数域的衰减曲线指导检测非稳态影响并以此实现对非稳态均衡与校正。与常规不同,本文利用对数域Gabor反褶积仅移除非稳态影响,而将分离震源子波和反射系数的任务交给具有更符合实际条件的稀疏约束反褶积处理,因此结合两种反褶积技术即可有效解决非稳态特征影响,又能避免反射系数和地震子波理想化假设的不利影响。海上地震资料的应用实际表明,校正非稳态影响有助于恢复更丰富的反射系数信息,使得与地质沉积和构造相关的细节特征得到更加清晰的展现。     

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

常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的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%,该时延主要是由地震同震效应引起的断裂带内应力调整和物性变化造成的.所使用的精密控制人工震源性能稳定,适于开展断层监测.这一结果对于发展断裂带主动震源动态监测技术具有重要意义.