一种加速微震数值模拟的方法及装置制造方法及图纸

本发明专利技术涉及一种微震数值模拟方法及装置，属于微震监测领域，具体是涉及一种加速微震数值模拟的方法及装置。本发明专利技术结合微震监测中震源与接收检波器数量特点，依据源检互换原理，在现有技术基础上，将原震源点作为接收点、原接收点作为震源点，计算地震波场，能够快速得到微震数值模拟中各个接收点的波场信息，与现有技术相比，大大提高了计算效率。

【技术实现步骤摘要】
一种加速微震数值模拟的方法及装置
本专利技术涉及一种微震数值模拟方法及装置，属于微震监测领域，具体是涉及一种加速微震数值模拟的方法及装置。
技术介绍
微震监测技术是近几年在石油、煤田领域应用较广的一项物探技术，被广泛应用于压裂监测、冲击地压监测等领域。微震监测技术主要涉及正演和反演两大技术问题。微震监测正演问题，即已知速度模型、观测系统、震源坐标等参数，求解微震波场(动力学)及地震波走时(运动学)等信息，为微震监测观测系统设计、定位方法正确性验证等提供依据，是微震监测中的基础问题，也是反演的基础。微震监测反演问题，是正演的逆过程，即已知观测系统、地震波场和走时等信息，求取速度模型、震源参数等信息的过程。微震正演方法主要包括：基于波动方程理论的有限差分法、有限元方法、伪谱法、波前快速推进法等，基于积分方程的边界积分法、边界元法等，基于射线追踪理论的试射法、弯曲法等等。常规微震监测正演数值模拟算法的基本流程是：1.给定模拟区域的三维网格速度模型、震源及检波点坐标等参数；2.计算每个震源点在三维模型空间的地震波场(运动学或者动力学)；3.在地震波场中提取每个检波点处的信息。由于微震监测正演是一个三维问题，模型空间较大，计算三维模型空间的地震波场是一个非常耗时的过程，其主要影响因素包括模型空间大小、网格划分尺寸等。若模型空间较大、网格划分较细，计算一个震源点的三维地震波场需要几十分钟甚至几百分钟时间。假设有M个微震震源，计算1个震源需要N分钟时间，完成所有震源点正演模拟的时间大致是M×N分钟。因此，上述第2步是制约微震正演模拟效率的主要因素。
技术实现思路
本专利技术主要是解决现有技术所存在的微震数值模拟(正演)速度较慢的技术问题，提供了一种加速微震数值模拟的方法及装置。该方法及装置结合微震监测中震源与检波点数量特点，依据源检互换原理，能够快速得到微震数值模拟中各个检波点的波场信息，与现有技术相比，大大提高了计算效率。本专利技术的上述技术问题主要是通过下述技术方案得以解决的：一种加速微震数值模拟的方法，具体包括以下步骤：正演模拟参数给定步骤，根据待模拟对象的地质条件，建立网格速度模型，并指定检波点坐标；波场计算步骤，根据震源点个数Ns和检波点个数Nr，确定地震波场计算方式。若Ns＞Nr，进行源检互换及波场计算，即将每个检波点当做震源点，分别计算每个检波点所在点作为震源的地震波场；若Ns＜Nr，则按照常规方法，计算每个震源点地震波场。微震监测中，一般属于前者。波场抽取步骤，基于上述计算得到的三维地震波场，抽取每个震源点的地震波场，正演结束。优选的，所述正演模拟参数给定步骤执行以下子步骤：根据待模拟区域建立XYZ三维直角坐标系；基于所述三维直角坐标系对待模拟区域进行网格划分；基于所述网格划分结果，对网格点进行速度赋值，建立速度模型；指定每个检波点的空间位置。其中，网格速度模型、检波点空间坐标及微震震源点分别记为：网格速度模型：vi,j,k(1≤i≤Nx,1≤j≤Ny,1≤k≤Nz)检波点坐标:(xri,yri,zrri)(1≤i≤Nr)震源点坐标：(xsi,ysi,zsri)(1≤i≤Ns)变量说明：vi,j,k为网格节点(i,j,k)对应的速度，Nx为模型X方向网格数，Ny为模型Y方向网格数，Nz为模型Z方向网格数，xri为第i个检波点的X坐标，yri为第i个检波点的Y坐标，zri为第i个检波点的Z坐标，Nr为检波点个数，xsi为第i个震源点的X坐标，ysi为第i个震源点的Y坐标，zsi为第i个震源点的Z坐标，Ns为震源点个数。优选的，所述波场计算模块执行以下子步骤：首先，根据震源点个数Ns和检波点个数Nr，确定地震波场计算方式。若Ns＜Nr，按常规遍历震源点，若Ns＞Nr，依次遍历各个检波点，将检波点当做震源点，原震源点标记为检波点。其次，计算当前震源点，在已知速度模型下的地震波场。以快速推进算法为例，求解下式二维直角坐标系(x-z)下的程函方程：其中t是地震波走时，t(x,z)表示在坐标点(x,z)的地震波走时，s是慢度，即速度的倒数，s2(x,z)表示在坐标点(x,z)的慢度值的平方，x、z代表二维直角坐标系的2个坐标轴。采用上风差分格式代替上式程函方程中的偏微分格式，即采用上风差分法对上式进行离散，其离散式可表达如下：其中D为差分算子，i、j分别为x方向和z方向的网格节点索引号，max为最大值函数，min为最小值函数。上式可简化为：表示走时t在点(i,j)处x方向上的向前、向后的差分算子，表示走时t在点(i,j)处z方向上的向前、向后的差分算子。为保证计算精度，采用二阶上风差分算子。即：式中Δx、Δz分别为X、Z方向的网格尺寸,ti,j表示在点(i,j)处的走时值。最终，利用目标点(i,j)周围八个节点(i,j±1)，(i,j±2)，(i±1,j)，(i±2,j)处的走时值，共有6种组合方式计算(i,j)点的走时值，6种组合方式的二阶差分计算公式如下：式中，t是地震波走时，s是地震波的慢度，h＝Δx＝Δz为正方形网格尺寸。利用上式计算出多个走时，按照窄带扩展基本原理，在求出的多值解中选取最小走时点，即为目标点走时，计算公式如下：ti,j＝min(t1i,j,t2i,j,t3i,j,t4i,j,t5i,j,t6i,j)遍历完所有节点，即可获得地震波时间场。优选的，所述波场抽取步骤执行以下子步骤：抽取每个检波点的地震波场，即在上步中求得的时间场t(x,z)中找到检波点位置对应的时间，就是原震源点到检波点的地震走时，正演结束。一种加速微震数值模拟的装置，具体包括以下步骤：正演模拟参数给定装置，根据待模拟对象的地质条件，建立网格速度模型，并指定检波点坐标。波场计算装置，根据震源点个数Ns和检波点个数Nr，确定地震波场计算方式。若Ns＞Nr，进行源检互换及波场计算，即将每个检波点当做震源点，分别计算每个检波点所在点作为震源的地震波场；若Ns＜Nr，则按照常规方法，计算每个震源点地震波场。微震监测中，一般属于前者。波场抽取装置，基于上述计算得到的三维地震波场，抽取每个检波点的地震波场，正演结束。优选的，所述正演模拟参数给定模块执行以下子步骤：根据待模拟区域建立XYZ三维直角坐标系；基于所述三维直角坐标系对待模拟区域进行网格划分，基于所述网格划分结果，对网格点进行速度赋值，建立速度模型；并指定各个检波点的空间位置。其中，网格速度模型、检波点空间坐标及潜在震源点分别记为：网格速度模型：vi,j,k(1≤i≤Nx,1≤j≤Ny,1≤k≤Nz)检波点坐标:(xri,yri,zrri)(1≤i≤Nr)震源点坐标：(xsi,ysi,zsri)(1≤i≤Ns)变量说明：vi,j,k为网格节点(i,j本文档来自技高网...

【技术保护点】
1.一种加速微震数值模拟的方法，其特征在于，包括：/n正演参数给定步骤，根据待模拟对象的地质条件，建立网格速度模型，并指定检波器坐标；/n时间场计算步骤，根据震源点个数Ns和检波点个数Nr，确定地震波场计算方式；若Ns＞Nr，进行源检互换及波场计算，即将每个检波点当做震源点，分别计算每个检波点所在点作为震源的地震波场；/n波场抽取步骤，抽取每个检波点的地震波场，正演结束。/n

【技术特征摘要】
1.一种加速微震数值模拟的方法，其特征在于，包括：
正演参数给定步骤，根据待模拟对象的地质条件，建立网格速度模型，并指定检波器坐标；
时间场计算步骤，根据震源点个数Ns和检波点个数Nr，确定地震波场计算方式；若Ns＞Nr，进行源检互换及波场计算，即将每个检波点当做震源点，分别计算每个检波点所在点作为震源的地震波场；
波场抽取步骤，抽取每个检波点的地震波场，正演结束。


2.根据权利要求1所述的一种加速微震数值模拟的方法，其特征在于，所述时间场计算步骤具体包括：
地震波场确定子步骤，计算当前震源点在已知速度模型下地震波场的程函方程；
程函方程计算子步骤，采用上风差分法对程函方程进行离散；
地震走时计算子步骤，利用目标点周围八个方向节点处的走时值计算其地震走时的二阶差分计算公式得到目标点走时的多值解；
目标走时确定子步骤，按照窄带扩展基本原理，在求出的目标点走时的多值解中选取最小走时点，即为目标点走时。


3.根据权利要求2所述的一种加速微震数值模拟的方法，其特征在于，所述地震波场确定子步骤中，基于下式确定地震波场的程函方程：


<...

【专利技术属性】
技术研发人员：王云宏，王保利，崔伟雄，赵朋朋，金丹，
申请(专利权)人：中煤科工集团西安研究院有限公司，
类型：发明
国别省市：陕西;61