CN112444865B - 勘探缝洞储层的方法、***、装置及存储介质 - Google Patents
勘探缝洞储层的方法、***、装置及存储介质 Download PDFInfo
- Publication number
- CN112444865B CN112444865B CN201910804363.3A CN201910804363A CN112444865B CN 112444865 B CN112444865 B CN 112444865B CN 201910804363 A CN201910804363 A CN 201910804363A CN 112444865 B CN112444865 B CN 112444865B
- Authority
- CN
- China
- Prior art keywords
- random
- function
- dimensional
- fracture
- coordinate system
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000005311 autocorrelation function Methods 0.000 claims abstract description 44
- 238000001228 spectrum Methods 0.000 claims abstract description 39
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 77
- 230000008569 process Effects 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 230000003595 spectral effect Effects 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004891 communication Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 230000005236 sound signal Effects 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000009918 complex formation Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 201000001098 delayed sleep phase syndrome Diseases 0.000 description 1
- 208000033921 delayed sleep phase type circadian rhythm sleep disease Diseases 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000019800 disodium phosphate Nutrition 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004549 pulsed laser deposition Methods 0.000 description 1
- 238000004088 simulation Methods 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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
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
技术领域
本公开涉及地震勘探技术领域,特别地涉及一种勘探缝洞储层的方法、装置及存储介质。
背景技术
油藏(或油储)地球物理学(Reservoir Geophysics)从上世纪80年代提出到现在,国内外地球物理学家对其广泛关注,成了一个热点研究领域。油藏地球物理学的基础为非均匀介质模型,主要内容就是精细研究油藏(或油储)特征、参数和油气储层,为油气勘探特别是油气开发服务。因为油藏地球物理学与复杂、精细的非均匀介质联系很大,所以用常规方法很难完整描述,因此要提出一种灵活、方便、能完整描述油藏非均匀性的正演模型,由此产生出了随机介质模型理论。
发明内容
本公开提供一种勘探缝洞储层的方法、装置及存储介质,以解决相关技术中地震勘探研究缺少具有方向性特征的三维随机介质建模的技术问题。
为实现上述目的,本公开实施例的第一方面,提供一种勘探缝洞储层的方法,所述方法包括:
在三维坐标系下构建指数型椭圆自相关函数;
通过对构建的所述指数型椭圆自相关函数进行傅里叶变换,获得随机扰动函数的功率谱函数;
根据所述功率谱函数,获得随机扰动函数并规范化;
根据规范化后的所述随机扰动函数和缝洞储层内的随机介质的背景介质参数,构建三维随机介质模型;
在构建的所述三维随机介质模型中输入待勘探的缝洞储层内的随机介质的背景介质参数,以获得所述待勘探的缝洞储层对应的随机介质模型;
利用所述随机介质模型进行缝洞储层的勘探。
可选地,所述在三维坐标系下构建指数型椭圆自相关函数,包括:
在三维全局坐标系中选取点的三个局部自相关长度;
根据三维局部坐标系与所述三维全局坐标系中的点与点之间的空间偏移量,选取与三个所述局部自相关长度对应的三个旋转角度;
根据三个所述局部自相关长度和三个所述旋转角度,构建指数型椭圆自相关函数。
可选地,设定a(x,y,z)、b(x,y,z)、c(x,y,z)分别为三维全局坐标系OXYZ中(x,y,z)点的三个局部自相关长度;根据三个所述局部自相关长度和三个所述旋转角度,构建的指数型椭圆自相关函数的表达式为:
式(1)中:a(x,y,z)、b(x,y,z)、c(x,y,z)分别为三维全局坐标系OXYZ点(x,y,z)的局部自相关函数和局部自相关长度;x′1,y′1,z′1为三维局部坐标系O'X1'Y1'Z1'中相对于点(x,y,z)的空间偏移量,三维局部坐标系O'X1'Y1'Z1'是由三维局部坐标系O'X1Y1 Z1沿逆时针方向旋转角度α,β,γ得到,所述空间偏移量的表达式为:
式(2)中:D为空间偏移量,α,β,γ为a(x,y,z)、b(x,y,z)、c(x,y,z)对应的旋转角度。
可选地,构建的所述指数型椭圆自相关函数进行傅里叶变换的公式为:
式(3)中:Γ表示空间傅里叶变换算子,φ(kx,ky,kz)表示随机扰动函数的功率谱函数,k为空间波数,kx为x坐标轴方向的空间波数,ky为y坐标轴方向的空间波数,kz为z坐标轴方向的空间波数。
可选地,所述根据所述功率谱函数,获得随机扰动函数,包括:
由功率谱函数,根据随机过程的方法产生对应的随机谱序列;
应用随机过程的谱展开公式,获得由构建的所述指数型椭圆自相关函数描述的随机扰动函数。
可选地,规范化所述随机扰动函数,包括:
按随机扰动均值和方差的要求,规范化所述随机扰动函数。
可选地,根据规范化后的所述随机扰动函数和随机介质的背景介质参数,构建的三维随机介质模型的表达式为:
式(3)中:ρ为密度,λ、μ为拉梅参数,ρ0、λ0、μ0为随机介质的背景介质参数;δρ、δλ、δμ为非均匀扰动量。
本公开实施例的第二方面,提供一种勘探缝洞储层的***,所述***包括:
第一构建模块,被配置成在三维坐标系下构建指数型椭圆自相关函数;
傅里叶变换模块,被配置成通过对构建的所述指数型椭圆自相关函数进行傅里叶变换,获得随机扰动函数的功率谱函数;
第一获得模块,被配置成根据所述功率谱函数,获得随机扰动函数并规范化;
第二构建模块,被配置成根据规范化后的所述随机扰动函数和缝洞储层内的随机介质的背景介质参数,构建三维随机介质模型;
第二获得模块,被配置成在构建的所述三维随机介质模型中输入待勘探的缝洞储层内的随机介质的背景介质参数,以获得所述待勘探的缝洞储层对应的随机介质模型,以利用所述随机介质模型进行缝洞储层的勘探。
本公开实施例的第三方面,提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述第一方面中任一项所述方法的步骤。
本公开实施例的第四方面,提供一种勘探缝洞储层的装置,包括:
存储器,其上存储有计算机程序;以及
处理器,用于执行所述存储器中的所述计算机程序,以实现上述第一方面中任一项所述方法的步骤。
采用上述技术方案,至少能够达到如下技术效果:
本公开基于随机介质建模理论,开展具有方向特征的三维随机介质模型构建,通过引入自相关角度来表征地层中复杂地层不同的角度,对缝洞型储层进行三维建模,使得构建的模型更加符合实际的缝洞储层地质情况,为缝洞地震响应模式及机理进行研究提供模型数据,更好地为野外实际深层微幅缝洞储层地震响应特征进行标定和识别,提高野外钻井的准确率,为油气增产提供支撑。
本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本公开的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本公开,但并不构成对本公开的限制。在附图中:
图1是本公开一示例性实施例示出的一种勘探缝洞储层的方法流程图。
图2是本公开一示例性实施例示出的一种勘探缝洞储层的方法包括的步骤中构建指数型椭圆自相关函数的流程图。
图3是本公开一示例性实施例示出的三维局部坐标系与三维全局坐标系之间的关系示意图。
图4是本公开一示例性实施例示出的两个三维局部坐标系之间的关系示意图。
图5是本公开一示例性实施例示出的一种勘探缝洞储层的方法包括的步骤中获得随机扰动函数的流程图。
图6是本公开一示例性实施例示出的三维幅角示意图。
图7是本公开一示例性实施例示出的一种三维随机介质模型的示意图。
图8是本公开一示例性实施例示出的一种勘探缝洞储层的***框图。
图9是本公开一示例性实施例示出的一种勘探缝洞储层的装置框图。
具体实施方式
以下将结合附图及实施例来详细说明本公开的实施方式,借此对本公开如何应用技术手段来解决技术问题,并达到相应技术效果的实现过程能充分理解并据以实施。本申请实施例以及实施例中的各个特征,在不相冲突前提下可以相互结合,所形成的技术方案均在本公开的保护范围之内。
本公开发明人经研究发现,相关技术中,随机介质建模主要是基于二维情况,要合理地表征缝洞储层及地震响应特征应该在三维条件下进行,而针对三维随机介质建模技术的研究比较少,且没有发现涉及到缝洞储层不同的倾角因素的三维随机介质建模的研究,因此,本公开提供了一种勘探缝洞储层的方法,为实际储层地震响应的研究打下基础。
实施例一
图1是本公开一示例性实施例示出的一种勘探缝洞储层的方法流程图,以解决相关技术中地震勘探研究缺少具有方向性特征的三维随机介质建模的技术问题。如图1所示,本实施例示出的勘探缝洞储层的方法可以包括以下步骤:
S11,在三维坐标系下构建指数型椭圆自相关函数。
S12,通过对构建的所述指数型椭圆自相关函数进行傅里叶变换,获得随机扰动函数的功率谱函数。
S13,根据所述功率谱函数,获得随机扰动函数并规范化。
S14,根据规范化后的所述随机扰动函数和缝洞储层内的随机介质的背景介质参数,构建三维随机介质模型。
S15,在构建的所述三维随机介质模型中输入待勘探的缝洞储层内的随机介质的背景介质参数,以获得所述待勘探的缝洞储层对应的随机介质模型;
S16,利用所述随机介质模型进行缝洞储层的勘探。
如图2所示,步骤S11,即在三维坐标系下构建指数型椭圆自相关函数,可以包括以下步骤:
S111,在三维全局坐标系中选取点的三个局部自相关长度。
S112,根据三维局部坐标系与所述三维全局坐标系中的点与点之间的空间偏移量,选取与三个所述局部自相关长度对应的三个旋转角度。
S113,根据三个所述局部自相关长度和三个所述旋转角度,构建指数型椭圆自相关函数。
在步骤S111中,假设三维全局坐标系OXYZ中点(x,y,z)的三个局部自相关长度为a(x,y,z)、b(x,y,z)、c(x,y,z)。x′1,y′1,z′1为三维局部坐标系O'X1'Y1'Z1'中相对于(x,y,z)点的空间偏移量,三维局部坐标系O'X1'Y1'Z1'与三维全局坐标系OXYZ之间的关系如图3所示,这二者的关系可以根据坐标转换公式(5)得到。
式(5)D为空间偏移量。进而可以通过坐标旋转方法来构造坐标转换矩阵,如式(2)所示。
其中,如图4所示,三维局部坐标系O'X1'Y1'Z1'是由三维局部坐标系O'X1Y1 Z1沿各自坐标逆时针方向旋转角度α,β,γ得到。
进而根据选取的自相关长度a(x,y,z)、b(x,y,z)、c(x,y,z)和旋转角度α(x,y,z)、β(x,y,z)、γ(x,y,z),构建的指数型椭圆自相关函数的表达式可以为:
式(1)中:a(x,y,z)、b(x,y,z)、c(x,y,z)分别为三维全局坐标系OXYZ点(x,y,z)的局部自相关函数和局部自相关长度。
构建好指数型椭圆自相关函数后,可以执行步骤S12,通过对构建的所述指数型椭圆自相关函数进行傅里叶变换,获得随机扰动函数的功率谱函数。
构建的所述指数型椭圆自相关函数进行傅里叶变换的公式可以为:
式(3)中:Γ表示空间傅里叶变换算子,φ(kx,ky,kz)表示随机扰动函数的功率谱函数,k为空间波数,kx为x坐标轴方向的空间波数,ky为y坐标轴方向的空间波数,kz为z坐标轴方向的空间波数。功率谱函数φ(kx,ky,kz)即为随机扰动σ(x,y,z)的功率谱。
获得功率谱函数后,则可以执行步骤S13,根据所述功率谱函数,获得随机扰动函数并规范化。如图5所示,所述根据所述功率谱函数,获得随机扰动函数,可以包括以下步骤:
S131,由功率谱函数φ(kx,ky,kz),根据随机过程的方法产生对应的随机谱序列。
S132,应用随机过程的谱展开公式,获得由构建的所述指数型椭圆自相关函数描述的随机扰动函数σ(x,y,z)。其中,随机谱序列和随机扰动函数σ(x,y,z)均为离散数据。
获得随机扰动函数σ(x,y,z)后,按随机扰动均值和方差的要求,将随机扰动函数σ(x,y,z)规范化。其中,规范随机扰动函数数值的均值等于0,并且其方差为一定的数值(一般为背景值扰动的百分数来取值,如0.01,那么设置的扰动值满足均值为零,方差为0.01即可),其规范化的过程为数学方面的公知技术,为了说明书的简洁,在此不具体展开。
规范化获得的随机扰动函数σ(x,y,z)后,执行步骤S14,根据规范化后的所述随机扰动函数和缝洞储层内的随机介质的背景介质参数,构建随机介质模型。
相关技术中,随机介质可以理解为在大尺度均匀背景下分布着小尺度非均匀扰动的理论模型,其构建理论可以描述为以下表达式:
m(X)=m0+δm(X)=m0(1+σ(X)) (6)
式(6)中:m0表示大尺度均匀介质;δm(X)表示小尺度非均匀扰动介质;σ(X)为介质的随机扰动特性;m(X)为随机介质;X=(x,y,z)为空间位置矢量。
各向同性弹性介质可以由密度ρ和拉梅参数λ、μ确定,当m(X)包含ρ,λ,μ时,式(6)可以表示为:
式(7)中:ρ0、λ0、μ0为随机介质的背景介质参数,假设为常数或随空间坐标(x,y,z)缓慢变化;δρ、δλ、δμ为加在上述背景上的非均匀扰动量,并假设其均值为零、方差为定值且和某一自相关函数的空间平稳随机过程。
由于式(6)和式(7)没有考虑到方向性特征,比如缝洞储层不同的倾角因素,因此对于复杂地层构造的随机介质模型并不准确。本公开在建模过程中加入了方向性特征进去,在规范化获得的随机扰动函数σ(x,y,z)后,可以得到所要构造的以作为自相关函数、具有指定均值及方差的随机介质模型,构建的三维随机介质模型的表达式为:
式(4)中:ρ为密度,λ、μ为拉梅参数,ρ0、λ0、μ0为缝洞储层内的随机介质的背景介质参数;δρ、δλ、δμ为非均匀扰动量。
三维随机介质模型构建完成后,执行步骤S15,在构建的所述三维随机介质模型中输入待勘探的缝洞储层内的随机介质的背景介质参数,以获得所述待勘探的缝洞储层对应的随机介质模型,进而,利用所述随机介质模型进行缝洞储层的勘探。由于实际介质情况是含有非均质体,导致地下介质速度表现为一定的随机性,并且非均质体在地下介质的空间展布各式各样,本公开增加了随机介质空间的方向性,更能合理构建不同情况下的随机介质模型,通过不同随机介质产状的模型进行波场模拟,分析剖面波场特征,利用室内分析结果和经验指导野外实际地震勘探,为野外实际深层微幅缝洞储层地震响应特征进行标定和识别,提高野外钻井的准确率,为油气增产提供数据支撑。
下面通过一个实例,把具有方向特征的三维随机模型整个建立过程进行展示。以下为实例中三维随机模型参数:
模型大小:Nx=100,Ny=100,Nz=100;
网格间距:dx=1m,dy=1m,dz=1m;
速度均值:v0=5500m/s2;
速度扰动标准差:10%;
自相关长度:a=3m,b=10m,c=20m;自相关角度:α=0°,β=0°,γ=75°;
指数型椭圆自相关函数采用式(1)进行,首先按照给出的模型范围,完成指数椭圆函数傅里叶变换需要的幅角,幅角范围在0-2π之间并且呈随机特征变化,截取的一个二维剖面展示如图6所示,再结合变换函数完成扰动功率谱的实现,根据扰动均值和方差值进一步完成扰动量的求取,结合背景速度场,根据式(4)完成最后的含有自相关角度的三维随机介质模型,如图7所示。利用图7构建的三维随机介质模型,就可以将待勘探的缝洞储层内的随机介质的背景介质参数输入该模型中,获得待勘探的缝洞储层对应的随机介质模型。
本公开基于随机介质建模理论,开展具有方向特征的三维随机介质模型构建,通过引入自相关角度来表征地层中缝洞型储层不同的角度,对缝洞型储层进行三维建模,使得构建的模型更加符合实际地质情况,为缝洞地震响应模式及机理进行研究提供模型数据,更好地为野外实际深层微幅缝洞储层地震响应特征进行标定和识别,提高野外钻井的准确率,为油气增产提供支撑。
值得说明的是,对于图1所示的方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本公开并不受所描述的动作顺序的限制。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作并不一定是本公开所必须的。
实施例二
图8是本公开一示例性实施例示出的一种勘探缝洞储层的***,如图8所示,所述勘探缝洞储层的***300包括:
第一构建模块310,被配置成在三维坐标系下构建指数型椭圆自相关函数;
傅里叶变换模块320,被配置成通过对构建的所述指数型椭圆自相关函数进行傅里叶变换,获得随机扰动函数的功率谱函数;
第一获得模块330,被配置成根据所述功率谱函数,获得随机扰动函数并规范化;
第二构建模块340,被配置成根据规范化后的所述随机扰动函数和缝洞储层内的随机介质的背景介质参数,构建随机介质模型;
第二获得模块350,被配置成在构建的所述三维随机介质模型中输入待勘探的缝洞储层内的随机介质的背景介质参数,以获得所述待勘探的缝洞储层对应的随机介质模型,以利用所述随机介质模型进行缝洞储层的勘探。
可选地,第一构建模块310还被配置成:在三维全局坐标系中选取点的三个局部自相关长度;根据三维局部坐标系与所述三维全局坐标系中的点与点之间的空间偏移量,选取与三个所述局部自相关长度对应的三个旋转角度;根据三个所述局部自相关长度和三个所述旋转角度,构建指数型椭圆自相关函数。
可选地,设定a(x,y,z)、b(x,y,z)、c(x,y,z)分别为三维全局坐标系OXYZ中(x,y,z)点的三个局部自相关长度;根据三个所述局部自相关长度和三个所述旋转角度,第一构建模块310构建的指数型椭圆自相关函数的表达式为:
式(1)中:a(x,y,z)、b(x,y,z)、c(x,y,z)分别为三维全局坐标系OXYZ点(x,y,z)的局部自相关函数和局部自相关长度;x′1,y′1,z′1为三维局部坐标系O1X′1Y′1Z′1中相对于点(x,y,z)的空间偏移量,三维局部坐标系O1X′1Y′1Z′1是由三维局部坐标系O1X1Y1Z1沿各自坐标逆时针方向旋转角度α,β,γ得到,所述空间偏移量的表达式为:
式(2)中:D为空间偏移量,α,β,γ为a(x,y,z)、b(x,y,z)、c(x,y,z)对应的旋转角度。
可选地,傅里叶变换模块320进行傅里叶变换的公式为:
式(3)中:Γ表示空间傅里叶变换算子,φ(kx,ky,kz)表示随机扰动函数的功率谱函数,k为空间波数,kx为x坐标轴方向的空间波数,ky为y坐标轴方向的空间波数,kz为z坐标轴方向的空间波数。
可选地,获得模块340还被配置成:由功率谱函数,根据随机过程的方法产生对应的随机谱序列;应用随机过程的谱展开公式,获得由构建的所述指数型椭圆自相关函数描述的随机扰动函数。
可选地,获得模块340还被配置成包括:按随机扰动均值和方差的要求,规范化所述随机扰动函数。
可选地,第二构建模块350构建的随机介质模型的表达式为:
式(3)中:ρ0、λ0、μ0为随机介质的背景介质参数;δρ、δλ、δμ为非均匀扰动量。
关于上述实施例中的***,其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述,此处将不做详细阐述说明。
实施例三
本公开还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述任一项可选实施例所述的勘探缝洞储层的方法步骤。
其中,在所述处理器上运行的勘探缝洞储层的方法的计算机程序被执行时所实现的方法可参照本公开勘探缝洞储层的方法的具体实施例,此处不再赘述。
所述处理器可以是一种集成电路芯片,具有信息处理能力。所述处理器可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、网络处理器(NetworkProcessor,NP)等。
实施例四
本公开还提供一种勘探缝洞储层的装置,包括:
存储器,其上存储有计算机程序;以及
处理器,用于执行所述存储器中的所述计算机程序,以实现上述任一项可选实施例所述的勘探缝洞储层的方法步骤。
图9是根据一示例性实施例示出的一种勘探缝洞储层的装置400的框图。如图9所示,该装置400可以包括:处理器401,存储器402,多媒体组件403,输入/输出(I/O)接口404,以及通信组件405。
其中,处理器401用于控制该装置400的整体操作,以完成上述的三维随机介质模型的构建方法中的全部或部分步骤。存储器402用于存储各种类型的数据以支持在该装置400的操作,这些数据例如可以包括用于在该装置400上操作的任何应用程序或方法的指令,以及应用程序相关的数据。该存储器402可以由任何类型的易失性或非易失性存储设备或者它们的组合实现,例如静态随机存取存储器(Static Random Access Memory,简称SRAM),电可擦除可编程只读存储器(Electrically Erasable Programmable Read-OnlyMemory,简称EEPROM),可擦除可编程只读存储器(Erasable Programmable Read-OnlyMemory,简称EPROM),可编程只读存储器(Programmable Read-Only Memory,简称PROM),只读存储器(Read-Only Memory,简称ROM),磁存储器,快闪存储器,磁盘或光盘。多媒体组件403可以包括屏幕和音频组件。其中屏幕例如可以是触摸屏,音频组件用于输出和/或输入音频信号。例如,音频组件可以包括一个麦克风,麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器402或通过通信组件405发送。音频组件还包括至少一个扬声器,用于输出音频信号。I/O接口404为处理器401和其他接口模块之间提供接口,上述其他接口模块可以是键盘,鼠标,按钮等。这些按钮可以是虚拟按钮或者实体按钮。通信组件405用于该装置400与其他设备之间进行有线或无线通信。无线通信,例如Wi-Fi,蓝牙,近场通信(Near Field Communication,简称NFC),2G、3G或4G,或它们中的一种或几种的组合,因此相应的该通信组件405可以包括:Wi-Fi模块,蓝牙模块,NFC模块。
在一示例性实施例中,装置400可以被一个或多个应用专用集成电路(Application Specific Integrated Circuit,简称ASIC)、数字信号处理器(DigitalSignal Processor,简称DSP)、数字信号处理设备(Digital Signal Processing Device,简称DSPD)、可编程逻辑器件(Programmable Logic Device,简称PLD)、现场可编程门阵列(Field Programmable Gate Array,简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现,用于执行上述的勘探缝洞储层的方法。
在另一示例性实施例中,还提供了一种包括程序指令的计算机可读存储介质,例如包括程序指令的存储器402,上述程序指令可由装置400的处理器401执行以完成上述的勘探缝洞储层的方法。
以上结合附图详细描述了本公开的优选实施方式,但是,本公开并不限于上述实施方式中的具体细节,在本公开的技术构思范围内,可以对本公开的技术方案进行多种简单变型,这些简单变型均属于本公开的保护范围。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合。为了避免不必要的重复,本公开对各种可能的组合方式不再另行说明。
此外,本公开的各种不同的实施方式之间也可以进行任意组合,只要其不违背本公开的思想,其同样应当视为本公开所公开的内容。
Claims (8)
1.一种勘探缝洞储层的方法,其特征在于,所述方法包括:
在三维坐标系下构建指数型椭圆自相关函数;具体包括:在三维全局坐标系中选取点的三个局部自相关长度;根据三维局部坐标系与所述三维全局坐标系中的点与点之间的空间偏移量,选取与三个所述局部自相关长度对应的三个旋转角度;根据三个所述局部自相关长度和三个所述旋转角度,构建指数型椭圆自相关函数;设定分别为三维全局坐标系/>中/>点的三个局部自相关长度;根据三个所述局部自相关长度和三个所述旋转角度,构建的指数型椭圆自相关函数的表达式为:
(1)
式(1)中:,/>分别为三维全局坐标系点/>的局部自相关函数和局部自相关长度;/>为三维局部坐标系中相对于点/>的空间偏移量,三维局部坐标系/>是由三维局部坐标系沿逆时针方向旋转角度/>得到,所述空间偏移量的表达式为:
(2)
式(2)中:D为空间偏移量,为/>对应的旋转角度;
通过对构建的所述指数型椭圆自相关函数进行傅里叶变换,获得随机扰动函数的功率谱函数;
根据所述功率谱函数,获得随机扰动函数并规范化;
根据规范化后的所述随机扰动函数和缝洞储层内的随机介质的背景介质参数,构建三维随机介质模型;
在构建的所述三维随机介质模型中输入待勘探的缝洞储层内的随机介质的背景介质参数,以获得所述待勘探的缝洞储层对应的随机介质模型;
利用所述随机介质模型进行缝洞储层的勘探。
2.根据权利要求1所述的方法,其特征在于,构建的所述指数型椭圆自相关函数进行傅里叶变换的公式为:
(3)
式(3)中:表示空间傅里叶变换算子,/>表示随机扰动函数的功率谱函数,k为空间波数,kx为x坐标轴方向的空间波数,ky为y坐标轴方向的空间波数,kz为z坐标轴方向的空间波数。
3.根据权利要求2所述的方法,其特征在于,所述根据所述功率谱函数,获得随机扰动函数,包括:
由功率谱函数,根据随机过程的方法产生对应的随机谱序列;
应用随机过程的谱展开公式,获得由构建的所述指数型椭圆自相关函数描述的随机扰动函数。
4.根据权利要求3所述的方法,其特征在于, 规范化所述随机扰动函数,包括:
按随机扰动均值和方差的要求,规范化所述随机扰动函数。
5.根据权利要求1所述的方法,其特征在于,根据规范化后的所述随机扰动函数和随机介质的背景介质参数,构建的三维随机介质模型的表达式为:
(4)
式(4)中:为密度 ,/>为拉梅参数,/>为缝洞储层内的随机介质的背景介质参数;/>为非均匀扰动量。
6.一种勘探缝洞储层的***,其特征在于,包括:
第一构建模块,被配置成在三维坐标系下构建指数型椭圆自相关函数;具体包括:在三维全局坐标系中选取点的三个局部自相关长度;根据三维局部坐标系与所述三维全局坐标系中的点与点之间的空间偏移量,选取与三个所述局部自相关长度对应的三个旋转角度;根据三个所述局部自相关长度和三个所述旋转角度,构建指数型椭圆自相关函数;设定分别为三维全局坐标系/>中/>点的三个局部自相关长度;根据三个所述局部自相关长度和三个所述旋转角度,构建的指数型椭圆自相关函数的表达式为:
(1)
式(1)中:,/>分别为三维全局坐标系点/>的局部自相关函数和局部自相关长度;/>为三维局部坐标系中相对于点/>的空间偏移量,三维局部坐标系/>是由三维局部坐标系沿逆时针方向旋转角度/>得到,所述空间偏移量的表达式为:
(2)
式(2)中:D为空间偏移量,为/>对应的旋转角度;
傅里叶变换模块,被配置成通过对构建的所述指数型椭圆自相关函数进行傅里叶变换,获得随机扰动函数的功率谱函数;
第一获得模块,被配置成根据所述功率谱函数,获得随机扰动函数并规范化;
第二构建模块,被配置成根据规范化后的所述随机扰动函数和缝洞储层内的随机介质的背景介质参数,构建三维随机介质模型;
第二获得模块,被配置成在构建的所述三维随机介质模型中输入待勘探的缝洞储层内的随机介质的背景介质参数,以获得所述待勘探的缝洞储层对应的随机介质模型,以利用所述随机介质模型进行缝洞储层的勘探。
7.一种勘探缝洞储层的装置,其特征在于,包括:
存储器,其上存储有计算机程序;以及
处理器,用于执行所述存储器中的所述计算机程序,以实现权利要求1至5中任一项所述方法的步骤。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现权利要求1至5中任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910804363.3A CN112444865B (zh) | 2019-08-28 | 2019-08-28 | 勘探缝洞储层的方法、***、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910804363.3A CN112444865B (zh) | 2019-08-28 | 2019-08-28 | 勘探缝洞储层的方法、***、装置及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112444865A CN112444865A (zh) | 2021-03-05 |
CN112444865B true CN112444865B (zh) | 2024-05-17 |
Family
ID=74741452
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910804363.3A Active CN112444865B (zh) | 2019-08-28 | 2019-08-28 | 勘探缝洞储层的方法、***、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112444865B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115267892A (zh) * | 2022-07-21 | 2022-11-01 | 北京化工大学 | 一种富有机质泥页岩地震响应数值模拟方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004090575A1 (en) * | 2003-04-09 | 2004-10-21 | Norsar | Method for simulating local prestack depth migrated seismic images |
EA200702257A2 (ru) * | 2007-11-15 | 2008-02-28 | Блаас Холдингс Лимитед | Способ радиолокационного зондирования земных недр и устройство для его реализации - комплекс георадиолокационной разведки |
CN102253415A (zh) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | 基于裂缝等效介质模型的地震响应模式建立方法 |
WO2014003597A1 (en) * | 2012-06-26 | 2014-01-03 | Schlumberger, Holding Limited | A method for determining pore volume characteristics and porous materials' matrix thermal conductivity |
CN103792573A (zh) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | 一种基于频谱融合的地震波阻抗反演方法 |
CN104077738A (zh) * | 2013-12-30 | 2014-10-01 | 辽宁师范大学 | 基于局部直方图特征的彩色图像水印方法 |
CN104977606A (zh) * | 2014-04-02 | 2015-10-14 | 中国石油化工股份有限公司 | 一种缝洞型储层地震数值模型建立方法 |
CN104992029A (zh) * | 2015-07-20 | 2015-10-21 | 中国科学院国家天文台 | 一种多尺度非均匀月壤层内离散随机介质建模方法 |
CN105354421A (zh) * | 2015-11-06 | 2016-02-24 | 吉林大学 | 一种随机导电媒质模型的大地电磁无网格数值模拟方法 |
CN107957594A (zh) * | 2017-11-15 | 2018-04-24 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据的椭圆校正方法、动校正方法及速度谱计算方法 |
CN109143362A (zh) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | 基于共散射角道集的散射波分离方法 |
CN109490965A (zh) * | 2018-10-15 | 2019-03-19 | 长江大学 | 一种定量评价地层非均匀性的方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2869421B1 (fr) * | 2004-04-27 | 2006-06-02 | Inst Francais Du Petrole | Methode de reconstruction d'un modele stochastique, representatif d'un milieu heterogene poreux, pour ameliorer son calage par les donnees de production |
US7619220B2 (en) * | 2005-11-30 | 2009-11-17 | Jeol Ltd. | Method of measuring aberrations and correcting aberrations using Ronchigram and electron microscope |
CN111615706A (zh) * | 2017-11-17 | 2020-09-01 | 脸谱公司 | 基于子流形稀疏卷积神经网络分析空间稀疏数据 |
-
2019
- 2019-08-28 CN CN201910804363.3A patent/CN112444865B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004090575A1 (en) * | 2003-04-09 | 2004-10-21 | Norsar | Method for simulating local prestack depth migrated seismic images |
EA200702257A2 (ru) * | 2007-11-15 | 2008-02-28 | Блаас Холдингс Лимитед | Способ радиолокационного зондирования земных недр и устройство для его реализации - комплекс георадиолокационной разведки |
CN102253415A (zh) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | 基于裂缝等效介质模型的地震响应模式建立方法 |
WO2014003597A1 (en) * | 2012-06-26 | 2014-01-03 | Schlumberger, Holding Limited | A method for determining pore volume characteristics and porous materials' matrix thermal conductivity |
CN103792573A (zh) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | 一种基于频谱融合的地震波阻抗反演方法 |
CN104077738A (zh) * | 2013-12-30 | 2014-10-01 | 辽宁师范大学 | 基于局部直方图特征的彩色图像水印方法 |
CN104977606A (zh) * | 2014-04-02 | 2015-10-14 | 中国石油化工股份有限公司 | 一种缝洞型储层地震数值模型建立方法 |
CN104992029A (zh) * | 2015-07-20 | 2015-10-21 | 中国科学院国家天文台 | 一种多尺度非均匀月壤层内离散随机介质建模方法 |
CN105354421A (zh) * | 2015-11-06 | 2016-02-24 | 吉林大学 | 一种随机导电媒质模型的大地电磁无网格数值模拟方法 |
CN109143362A (zh) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | 基于共散射角道集的散射波分离方法 |
CN107957594A (zh) * | 2017-11-15 | 2018-04-24 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据的椭圆校正方法、动校正方法及速度谱计算方法 |
CN109490965A (zh) * | 2018-10-15 | 2019-03-19 | 长江大学 | 一种定量评价地层非均匀性的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112444865A (zh) | 2021-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2879773C (en) | Multi-level reservoir history matching | |
US20200019882A1 (en) | Systems and Methods for Generating, Deploying, Discovering, and Managing Machine Learning Model Packages | |
MX2015004001A (es) | Propagacion de actualizaciones del plano de fractura. | |
JP2014535044A5 (zh) | ||
KR20190017454A (ko) | 인공 신경망을 이용한 자기장 기반 위치 추정 모델 생성 장치 및 방법 | |
CN111967169B (zh) | 二度体重力异常积分解数值模拟方法和装置 | |
US9581707B2 (en) | System for automated identification of surfaces for building of geologic hydrodynamic model of oil and gas deposit by seismic data | |
CN112444865B (zh) | 勘探缝洞储层的方法、***、装置及存储介质 | |
EP2975438A1 (en) | Multiscale method for reservoir models | |
Liang et al. | Uncertainty quantification of geologic model parameters in 3D gravity inversion by Hessian-informed Markov chain Monte Carlo | |
CN107870361A (zh) | 一种地震回折波层析成像方法、装置及终端设备 | |
CN110083851B (zh) | 气井井底压力的确定方法、装置及存储介质 | |
Fischer et al. | A cellular automata‐based deterministic inversion algorithm for the characterization of linear structural heterogeneities | |
CN113221409B (zh) | 一种有限元和边界元耦合的声波二维数值模拟方法和装置 | |
CN113608262B (zh) | 利用平动分量计算旋转分量的地震数据处理方法及装置 | |
KR20190042540A (ko) | 인공 신경망을 이용한 자기장 기반 위치 추정 모델 생성 장치 및 방법 | |
CN115809570A (zh) | 基于云平台的断面水深流速分布模型的建立方法及*** | |
Wu et al. | Micro-seismic source location determined by a modified objective function | |
CA3111539C (en) | Distributed sequencial gaussian simulation | |
US20150205003A1 (en) | Method of determining channelway trajectories | |
CN114254505A (zh) | 渗透率模型的构建方法、装置、电子设备及存储介质 | |
CN114063161A (zh) | 一种地震反演方法、装置及*** | |
US8494778B2 (en) | Variable grid for finite difference computation | |
WO2024125606A1 (zh) | 基于多源多任务深度学习的地震构造解释方法及装置 | |
US11994639B2 (en) | Seismic velocity modeling |
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 |