CN109143340A - 一种基于常q模型的粘弹介质地震波模拟方法及系统 - Google Patents
一种基于常q模型的粘弹介质地震波模拟方法及系统 Download PDFInfo
- Publication number
- CN109143340A CN109143340A CN201810945943.XA CN201810945943A CN109143340A CN 109143340 A CN109143340 A CN 109143340A CN 201810945943 A CN201810945943 A CN 201810945943A CN 109143340 A CN109143340 A CN 109143340A
- Authority
- CN
- China
- Prior art keywords
- wave
- hermite
- speed
- underground medium
- discretization
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000002245 particle Substances 0.000 claims description 18
- 230000009466 transformation Effects 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 6
- 241000209094 Oryza Species 0.000 claims description 5
- 235000007164 Oryza sativa Nutrition 0.000 claims description 5
- 235000009566 rice Nutrition 0.000 claims description 5
- 238000004088 simulation Methods 0.000 abstract description 11
- 230000000694 effects Effects 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 235000003332 Ilex aquifolium Nutrition 0.000 description 1
- 241000209027 Ilex aquifolium Species 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000750 progressive effect Effects 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
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
本发明公开了一种基于常Q模型的粘弹介质地震波模拟方法及系统。该方法包括:确定地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度;选取地震子波;确定离散化的埃尔米特分布近似函数;确定所述埃尔米特分布近似函数的导数求解模型;根据地下介质的密度、品质因子、横波在地下介质中的传播速度、纵波在地下介质中的传播速度、地震子波以及离散化的埃尔米特分布近似函数的导数求解模型,确定地震波的粘弹速度‑应力波动参数;根据所述地震波的粘弹速度‑应力波动参数,模拟所述地震波。本发明提供的基于常Q模型的粘弹介质地震波模拟方法及系统具有灵活性强、模拟精度高的特点。
Description
技术领域
本发明涉及石油勘探领域,特别是涉及一种基于常Q模型的粘弹介质地震波模拟方法及系统。
背景技术
地震波在实际地下介质中的传播具有粘滞性的特性。介质对地震波的衰减造成能量衰减和相位畸变。为了更准确描述波在实际介质中的传播规律,需要在波场模拟中考虑介质的吸收衰减作用。Kjartansson从Strick-Azimi Q模型出发,提出了一种由相速度和品质因子Q表示的常Q模型。Carcione et al.和Carcione基于该模型推导了声波和弹性波波动方程,采用Grünwald-Letnikov近似法求解。Chen and Holm提出使用分数阶拉普拉斯算子来求解均匀介质中的波动方程。Carcione给出了均匀介质中使用分数阶拉普拉斯算子的常Q粘声波动方程,在空间-频率域中使用广义傅里叶方法求解。Zhu et al.推导了由二阶时间导数和两个拉普拉斯算子表示的解耦的常Q波动方程并推广到非均匀介质中。但是,现有方法中求解分数阶拉普拉斯算子的傅里叶伪谱法使用全局傅里叶算子,在处理Q值变化剧烈的复杂地质模型时,会面临导致波场畸变的吉普斯效应,进而使计算精度降低。
发明内容
本发明的目的是提供基于常Q模型的粘弹介质地震波模拟方法及系统,具有灵活性强、模拟精度高的特点。
为实现上述目的,本发明提供了如下方案:
一种基于常Q模型的粘弹介质地震波模拟方法,所述方法包括:
确定地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度;
选取地震子波;
确定离散化的埃尔米特分布近似函数;
确定所述埃尔米特分布近似函数的导数求解模型;
根据地下介质的密度、品质因子、横波在地下介质中的传播速度、纵波在地下介质中的传播速度、地震子波以及离散化的埃尔米特分布近似函数的导数求解模型,确定地震波的粘弹速度-应力波动参数;
根据所述地震波的粘弹速度-应力波动参数,模拟所述地震波。
可选的,所述确定离散化的埃尔米特分布近似函数,具体包括:
确定所述离散化的埃尔米特分布近似函数中的参数:高斯窗函数的方差、埃尔米特多项式阶数以及埃尔米特分布近似函数网格离散化算子长度;
根据所述参数确定离散化的埃尔米特分布近似函数其中,δM为埃尔米特分布近似泛函对狄拉克函数的近似离散网格化算子,x为网格中心点,x'为网格的空间坐标,x-x'为网格上每一点到中心点的空间距离,满足|x-x'|/Δx<(W-1)/2,σ为高斯窗函数的方差,M为埃尔米特多项式阶数,W为及埃尔米特分布近似函数网格离散化算子长度,Hi(x)为i阶埃尔米特多项式,
可选的,所述确定所述埃尔米特分布近似函数的导数求解模型,具体包括:
确定所述埃尔米特分布近似函数的整数阶导数求解模型为其中,α为导数阶数;
确定所述埃尔米特分布近似函数的分数阶导数求解模型为其中,F和F-1分别表示傅里叶变换和傅里叶反变换,为虚数单位,k为波数。
可选的,所述根据地下介质的密度、品质因子、横波在地下介质中的传播速度、纵波在地下介质中的传播速度、地震子波以及离散化的埃尔米特分布近似函数的导数求解模型,确定地震波的粘弹速度-应力波动参数,具体包括:
获取地震波的粘弹速度-应力波动方程;
对所述粘弹速度-应力波动方程中的速度、应变分量的时间偏导数做差分离散,得到粘弹速度-应力波动参数的离散形式
其中,分别为n时刻对应的x方向、z方向的地下介质质点的振动速度,分别为n时刻对应的x方向、xz方向、z方向的地下介质质点的应力,分别为n时刻对应的x方向、xz方向、z方向地下介质质点的应变,fx、fz分别为n时刻对应的x方向、z方向所述地震子波的体力分量,ρ为地下介质的密度,Δt为时间采样间隔, 为中间变量,ω0为设定参考频率,γp、γs为空变分数阶阶数,Qp、Qs分别为纵波品质因子和横波品质因子,为体积弹性模量,为剪切弹性模量,和分别为所述设定频率下对应的纵波和横波的速度,(-▽2)γ为γ阶拉普拉斯算子,Dx、Dz为空间偏导数算子,DxU表示对波场U应用空间偏导数算子,,DzU表示对波场U应用空间偏导数算子,Δx、Δz为空间网格间距,为差分系数,N为差分阶数,Δx,Δz分别为空间网格间距,W1、W2分别为x方向、z方向上的网格点数,u为应用所述拉普拉斯算子的波场,xi,zj为应用所述拉普拉斯算子所需的空间褶积对应的波场空间坐标,xm,zn为应用所述拉普拉斯算子后波场空间坐标,σ1和σ2分别为x方向、z方向上的高斯窗函数的方差。
本发明还提供了一种基于常Q模型的粘弹介质地震波模拟系统,所述系统包括:
参数确定模块,用于确定地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度;
地震子波选取模块,用于选取地震子波;
离散化函数确定模块,用于确定离散化的埃尔米特分布近似函数;
导数求解模型确定模块,用于确定所述埃尔米特分布近似函数的导数求解模型;
地震波波动参数确定模块,用于根据地下介质的密度、品质因子、横波在地下介质中的传播速度、纵波在地下介质中的传播速度、地震子波以及离散化的埃尔米特分布近似函数的导数求解模型,确定地震波的粘弹速度-应力波动参数;
地震波模拟模块,用于根据所述地震波的粘弹速度-应力波动参数,模拟所述地震波。
可选的,所述离散化函数确定模块,具体包括:
参数确定单元,用于确定所述离散化的埃尔米特分布近似函数中的参数:高斯窗函数的方差、埃尔米特多项式阶数以及埃尔米特分布近似函数网格离散化算子长度;
离散化函数确定单元,用于根据所述参数确定离散化的埃尔米特分布近似函数其中,δM为埃尔米特分布近似泛函对狄拉克函数的近似离散网格化算子,x为网格中心点,x'为网格的空间坐标,x-x'为网格上每一点到中心点的空间距离,满足|x-x'|/Δx<(W-1)/2,σ为高斯窗函数的方差,M为埃尔米特多项式阶数,W为及埃尔米特分布近似函数网格离散化算子长度,Hi(x)为i阶埃尔米特多项式,
可选的,所述导数求解模型确定模块,具体包括:
整数阶导数求解模型确定单元,用于确定所述埃尔米特分布近似函数的整数阶导数求解模型为其中,α为导数阶数;
分数阶导数求解模型确定单元,用于确定所述埃尔米特分布近似函数的分数阶导数求解模型为其中,F和F-1分别表示傅里叶变换和傅里叶反变换,为虚数单位,k为波数。
可选的,所述地震波波动参数确定模块,具体包括:
波动方程获取单元,用于获取地震波的粘弹速度-应力波动方程;
波动参数离散化单元,用于对所述粘弹速度-应力波动方程中的速度、应变分量的时间偏导数做差分离散,得到粘弹速度-应力波动参数的离散形式
其中,分别为n时刻对应的x方向、z方向的地下介质质点的振动速度,分别为n时刻对应的x方向、xz方向、z方向的地下介质质点的应力,分别为n时刻对应的x方向、xz方向、z方向地下介质质点的应变,fx、fz分别为n时刻对应的x方向、z方向所述地震子波的体力分量,ρ为地下介质的密度,Δt为时间采样间隔, 为中间变量,ω0为设定参考频率,γp、γs为空变分数阶阶数,Qp、Qs分别为纵波品质因子和横波品质因子,为体积弹性模量,为剪切弹性模量,和分别为所述设定频率下对应的纵波和横波的速度,(-▽2)γ为γ阶拉普拉斯算子,Dx、Dz为空间偏导数算子,DxU表示对波场U应用空间偏导数算子,DzU表示对波场U应用空间偏导数算子,Δx、Δz为空间网格间距,为差分系数,N为差分阶数,Δx,Δz分别为空间网格间距,W1、W2分别为x方向、z方向上的网格点数,u为应用所述拉普拉斯算子的波场,xi,zj为应用所述拉普拉斯算子所需的空间褶积对应的波场空间坐标,xm,zn为应用所述拉普拉斯算子后波场空间坐标,σ1和σ2分别为x方向、z方向上的高斯窗函数的方差。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供的基于常Q模型的粘弹介质地震波模拟方法及系统根据工区地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度、选取的地震子波以及通过埃尔米特分布近似泛函方法计算常Q模型粘弹性波方程中出现的空变分数阶拉普拉斯算子,在时间方向上迭代,对地震波场进行模拟,避免了传统方法中对波场进行的全局傅里叶变换和反变换,实现了对真实地下介质中的吸收衰减和频散效应的模拟,具有灵活性强、模拟精度高的特点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例基于常Q模型的粘弹介质地震波模拟方法流程图;
图2为本发明实施例基于常Q模型的粘弹介质地震波模拟系统结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于基于常Q模型的粘弹介质地震波模拟方法及系统,具有灵活性强、模拟精度高的特点。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例基于常Q模型的粘弹介质地震波模拟方法流程图,如图1所示,本发明提供的基于常Q模型的粘弹介质地震波模拟方法流程图包括以下步骤:
步骤101:确定地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度;
步骤102:选取地震子波;
步骤103:确定离散化的埃尔米特分布近似函数;
步骤104:确定所述埃尔米特分布近似函数的导数求解模型;
步骤105:根据地下介质的密度、品质因子、横波在地下介质中的传播速度、纵波在地下介质中的传播速度、地震子波以及离散化的埃尔米特分布近似函数的导数求解模型,确定地震波的粘弹速度-应力波动参数;
步骤106:根据所述地震波的粘弹速度-应力波动参数,模拟所述地震波。
其中,步骤103中使用如下公式求取离散化的埃尔米特分布近似函数:
所述方程(1)中,δM为埃尔米特分布近似泛函对狄拉克函数的近似离散网格化算子,该算子为一维离散网格,x为网格中心点,x'为网格的空间坐标,x-x'为网格上每一点到中心点的空间距离,满足|x-x'|/Δx<(W-1)/2。其中,Δx为空间间隔,W为埃尔米特分布近似函数网格离散化算子长度,σ和M为狄拉克函数参数,σ为高斯窗函数的方差,M为埃尔米特多项式阶数,当σ不变,M趋于无穷大或M不变,σ趋于0时,δM趋近于精确的狄拉克函数。参数W和M的取值越大,算子精度越高,计算量越大,实际取值时应结合允的误差范围和计算效率进行参数选取。Hi(x)为i阶埃尔米特多项式:
所述方程(1)中埃尔米特多项式n阶项通过如下公式递推求取:
步骤104中包括:所述埃尔米特分布近似函数分数阶导数通过如下公式求取:
所述方程(4)中α为导数阶数,当α为整数时,使用上述公式求取;当α为非整数时,使用如下公式:
方程(5)中F和F-1分别表示傅里叶变换和傅里叶反变换,为虚数单位,k为波数。
步骤105中包括:
获取二维黏弹速度-应力波动方程:
对方程(6)中速度、应变分量的时间偏导数做差分离散,得到二维粘弹速度-应力方程递推公式:
方程(6)、(7)中,分别为n时刻对应的x方向、z方向的地下介质质点的振动速度,分别为n时刻对应的x方向、xz方向、z方向的地下介质质点的应力,分别为n时刻对应的x方向、xz方向、z方向地下介质质点的应变,fx、fz分别为n时刻对应的x方向、z方向所述地震子波的体力分量,ρ为地下介质的密度,Δt为时间采样间隔,为中间变量,ω0为设定参考频率,γp、γs为空变分数阶阶数,Qp、Qs分别为纵波品质因子和横波品质因子,为体积弹性模量,为剪切弹性模量,和分别为所述设定频率下对应的纵波和横波的速度,(-▽2)γ为γ阶拉普拉斯算子,Dx、Dz为空间偏导数算子。
DxU、DzU表示对波场U应用空间偏导数算子,Δx、Δz为空间网格间距,为差分系数,N为差分阶数。
分数阶拉普拉斯算子通过如下公式计算:
Δx,Δz分别为空间网格间距,W1、W2分别为x方向、z方向上的网格点数,u为应用所述拉普拉斯算子的波场,xi,zj为应用所述拉普拉斯算子所需的空间褶积对应的波场空间坐标,xm,zn为应用所述拉普拉斯算子后波场空间坐标,σ1和σ2分别为x方向、z方向上的高斯窗函数的方差。
本发明提供的基于常Q模型的粘弹介质地震波模拟方法根据工区地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度、选取的地震子波以及通过埃尔米特分布近似泛函方法计算常Q模型粘弹性波方程中出现的空变分数阶拉普拉斯算子,在时间方向上迭代,对地震波场进行模拟,避免了传统方法中对波场进行的全局傅里叶变换和反变换,实现了对真实地下介质中的吸收衰减和频散效应的模拟,具有灵活性强、模拟精度高的特点。
本发明还提供了一种基于常Q模型的粘弹介质地震波模拟系统,图2为本发明实施例基于常Q模型的粘弹介质地震波模拟系统结构图,如图2所示,本发明提供的地震波模拟系统包括:
参数确定模块201,用于确定地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度;
地震子波选取模块202,用于选取地震子波;
离散化函数确定模块203,用于确定离散化的埃尔米特分布近似函数;
导数求解模型确定模块204,用于确定所述埃尔米特分布近似函数的导数求解模型;
地震波波动参数确定模块205,用于根据地下介质的密度、品质因子、横波在地下介质中的传播速度、纵波在地下介质中的传播速度、地震子波以及离散化的埃尔米特分布近似函数的导数求解模型,确定地震波的粘弹速度-应力波动参数;
地震波模拟模块206,用于根据所述地震波的粘弹速度-应力波动参数,模拟所述地震波。
其中,离散化函数确定模块203,具体包括:
参数确定单元,用于确定所述离散化的埃尔米特分布近似函数中的参数:高斯窗函数的方差、埃尔米特多项式阶数以及埃尔米特分布近似函数网格离散化算子长度;
离散化函数确定单元,用于根据所述参数确定离散化的埃尔米特分布近似函数其中,δM为埃尔米特分布近似泛函对狄拉克函数的近似离散网格化算子,x为网格中心点,x'为网格的空间坐标,x-x'为网格上每一点到中心点的空间距离,满足|x-x'|/Δx<(W-1)/2,σ为高斯窗函数的方差,M为埃尔米特多项式阶数,W为及埃尔米特分布近似函数网格离散化算子长度,Hi(x)为i阶埃尔米特多项式,
导数求解模型确定模块204,具体包括:
整数阶导数求解模型确定单元,用于确定所述埃尔米特分布近似函数的整数阶导数求解模型为其中,α为导数阶数;
分数阶导数求解模型确定单元,用于确定所述埃尔米特分布近似函数的分数阶导数求解模型为其中,F和F-1分别表示傅里叶变换和傅里叶反变换,为虚数单位,k为波数。
地震波波动参数确定模块205,具体包括:
波动方程获取单元,用于获取地震波的粘弹速度-应力波动方程;
波动参数离散化单元,用于对所述粘弹速度-应力波动方程中的速度、应变分量的时间偏导数做差分离散,得到粘弹速度-应力波动参数的离散形式
其中,分别为n时刻对应的x方向、z方向的地下介质质点的振动速度,分别为n时刻对应的x方向、xz方向、z方向的地下介质质点的应力,分别为n时刻对应的x方向、xz方向、z方向地下介质质点的应变,fx、fz分别为n时刻对应的x方向、z方向所述地震子波的体力分量,ρ为地下介质的密度,Δt为时间采样间隔, 为中间变量,ω0为设定参考频率,γp、γs为空变分数阶阶数,Qp、Qs分别为纵波品质因子和横波品质因子,为体积弹性模量,为剪切弹性模量,和分别为所述设定频率下对应的纵波和横波的速度,(-▽2)γ为γ阶拉普拉斯算子,Dx、Dz为空间偏导数算子,DxU表示对波场U应用空间偏导数算子,DzU表示对波场U应用空间偏导数算子,Δx、Δz为空间网格间距,为差分系数,N为差分阶数,Δx,Δz分别为空间网格间距,W1、W2分别为x方向、z方向上的网格点数,u为应用所述拉普拉斯算子的波场,xi,zj为应用所述拉普拉斯算子所需的空间褶积对应的波场空间坐标,xm,zn为应用所述拉普拉斯算子后波场空间坐标,σ1和σ2分别为x方向、z方向上的高斯窗函数的方差。
本发明提供的基于常Q模型的粘弹介质地震波模拟系统根据工区地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度、选取的地震子波以及通过埃尔米特分布近似泛函方法计算常Q模型粘弹性波方程中出现的空变分数阶拉普拉斯算子,在时间方向上迭代,对地震波场进行模拟,避免了传统方法中对波场进行的全局傅里叶变换和反变换,实现了对真实地下介质中的吸收衰减和频散效应的模拟,具有灵活性强、模拟精度高的特点。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种基于常Q模型的粘弹介质地震波模拟方法,其特征在于,所述方法包括:
确定地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度;
选取地震子波;
确定离散化的埃尔米特分布近似函数;
确定所述埃尔米特分布近似函数的导数求解模型;
根据地下介质的密度、品质因子、横波在地下介质中的传播速度、纵波在地下介质中的传播速度、地震子波以及离散化的埃尔米特分布近似函数的导数求解模型,确定地震波的粘弹速度-应力波动参数;
根据所述地震波的粘弹速度-应力波动参数,模拟所述地震波。
2.根据权利要求1所述的基于常Q模型的粘弹介质地震波模拟方法,其特征在于,所述确定离散化的埃尔米特分布近似函数,具体包括:
确定所述离散化的埃尔米特分布近似函数中的参数:高斯窗函数的方差、埃尔米特多项式阶数以及埃尔米特分布近似函数网格离散化算子长度;
根据所述参数确定离散化的埃尔米特分布近似函数其中,δM为埃尔米特分布近似泛函对狄拉克函数的近似离散网格化算子,x为网格中心点,x'为网格的空间坐标,x-x'为网格上每一点到中心点的空间距离,满足|x-x'|/Δx<(W-1)/2,σ为高斯窗函数的方差,M为埃尔米特多项式阶数,W为及埃尔米特分布近似函数网格离散化算子长度,Hi(x)为i阶埃尔米特多项式,
3.根据权利要求1所述的基于常Q模型的粘弹介质地震波模拟方法,其特征在于,所述确定所述埃尔米特分布近似函数的导数求解模型,具体包括:
确定所述埃尔米特分布近似函数的整数阶导数求解模型为其中,α为导数阶数;
确定所述埃尔米特分布近似函数的分数阶导数求解模型为其中,F和F-1分别表示傅里叶变换和傅里叶反变换,为虚数单位,k为波数。
4.根据权利要求1所述的基于常Q模型的粘弹介质地震波模拟方法,其特征在于,所述根据地下介质的密度、品质因子、横波在地下介质中的传播速度、纵波在地下介质中的传播速度、地震子波以及离散化的埃尔米特分布近似函数的导数求解模型,确定地震波的粘弹速度-应力波动参数,具体包括:
获取地震波的粘弹速度-应力波动方程;
对所述粘弹速度-应力波动方程中的速度、应变分量的时间偏导数做差分离散,得到粘弹速度-应力波动参数的离散形式
其中,分别为n时刻对应的x方向、z方向的地下介质质点的振动速度,分别为n时刻对应的x方向、xz方向、z方向的地下介质质点的应力,分别为n时刻对应的x方向、xz方向、z方向地下介质质点的应变,fx、fz分别为n时刻对应的x方向、z方向所述地震子波的体力分量,ρ为地下介质的密度,Δt为时间采样间隔, 为中间变量,ω0为设定参考频率,γp、γs为空变分数阶阶数,Qp、Qs分别为纵波品质因子和横波品质因子,为体积弹性模量,为剪切弹性模量,和分别为所述设定频率下对应的纵波和横波的速度,(-▽2)γ为γ阶拉普拉斯算子,Dx、Dz为空间偏导数算子,DxU表示对波场U应用空间偏导数算子,DzU表示对波场U应用空间偏导数算子,Δx、Δz为空间网格间距,为差分系数,N为差分阶数,Δx,Δz分别为空间网格间距,W1、W2分别为x方向、z方向上的网格点数,u为应用所述拉普拉斯算子的波场,xi,zj为应用所述拉普拉斯算子所需的空间褶积对应的波场空间坐标,xm,zn为应用所述拉普拉斯算子后波场空间坐标,σ1和σ2分别为x方向、z方向上的高斯窗函数的方差。
5.一种基于常Q模型的粘弹介质地震波模拟系统,其特征在于,所述系统包括:
参数确定模块,用于确定地下介质的密度、地下介质的品质因子、横波在所述地下介质中的传播速度和纵波在所述地下介质中的传播速度;
地震子波选取模块,用于选取地震子波;
离散化函数确定模块,用于确定离散化的埃尔米特分布近似函数;
导数求解模型确定模块,用于确定所述埃尔米特分布近似函数的导数求解模型;
地震波波动参数确定模块,用于根据地下介质的密度、品质因子、横波在地下介质中的传播速度、纵波在地下介质中的传播速度、地震子波以及离散化的埃尔米特分布近似函数的导数求解模型,确定地震波的粘弹速度-应力波动参数;
地震波模拟模块,用于根据所述地震波的粘弹速度-应力波动参数,模拟所述地震波。
6.根据权利要求5所述的基于常Q模型的粘弹介质地震波模拟系统,其特征在于,所述离散化函数确定模块,具体包括:
参数确定单元,用于确定所述离散化的埃尔米特分布近似函数中的参数:高斯窗函数的方差、埃尔米特多项式阶数以及埃尔米特分布近似函数网格离散化算子长度;
离散化函数确定单元,用于根据所述参数确定离散化的埃尔米特分布近似函数其中,δM为埃尔米特分布近似泛函对狄拉克函数的近似离散网格化算子,x为网格中心点,x'为网格的空间坐标,x-x'为网格上每一点到中心点的空间距离,满足|x-x'|/Δx<(W-1)/2,σ为高斯窗函数的方差,M为埃尔米特多项式阶数,W为及埃尔米特分布近似函数网格离散化算子长度,Hi(x)为i阶埃尔米特多项式,
7.根据权利要求5所述的基于常Q模型的粘弹介质地震波模拟系统,其特征在于,所述导数求解模型确定模块,具体包括:
整数阶导数求解模型确定单元,用于确定所述埃尔米特分布近似函数的整数阶导数求解模型为其中,α为导数阶数;
分数阶导数求解模型确定单元,用于确定所述埃尔米特分布近似函数的分数阶导数求解模型为其中,F和F-1分别表示傅里叶变换和傅里叶反变换,为虚数单位,k为波数。
8.根据权利要求5所述的基于常Q模型的粘弹介质地震波模拟系统,其特征在于,所述地震波波动参数确定模块,具体包括:
波动方程获取单元,用于获取地震波的粘弹速度-应力波动方程;
波动参数离散化单元,用于对所述粘弹速度-应力波动方程中的速度、应变分量的时间偏导数做差分离散,得到粘弹速度-应力波动参数的离散形式
其中,分别为n时刻对应的x方向、z方向的地下介质质点的振动速度,分别为n时刻对应的x方向、xz方向、z方向的地下介质质点的应力,分别为n时刻对应的x方向、xz方向、z方向地下介质质点的应变,fx、fz分别为n时刻对应的x方向、z方向所述地震子波的体力分量,ρ为地下介质的密度,Δt为时间采样间隔, 为中间变量,ω0为设定参考频率,γp、γs为空变分数阶阶数,Qp、Qs分别为纵波品质因子和横波品质因子,为体积弹性模量,为剪切弹性模量,和分别为所述设定频率下对应的纵波和横波的速度,(-▽2)γ为γ阶拉普拉斯算子,Dx、Dz为空间偏导数算子,DxU表示对波场U应用空间偏导数算子,DzU表示对波场U应用空间偏导数算子,Δx、Δz为空间网格间距,为差分系数,N为差分阶数,Δx,Δz分别为空间网格间距,W1、W2分别为x方向、z方向上的网格点数,u为应用所述拉普拉斯算子的波场,xi,zj为应用所述拉普拉斯算子所需的空间褶积对应的波场空间坐标,xm,zn为应用所述拉普拉斯算子后波场空间坐标,σ1和σ2分别为x方向、z方向上的高斯窗函数的方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810945943.XA CN109143340B (zh) | 2018-08-20 | 2018-08-20 | 一种基于常q模型的粘弹介质地震波模拟方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810945943.XA CN109143340B (zh) | 2018-08-20 | 2018-08-20 | 一种基于常q模型的粘弹介质地震波模拟方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109143340A true CN109143340A (zh) | 2019-01-04 |
CN109143340B CN109143340B (zh) | 2020-03-10 |
Family
ID=64790178
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810945943.XA Expired - Fee Related CN109143340B (zh) | 2018-08-20 | 2018-08-20 | 一种基于常q模型的粘弹介质地震波模拟方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109143340B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109375251A (zh) * | 2018-09-29 | 2019-02-22 | 山东大学 | 利用城市既有地下空间与地表的探测方法及系统 |
CN111009087A (zh) * | 2020-01-08 | 2020-04-14 | 桂林赛普电子科技有限公司 | 一种矿区防越界盗采监测预警方法及系统 |
CN111025396A (zh) * | 2020-01-06 | 2020-04-17 | 中国石油化工股份有限公司 | 基于人工智能算法的油藏物性参数地震预测方法 |
CN112327358A (zh) * | 2020-10-29 | 2021-02-05 | 中国石油大学(华东) | 一种粘滞性介质中声波地震数据正演模拟方法 |
CN112394391A (zh) * | 2019-12-23 | 2021-02-23 | 中国海洋石油集团有限公司 | 雷克子波分数阶导数的地震相定量表征方法 |
CN113341455A (zh) * | 2021-06-24 | 2021-09-03 | 中国石油大学(北京) | 一种粘滞各向异性介质地震波数值模拟方法、装置及设备 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103823239A (zh) * | 2013-10-13 | 2014-05-28 | 中国石油集团西北地质研究所 | 频率域优化混合交错网格有限差分正演模拟方法 |
CN104732093A (zh) * | 2015-03-30 | 2015-06-24 | 西安交通大学 | 一种基于弥散黏滞性波动方程的fct-fdm正演模拟方法 |
US20160091623A1 (en) * | 2013-12-30 | 2016-03-31 | Petrochina Company Limited | Method and Device for Estimating Quality Factor Based on Zero Offset Vertical Seismic Profile Data |
US20160209537A1 (en) * | 2015-01-21 | 2016-07-21 | Baker Hughes Incorporated | Processing seismic data to remove noise |
JP2016133497A (ja) * | 2015-01-22 | 2016-07-25 | 株式会社セオコンプ | 地下構造探査システム及び地下構造探査方法 |
CN106842306A (zh) * | 2017-04-18 | 2017-06-13 | 中国科学院地质与地球物理研究所 | 一种全局优化的交错网格有限差分正演模拟方法和装置 |
CN106959469A (zh) * | 2017-04-14 | 2017-07-18 | 中国石油天然气股份有限公司 | 地震波的速度及衰减模拟分析方法及装置 |
US20180024263A1 (en) * | 2016-07-25 | 2018-01-25 | Schlumberger Technology Corporation | Seismic Spectral Balancing |
CN107894612A (zh) * | 2017-10-23 | 2018-04-10 | 中国地质大学(武汉) | 一种q吸收衰减补偿的声波阻抗反演方法及系统 |
-
2018
- 2018-08-20 CN CN201810945943.XA patent/CN109143340B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103823239A (zh) * | 2013-10-13 | 2014-05-28 | 中国石油集团西北地质研究所 | 频率域优化混合交错网格有限差分正演模拟方法 |
US20160091623A1 (en) * | 2013-12-30 | 2016-03-31 | Petrochina Company Limited | Method and Device for Estimating Quality Factor Based on Zero Offset Vertical Seismic Profile Data |
US20160209537A1 (en) * | 2015-01-21 | 2016-07-21 | Baker Hughes Incorporated | Processing seismic data to remove noise |
JP2016133497A (ja) * | 2015-01-22 | 2016-07-25 | 株式会社セオコンプ | 地下構造探査システム及び地下構造探査方法 |
CN104732093A (zh) * | 2015-03-30 | 2015-06-24 | 西安交通大学 | 一种基于弥散黏滞性波动方程的fct-fdm正演模拟方法 |
US20180024263A1 (en) * | 2016-07-25 | 2018-01-25 | Schlumberger Technology Corporation | Seismic Spectral Balancing |
CN106959469A (zh) * | 2017-04-14 | 2017-07-18 | 中国石油天然气股份有限公司 | 地震波的速度及衰减模拟分析方法及装置 |
CN106842306A (zh) * | 2017-04-18 | 2017-06-13 | 中国科学院地质与地球物理研究所 | 一种全局优化的交错网格有限差分正演模拟方法和装置 |
CN107894612A (zh) * | 2017-10-23 | 2018-04-10 | 中国地质大学(武汉) | 一种q吸收衰减补偿的声波阻抗反演方法及系统 |
Non-Patent Citations (1)
Title |
---|
罗伟平 等: "谐振Q模型介质地震波场数值模拟", 《地球物理学进展》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109375251A (zh) * | 2018-09-29 | 2019-02-22 | 山东大学 | 利用城市既有地下空间与地表的探测方法及系统 |
CN109375251B (zh) * | 2018-09-29 | 2021-04-13 | 山东大学 | 利用城市既有地下空间与地表的探测方法及系统 |
CN112394391A (zh) * | 2019-12-23 | 2021-02-23 | 中国海洋石油集团有限公司 | 雷克子波分数阶导数的地震相定量表征方法 |
CN111025396A (zh) * | 2020-01-06 | 2020-04-17 | 中国石油化工股份有限公司 | 基于人工智能算法的油藏物性参数地震预测方法 |
CN111025396B (zh) * | 2020-01-06 | 2021-11-05 | 中国石油化工股份有限公司 | 基于人工智能算法的油藏物性参数地震预测方法 |
CN111009087A (zh) * | 2020-01-08 | 2020-04-14 | 桂林赛普电子科技有限公司 | 一种矿区防越界盗采监测预警方法及系统 |
CN112327358A (zh) * | 2020-10-29 | 2021-02-05 | 中国石油大学(华东) | 一种粘滞性介质中声波地震数据正演模拟方法 |
CN112327358B (zh) * | 2020-10-29 | 2022-05-24 | 中国石油大学(华东) | 一种粘滞性介质中声波地震数据正演模拟方法 |
CN113341455A (zh) * | 2021-06-24 | 2021-09-03 | 中国石油大学(北京) | 一种粘滞各向异性介质地震波数值模拟方法、装置及设备 |
CN113341455B (zh) * | 2021-06-24 | 2024-02-09 | 中国石油大学(北京) | 一种粘滞各向异性介质地震波数值模拟方法、装置及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN109143340B (zh) | 2020-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109143340A (zh) | 一种基于常q模型的粘弹介质地震波模拟方法及系统 | |
CN103630933B (zh) | 基于非线性优化的时空域交错网格有限差分方法和装置 | |
CN104965223B (zh) | 粘声波全波形反演方法及装置 | |
CN105425289B (zh) | 确定低频波阻抗的方法和装置 | |
CN108549100B (zh) | 基于非线性高次拓频的时间域多尺度全波形反演方法 | |
CN104133241B (zh) | 波场分离方法和装置 | |
CN101425152B (zh) | 一种基于变学习率神经网络的fir滤波器的设计方法 | |
CN107798156B (zh) | 一种频率域2.5维粘弹性波数值模拟方法及装置 | |
CN108873063B (zh) | 一种微地震矩张量反演的方法及装置 | |
CN106033124A (zh) | 一种基于随机最优化的多震源粘声最小二乘逆时偏移方法 | |
CN113341455B (zh) | 一种粘滞各向异性介质地震波数值模拟方法、装置及设备 | |
CN109946742B (zh) | 一种TTI介质中纯qP波地震数据模拟方法 | |
CN110542928A (zh) | 基于vti各向异性传播矩阵的地震响应模拟方法 | |
CN106324675A (zh) | 一种宽频地震波阻抗低频信息预测方法及系统 | |
CN109212599A (zh) | 一种地震数据的叠前同步反演方法 | |
CN106569262A (zh) | 低频地震数据缺失下的背景速度模型重构方法 | |
CN107589452B (zh) | 纵波和转换波的数据匹配方法和装置 | |
CN103744114B (zh) | 基于零偏垂直地震剖面数据估计品质因子的方法和装置 | |
CN107229075B (zh) | 深度域地震子波的确定方法和装置 | |
CN111257930B (zh) | 一种黏弹各向异性双相介质区域变网格求解算子 | |
Lu et al. | Numerical modelling method for wave propagation in a linear viscoelastic medium with singular memory | |
CN104732093B (zh) | 一种基于弥散黏滞性波动方程的fct‑fdm正演模拟方法 | |
CN110658558A (zh) | 吸收衰减介质叠前深度逆时偏移成像方法及系统 | |
CN109425892A (zh) | 地震子波的估计方法及系统 | |
CN113866823A (zh) | 一种粘声各向异性介质中的正演成像方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200310 Termination date: 20200820 |