基于等效介质模型和频变AVO反演的裂缝储层参数估算方法

由于裂缝性油气藏具有突出的资源潜力和经济效益,利用地震方法对裂缝储层进行精细的定量描述逐渐成为勘探地球物理的关键任务之一。为了克服以往数据驱动类反演方法无法直接获得裂缝参数、而基于静态等效介质模型驱动的反演方法无法描述孔隙内部结构和流体信息的缺点,笔者提出一种基于动态等效介质模型的储层定量描述新方法。该方法通过频变AVO(amplitude variation with offset)理论建立目标函数并使用全局最优化算法反演裂缝参数。一维和二维模型测试证实,由于充分利用了反射系数频变响应对裂缝密度和时间尺度因子的敏感性,反演方法可以对裂缝储层实现有效描述。     

地震波正演问题的物理信息神经网络算法研究进展

偏微分方程作为一种描述客观物理规律的重要工具,可应用于解决大多数科学问题和工程问题.在地球物理学中,求解波动方程是地震波场数值模拟的重要步骤,对于偏移成像和反演起着至关重要的作用.随着勘探目标复杂化,地震波方程的数学描述也更加复杂,现有的经典地震波数值模拟方法存在一些待完善之处,需要一些新型正演算法的应用.物理信息神经网络(Physics-Informed Neural Network,简称PINN)是一种利用神经网络求解偏微分方程近似解的有效手段,通过自动微分将偏微分方程及其定解条件引入损失函数中,在最小化数据驱动误差的同时,使其满足方程的物理约束条件.本文在分析PINN原理并总结其研究现状的基础上,对其在波动方程正演中已有成果进行归纳和分析,展望了其未来可能的研究方向,并对PINN在一维声波方程正演中的应用进行了探究.     

改进的地震AVAZ裂缝弱度反演方法及其在致密砂岩储层中的应用

致密砂岩气储层具有低孔、低渗的特征,裂缝的存在可以提高储层的渗透率,同时裂缝是油气重要的储存空间和运移通道,裂缝发育也有利于水力压裂过程中裂缝网络的形成,裂缝预测可为致密砂岩气储层的开发和部署提供重要依据.地震振幅随方位角的变化可以提供储层中垂直裂缝的信息,本文针对HTI介质提出了一种改进的方位振幅差异反演方法,并结合岩石物理理论预测表征裂缝性质的裂缝弱度参数.常规的反演方法一般同时反演弹性参数和裂缝参数,改进的方位振幅差异方法引入一个参考方位,构建消除各向同性背景的方位振幅差异道集,仅反演与各向异性项相关的裂缝弱度参数,充分利用方位各向异性响应,提高裂缝识别的敏感性与裂缝参数反演的准确性.实际数据应用表明,方位振幅差异反演方法对裂缝参数预测的敏感性较常规方法有所提高,预测的裂缝弱度与测井渗透率曲线相吻合,并且与致密砂岩气储层产气性具有明显的相关性.因此,利用方位振幅差异方法预测裂缝分布及其发育程度可为致密砂岩储层含气有利区的识别与开发提供可靠的指标.     

地震数据Bregman整形迭代插值方法

人工地震方法由于受到野外观测系统和经济因素等的限制,采集的数据在空间方向总是不规则分布.但是,许多地震数据处理技术的应用（如：多次波衰减,偏移和时移地震）都基于空间规则分布条件下的地震数据体.因此,数据插值技术是地震数据处理流程中关键环节之一.失败的插值方法往往会引入虚假信息,给后续处理环节带来严重的影响.迭代插值方法是目前广泛应用的地震数据重建思路,但是常规的迭代插值方法往往很难保证插值精度,并且迭代收敛速度较慢,尤其存在随机噪声的情况下,插值地震道与原始地震道之间存在较大的信噪比差异.因此开发快速的、有效的迭代数据插值方法具有重要的工业价值.本文将地震数据插值归纳为数学基追踪问题,在压缩感知理论框架下,提出新的非线性Bregman整形迭代算法来求解约束最小化问题,同时在迭代过程中提出两种匹配的迭代控制准则,通过有效的稀疏变换对缺失数据进行重建.通过理论模型和实际数据测试本文方法,并且与常规迭代插值算法进行比较,结果表明Bregman整形迭代插值方法能够更加有效地恢复含有随机噪声的缺失地震信息.     

