一种中深层干热岩体的自适应变网格正演模拟方法技术

本说明书实施例公开了一种中深层干热岩体的自适应变网格正演模拟方法。通过输入初始地下介质纵波速度V

【技术实现步骤摘要】
一种中深层干热岩体的自适应变网格正演模拟方法


[0001]本说明书涉及地热资源勘探
，尤其涉及一种中深层干热岩体的自适应变网格正演模拟方法。

技术介绍

[0002]传统的常规网格正演模拟方法在横、纵向网格间距设为特定值，若网格间距选取过大，则低速区域易出现数值频散，波场模拟精度降低；若网格间距过小，则高速区域易过采样，造成计算资源浪费，计算成本大。此外，传统的变网格正演模拟方法通常在目标区域采用细网格，其他区域采用粗网格，当模拟参数选取不合理时，在粗、细网格的过渡区域会产生强烈的虚假反射，从而降低模拟波场的精度。
[0003]基于此，需要一种计算速度更快的中深层干热岩体的自适应变网格正演模拟方案。

技术实现思路

[0004]本说明书实施例提供一种中深层干热岩体的自适应变网格正演模拟方法，用以解决如下技术问题：需要一种计算速度更快的中深层干热岩体的自适应变网格正演模拟方案。
[0005]为解决上述技术问题，本说明书一个或多个实施例是这样实现的：
[0006]在第一方面，本说明书实施例提供一种中深层干热岩体的自适应变网格正演模拟方法，包括：
[0007]输入初始地下介质纵波速度V
p
、横波速度V
s
和介质密度ρ模型，预先设置正演模拟参数，所述参数包括：横向网格点数nx，纵向网格点数nz，吸收边界大小d，横向网格间距dx，纵向网格间距dz，时间切片数nt，时间采样率dt，子波主频f
m
，初始炮点位置(S
x<br/>，S
z
),炮间距ds，总炮数ns；
[0008]对所述速度模型进行自适应变网格处理，确定坐标系之间的映射关系，基于所述映射关系确定在变网格下的波动方程：
[0009][0010]其中，ν
x
表示总波场的横向速度、ν
xp
表示纵波的横向速度、ν
xs
表示横波的横向速度，ν
z
表示总波场的纵向速度、ν
zp
表示纵波的纵向速度、ν
zs
表示横波的纵向速度，v
p
表示纵波速度、v
s
表示横波速度，ρ表示介质密度大小；τ
xx
表示介质受到的横向正应力、τ
zz
表示介质受到的纵向正应力和τ
xz
表示介质受到的剪切应力；是自适应变网格坐标系和常规网格坐标系的映射关系，为选取的自适应变网格的垂向网格间距的拟合函数；
[0011]根据所述变网格后的地下介质纵波速度、横波速度和密度，采用所述动方程进行正演模拟，获得在变网格后的特定时刻的各方向波场信息，根据所述变网格后的特定时刻的各方向波场信息，逆向使用映射关系生成最终的规则网格水平方向纵波波场、水平方向横波波场、垂直方向纵波波场和垂直方向横波波场信息。
[0012]本说明书一个或多个实施例采用的上述至少一个技术方案能够达到以下有益效果：通过输入初始地下介质纵波速度V
p
、横波速度V
s
和介质密度ρ模型，预先设置正演模拟参数；对所述速度模型进行自适应变网格处理，确定坐标系之间的映射关系，基于所述映射关系确定在变网格下的波动方程；对所述速度模型进行自适应变网格处理，确定坐标系之间的映射关系，获取变网格后的地下介质纵波速度、横波速度和密度，基于所述映射关系确定在变网格下的波动方程，；根据所述变网格后的特定时刻的各方向波场信息，逆向使用映射关系生成最终的规则网格水平方向纵波波场、水平方向横波波场、垂直方向纵波波场和垂直方向横波波场信息。从而实现在不同尺寸网格间无需设置过渡区域，改变差分格式，或进行插值处理，即可有效避免虚假反射问题，且实现大幅度地降低计算量和计算机内存占用
量，计算速度快。为中深层干热岩体区域的地震数据处理提供了高效率成像保障。
附图说明
[0013]为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0014]图1为本说明书实施例所提供的流程示意图；
[0015]图2为本说明书实施例所提供的自适应变网格重采样方法图；
[0016]图3为本说明书实施例所提供的初始模型纵波速度场图；
[0017]图4为本说明书实施例所提供的初始模型横波速度场图；
[0018]图5为本说明书实施例所提供的初始模型介质密度分布图；
[0019]图6为本说明书实施例所提供的经自适应变网格重采样后的纵波速度场图；
[0020]图7为本说明书实施例所提供的经自适应变网格重采样后的横波速度场图；
[0021]图8为本说明书实施例所提供的经自适应变网格重采样后的介质密度分布图；
[0022]图9为本说明书实施例所提供的本专利技术方案得到的水平方向纵波波场信息图；
[0023]图10为本说明书实施例所提供的本专利技术方案得到的水平方向横波波场信息图；
[0024]图11为本说明书实施例所提供的本专利技术方案得到的垂直方向纵波波场信息图；
[0025]图12为本说明书实施例所提供的本专利技术方案得到的垂直方向横波波场信息图；
[0026]图13为本说明书实施例所提供的常规网格和自适应变网格弹性波正演模拟耗时对比图。
具体实施方式
[0027]本说明书实施例提供一种中深层干热岩体的自适应变网格正演模拟方法。
[0028]为了使本
的人员更好地理解本说明书中的技术方案，下面将结合本说明书实施例中的附图，对本说明书实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本说明书实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
[0029]如图1所示，图1为本说明书实施例提供的一种中深层干热岩体的自适应变网格正演模拟方法的流程示意图。
[0030]图1中的流程可以包括以下步骤：
[0031]S1，输入初始地下介质纵波速度V
p
、横波速度V
s
和介质密度ρ模型，预先设置正演模拟参数。
[0032]所述参数包括：横向网格点数nx，纵向网格点数nz，吸收边界大小d，横向网格间距dx，纵向网格间距dz，时间切片数nt，时间采样率dt，子波主频f
m
，初始炮点位置(S
x
，S
z
),炮间距ds，总炮数ns等。
[0033]S2，对所述速度模型进行自适应变网格处理，确定坐标系之间的映射关系，获取变网格后的地下介质纵波速度、横波速度和密度，基于所述映射关系确定在变网格下的波动本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种中深层干热岩体的自适应变网格正演模拟方法，包括：输入初始地下介质纵波速度V
p
、横波速度V
s
和介质密度ρ模型，预先设置正演模拟参数，所述参数包括：横向网格点数nx，纵向网格点数nz，吸收边界大小d，横向网格间距dx，纵向网格间距dz，时间切片数nt，时间采样率dt，子波主频f
m
，初始炮点位置(S
x
，S
z
),炮间距ds，总炮数ns；对所述速度模型进行自适应变网格处理，确定坐标系之间的映射关系，获取变网格后的地下介质纵波速度、横波速度和密度，基于所述映射关系确定在变网格下的波动方程：其中，ν
x
表示总波场的横向速度、ν
xp
表示纵波的横向速度、ν
xs
表示横波的横向速度，ν
z
表示总波场的纵向速度、ν
zp
表示纵波的纵向速度、ν
zs
表示横波的纵向速度，v
p
表示纵波速度、v
s
表示横波速度，ρ表示介质密度大小；τ
xx
表示介质受到的横向正...

【专利技术属性】
技术研发人员：杨秀金，黄建平，杨继东，张东林，王伟奇，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：