基于方位各向异性的裂缝密度反演方法及应用

裂缝预测是页岩气甜点预测的重要内容.页岩气储层高角度裂缝发育,具有明显的HTI介质各向异性特征,基于各向异性参数与裂缝密度的相关性,提出了一种裂缝密度反演的新方法.首先推导了一种能反映裂缝密度的方位AVO近似公式及其方位弹性阻抗方程（EVAAz方程）;然后应用贝叶斯理论,形成了组稀疏方位弹性阻抗反演目标泛函,在此基础上通过阻尼最小二乘算法求解弹性矩阵,稳健得到了剪切模量、法向弱度与切向弱度;最后应用裂缝弱度与裂缝密度的关系,实现了裂缝密度的定量预测.经实际数据试验,预测结果与页岩压裂参数吻合度高,证实该方法是有效的.

     

基于深度卷积神经网络的2021年5月21-22日云南漾濞地震和青海玛多地震震级估算

以探索深度卷积神经网络震级估算模型对2021年5月21日云南漾濞和5月22日青海玛多地震预警震级估算的可行性为目标,本文使用P波到达后不同时间窗下的特征参数作为深度卷积神经网络输入进行训练和验证,构建了1~40 s内不同时间窗下的深度卷积神经网络单台震级估算模型,并采用多台加权平均方法对玛多地震主震以及漾濞地震主震、前震和余震共11次地震事件进行实时震级估算.结果表明：在P波到达后1~40 s的不同时间窗下,随着时间窗的增加,深度卷积神经网络单台震级估算模型对训练集和验证集的震级估算误差和标准差逐渐减小,且逐渐趋于稳定;对于云南漾濞M_S6.4主震与青海玛多M_S7.4主震,大部分单台震级估算误差主要分布在±0.5震级单位内;通过多台加权平均方法计算预测震级,漾濞主震与玛多主震在首台触发后2 s时,震级估算结果分别是M6.0和M7.0;同时,随着首台触发后时间的增加,漾濞与玛多主震的震级估算逐渐接近实际震级,且在首台触发后9 s时,漾濞主震和玛多主震的震级估算分别是M6.2和M7.3;在首台触发后1 s时,对于玛多与漾濞地震的前震和余震,多台加权平均震级估算结果也主要分布在±0.5震级单位误差范围内.该方法为中国地震预警系统的震级估算方案提供了潜在可能.

     

基于Metropolis抽样的弹性阻抗随机反演

笔者提出了基于Metropolis抽样的弹性阻抗随机反演方法,该方法是一种基于蒙特卡洛的非线性反演方法,能够有效地融合测井资料中的高频信息,提高反演结果的分辨率。应用贝叶斯理论框架,首先通过快速傅里叶滑动平均模拟算法( fast Fourier transform?moving average,FFT?MA)和逐渐变形算法( gradual deformation method,GDM)得到基于地质统计学的先验信息,然后用Metropolis抽样算法对后验概率密度进行抽样,得到反演问题的解。其中FFT?MA模拟作为一种高效的频率域模拟方法,结合GDM后,在保持模拟空间结构不变的前提下,可以连续修改储层模型,直至满足实际观测地震记录。数值试验表明:FFT?MA模拟有效地提高了计算效率,融入GDM更新算法,可以保证反演结果有效地收敛,并且反演结果与理论模型吻合较好,具有较高的分辨率,该方法采用两步法反演弹性参数,运算速度较快。     

大地电磁资料精细处理和二维反演解释技术研究（四）——阻抗张量分解的多测点-多频点统计成像分析

基于大地电磁阻抗张量分解技术,本文提出了两种电性主轴方位的统计描述图像：随频率变化的统计分布成像（频率分布云图）和随测点序列变化的统计分布成像（测点分布云图）.这两种图像与传统的统计玫瑰图一起,较全面地描述了最佳主轴的分布特征.在进行构造维性分析过程中,通过定义二维有效因子e2d,来压制一维结构和三维结构、突出纯二维结构的影响.e2d被用于电性主轴的统计加权,有效地起到了滤波的作用;同时,统计成像中还考虑了数据质量的影响.为了得到稳定、高质量的区域阻抗张量数据,提出并实现了共主轴的多测点-多频点阻抗张量分解新算法.最终,完成了以上各项处理手段的可视化实现.本文通过两个理论模型和一个实测算例,以共轭阻抗法（CCZ法）为基础,展示了这一新技术的有效性.     