“松科二井”邻域岩石圈精细结构特征及动力学环境——深地震反射剖面的揭示

为揭示"松科二井"邻域岩石圈精细结构特征,布设了一条过井近南北向的深地震反射剖面,采用多尺度药量结合的激发技术,长排列、重点段加密接收高次覆盖的采集方式,通过高保真、高保幅的处理流程获得高分辨率、高信噪比的叠前时间偏移剖面.结果显示:T4反射轴之下沉积岩层状反射与火成岩杂乱反射相间,有两处疑似上古生界地层,分别位于松科二井下方(双程走时,Two way travel time,TWT,3.5～4s)和任民—永安断隆下方(CDP,3500～4500,TWT,3～4s);中下地壳可观察到近平行的北倾反射,中和断陷下方存在上下关联的透镜状强反射,整体轮廓呈蘑菇云状,解释其为岩石圈伸展构造中发育的深部热流底辟体;北部莫霍面呈近水平连续强反射,南部徐家围子断陷区域莫霍面反射较弱,同时剖面上存在3种明显的岩石圈上地幔反射,包括倾斜地幔反射、近水平地幔反射和超深地幔反射,推断其分别为早期俯冲遗迹、早期增厚地壳底界面以及现今岩石圈底界面.本文利用过松科二井地震剖面的最新成果,揭示出古亚洲洋、蒙古—鄂霍茨克洋和古太平洋三大构造域先后作用下,保留在松嫩地块岩石圈内的结构特征,为探讨松辽盆地形成原因、构造背景及动力学因素提供新视野.     

基于岩石物理参数转换模型的页岩气含气量评价新方法研究

页岩气在储层中主要以游离气和吸附气的形式存在,其含气量主要通过现场解吸和测井方法计算.由于页岩储层十分致密,并且常伴有微裂缝发育,这就造成了页岩储层的储集空间从纳米级孔隙到微米-毫米级的微裂缝并存的特殊情况.针对微孔和介孔段(50nm以内)的等温吸附方法,以及针对宏孔(50nm以上)的高压压汞方法在表征页岩孔隙结构上都存在着局限性.而核磁(NMR)T2分布本身虽然能够表征孔径分布,但是需要精确刻度横向弛豫时间与孔径的关系.本文在现有方法的基础上,结合上述方法的优势,采用拼接的方式表征页岩整体孔径分布,并将核磁T2分布刻度为孔径分布.基于横向弛豫时间和孔径的关系,结合T2-I转换模型,采用动态参数计算页岩储层含水饱和度和含气量,这为页岩储层含气量评价提供了新的可靠的方法.     

2017年9月23日朝鲜ML3.4地震矩张量反演

利用区域地震台网数字波形资料,对2017年9月23日朝鲜M_L3.4地震进行地震矩张量反演计算与参数稳定性评估,获得了此次地震的震源机制解.结果表明,地震矩心深度为3 km,标量地震矩为1.34×10¹⁴ N·m,矩震级为M_W3.4.地震矩张量结果分解后,双力偶分量（DC）为96.4%,补偿线性矢量偶极分量（CLVD）为-0.8%,震源体体积变化的各向同性分量（ISO）为-2.8%.主压应力P轴方位角为144°,倾角为74°,主张应力T轴方位角为341°,倾角为15°.其中一个节面的参数为：走向248°,倾向60°,滑动角-94°.地震震源体积变化分量很小,震源机制类型属于典型的由断层剪切位错引起的正断层型地震事件,且主张应力T轴方向与区域近地表应变率场方向一致.由于朝鲜2017年9月3日核试验释放的能量对局部区域应力场进行了扰动,致使核试验场附近地壳岩体处于破裂的临界状态,2017年9月23日朝鲜M_L3.4地震事件可能是区域应力场作用下的一次山体滑动事件.