模拟微地震的弹性波正演方法技术

本发明专利技术提供一种模拟微地震的弹性波正演方法，包括：步骤1，根据地质勘测的实际资料，分析岩石的性质，确定岩石破裂的模式，生成表达微地震源的矩张量矩阵；步骤2，建立所需勘探区域的地质速度场模型；步骤3，波动理论下，建立适用于描述微地震波场传播的弹性波波动方程；步骤4，在波动理论下，求解微地震弹性波波动方程，实现微地震的波场传播模拟；步骤5，模拟完成，输出微地震地面、井中多分量记录和纵横波的能量辐射模式。该模拟微地震的弹性波正演方法能够快速的得到较射线理论更精确的模拟结果，且信息量更为丰富，并为后续的特征规律分析、定位方法研究以及震源参数反演提供更好的技术手段和工具。

【技术实现步骤摘要】

本专利技术涉及油田开发
，特别是涉及到一种模拟微地震的弹性波正演方法。
技术介绍
进入21世纪以来，非常规油气资源表现出巨大的开发潜力。在众多非常规油气的勘探和开发方法中，微地震监测技术起着举足轻重的作用。为了充分掌握微地震技术，需要对微地震的正问题、震源机制进行深入的了解和研究，而实现这一目标最便捷的途径便是地震波数值模拟技术。常规的地震正演技术使用地震子波作为震源，仅适合常规地震勘探的爆炸震源或可控震源，而不能很好的对微地震震源进行描述。目前微地震正演模拟的一种主要方式是基于射线理论，在给点介质模型和震源位置前提下通过射线追踪方式得到微地震记录。但射线理论下模拟的微地震记录无法再现不同频率成分对介质的不同的响应，记录中缺少首波信息，且只能适用于光滑介质。不仅如此，实际中观测到的微地震源具有很宽的频带范围且依赖于介质的变化情况，因此高频射线的假设严重影响了微地震数值模拟的适用范围。另一种方式基于波动理论，使用天然地震的矩张量对岩石破裂模式进行刻画，并使用Green函数法进行弹性波传播数值模拟。但Green函数法求取过程中需要多次求取格林函数，计算量较大，且求取格林函数需要对其进行近似展开，这样就降低了模拟的精度。为此我们专利技术了一种新的模拟微地震的弹性波正演方法，解决了以上技术问题。
技术实现思路
本专利技术的目的是提供一种引入天然地震中描述刻画震源的矩张量概念，实现对压裂震源的刻画的模拟微地震的弹性波正演方法。本专利技术的目的可通过如下技术措施来实现：模拟微地震的弹性波正演方法，该模拟微地震的弹性波正演方法包括：步骤1，根据地质勘测的实际资料，分析岩石的性质，确定岩石破裂的模式，生成表达微地震源的矩张量矩阵；步骤2，建立所需勘探区域的地质速度场模型；步骤3，波动理论下，建立适用于描述微地震波场传播的弹性波波动方程；步骤4，在波动理论下，求解微地震弹性波波动方程，实现微地震的波场传播模拟；步骤5，模拟完成，输出微地震地面、井中多分量记录和纵横波的能量辐射模式。本专利技术的目的还可通过如下技术措施来实现：在步骤1中，通过对水力压裂过程中接收到的实际微地震数据的分析，获得微地震的量级，从而生成表达微地震源的矩张量矩阵 M = M 11 M 12 M 13 M 21 M 22 M 23 M 31 M 32 M 33 ]]>其中Mij为岩石破裂量级，矩张量中的每个元素代表一对方向相反的力的作用。在步骤2中，读取预先建立的所需勘探区域的地质速度场模型信息，建立所需勘探区域的地质速度场模型，即完成区域速度场建模。在步骤3中，建立的微地震弹性波波动方程为： ρ ∂ 2 u i ( x , t ) ∂ t 2 - c i j k l ∂ 2 u k ( x , t ) ∂ x j ∂ x l = M i j ∂ [ δ ( x - x s ) ] ∂ x j S ( t ) , ]]>其中ui(x,t)和uk(x,t)为波场位移，cijkl为刚度系数，ρ为介质密度，x为空间点坐标，t为时间，三维情况下，i,j,k,l＝x,y,z，xs为震源点位置，S(t)为震源子波，M为矩张量。在步骤3中，利用步骤1获得的矩张量，代入微地震波动方程，使用等效体力法对方程右边的震源项进行等效离散。在步骤4中，利用交错网格高阶有限差分技术、伪谱法、有限元法、谱元法或迭代求解法求解微地震弹性波波动方程。本专利技术中的模拟微地震的弹性波正演方法，克服了射线理论高频假设的缺陷，避免Green函数法计算量大、模拟精度低的缺点，在充分考虑不同压裂本文档来自技高网...

【技术保护点】
模拟微地震的弹性波正演方法，其特征在于，该模拟微地震的弹性波正演方法包括：步骤1，根据地质勘测的实际资料，分析岩石的性质，确定岩石破裂的模式，生成表达微地震源的矩张量矩阵；步骤2，建立所需勘探区域的地质速度场模型；步骤3，波动理论下，建立适用于描述微地震波场传播的弹性波波动方程；步骤4，在波动理论下，求解微地震弹性波波动方程，实现微地震的波场传播模拟；步骤5，模拟完成，输出微地震地面、井中多分量记录和纵横波的能量辐射模式。

【技术特征摘要】
1.模拟微地震的弹性波正演方法，其特征在于，该模拟微地震的弹性波正演方法包括：步骤1，根据地质勘测的实际资料，分析岩石的性质，确定岩石破裂的模式，生成表达微地震源的矩张量矩阵；步骤2，建立所需勘探区域的地质速度场模型；步骤3，波动理论下，建立适用于描述微地震波场传播的弹性波波动方程；步骤4，在波动理论下，求解微地震弹性波波动方程，实现微地震的波场传播模拟；步骤5，模拟完成，输出微地震地面、井中多分量记录和纵横波的能量辐射模式。2.根据权利要求1所述的模拟微地震的弹性波正演方法，其特征在于，在步骤1中，通过对水力压裂过程中接收到的实际微地震数据的分析，获得微地震的量级，从而生成表达微地震源的矩张量矩阵 M = M 11 M 12 M 13 M 21 M 22 ...

【专利技术属性】
技术研发人员：隆文韬，廉西猛，王延光，黄超，王修敏，陈云峰，张猛，王修银，王贤真，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：山东;37