基于弹性阻抗的储层物性参数预测方法

储层物性参数是储层描述的重要参数,常规的基于贝叶斯理论的储层物性参数反演方法大多是通过反演获得的弹性参数进一步转换而获得物性参数,本文提出一种基于弹性阻抗数据预测储层物性参数的反演方法.该方法主要通过建立可以表征弹性阻抗与储层物性参数之间关系的统计岩石物理模型,联合蒙特卡罗仿真模拟技术,在贝叶斯理论框架的指导下,应用期望最大化算法估计物性参数的后验概率分布,最终实现储层物性参数反演.经过模型测试和实际资料的处理,其结果表明本文提出的方法具有预测精度高,稳定性强,横向连续性好等优点.     

资源储量估算过程中矿业权边界的确定

在垂直纵投影图上如何确定矿业权边界投影位置, 对准确圈定资源储量估算范围具有重要的现实意义.在充分研究资源储量估算过程的基础上, 结合MapGIS、AutoCAD等常用地质绘图软件, 通过构建矿体三维数学模型, 进行三角函数公式推导, 得出资源储量垂直纵投影图上确定矿业权边界投影位置的公式, 可快速、准确地在资源储量垂直纵投影图上确定矿业权边界的投影位置.这一方法既能为合理开采、利用矿业权范围内资源储量提供数据支持, 也能为矿政管理部门核准、监管矿业权范围内的资源储量提供技术支持.     

强震区岷江沿岸泥石流物源体积估算模型与演变特征

2008年汶川大地震后,泥石流灾害表现活跃,至今经历了多个暴雨过程,那么如何才能较为准确地分析泥石流源区演变特征。首要是确定泥石流源区崩滑体体积。本文选择映秀镇附近,岷江沿岸的5条泥石流沟为研究对象,针对泥石流源区崩滑体体积难以确定的问题,野外调查研究区内167个崩滑体点,建立崩滑体体积与面积间的关系模型,并通过野外调查研究区内48个崩滑点进行验证。定量解译研究区3期航空影像的物源面积,利用崩滑体体积估算模型,进一步分析研究区泥石流物源的演变特征。     

大地电磁资料精细处理和二维反演解释技术研究(五)--利用阻抗张量成像识别大地线性构造

大地线性构造往往与块体或物性的边界、断裂的几何延伸相关联,是地学研究的重要目标构造.在大地电磁法中,传统的方法是通过反演得到的电性结构的差异来解译线性构造,是一种间接分析方法.由于线性构造发育处走向和倾向明显,二维性较其他地方一般会更强,故而在阻抗张量成像分析中,相应的二维有效因子就会比较大,最佳主轴统计分布会更加明显.据此,本文利用大地电磁阻抗张量成像（包括最佳主轴统计成像、构造维性成像）来探测识别大地线性构造.本文首先通过三维理论模型的研究结果来验证这一技术在线性结构探测中的有效性,然后将其应用于郯庐断裂带南段实测剖面的具体分析中,展示了这一技术在线性构造探测中的实用性及其数据处理与分析过程.本文研究表明,大地电磁阻抗张量成像技术可独立于反演为地质解释尤其是线性构造的解释提供诸多参考信息.

     

基于EEI反演的储层油气识别研究

埕岛油田是我国极浅海地区第一个投入开发的大油田,其馆陶组上段属河流相砂泥岩沉积,油气富集程度高。该层段分为五套砂组,其中三砂组和四砂组是主要的含油储层。此次研究的目标工区内主要的烃类是重质油,含油储层与含水储层的物性差异较小,通过反演常规的弹性参数无法准确指示油气。基于测井曲线和岩石物理交汇分析,本文构建了敏感的油气指示因子,通过扩展的弹性阻抗反演技术反演该指示因子从而进行油气展布预测。该方法在工区得到了良好的应用效果。