CN114442154A - 变热物性地震波传播模拟方法、系统、设备 - Google Patents

变热物性地震波传播模拟方法、系统、设备 Download PDF

Info

Publication number
CN114442154A
CN114442154A CN202210370973.9A CN202210370973A CN114442154A CN 114442154 A CN114442154 A CN 114442154A CN 202210370973 A CN202210370973 A CN 202210370973A CN 114442154 A CN114442154 A CN 114442154A
Authority
CN
China
Prior art keywords
temperature
thermophysical
medium
thermal
seismic
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.)
Granted
Application number
CN202210370973.9A
Other languages
English (en)
Other versions
CN114442154B (zh
Inventor
王志伟
符力耘
王恩江
薛章纳
邓武兵
侯婉婷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202210370973.9A priority Critical patent/CN114442154B/zh
Publication of CN114442154A publication Critical patent/CN114442154A/zh
Application granted granted Critical
Publication of CN114442154B publication Critical patent/CN114442154B/zh
Priority to JP2022159420A priority patent/JP7199132B1/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis

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

本发明属于地球物理勘探、深水/深层油气勘探领域，具体涉及一种变热物性地震波传播模拟方法、系统、设备，旨在解决在高温地层的地下构造成像时，现有技术无法对弹性波在高温介质内的波场快照进行模拟，进而导致地震成像不准确的问题。本方法包括：获取在地震勘探过程中采集的震源参数、变热物性介质参数以及速度模型参数；通过热介质的热导率、比热随温度变化的线性公式计算斜率值；构建变热物性热弹性动力学方程；结合变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照；获取弹性波在变热物性介质中的真实波场快照。本发明完成了弹性波在不同地层温度介质内的波场快照模拟，进而提升地震成像质量。

Description

变热物性地震波传播模拟方法、系统、设备
技术领域
本发明属于地球物理勘探、深水（海洋地质勘察）/深层（陆地地质勘察）油气勘探领域，具体涉及一种变热物性地震波传播模拟方法、系统、设备。
背景技术
我国油气勘探开发逐渐从常规油气资源向深层/深水油气资源发展，深层/深水油气的勘探主要依赖于地球物理技术，而岩石物理研究是油气测井评价和地震预测的基础，虽然我国深层/深水油气勘探取得局部突破，但深层/深水高温高压条件下的岩石物理规律认识并不清楚，高温岩石物理的理论和实验研究才刚刚起步，深层/深水环境下的高温对储层岩石的地球物理响应特征的研究几为空白。对于地下真实的高陡构造，只考虑常温情况，能够利用数值方法模拟弹性波在介质内的波场快照，也能用高精度偏移算法进行成像，恢复其高陡构造形态。但是如果对于一个高陡构造，不同地层的温度不同，这种假设更符合实际地层，弹性波在该介质中的传播过程是怎样，它的波场特征、相速度、走时是怎样；采集的信号是否衰减；利用高精度的偏移算法成像，其界面是否准确；在高温介质中，目前还不清楚，但若想回答以上问题，首先需要研究弹性波在高温介质中的传播规律。基于此，本发明提出了一种变热物性地震波传播模拟方法。
发明内容
为了解决现有技术中的上述问题，即为了解决在高温地层的地下构造成像时，现有技术无法对弹性波在高温介质内的波场快照进行模拟，进而导致地震成像不准确的问题，本发明第一方面，提出了一种变热物性地震波传播模拟方法，该方法包括以下步骤：
S100，获取在地震勘探过程中采集的震源参数、变热物性介质参数以及速度模型参数；所述震源参数包括震源类型、震源主频；所述变热物性介质参数包括变热物性介质的密度、热导率、比热、热膨胀系数、参考温度、拉梅常数；所述速度模型参数包括速度模型的大小、网格间距、网格数量、时间步长以及吸收边界厚度；所述变热物性介质为岩石；
S200，基于所述变热物性介质参数，通过热介质的热导率、比热随温度变化的线性公式，计算热导率参数随温度变化的快慢程度，即斜率值；
S300，利用Kirchhoff变换，将随温度变化的热导率和比热引入均匀各向同性的L-S热弹性动力学方程，构建变热物性热弹性动力学方程；
S400，结合所述变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照；
S500，基于步骤S400获取的波场快照，结合所述斜率值，通过参考温度的温度增量 与预设的变量
Figure 977505DEST_PATH_IMAGE001
的换算关系，得到弹性波在变热物性介质中的真实波场快照。
在一些优选的实施方式中，热介质的热导率、比热随温度变化的线性公式为：
Figure 301170DEST_PATH_IMAGE002
Figure 612066DEST_PATH_IMAGE003
其中，
Figure 413800DEST_PATH_IMAGE004
为参考温度为
Figure 962593DEST_PATH_IMAGE005
的温度增量，
Figure 316213DEST_PATH_IMAGE006
是随温度变化的热导率，
Figure 488307DEST_PATH_IMAGE007
为岩石热导 率，
Figure 483944DEST_PATH_IMAGE008
为常温下的岩石的热导率，
Figure 559348DEST_PATH_IMAGE009
为系数，表示热导率参数随温度变化的快慢程度，
Figure 287132DEST_PATH_IMAGE010
是 比热，
Figure 572620DEST_PATH_IMAGE011
是随温度变化的比热，
Figure 981736DEST_PATH_IMAGE012
为热扩散率，
Figure 36280DEST_PATH_IMAGE013
为密度。
在一些优选的实施方式中，所述均匀各向同性的L-S热弹性动力学方程为：
Figure 872649DEST_PATH_IMAGE014
Figure 848695DEST_PATH_IMAGE015
其中，
Figure 654977DEST_PATH_IMAGE016
为拉普拉斯算子，
Figure 469087DEST_PATH_IMAGE017
为驰豫时间，
Figure 335412DEST_PATH_IMAGE018
为应力系数，
Figure 798754DEST_PATH_IMAGE019
，
Figure 284093DEST_PATH_IMAGE020
为热膨 胀系数，
Figure 578808DEST_PATH_IMAGE021
和
Figure 756980DEST_PATH_IMAGE022
合称为拉梅常数，
Figure 442039DEST_PATH_IMAGE023
、
Figure 855703DEST_PATH_IMAGE024
分别表示
Figure 880291DEST_PATH_IMAGE025
、
Figure 353997DEST_PATH_IMAGE026
对时间求一阶导数，
Figure 526353DEST_PATH_IMAGE027
、
Figure 648767DEST_PATH_IMAGE028
、
Figure 855758DEST_PATH_IMAGE029
分别 表示
Figure 641311DEST_PATH_IMAGE025
、
Figure 97700DEST_PATH_IMAGE030
、
Figure 728533DEST_PATH_IMAGE031
对时间求二阶导数，
Figure 524450DEST_PATH_IMAGE031
和
Figure 339960DEST_PATH_IMAGE030
表示位移，
Figure 424590DEST_PATH_IMAGE032
代表x方向的位移，
Figure 983748DEST_PATH_IMAGE033
代表
Figure 899751DEST_PATH_IMAGE034
方向的位移，右下角标
Figure 17921DEST_PATH_IMAGE035
或者
Figure 714482DEST_PATH_IMAGE036
表示对x方向或者z方向求一阶导数，右下角标
Figure 687117DEST_PATH_IMAGE037
或者
Figure 723206DEST_PATH_IMAGE038
表示对x方向或者z方向求二阶导。
在一些优选的实施方式中，变量
Figure 614939DEST_PATH_IMAGE039
为：
Figure 674162DEST_PATH_IMAGE040
。
在一些优选的实施方式中，所述变热物性的热弹性动力学方程为：
Figure 575122DEST_PATH_IMAGE041
Figure 200138DEST_PATH_IMAGE042
其中，
Figure 403717DEST_PATH_IMAGE043
表示
Figure 809291DEST_PATH_IMAGE044
对时间求一阶导数，
Figure 887843DEST_PATH_IMAGE045
为
Figure 632945DEST_PATH_IMAGE046
对位移求二阶导数，
Figure 132060DEST_PATH_IMAGE047
，分别表示 对x方向或者z方向求二阶导数，
Figure 900296DEST_PATH_IMAGE048
表示对时间
Figure 143058DEST_PATH_IMAGE049
求导数，
Figure 414771DEST_PATH_IMAGE005
为参考温度。
在一些优选的实施方式中，结合所述变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照，其方法为：
将所述变热物性的热弹性动力学方程按照速度-应力速度格式展开，空间导数采用傅里叶变换，时间导数采用中心差分，利用伪谱法模拟弹性波在每一个时间点和每一个位置的振幅值，即弹性波在热物性参数不随温度变化的介质中的波场快照。
在一些优选的实施方式中，所述参考温度的温度增量T与变量ϑ的换算关系为：
Figure 288049DEST_PATH_IMAGE050
。
本发明的第二方面，提出了一种变热物性地震波传播模拟系统，该系统包括：参数获取模块、斜率值计算模块、方程构建模块、伪谱法模拟模块、真实波场快照获取模块；
所述参数获取模块，配置为获取在地震勘探过程中采集的震源参数、变热物性介质参数以及速度模型参数；所述震源参数包括震源类型、震源主频；所述变热物性介质参数包括变热物性介质的密度、热导率、比热、热膨胀系数、参考温度、拉梅常数；所述速度模型参数包括速度模型的大小、网格间距、网格数量、时间步长以及吸收边界厚度；所述变热物性介质为岩石；
所述斜率值计算模块，配置为基于所述变热物性介质参数，通过热介质的热导率、比热随温度变化的线性公式，计算热导率参数随温度变化的快慢程度，即斜率值；
所述方程构建模块，配置为利用Kirchhoff变换，将随温度变化的热导率和比热引入均匀各向同性的L-S热弹性动力学方程，构建变热物性热弹性动力学方程；
所述伪谱法模拟模块，配置为结合所述变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照；
所述真实波场快照获取模块，配置为基于伪谱法模拟模块获取的波场快照，结合 所述斜率值，通过参考温度的温度增量与预设的变量
Figure 402635DEST_PATH_IMAGE001
的换算关系，得到弹性波在变热物 性介质中的真实波场快照。
本发明的第三方面，提出了一种设备，至少一个处理器；以及与至少一个所述处理器通信连接的存储器；其中，所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的变热物性地震波传播模拟方法。
本发明的第四方面，提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的变热物性地震波传播模拟方法。
本发明的有益效果：
本发明完成了弹性波在不同地层温度介质内的波场快照模拟，进而提升地震成像质量。
1）本发明建立耦合温度场和位移场的热弹性动力学方程，研究高温地层介质条件下弹性波的波场特征，实现了在不同地层温度介质内的波场快照模拟；
2）本发明通过精度高的伪谱法模拟弹性波在热物性参数随温度变化的介质中的波场快照，使得模拟的波场快照结果可信度高；
3）本发明考虑了热导率和比热随温度变化的情况，与真实地层中的热物性参数随温度变化的规律相符合，更加真实的反映地下弹性波在地层传播的特征，并提升了地震成像质量，使得本发明可以广泛应用于深层/深水环境下的陆地/海洋地质勘查业相关的地球物理领域。
附图说明
通过阅读参照以下附图所做的对非限制性实施例所做的详细描述，本申请的其他特征、目的和优点将会变得更明显。
图1是本发明一种实施例的变热物性地震波传播模拟方法的流程示意图；
图2是本发明一种实施例的变热物性地震波传播模拟系统的框架示意图；
图3是本发明一种实施例的变热物性地震波传播模拟方法的简略示意图；
图4是本发明一种实施例的各类岩石的热导率随温度变化的曲线的示意图；
图5是本发明一种实施例的未考虑热物性参数变化的热波波场模拟快照的示意图；
图6是本发明一种实施例的考虑热物性参数变化的热波波场快照的示意图；
图7是本发明一种实施例的在未考虑变热物性的道集设定位置提取的曲线特征的示意图；
图8是本发明一种实施例的在考虑变热物性的道集设定位置提取的曲线特征的示意图；
图9是本发明一种实施例的适于用来实现本申请实施例的电子设备的计算机系统的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅用于解释相关发明，而非对该发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与有关发明相关的部分。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
本发明的一种变热物性地震波传播模拟方法，如图1所示，该方法包括以下步骤：
S100，获取在地震勘探过程中采集的震源参数、变热物性介质参数以及速度模型参数；所述震源参数包括震源类型、震源主频；所述变热物性介质参数包括变热物性介质的密度、热导率、比热、热膨胀系数、参考温度、拉梅常数；所述速度模型参数包括速度模型的大小、网格间距、网格数量、时间步长以及吸收边界厚度；所述变热物性介质为岩石；
S200，基于所述变热物性介质参数，通过热介质的热导率、比热随温度变化的线性公式，计算热导率参数随温度变化的快慢程度，即斜率值；
S300，利用Kirchhoff变换，将随温度变化的热导率和比热引入均匀各向同性的L-S热弹性动力学方程，构建变热物性热弹性动力学方程；
S400，结合所述变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照；
S500，基于步骤S400获取的波场快照，结合所述斜率值，通过参考温度的温度增量 与预设的变量
Figure 324455DEST_PATH_IMAGE001
的换算关系，得到弹性波在变热物性介质中的真实波场快照。
为了更清晰地对本发明变热物性地震波传播模拟方法进行说明，下面结合附图3对本发明方法一种实施例中各步骤进行展开详述。
本发明提出一种新的变热物性地震波传播模拟技术，首先建立了热弹性动力学方程，在此基础上，考虑了真实情况下岩石的比热和热导率随温度发生变化的情况，将两种变量（比热和热导率）引入 ，利用伪谱法模拟变热物性介质中地震波传播过程。对比考虑和未考虑热导率和比热性质随温度变化的波场快照，从而说明实际模拟过程中考虑热导率和比热性质随温度变化是非常有必要的。具体过程如下：
S100，获取在地震勘探过程中采集的震源参数、变热物性介质参数以及速度模型参数；所述震源参数包括震源类型、震源主频；所述变热物性介质参数包括变热物性介质的密度、热导率、比热、热膨胀系数、参考温度、拉梅常数；所述速度模型参数包括速度模型的大小、网格间距、网格数量、时间步长以及吸收边界厚度；所述变热物性介质为岩石；
在本实施例中，首先获取震源参数、变热物性介质参数以及速度模型参数。为了更好区分，将震源参数、变热物性介质参数作为方程参数，将速度模型参数作为模型参数。
其中，在本发明中，震源类型、震源主频优选设置为子波为主频3.5 MHz的雷克子波，密度优选设置为2600 kg/m3，热导率优选设置为2 Wm-1K-1、比热优选设置为104 m/(s2°K)，即本发明优选的变热物性介质为石榴斜长角闪岩、热膨胀系数优选设置为4.09×10-6K-1、参考温度优选设置为298K、拉梅常数优选设置为λ = 4.0×109 Pa和μ = 6.0×109 Pa，横向和纵向网格间距优选设置为100 μm，网格数量优选设置为231个网格，时间步长优选设置为100 μm、吸收边界厚度优选设置为20个网格。
S200，基于所述变热物性介质参数，通过热介质的热导率和比热随温度变化的线性公式，计算热导率参数随温度变化的快慢程度，即斜率值；
在本实施例中，通过热介质的热导率、比热随温度变化的线性公式，计算热导率参 数随温度变化的快慢程度，即斜率值
Figure 840887DEST_PATH_IMAGE009
，具体如同时（1）（2）所示：
Figure 885066DEST_PATH_IMAGE002
（1）
Figure 860850DEST_PATH_IMAGE051
（2）
其中，
Figure 445415DEST_PATH_IMAGE004
为参考温度的温度增量，
Figure 957299DEST_PATH_IMAGE006
是随温度变化的热导率，
Figure 969118DEST_PATH_IMAGE007
为岩石热导率，
Figure 995979DEST_PATH_IMAGE008
为常温下的岩石的热导率，
Figure 259602DEST_PATH_IMAGE009
为系数，表示热导率参数随温度变化的快慢程度，
Figure 750626DEST_PATH_IMAGE052
是比 热，
Figure 808712DEST_PATH_IMAGE011
是随温度变化的比热，
Figure 588449DEST_PATH_IMAGE012
为热扩散率，
Figure 514817DEST_PATH_IMAGE013
为密度。
其中，各类岩石（即热介质）的热导率随温度变化的曲线，如图4所示。
S300，利用Kirchhoff变换，将随温度变化的热导率和比热引入均匀各向同性的L-S热弹性动力学方程，构建变热物性热弹性动力学方程；
在本实施例中，均匀各向同性的L-S热弹性动力学方程（即热弹性波动方程和热传导方程）为：
Figure 234249DEST_PATH_IMAGE053
（3）
Figure 853449DEST_PATH_IMAGE054
（4）
其中，
Figure 854903DEST_PATH_IMAGE055
为拉普拉斯算子，
Figure 460328DEST_PATH_IMAGE056
为驰豫时间，
Figure 394786DEST_PATH_IMAGE057
为应力系数，
Figure 325832DEST_PATH_IMAGE058
，
Figure 611320DEST_PATH_IMAGE059
为热膨胀 系数，
Figure 82753DEST_PATH_IMAGE060
和
Figure 12663DEST_PATH_IMAGE061
合称为拉梅常数，
Figure 239245DEST_PATH_IMAGE062
、
Figure 120351DEST_PATH_IMAGE063
分别表示
Figure 395474DEST_PATH_IMAGE064
、
Figure 570104DEST_PATH_IMAGE026
对时间求一阶导数，
Figure 577374DEST_PATH_IMAGE065
、
Figure 837454DEST_PATH_IMAGE066
、
Figure 650689DEST_PATH_IMAGE067
分别表 示
Figure 555191DEST_PATH_IMAGE064
、
Figure 592418DEST_PATH_IMAGE068
、
Figure 480739DEST_PATH_IMAGE069
对时间求二阶导数，
Figure 894403DEST_PATH_IMAGE069
和
Figure 981308DEST_PATH_IMAGE068
表示位移，
Figure 828916DEST_PATH_IMAGE070
，，
Figure 798009DEST_PATH_IMAGE071
代表x方向的位移，
Figure 890730DEST_PATH_IMAGE072
代 表
Figure 894458DEST_PATH_IMAGE073
方向的位移，右下角标
Figure 476749DEST_PATH_IMAGE074
或者
Figure 74083DEST_PATH_IMAGE075
表示对x方向或者z方向求一阶导数，右下角标
Figure 829550DEST_PATH_IMAGE076
或者
Figure 891047DEST_PATH_IMAGE077
表示对x方向或者z方向求二阶导。以上方程满足爱因斯坦求和公式。
利用Kirchhoff变换公式，引入变量
Figure 581922DEST_PATH_IMAGE001
：
Figure 791187DEST_PATH_IMAGE078
（5）
进而导出变热物性的热弹性动力学方程，如公式（6）（7）：
Figure 458666DEST_PATH_IMAGE079
（6）
Figure 374669DEST_PATH_IMAGE080
（7）
其中，
Figure 361080DEST_PATH_IMAGE081
表示
Figure 933007DEST_PATH_IMAGE082
对时间求一阶导数，
Figure 30276DEST_PATH_IMAGE083
为
Figure 800786DEST_PATH_IMAGE001
对位移求二阶导数，
Figure 833464DEST_PATH_IMAGE084
，分别表示 对x方向或者z方向求二阶导数，
Figure 17320DEST_PATH_IMAGE085
表示对时间
Figure 528067DEST_PATH_IMAGE086
求导数，
Figure 215401DEST_PATH_IMAGE087
为参考温度。
S400，结合所述变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照；
在本实施中，具体来说，将上述的变热物性热弹性动力学方程按照速度-应力速度格式展开，空间导数采用傅里叶变换，时间导数利用中心差分，利用伪谱法可以模拟弹性波在每一个时间点和每一个位置的振幅值，即弹性波在介质中的波场快照。
其中，在热物性参数不随温度变化的介质中的波场快照，即未考虑变热物性参数变化的热波波场模拟快照，如图5所示。
S500，基于步骤S400获取的波场快照，结合所述斜率值，通过参考温度的温度增量 与预设的变量
Figure 746876DEST_PATH_IMAGE001
的换算关系，得到弹性波在变热物性介质中的真实波场快照。
在本实施例中，参考温度的温度增量与变量
Figure 260772DEST_PATH_IMAGE001
的换算关系为：
Figure 231002DEST_PATH_IMAGE088
（8）
得到的弹性波在变热物性介质中的真实波场快照，即考虑热物性参数随温度变化的热波波场快照，如图6所示。
基于得到的弹性波在变热物性介质中的真实波场快照，可以进一步进行地震成像，实现高效率、高质量的地质勘探。基于波场快照进行地震成像，属于现有技术，此处不再一一阐述。
另外，为了说明实际模拟过程中考虑热导率和比热性质随温度变化是非常有必要的，本发明将考虑和未考虑热导率和比热性质随温度变化的波场快照进行了对比，即分别提取变热物性和未变热物性的一个道集（即设定位置提取的道集），对比考虑热物性随温度变化前后的曲线特征，对比结果如图7、8所示。
本发明第二实施例的一种变热物性地震波传播模拟系统，如图2所示，具体包括以下模块：参数获取模块100、斜率值计算模块200、方程构建模块300、伪谱法模拟模块400、真实波场快照获取模块500；
所述参数获取模块100，配置为获取在地震勘探过程中采集的震源参数、变热物性介质参数以及速度模型参数；所述震源参数包括震源类型、震源主频；所述变热物性介质参数包括变热物性介质的密度、热导率、比热、热膨胀系数、参考温度、拉梅常数；所述速度模型参数包括速度模型的大小、网格间距、网格数量、时间步长以及吸收边界厚度；所述变热物性介质为岩石；
所述斜率值计算模块200，配置为基于所述变热物性介质参数，通过热介质的热导率、比热随温度变化的线性公式，计算热导率参数随温度变化的快慢程度，即斜率值；
所述方程构建模块300，配置为利用Kirchhoff变换，将随温度变化的热导率和比热引入均匀各向同性的L-S热弹性动力学方程，构建变热物性热弹性动力学方程；
所述伪谱法模拟模块400，配置为结合所述变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照；
所述真实波场快照获取模块500，配置为基于伪谱法模拟模块400获取的波场快 照，结合所述斜率值，通过参考温度的温度增量与预设的变量
Figure 648208DEST_PATH_IMAGE046
的换算关系，得到弹性波在 变热物性介质中的真实波场快照。
所述技术领域的技术人员可以清楚的了解到，为描述的方便和简洁，上述描述的系统具体的工作过程及有关说明，可以参考前述方法实施例中的对应过程，在此不再赘述。
需要说明的是，上述实施例提供的变热物性地震波传播模拟系统，仅以上述各功能模块的划分进行举例说明，在实际应用中，可以根据需要而将上述功能分配由不同的功能模块来完成，即将本发明实施例中的模块或者步骤再分解或者组合，例如，上述实施例的模块可以合并为一个模块，也可以进一步拆分成多个子模块，以完成以上描述的全部或者部分功能。对于本发明实施例中涉及的模块、步骤的名称，仅仅是为了区分各个模块或者步骤，不视为对本发明的不当限定。
本发明第三实施例的一种设备，至少一个处理器；以及与至少一个所述处理器通信连接的存储器；其中，所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的变热物性地震波传播模拟方法。
本发明第四实施例的一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的变热物性地震波传播模拟方法。
所述技术领域的技术人员可以清楚的了解到，未描述的方便和简洁，上述描述的设备、计算机可读存储介质的具体工作过程及有关说明，可以参考前述方法实例中的对应过程，在此不再赘述。
下面参考图9，其示出了适于用来实现本申请方法、系统、装置实施例的服务器的计算机系统的结构示意图。图9示出的服务器仅仅是一个示例，不应对本申请实施例的功能和使用范围带来任何限制。
如图9所示，计算机系统包括中央处理单元(CPU，Central Processing Unit)901，其可以根据存储在只读存储器(ROM，Read Only Memory)902中的程序或者从存储部分908加载到随机访问存储器(RAM，Random Access Memory)903中的程序而执行各种适当的动作和处理。在RAM903 中，还存储有系统操作所需的各种程序和数据。CPU 901、ROM 902以及RAM 903通过总线904彼此相连。输入/输出(I/O，Input/Output)接口905也连接至总线904。
以下部件连接至I/O接口905：包括键盘、鼠标等的输入部分906；包括诸如阴极射线管(CRT，Cathode Ray Tube)、液晶显示器(LCD，Liquid Crystal Display)等以及扬声器等的输出部分907；包括硬盘等的存储部分908；以及包括诸如LAN(局域网，Local AreaNetwork)卡、调制解调器等的网络接口卡的通信部分909。通信部分909经由诸如因特网的网络执行通信处理。驱动器910也根据需要连接至I/O接口905。可拆卸介质911，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器910上，以便于从其上读出的计算机程序根据需要被安装入存储部分908。
特别地，根据本公开的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本公开的实施例包括一种计算机程序产品，其包括承载在计算机可读介质上的计算机程序，该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分909从网络上被下载和安装，和/或从可拆卸介质911被安装。在该计算机程序被中央处理单元(CPU)901执行时，执行本申请的方法中限定的上述功能。需要说明的是，本申请上述的计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意组合。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本申请中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。而在本申请中，计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于：无线、电线、光缆、RF等等，或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本申请的操作的计算机程序代码，上述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++，还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
附图中的流程图和框图，图示了按照本申请各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分，该模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个接连地表示的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
术语“第一”、 “第二”等是用于区别类似的对象，而不是用于描述或表示特定的顺序或先后次序。
术语“包括”或者任何其它类似用语旨在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备/装置不仅包括那些要素，而且还包括没有明确列出的其它要素，或者还包括这些过程、方法、物品或者设备/装置所固有的要素。
至此，已经结合附图所示的优选实施方式描述了本发明的技术方案，但是，本领域技术人员容易理解的是，本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下，本领域技术人员可以对相关技术特征作出等同的更改或替换，这些更改或替换之后的技术方案都将落入本发明的保护范围之内。

Claims (10)

1.一种变热物性地震波传播模拟方法，其特征在于，该方法包括以下步骤：
S100，获取在地震勘探过程中采集的震源参数、变热物性介质参数以及速度模型参数；所述震源参数包括震源类型、震源主频；所述变热物性介质参数包括变热物性介质的密度、热导率、比热、热膨胀系数、参考温度、拉梅常数；所述速度模型参数包括速度模型的大小、网格间距、网格数量、时间步长以及吸收边界厚度；所述变热物性介质为岩石；
S200，基于所述变热物性介质参数，通过热介质的热导率、比热随温度变化的线性公式，计算热导率参数随温度变化的快慢程度，即斜率值；
S300，利用Kirchhoff变换，将随温度变化的热导率和比热引入均匀各向同性的L-S热弹性动力学方程，构建变热物性热弹性动力学方程；
S400，结合所述变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照；
S500，基于步骤S400获取的波场快照，结合所述斜率值，通过参考温度的温度增量与预 设的变量
Figure 364215DEST_PATH_IMAGE001
的换算关系，得到弹性波在变热物性介质中的真实波场快照。
2.根据权利要求1所述的变热物性地震波传播模拟方法，其特征在于，变热物性介质的热导率、比热随温度变化的线性公式为：
Figure 981142DEST_PATH_IMAGE002
Figure 864784DEST_PATH_IMAGE003
其中，
Figure 604070DEST_PATH_IMAGE004
为参考温度的温度增量，
Figure 42004DEST_PATH_IMAGE005
是随温度变化的热导率，
Figure 462621DEST_PATH_IMAGE006
为岩石热导率，
Figure 466350DEST_PATH_IMAGE007
为 常温下的岩石的热导率，
Figure 314220DEST_PATH_IMAGE008
为系数，表示热导率参数随温度变化的快慢程度，
Figure 36188DEST_PATH_IMAGE009
是比热，
Figure 994917DEST_PATH_IMAGE010
是随温度变化的比热，
Figure 321993DEST_PATH_IMAGE011
为热扩散率，
Figure 137502DEST_PATH_IMAGE012
为密度。
3.根据权利要求2 所述的变热物性地震波传播模拟方法，其特征在于，所述均匀各向同性的L-S热弹性动力学方程为：
Figure 284450DEST_PATH_IMAGE013
Figure 109186DEST_PATH_IMAGE014
其中，
Figure 25190DEST_PATH_IMAGE015
为拉普拉斯算子，τ为驰豫时间，
Figure 277180DEST_PATH_IMAGE016
为应力系数，
Figure 645844DEST_PATH_IMAGE017
，
Figure 211955DEST_PATH_IMAGE018
为热膨胀系 数，
Figure 310361DEST_PATH_IMAGE019
和
Figure 670935DEST_PATH_IMAGE020
合称为拉梅常数，
Figure 589212DEST_PATH_IMAGE021
、
Figure 693434DEST_PATH_IMAGE022
分别表示
Figure 646347DEST_PATH_IMAGE004
、
Figure 177822DEST_PATH_IMAGE024
对时间求一阶导数，
Figure 583396DEST_PATH_IMAGE025
、
Figure 491309DEST_PATH_IMAGE026
、
Figure 236411DEST_PATH_IMAGE027
分别表示
Figure 1105DEST_PATH_IMAGE004
、
Figure 831658DEST_PATH_IMAGE028
、
Figure 339999DEST_PATH_IMAGE029
对时间求二阶导数，
Figure 205187DEST_PATH_IMAGE029
和
Figure 78465DEST_PATH_IMAGE028
表示位移，
Figure 458631DEST_PATH_IMAGE030
代表x方向的位移，
Figure 708347DEST_PATH_IMAGE031
代表
Figure 224779DEST_PATH_IMAGE032
方向的位移，右下角标
Figure 268958DEST_PATH_IMAGE033
或者
Figure 339682DEST_PATH_IMAGE034
表示对x方向或者z方向求一阶导数，右下角标
Figure 189827DEST_PATH_IMAGE035
或者
Figure 764028DEST_PATH_IMAGE036
表示对x方向或者z方向求二阶导。
4.根据权利要求3所述的变热物性地震波传播模拟方法，其特征在于，变量
Figure 41425DEST_PATH_IMAGE001
为：
Figure 333866DEST_PATH_IMAGE037
。
5.根据权利要求4所述的变热物性地震波传播模拟方法，其特征在于，所述变热物性的热弹性动力学方程为：
Figure 722122DEST_PATH_IMAGE038
Figure 150829DEST_PATH_IMAGE039
其中，
Figure 864708DEST_PATH_IMAGE040
表示
Figure 378866DEST_PATH_IMAGE041
对时间求一阶导数，
Figure 570812DEST_PATH_IMAGE042
为
Figure 385185DEST_PATH_IMAGE043
对位移求二阶导数，
Figure 942068DEST_PATH_IMAGE044
，分别表示对x方 向或者z方向求二阶导数，
Figure 5839DEST_PATH_IMAGE045
表示对时间
Figure 939160DEST_PATH_IMAGE046
求导数，
Figure 404776DEST_PATH_IMAGE047
为参考温度。
6.根据权利要求5所述的变热物性地震波传播模拟方法，其特征在于，结合所述变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照，其方法为：
将所述变热物性的热弹性动力学方程按照速度-应力速度格式展开，空间导数采用傅里叶变换，时间导数采用中心差分，利用伪谱法模拟弹性波在每一个时间点和每一个位置的振幅值，即弹性波在热物性参数不随温度变化的介质中的波场快照。
7.根据权利要求5所述的变热物性地震波传播模拟方法，其特征在于，所述参考温度的 温度增量T与预设的变量ϑ的换算关系为：
Figure 132561DEST_PATH_IMAGE048
。
8.一种变热物性地震波传播模拟系统，其特征在于，该系统包括：参数获取模块、斜率值计算模块、方程构建模块、伪谱法模拟模块、真实波场快照获取模块；
所述参数获取模块，配置为获取在地震勘探过程中采集的震源参数、变热物性介质参数以及速度模型参数；所述震源参数包括震源类型、震源主频；所述变热物性介质参数包括变热物性介质的密度、热导率、比热、热膨胀系数、参考温度、拉梅常数；所述速度模型参数包括速度模型的大小、网格间距、网格数量、时间步长以及吸收边界厚度；所述变热物性介质为岩石；
所述斜率值计算模块，配置为基于所述变热物性介质参数，通过热介质的热导率、比热随温度变化的线性公式，计算热导率参数随温度变化的快慢程度，即斜率值；
所述方程构建模块，配置为利用Kirchhoff变换，将随温度变化的热导率和比热引入均匀各向同性的L-S热弹性动力学方程，构建变热物性热弹性动力学方程；
所述伪谱法模拟模块，配置为结合所述变热物性热弹性动力学方程，利用伪谱法模拟弹性波在热物性参数不随温度变化的介质中的波场快照；
所述真实波场快照获取模块，配置为基于伪谱法模拟模块获取的波场快照，结合所述 斜率值，通过参考温度的温度增量与预设的变量
Figure 886890DEST_PATH_IMAGE001
的换算关系，得到弹性波在变热物性介 质中的真实波场快照。
9.一种设备，其特征在于，包括：
至少一个处理器；以及
与至少一个所述处理器通信连接的存储器；其中，
所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现权利要求1-7任一项所述的变热物性地震波传播模拟方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现权利要求1-7任一项所述的变热物性地震波传播模拟方法。
CN202210370973.9A 2022-04-11 2022-04-11 变热物性地震波传播模拟方法、系统、设备 Active CN114442154B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210370973.9A CN114442154B (zh) 2022-04-11 2022-04-11 变热物性地震波传播模拟方法、系统、设备
JP2022159420A JP7199132B1 (ja) 2022-04-11 2022-10-03 変化熱物性地震波伝播模擬方法、システム、及び機器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210370973.9A CN114442154B (zh) 2022-04-11 2022-04-11 变热物性地震波传播模拟方法、系统、设备

Publications (2)

Publication Number Publication Date
CN114442154A true CN114442154A (zh) 2022-05-06
CN114442154B CN114442154B (zh) 2022-06-28

Family

ID=81359175

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210370973.9A Active CN114442154B (zh) 2022-04-11 2022-04-11 变热物性地震波传播模拟方法、系统、设备

Country Status (2)

Country Link
JP (1) JP7199132B1 (zh)
CN (1) CN114442154B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008062199A2 (en) * 2006-11-24 2008-05-29 Ucl Business Plc Imaging method and apparatus
US20130189867A1 (en) * 2012-01-25 2013-07-25 INOVA, Ltd. Sealing feature for use with connectors
CN105044771A (zh) * 2015-08-05 2015-11-11 北京多分量地震技术研究院 基于有限差分法的三维tti双相介质地震波场数值模拟方法
CN105759310A (zh) * 2014-12-16 2016-07-13 中国石油化工股份有限公司 复杂非均质储层介质中地震波衰减及速度频散预测方法
CN109490947A (zh) * 2018-10-16 2019-03-19 中国科学院地质与地球物理研究所 一种高温介质地震波传播模拟方法
CN110632677A (zh) * 2019-09-26 2019-12-31 中石化石油工程技术服务有限公司 一种碳酸盐岩地热储层分类方法
CN112329311A (zh) * 2020-11-10 2021-02-05 中国石油大学(华东) 地震波传播有限差分模拟方法、装置及计算机存储介质

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4676664A (en) * 1983-07-15 1987-06-30 The Trustees Of Columbia University In The City Of New York Exploring for subsurface hydrocarbons by sea floor temperature gradients preferably using a multiplexed thermistor probe
JP4332643B2 (ja) 2005-04-05 2009-09-16 独立行政法人産業技術総合研究所 非分極性電極の岩盤或いは地盤への設置方法と、それを用いた電気探査方法若しくは電磁探査方法
BR112015003662B1 (pt) 2012-08-20 2020-11-17 Equinor Energy As MODELAGEM DE TEMPERATURA RESTRITA A DADOS GEOFfSICOS E RESTAURAQAO CINEMATICA
CN111665574B (zh) 2019-03-05 2023-04-07 中国石油化工集团公司 岩石热学参数测井解释方法及系统
CN112485096B (zh) 2020-11-16 2022-01-07 中国石油大学(北京) 用于地热试验的具有非均质岩样的模拟地层的制作方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008062199A2 (en) * 2006-11-24 2008-05-29 Ucl Business Plc Imaging method and apparatus
US20130189867A1 (en) * 2012-01-25 2013-07-25 INOVA, Ltd. Sealing feature for use with connectors
CN105759310A (zh) * 2014-12-16 2016-07-13 中国石油化工股份有限公司 复杂非均质储层介质中地震波衰减及速度频散预测方法
CN105044771A (zh) * 2015-08-05 2015-11-11 北京多分量地震技术研究院 基于有限差分法的三维tti双相介质地震波场数值模拟方法
CN109490947A (zh) * 2018-10-16 2019-03-19 中国科学院地质与地球物理研究所 一种高温介质地震波传播模拟方法
CN110632677A (zh) * 2019-09-26 2019-12-31 中石化石油工程技术服务有限公司 一种碳酸盐岩地热储层分类方法
CN112329311A (zh) * 2020-11-10 2021-02-05 中国石油大学(华东) 地震波传播有限差分模拟方法、装置及计算机存储介质

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
JIA WEI ET AL.: "Green’s function of the Lord–Shulman thermo-poroelasticity theory", 《GEOPHYS. J. INT》 *
JOS E M. CARCIONE ET AL.: "Thermoelasticity and P-wave simulation based on the Cole-Cole model", 《JOURNAL OF THERMAL STRESSES》 *
JOSÉ M. CARCIONE ET AL.: "Physics and Simulation of Wave Propagation in Linear Thermoporoelastic Media", 《JOURNAL OF GEOPHYSICAL RESEARCH: SOLID EARTH》 *
WANTING HOU ET AL.: "Simulation of thermoelastic waves based on the Lord-Shulman theory", 《GEOPHYSICS》 *
侯婉婷等: "热弹性介质中波传播特征", 《地球物理学报》 *
李元燮等: "基于旋转交错网格伪谱法和L-S理论的热弹性介质波场模拟与特征分析", 《中国科学: 地球科学》 *
王志伟等: "岩石热弹性理论及其在地球物理中的应用", 《地球与行星物理论评》 *

Also Published As

Publication number Publication date
JP2023155874A (ja) 2023-10-23
CN114442154B (zh) 2022-06-28
JP7199132B1 (ja) 2023-01-05

Similar Documents

Publication Publication Date Title
Yang et al. Wavefield reconstruction in attenuating media: A checkpointing-assisted reverse-forward simulation method
Ladopoulos Petroleum and gas reserves exploration by real-time expert seismology and non-linear seismic wave motion'
Kulkarni et al. A streamline approach for integrating transient pressure data into high-resolution reservoir models
Zhang et al. Fast-marching methods for complex grids and anisotropic permeabilities: Application to unconventional reservoirs
Dupuy et al. Wave propagation in heterogeneous porous media formulated in the frequency-space domain using a discontinuous Galerkin method
CN105308479A (zh) 通过与偏移距相关的弹性fwi的多参数反演
WO2012170201A2 (en) System and method for seismic data inversion by non-linear model update
Withers et al. Memory-efficient simulation of frequency-dependent Q
Zhang et al. From streamlines to fast marching: rapid simulation and performance assessment of shale gas reservoirs using diffusive time of flight as a spatial coordinate
CN114428324A (zh) 叠前高角度快速傅里叶变换地震成像方法、系统、设备
CN109490947B (zh) 一种高温介质地震波传播模拟方法
CN114442154B (zh) 变热物性地震波传播模拟方法、系统、设备
KR101459388B1 (ko) 지하 속도정보 도출 방법
CN111239809B (zh) 一种基于属性融合原理的优势储层确定方法及系统
CN116738794A (zh) 孔裂隙介质的岩石物理数值模拟方法、装置、设备及介质
KR101695080B1 (ko) 단일채널 탄성파 탐사자료로부터 지하 속도정보 도출을 위한 탄성파 파형역산 방법 및 장치
CN113311484B (zh) 利用全波形反演获取粘弹性介质弹性参数的方法及装置
Fontara et al. Numerical simulation of seismic wave field in graded geological media containing multiple cavities
Maghoul et al. Wave propagation in unsaturated poroelastic media: Boundary integral formulation and three-dimensional fundamental solution
Ma et al. A model expansion criterion for treating surface topography in ray path calculations using the eikonal equation
CN113051697B (zh) 古气藏地层孔隙压力演化史曲线生成方法及装置
Slimani et al. Identification of the source function for a seawater intrusion problem in unconfined aquifer
da Silva et al. Target-oriented inversion using the patched Green’s function method
Sun et al. An efficient 3D cell-based discrete fracture-matrix flow model for digitally captured fracture networks
Li et al. Rapid Inference of Reservoir Permeability From Inversion of Travel Time Data Under a Fast Marching Method Based Deep Learning Framework

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
GR01 Patent grant
GR01 Patent grant