强震震源破裂过程浅析

通过了解强震震源区域周边地带地震破裂过程及破裂长度,根据各地震台网记录的地震数据设定计算模型,对地表或深井中的破裂效应进行研究,得出微破裂或破裂链形成的动力过程,可对地震进行一定预测。     

添加正则化项的黏声逆时偏移成像方法研究

在实际地球介质中传播的地震波会产生衰减和频散现象,因此其更接近于黏弹性介质,在地震处理中补偿这些黏性影响是十分必要的。基于波动方程的叠前深度偏移中进行吸收衰减补偿更准确,也更有物理意义,直接求解双程波动方程的逆时偏移(RTM)能够成像大倾角复杂构造,具有诸多优势。然而当考虑吸收衰减补偿时通常会产生不稳定现象,大部分研究都是在逆时偏移的波场延拓中进行波数域的低通滤波来解决这个问题。本文采用广义标准线性固体的黏声波动方程进行吸收衰减补偿的Q--RTM方法,通过添加正则化项的方式来稳定延拓过程。添加正则化项本质上是低通滤波,滤波窗口是指数形式,在时空域有明确的形式,可以阻止发生高频不稳定。与直接滤波相比,正则化参数可以是空变的,因此比较适合剧烈变化的区域,我们还发现震源归一化互相关成像条件更适合Q--RTM方法。      

利用重复地震研究盖州震群地区介质变化

盖州震群自2012年2月2日开始活动,截至2016年12月已活动了近5年的时间,其中持续活动的青石岭序列精定位后的空间展布在5km以内,深度也较为集中,可作为重复地震进行研究。通过计算盖州地区2012年2月—2016年12月M_L≥2.0重复地震序列经射线追踪法校正后的P波走时差,得到了盖州震群区域至岫岩台、辽阳台、高升台、瓦房店台四象限的介质变化情况,结果介质短期的变化受震群本身显著事件活动影响较为明显。中强地震前后周边地区介质变化显著。     

含流体致密砂岩的纵波频散及衰减:基于双重双重孔隙结构模型描述的特征分析

致密砂岩普遍具有低孔、低渗及微裂缝发育的地质特征,并且呈现出很强的非均匀性.致密砂岩储层与常规砂岩储层比较,具有明显的岩石物理性质、渗流力学性质方面的差异.致密砂岩内部的非均匀性对弹性波频散、衰减有显著影响,其中包括孔隙结构的非均匀性,即岩石内部孔隙参数的不均一性,以及孔隙内部不相混溶流体的非均匀分布;此外,非均匀性的尺度也决定了波出现显著频散与衰减的频段.综合考虑致密砂岩孔隙结构非均匀性及流体斑块状饱和的非均匀性,本文采用双双重孔隙介质结构模拟了致密砂岩的弹性波响应,分析了同时具备两类非均质性岩石中的波传播特征.调查分析了两组分别来自中国鄂尔多斯盆地苏里格气田及四川盆地广安气田的不同类型致密砂岩储层的岩芯超声波实验数据,给出了岩石样本的弹性波速度频散与衰减曲线.结果显示理论模型预测结果与完全饱和、部分饱和岩石的实验数据吻合良好.对两个地区致密砂岩岩芯数据进行对比分析,苏里格致密砂岩样本总体上比广安致密砂岩渗透率高,在各孔隙度范围内,特征模拟显示苏里格样本的裂隙尺寸明显大于广安样本.广安致密砂岩在低孔隙度范围内发育了更多、更小的颗粒裂隙/接触.致密砂岩的速度频散与衰减结果受流体黏度、晶体破裂及流体斑块状饱和的共同影响.此外,孔隙度越大,部分饱和岩石中斑块状饱和机制对总衰减的贡献越低,与之相对,结构非均质性所占的比重则有所增强.     

含孔隙介质的分层半空间表面瑞利波的衰减特性

以分层半空间内部含有一层孔隙介质为物理模型进行数值计算,研究半空间表面瑞利波的传播和衰减特性.为更加接近实际,结合瑞利波的激发特性,确定了瑞利波的主衰减曲线,并主要以此进行规律分析.针对速度递增和含低速层这两种典型的地质模型,讨论了瑞利波的传播衰减在不同地质模型下的特性,并分析了各自的规律.结果表明,在这两种模型下瑞利波的主衰减曲线都受孔隙介质所处空间位置影响产生比较明显的变化,但衰减系数极大值对应的波长与模型的表层厚度存在较明显的线性对应关系,利用这一关系,可以在实际勘探中快速得到表层介质厚度.另外,通过对比分析还可以看到,瑞利波主衰减曲线随孔隙介质的孔隙度和渗透率的变化都强于主频散曲线的变化,表明衰减曲线对孔隙度和渗透率的变化更加敏感,理论上更加适合进行介质参数反演工作.综合对比结果,我们认为瑞利波主衰减曲线中包含了更丰富的介质参数信息,如果能够有效利用,将可以提高瑞利波勘探的准确性和应用范围.     

基于FFT-MA的随机介质建模方法

