CN107526104A - 基于多机多核的裂缝介质地震波场数值模拟方法 - Google Patents
基于多机多核的裂缝介质地震波场数值模拟方法 Download PDFInfo
- Publication number
- CN107526104A CN107526104A CN201710749074.9A CN201710749074A CN107526104A CN 107526104 A CN107526104 A CN 107526104A CN 201710749074 A CN201710749074 A CN 201710749074A CN 107526104 A CN107526104 A CN 107526104A
- Authority
- CN
- China
- Prior art keywords
- numerical simulation
- mpi
- seismic wave
- wave field
- parallel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000004088 simulation Methods 0.000 title claims abstract description 32
- 230000008569 process Effects 0.000 claims description 33
- 230000010363 phase shift Effects 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 8
- 238000009826 distribution Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 4
- 238000003825 pressing Methods 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 208000010392 Bone Fractures Diseases 0.000 description 18
- 206010017076 Fracture Diseases 0.000 description 18
- 230000004044 response Effects 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种基于多机多核的裂缝介质地震波场数值模拟方法,应用于裂缝性地震波场领域,通过在波动方程数值模拟中引入集群处理,通过MPI消息传递接口在多机上面工作,使得计算机集群应用在了细网格化的裂缝介质地震波场数值模拟中,极大得提高了运算速度;并且结合多线程的应用,充分利用了单个节点的多核性质,充分利用了计算机的内存和闲置资源,又极大一步提升了算法的加速比。
Description
技术领域
本发明属于裂缝性地震波场领域,特别涉及一种裂缝性储层地质模型。
背景技术
裂缝油气藏已经成为全球石油、天然气储量的重要组成部分和增储重点领域。利用地 面地震资料预测裂缝油气藏一直被地球科学家和油藏工程师视为颇具挑战的研究内容,其 中涉及的关键问题之一是不同形状、密度、大小尺度裂缝的地震波场响应特征难以识别。 在当今国际低油价的环境下,清晰地认识裂缝产生的地震波场响应特征对提高裂缝型油气 藏预测精度和裂缝油气藏钻井成功率至关重要。地震波场数值模拟是人们正确了解地震波 的传播规律和地下介质构造、结构和岩性特征的有效手段,具有正确的理论和实际意义, 是利用地震勘探方法识别和描述油气储层的前提和基础。
常用的地震波场数值模拟方法有射线追踪和波动方程法。波动方程法能够全面地反映 地震波的运动学和动力学两方面的特征,并且采用合理的算法,能够获得高精度的地震波 场正演模拟记录,尤其是针对裂缝介质。众所周知,地层中的裂缝具有多尺度性,大小介 于几微米至几十米之间。按照常规方法只能模拟大尺度的裂缝,对于小尺度的裂缝无能为 力。通常裂缝介质的数值模拟均是基于等效介质理论进行完成,即需要根据裂缝介质的性 质(密度、充填物等)将其划分为若干裂缝发育带,并假定位于同一裂缝发育带内的所有裂缝和裂缝之间的围体介质具有相同的性质。尽管采用等效介质理论简化了裂缝介质地震波场数值模拟的复杂度,但是数值模拟结果只能获得裂缝发育地质体边界处的地震波场响应特征,不能有效地分析不同形状、大小裂缝对地震波场的响应特征,不能有效地指导后续裂缝储层的精细预测和描述。
发明内容
为解决上述技术问题,本申请提出了一种基于多机多核的裂缝介质地震波场数值模拟 方法,通过在波动方程数值模拟中引入集群处理,通过MPI消息传递接口在多机上面工作, 使得计算机集群应用在了细网格化的裂缝介质地震波场数值模拟中,极大得提高了运算速 度。
本申请采用的技术方案为:基于多机多核的裂缝介质地震波场数值模拟方法,包括:
S1、根据裂缝的方向、形状、体积大小,找出最小裂缝单元,采用网格细分技术构建地质模型;
S2、采用相移校正和相移插值结合对构建的地质模型进行地震波场数值模拟;
S3、在波动方程频率域采用MPI与OpenMP并行的方式,实现地震波场数值模拟,具体 为:在集群的计算节点间采用MPI的粗粒度的多进程并行模式,在计算节点内采用OpenMP 的按频率细分的细粒度的多线程并行模式。
进一步地,步骤S3所述在集群的计算节点间采用MPI的粗粒度的多进程并行模式具体 为:
A1、调用MPI_INIT_THREAD启动并行环境;
A2、根据计算节点的数量,调用MPI_Comm_size命令确定进程个数,同时调用 MPI_Comm_rank命令获得各进程的编号;
A3、其他进程向主进程发送请求需要处理的频率块的范围,主进程接收到请求后向其 他进程分发需要处理的频率块范围;
A4、各进程执行波动方程正演数值模拟并行计算;
A5、将成像结果文件发送到主进程,由主进程进行合并;
A6、调用MPI_Finalize,结束并行环境。
进一步地,步骤S3所述在计算节点内采用OpenMP的按频率细分的细粒度多线程并行模式,具体为:
B1、根据当前结点的CPU核数,调用int_omp_set_num_threads函数设置当前节点要 启动的线程数目;
B2、在频率循环外层调用#pragma omp parallel num_threads()语句,将频率波场数值模 拟设定为并行计算区域;
B3、线程私有变量申请;
B4、在频率循环外层通过调用omp_get thread_num函数,得到当前的线程号,为每个 线程分配起始和终止频率。
本发明的有益效果:本发明的基于多机多核的裂缝介质地震波场数值模拟方法,通过 MPI+OpenMP的混合并行模式,使得计算机内存和核数都得到充分利用来进行加速。通过 运用集群来实现MPI+OpenMP混合模型,使得计算效率得到了极大提升,可以达到以下有益效果:
1)通过在波动方程数值模拟中引入集群处理,通过MPI消息传递接口在多机上面工 作,使得计算机集群应用在了细网格化的裂缝介质地震波场数值模拟中,极大得提高了运 算速度;
2)再次应用多线程在上述过程中,充分利用了单个节点的多核性质,充分利用了计 算机的内存和闲置资源,又极大一步提升了算法的加速比;
3)从物理意义上讲,本发明该方法能够考虑不同裂缝方向、形状、体积大小等对地震波场的影响,得到的数值模拟结果更加符合实际情况,为裂缝储层的地震波场特征研究提供了一种的思路和方法。
附图说明
图1为本申请的方案流程图;
图2为本申请实施例提供的等效介质理论模拟与本申请方法的效果对比图;
其中,图2(a)为等效介质地震记录;图2(b)为本申请方法的地震记录;
图3为本申请实施例提供的偏移剖面;
图4为本申请实施例提供的各平台运行时间对比。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐 释。
如图1所示为本申请的方案流程图,本申请的技术方案为:基于多机多核的裂缝介质 地震波场数值模拟方法,包括:
S1、根据裂缝的方向、形状、体积大小,找出最小裂缝单元,采用网格细分技术构建地质模型;
S2、综合考虑裂缝发育区内裂缝与围岩介质的岩石物理性质,根据适合于速度横向变 化大的波动方程数值模拟理论,提出采用相移校正和相移插值结合对构建的地质模型进行 地震波场数值模拟;
S3、在程序的实现过程中,考虑到对采用网格细分后的地质模型进行波动方程数值模 拟计算量巨大,在波动方程频率域采用MPI与OpenMP并行的方式,实现地震波场数值模拟, 具体为:在集群的计算节点间采用MPI的粗粒度的多进程并行模式,在计算节点内采用 OpenMP的按频率细分的细粒度多线程并行模式。
步骤S2中相移校正和相移插值结合的波场延拓算子具体为:
从相移校正方法出发,三维零炮检距地震数据的波动方程如下:
当介质为非均匀的情况下,慢度场可分解成两部分:
在(2)式中,u0(z)为参考慢度,而慢度的变化都归入Δu(x,y,z)分量内。把(2)式代入(1)式并对时间t作傅立叶变换,得到:
或
式(3)中:4w2(2u0Δu+Δu2)p=s(x,y,z)(即为非均匀介质的波动方程)称为源函数。
这时忽略Δu2的作用来求解方程(4),得到另一个相移项eiwΔu(x,y,z)Δz,该相移项就是近 似反映横向慢度的扰动源项。相移加校正的波长外推方法可简单地概括为(5)式描述:
p(x,y,zn+1,w)=eiwΔuΔz2DDFT-1(2DFFT(P(x,y,znw))eiθΔz) (5)
在(5)式中:
为了适应速度波场横向剧烈的变化,对(5)式采用插值的方法,设在Δz的层中最小 速度为v1,最大速度为v2,则有慢度:
在Δz层段中用u1和u2进行延拓, u(x,y,z)是Δz间隔内的实际慢度。
其中,为空间域中的相移校正项。
实际波场由以下插值公式得到:
最后波场成像结果为:
尽管根据以上公式能够实现裂缝介质的地震波场数值模拟,但是该方法需要大量的计 算资源才能实现。消耗大量的计算资源主要体现在:(1)对每一频率进行地震波场数值 模拟时,需要在时间域和空间反复调用快速傅里叶变换;(2)对于裂缝与围岩速度横向的快 速变化,需要采用多个速度进行插值获得精度更高的波场;(3)针对裂缝介质,采用网格细分的方法构建的地质模型也显著地增加了计算工作量。为此本申请采用了多机多核的并行框架来进行加速。
基于MPI的相移校正和相移插值结合的波动方程数值模拟的节点间进程级并行计算实 现步骤如下:
A1、调用MPI_INIT_THREAD启动并行环境;
A2、根据计算节点的数量,调用MPI_Comm_size命令确定进程个数(一个节点执行一 个进程),同时调用MPI_Comm_rank命令获得各进程的编号;
A3、其他进程向主进程发送请求需要处理的频率块的范围,0号进程为主进程,接收 到请求后向其他分发需要处理的频率块范围;
A4、各进程执行波动方程正演数值模拟并行计算;波动方程表达式见式(1);
A5、将成像结果文件发送到0号进程,0号进程进行合并;
A6、调用MPI_Finalize,结束并行环境。
根据上述步骤实现的波动方程数值模拟并行计算,在每个节点上只启动一个进程,也 即每个节点上只调用一个CPU核参与计算,因此其可为各个进程提供足够大的内存。当然, 该并行模式下每个节点只有一个CPU核参与运算,也造成了计算资源的浪费,为充分利用 各个节点的计算资源,本发明在实现基于MPI的节点间波动方程数值模拟并行计算的基础 上,在每个节点内部又根据各个节点的CPU核数,基于OpenMP将每个节点的单进程分解为多个线程,每个CPU核执行一个线程,以充分利用计算资源,提高波动方程正演的并 行计算效率。
基于OpenMP的节点内波动方程数值模拟线程级并行计算实现步骤如下:
B1、根据当前结点的CPU核数,调用int_omp_set_num_threads函数设置当前节点要 启动的线程数目;
B2、在频率循环外层调用#pragma omp parallel num_threads()语句,将频率波场数值模 拟设定为并行计算区域;
B3、线程私有变量申请,例如:#pragma omp parallel for firstprivate(v,real,imag,vv,re,im) lastprivate(re,im);为了确保线程安全,因为每个线程中的变量名都是一样的,变量在线程 中发生改变时,容易发生冲突,造成线程不安全。
B4、在频率循环外层调用omp_get thread_num函数,得到当前的线程号,为每个线程 分配起始和终止频率。
根据等效介质理论模拟的地震记录如图2(a)所示。在图2(a)中仅能观察到裂缝发育带边 界处的各种绕射波,裂缝发育带内部无响应;图2(b)为采用本专利方法获得的计算地震记 录;图2(b)中可以清楚地识别各种绕射波,不仅能观察到裂缝带边界处的地震响应特征, 也能展示裂缝带内部的反射特征,因此可在实际的地震剖面上识别这些裂缝发育区。由于 实际地震资料偏移处理时,不能将裂缝与围体的速度加以区分,我们采用等效介质获得速 度对图2(b)的进行偏移,获得图3所示的偏移剖面。从图3中我们能够清晰地识别裂缝产生 的地震响应特征,也从另一侧面证实了提出方法的正确性和有效性。
采用集群实现,本专利针对构建的地质模型,速度体大小4096*4096,开启了6个运行 节点,多线程个数为3,可以看出各个运行平台的时间对比图如图4所示,基于MPI+OpenMP 的时间比单机单核的基本上提高了23倍。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的 原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的 技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任 何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (3)
1.基于多机多核的裂缝介质地震波场数值模拟方法,其特征在于,包括:
S1、根据裂缝的方向、形状、体积大小,找出最小裂缝单元,采用网格细分技术构建地质模型;
S2、采用相移校正和相移插值结合对构建的地质模型进行地震波场数值模拟;
S3、在波动方程频率域采用MPI与OpenMP并行的方式,实现地震波场数值模拟,具体为:在集群的计算节点间采用MPI的粗粒度的多进程并行模式,在计算节点内采用OpenMP的按频率细分的细粒度的多线程并行模式。
2.根据权利要求1所述的基于多机多核的裂缝介质地震波场数值模拟方法,其特征在于,步骤S3所述在集群的计算节点间采用MPI的粗粒度的多进程并行模式具体为:
A1、调用MPI_INIT_THREAD启动并行环境;
A2、根据计算节点的数量,调用MPI_Comm_size命令确定进程个数,同时调用MPI_Comm_rank命令获得各进程的编号;
A3、其他进程向主进程发送请求需要处理的频率块的范围,主进程接收到请求后向其他进程分发需要处理的频率块范围;
A4、各进程执行波动方程正演数值模拟并行计算;
A5、将成像结果文件发送到主进程,由主进程进行合并;
A6、调用MPI_Finalize,结束并行环境。
3.根据权利要求1所述的基于多机多核的裂缝介质地震波场数值模拟方法,其特征在于,步骤S3所述在计算节点内采用OpenMP的按频率细分的细粒度多线程并行模式,具体为:
B1、根据当前结点的CPU核数,调用int_omp_set_num_threads函数设置当前节点要启动的线程数目;
B2、在频率循环外层调用#pragma omp parallel num_threads()语句,将频率波场数值模拟设定为并行计算区域;
B3、线程私有变量申请;
B4、在频率循环外层通过调用omp_get thread_num函数,得到当前的线程号,为每个线程分配起始和终止频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710749074.9A CN107526104A (zh) | 2017-08-28 | 2017-08-28 | 基于多机多核的裂缝介质地震波场数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710749074.9A CN107526104A (zh) | 2017-08-28 | 2017-08-28 | 基于多机多核的裂缝介质地震波场数值模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107526104A true CN107526104A (zh) | 2017-12-29 |
Family
ID=60682718
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710749074.9A Pending CN107526104A (zh) | 2017-08-28 | 2017-08-28 | 基于多机多核的裂缝介质地震波场数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107526104A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110542920A (zh) * | 2019-09-03 | 2019-12-06 | 北京云庐科技有限公司 | 地震数据处理方法及其*** |
CN110927784A (zh) * | 2019-11-29 | 2020-03-27 | 中国地质大学(武汉) | 一种基于OpenMP的频率域高斯束格林函数计算方法 |
CN112632005A (zh) * | 2019-10-08 | 2021-04-09 | 中国石油化工股份有限公司 | 基于mpi的地震数据计算方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104461466A (zh) * | 2013-09-25 | 2015-03-25 | 广州中国科学院软件应用技术研究所 | 基于MPI和OpenMP混合编程模型并行计算提高计算速度的方法 |
CN102239507B (zh) * | 2008-12-03 | 2016-02-24 | 雪佛龙美国公司 | 用于离散裂隙建模的网格生成***和方法 |
-
2017
- 2017-08-28 CN CN201710749074.9A patent/CN107526104A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102239507B (zh) * | 2008-12-03 | 2016-02-24 | 雪佛龙美国公司 | 用于离散裂隙建模的网格生成***和方法 |
CN104461466A (zh) * | 2013-09-25 | 2015-03-25 | 广州中国科学院软件应用技术研究所 | 基于MPI和OpenMP混合编程模型并行计算提高计算速度的方法 |
Non-Patent Citations (3)
Title |
---|
宋鹏 等: "基于MPI+OpenMP的三维声波方程正演模拟", 《中国海洋大学学报》 * |
张选朋 等: "相移及相移内插(PSPI)偏移方法研究", 《城市建设理论研究(电子版)》 * |
张邦: "三维复杂介质弹性波场有限差分数值模拟方法研究", 《中国优秀硕士学位论文全文数据库·基础科学辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110542920A (zh) * | 2019-09-03 | 2019-12-06 | 北京云庐科技有限公司 | 地震数据处理方法及其*** |
CN110542920B (zh) * | 2019-09-03 | 2021-06-22 | 北京云庐科技有限公司 | 地震数据处理方法及其*** |
CN112632005A (zh) * | 2019-10-08 | 2021-04-09 | 中国石油化工股份有限公司 | 基于mpi的地震数据计算方法及*** |
CN112632005B (zh) * | 2019-10-08 | 2024-01-23 | 中国石油化工股份有限公司 | 基于mpi的地震数据计算方法及*** |
CN110927784A (zh) * | 2019-11-29 | 2020-03-27 | 中国地质大学(武汉) | 一种基于OpenMP的频率域高斯束格林函数计算方法 |
CN110927784B (zh) * | 2019-11-29 | 2020-09-29 | 中国地质大学(武汉) | 一种基于OpenMP的频率域高斯束格林函数计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ahmed et al. | Control-volume distributed multi-point flux approximation coupled with a lower-dimensional fracture model | |
US10995592B2 (en) | Method and system for analyzing the uncertainty of subsurface model | |
US10318663B2 (en) | Method of reservoir compartment analysis using topological structure in 3D earth model | |
AU2018277657A1 (en) | Constructing structural models of the subsurface | |
AU2011356658A1 (en) | Method of reservoir compartment analysis using topological structure in 3D earth model | |
Wheeler et al. | A parallel multiblock/multidomain approach for reservoir simulation | |
CN104297785A (zh) | 岩相约束储层物性参数反演方法及装置 | |
Zhang et al. | Determination of RVE with consideration of the spatial effect | |
Xu et al. | Building a three dimensional sealed geological model to use in numerical stress analysis software: a case study for a dam site | |
CN107526104A (zh) | 基于多机多核的裂缝介质地震波场数值模拟方法 | |
Chen et al. | Scalable generation of large-scale unstructured meshes by a novel domain decomposition approach | |
Wang et al. | Accurate modeling and simulation of seepage in 3D heterogeneous fractured porous media with complex structures | |
CN114895351A (zh) | 模拟地震波在任意不连续界面处传播的介质模型化的方法及装置 | |
Qian et al. | GAN inversion method of an initial in situ stress field based on the lateral stress coefficient | |
CN106597542A (zh) | 一种储层特征参数的预测方法及装置 | |
Mohajerani et al. | An efficient algorithm for generation of conforming mesh for three-dimensional discrete fracture networks | |
Liu et al. | Flat voxel-based modelling, assessment and visualization of urban underground space resource quality | |
Zhang et al. | A Low Overhead Heterogeneous Parallel Optimization Method Based on Three-Dimensional Elastic Wave Numerical Simulation | |
Zhang et al. | Conditional stochastic simulation of fluvial reservoirs using multi-scale concurrent generative adversarial networks | |
Gao | Inverse back analysis based on evolutionary neural networks for underground engineering | |
Tran et al. | Characterizing and modelling of fractured reservoirs with object-oriented global optimization | |
Hu et al. | An efficient computational framework for height-contained growing and intersecting hydraulic fracturing simulation via SGBEM–FEM | |
Käser et al. | SeisSol–a software for seismic wave propagation simulations | |
CN110162804A (zh) | 基于cpu加速的波场正演模拟优化方法 | |
Liang et al. | A novel workflow for fracture reconstruction and uncertainty analysis for unconventional reservoir development |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171229 |