随机介质是描述地球介质小尺度非均匀性的有效模型,在地震散射波场分析、储层描述等领域具有广泛应用,快速准确的随机介质建模方法是开展相关研究的基础与前提.本文首次将FFT-MA（Fast Fourier Transform Moving Average）算法引入到随机介质建模研究中,分析了该方法与传统随机介质建模方法相比具有的优势,并提出了基于FFT-MA的非平稳随机介质建模方法.建模实验表明,与传统的基于谱分解定理的随机介质建模方法相比,基于FFT-MA的方法在空间域产生随机数序列,使随机数序列与结构参数分离,因此,随机数序列与所建模型存在空间上的对应关系,可以分区域建模和局部修改模型.在非平稳随机介质模型建模时,滑动窗口的大小可以根据自相关长度变化而变化,避免了每个采样点都建立一次完整大小的模型,提高了建模效率.因此,FFT-MA随机介质建模方法能准确构建满足自相关函数要求的平稳及非平稳随机介质模型,具有建模效率高、灵活、实用的优点.     

用柱状成层各向异性介质的并矢Green函数模拟偏心对多分量感应测井响应的影响

采用偏心状态下柱状成层各向异性(横向各向同性)介质中并矢Green函数的解析表达式高效模拟多分量感应测井仪器在井眼中偏心时的响应.为提高精度,在模拟时考虑到了金属心轴、绝缘保护层的存在以及各分量线圈系的具体形状.数值模拟结果表明,当井眼内钻井液电导率相对较高、地层电导率相对较低时,偏心对仪器响应的影响较大,尤其是对短线圈距线圈系的影响更为明显,必须进行偏心效应校正.当钻井液电导率相对较低时,偏心对线圈系响应的影响可忽略不计.对位于相对低电导率井眼中的线圈系而言仪器方位角的影响可忽略不计,而当线圈系位于相对高电导率井眼中时仪器方位角的影响极为明显.当仪器偏心率较小时线圈系的响应随仪器方位角的变化较小,仪器偏心率越大线圈系的响应随仪器方位角的变化越明显.     

巴颜喀拉块体地壳结构多样性探测

青藏高原内部地壳岩性的改造、岩性随深度变化及形变构造是探索研究地壳增厚、物质运动问题的关键.巴颜喀拉块体位于青藏高原中北部,地域广袤,通过对块体内中、东部不同区域的深地震广角反射/折射震相的综合分析,利用反射率理论地震图方法对不同性质震相走时及振幅特征进行细致的模拟计算,进一步研究巴颜喀拉块体内部不同区域地壳精细结构.结果显示:巴颜喀拉块体地壳厚度50~60 km、整体向西逐渐增厚,结晶地壳平均速度6.07~6.18 km·s~(-1)、岩层速度大幅降低,壳内多强反射界面结构、但不同区域差异明显;东部若尔盖盆地地壳介质速度整体低速、壳幔边界较为清晰;中部玉树一玛多段下部地壳发现约6.8 km·s~(-1)的"高速度"介质结构,壳幔边界不清、被改造为2~4 km厚的高速度梯度层,显示了巴颜喀拉块体内部地壳增厚、介质岩性结构被改造的差异性.地壳内部多组强反射、低视速度走时震相揭示了介质岩性的低速破碎、弱化蠕变以及可能的壳内解耦构造.局部地区下地壳的高视速度震相特征显示了青藏高原地壳改造增厚大背景下可能存在稳定的"原始地壳"结构残留或是与上地幔物质的浸入交流.巴颜喀拉块体内不同区域地壳增厚、岩性结构、结晶基底及壳幔边界性质被改造的多样性为深入认识青藏高原地壳形变及动力学过程带来新的启迪.     

基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移

近年来,面向实际应用的TI介质准P波正演模拟与逆时偏移成像技术受到空前的关注.基于常规耦合型传播方程的正演模拟方法不仅存在伪横波及频散假象干扰,而且还遭受模型参数限制(η0)和不稳定影响;而纯qP波方程的推导繁琐,且由于方程中包含拟微分算子造成求解难度大且精度有限.为此,本文首先构建了一种适用于任意TI介质的纯qP波传播算子,然后借助Low-rank分解求取该算子中的空间-波数域矩阵,同时引入Cerjan衰减边界条件来压制边界反射干扰,最终实现了一种间接的纯qP波波场外推方案,并将其成功应用于复杂TI介质正演模拟与逆时偏移成像中.通过开展数值模拟,并与其他方法对比表明:①该方法既避免了纯qP波方程的繁琐推导,又克服了耦合型方程对模型参数的限制;②还彻底消除了残余伪横波噪音及数值频散;③且能适应较大时间或空间步长及高频震源,是一种相对准确且稳定的各向异性纵波正演与成像策略.     

大同地震窗P波频散衰减的异常变化特征

通过测量介质品质因子来开展大同窗附近测震台站记录P波衰减频散特征的研究,并分析大同窗近场区和外围区震例前后的变化。对山自皂台记录的2001年以来大同窗地震序列的研究结果表明,近场区3次和外围区6次地震前,介质品质因子显示超出1倍均方差的高值变化特征,震后逐渐恢复到背景起伏状态;与大同窗传统的频次指标相比,介质品质因子的映震效果更好;阳原台的结果表明,仅近场区1次和外围区2次地震前,介质品质因子显示超出1倍均方差的高值变化特征,且有虚报、漏报。所处地质构造位置不同、传播路径介质差异可能是引起上述两个台站附近介质品质因子的地震对应效果明显不同的原因。