CN108469633B - 一种地层品质因子的计算方法及装置 - Google Patents
一种地层品质因子的计算方法及装置 Download PDFInfo
- Publication number
- CN108469633B CN108469633B CN201810123034.8A CN201810123034A CN108469633B CN 108469633 B CN108469633 B CN 108469633B CN 201810123034 A CN201810123034 A CN 201810123034A CN 108469633 B CN108469633 B CN 108469633B
- Authority
- CN
- China
- Prior art keywords
- amplitude spectrum
- low
- frequency
- bandwidth
- cutoff frequency
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 32
- 238000001228 spectrum Methods 0.000 claims abstract description 340
- 238000000034 method Methods 0.000 claims abstract description 120
- 238000001914 filtration Methods 0.000 claims abstract description 37
- 230000008569 process Effects 0.000 claims abstract description 23
- 238000003860 storage Methods 0.000 claims description 23
- 101000832225 Homo sapiens Stabilin-1 Proteins 0.000 claims description 8
- 102100024471 Stabilin-1 Human genes 0.000 claims description 8
- 101150109354 FEX2 gene Proteins 0.000 claims description 7
- 101100534471 Homo sapiens STAB2 gene Proteins 0.000 claims description 7
- 102100024470 Stabilin-2 Human genes 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 description 38
- 238000010586 diagram Methods 0.000 description 18
- 230000006870 function Effects 0.000 description 14
- 230000006872 improvement Effects 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 7
- 230000006854 communication Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 239000000047 product Substances 0.000 description 5
- 238000009499 grossing Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 229910021389 graphene Inorganic materials 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 230000037431 insertion Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 101100129496 Arabidopsis thaliana CYP711A1 gene Proteins 0.000 description 1
- 101100083446 Danio rerio plekhh1 gene Proteins 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
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
技术领域
本申请属于地球物理勘探处理技术领域,尤其涉及一种地层品质因子的计算方法及装置。
背景技术
在有吸收的介质中,地震波每传播一个波长的距离,能量损失的程度可认为是固定的,它与介质的一种物理性质有关,这种物理性质可以称为品质因子,可以用符号Q表示。品质因子Q反映了地下介质对地震波吸收衰减的强弱,是油气检测和储层描述的重要参数。
现有技术中,通常可以使用谱比法进行品质因子的计算。现有技术在使用谱比法计算品质因子时,通常需要求取地震子波振幅谱,但是实际的地震数据,地震子波是未知的,地震子波的求取过程需要复杂的数学计算。如:采用多项式拟合法谱模拟法求取地震子波振幅谱,预设的多项式限制了子波振幅谱的形态,拟合出来的地震子波振幅谱,始终是具有雷克子波特征的单峰的光滑曲线,拟合参数及效果不易控制,阶数低时拟合误差大、拟合效果差,阶数高则计算容易溢出。可以看出,现有技术中计算品质因子的方法,在求取地震子波时处理过程繁琐、计算结果精度低,在生产中不能被广泛应用。有些情况下,品质因子场的建立只能靠处理员凭经验填写品质因子或通过反复试验估计,工作量大,应用结果也不稳定。
因此,提供一种减少品质因子计算时的数据处理的复杂程度,进一步提高品质因子计算结果的准确性的技术方案,是亟待解决的技术问题。
发明内容
本申请目的在于提供一种地层品质因子的计算方法及装置,减少了现有技术中谱比法估算品质因子时,繁琐的参数设置以及复杂的计算过程,提高了地层品质因子计算结果的准确性。
一方面本申请提供了一种地层品质因子的计算方法,包括:
根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱;
对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值;;
根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子。
进一步地,所述方法的另一个实施例中,所述对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,包括:
计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽;
根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率。
进一步地,所述方法的另一个实施例中,所述计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽,包括:
利用第三低通滤波器对所述第一振幅谱进行平滑处理,获得第一光滑振幅谱;
根据所述第一光滑振幅谱计算所述第一频宽;
利用第四低通滤波器对所述第二振幅谱进行平滑处理,获得第二光滑振幅谱;
根据所述第二光滑振幅谱计算所述第二频宽。
进一步地,所述方法的另一个实施例中,所述根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率,包括:
分别设置频宽阈值一对应的第一截止频率,和频宽阈值二对应的第二截止频率,所述频宽阈值一小于所述频宽阈值二,利用如下公式计算所述第一低通滤波器的截止频率,和所述第二低通滤波器的截止频率:
式中,FEX1表示所述第一低通滤波器的截止频率,FEX2表示所述第二低通滤波器的截止频率,f1表示所述第一截止频率,f2表示所述第二截止频率,BW1表示所述第一频宽,BW2表示所述第二频宽,B1表示所述频宽阈值一,B2表示所述频宽阈值二。
进一步地,所述方法的另一个实施例中,所述根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率,包括:
若所述第一频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第二低通滤波器的截止频率;
若所述第一频宽大于所述频宽阈值二,则将所述第二截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽大于所述频宽阈值一,则将所述第二截止频率作为所述第一低通滤波器的截止频率。
进一步地,所述方法的另一个实施例中,所述根据所述第一估计值和所述第二估计值,采用谱比法计算地层的品质因子,包括:
根据所述第一估计值和所述第二估计值,计算所述第一时刻和所述第二时刻对应的地震子波的振幅谱比;
根据所述振幅谱比获得地震子波的振幅谱比的自然对数曲线;
根据所述第一估计值或所述第二估计值,确定所述自然对数曲线中的有效范围;
在所述有效范围内对所述自然对数曲线进行最小平方直线拟合,获得拟合后直线的斜率;
根据所述斜率,计算获得所述地层品质因子。
进一步地,所述方法的另一个实施例中,所述根据所述第一估计值或所述第二估计值,确定所述自然对数曲线中的有效范围,包括:
获取所述第一估计值对应的第一最大值,或者获取所述第二估计值对应的第二最大值;
在所述自然对数曲线中以所述第一最大值对应的第一频率或以所述第二最大值对应的第二频率作为起始点,分别向低频端和高频端搜索,确定低频起始点和高频终止点;
将所述低频起始点和所述高频终止点之间的范围作为所述有效范围。
进一步地,所述方法的另一个实施例中,所述分别向低频端和高频端搜索,确定低频起始点和高频终止点,包括:
在所述自然对数曲线中以所述第一频率或所述第二频率作为起始点,向低频端进行搜索,当高频点对应的函数值大于或等于相邻的低频点对应的函数值时,将所述高频点对应的频率作为所述低频起始点;
在所述自然对数曲线中以所述第一最大值或所述第二最大值作为起始点,向高频端进行搜索,当高频点对应的函数值大于或等于相邻低频点对应的函数值时,将所述低频点对应的频率作为所述高频终止点。
进一步地,所述方法的另一个实施例中,所述将所述低频起始点和所述高频终止点之间的范围作为所述有效范围,包括:
采用如下公式调整所述低频起始点和高频终止点:
式中,F'1表示调整后的低频起始点,F'2表示调整后的高频终止点,FMAX表示所述第一频率或所述第二频率,F1表示所述低频起始点,F2表示所述高频终止点,X1、X2表示调整系数,X1、X2是(0,1)之间的实数;
将调整后的低频起始点和调整后的高频终止点之间的范围作为所述有效范围。
进一步地,所述方法的另一个实施例中,所述方法还包括:
根据不同时间段、不同位置的地震数据计算所述地层品质因子;
根据不同时间段、不同位置对应的地层品质因子,构建品质因子模型。
另一方面,本申请提供了一种地层品质因子的计算装置,包括:
振幅谱获取模块,用于根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱;
子波振幅谱估算模块,用于对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值;
品质因子计算模块,用于根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子。
进一步地,所述装置的另一实施例中,所述子波振幅谱估算模块包括:
频宽计算模块,用于计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽;
滤波器设置模块,用于根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率。
进一步地,所述装置的另一个实施例中,所述频宽计算模块包括:
第一光滑振幅谱获取单元,用于利用第三低通滤波器对所述第一振幅谱进行平滑处理,获得第一光滑振幅谱;
第一频宽获取单元,用于根据所述第一光滑振幅谱计算所述第一频宽;
第二光滑振幅谱获取单元,用于利用第四低通滤波器对所述第二振幅谱进行平滑处理,获得第二光滑振幅谱;
第二频宽获取单元,用于根据所述第二光滑振幅谱计算所述第二频宽。
进一步地,所述装置的另一个实施例中,所述滤波器设置模块具体用于:
分别设置频宽阈值一对应的第一截止频率,和频宽阈值二对应的第二截止频率,所述频宽阈值一小于所述频宽阈值二,利用如下公式计算所述第一低通滤波器的截止频率,和所述第二低通滤波器的截止频率:
式中,FEX1表示所述第一低通滤波器的截止频率,FEX2表示所述第二低通滤波器的截止频率,f1表示所述第一截止频率,f2表示所述第二截止频率,BW1表示所述第一频宽,BW2表示所述第二频宽,B1表示所述频宽阈值一,B2表示所述频宽阈值二。
进一步地,所述装置的另一个实施例中,所述滤波器设置模块还用于:
若所述第一频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第二低通滤波器的截止频率;
若所述第一频宽大于所述频宽阈值二,则将所述第二截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽大于所述频宽阈值一,则将所述第二截止频率作为所述第一低通滤波器的截止频率。
进一步地,所述装置的另一个实施例中,所述品质因子计算模块包括:
振幅谱比计算单元,用于根据所述第一估计值和所述第二估计值,计算所述第一时刻和所述第二时刻对应的地震子波的振幅谱比;
对数处理单元,用于根据所述振幅谱比获得地震子波的振幅谱比的自然对数曲线;
有效范围获取单元,用于根据所述第一估计值或所述第二估计值,确定所述自然对数曲线中的有效范围;
直线拟合单元,用于在所述有效范围内对所述自然对数曲线进行最小平方直线拟合,获得拟合后直线的斜率;
品质因子计算单元,用于根据所述斜率,计算获得所述地层品质因子。
进一步地,所述装置的另一个实施例中,所述有效范围获取单元包括:
最大值获取子单元,用于获取所述第一估计值对应的第一最大值,或者获取所述第二估计值对应的第二最大值;
有效范围搜索子单元,用于在所述自然对数曲线中以所述第一最大值对应的第一频率或以所述第二最大值对应的第二频率作为起始点,分别向低频端和高频端搜索,确定低频起始点和高频终止点;
有效范围确定子单元,用于将所述低频起始点和所述高频终止点之间的范围作为所述有效范围。
进一步地,所述装置的另一个实施例中,所述有效范围搜索子单元具体用于:
在所述自然对数曲线中以所述第一频率或所述第二频率作为起始点,向低频端进行搜索,当高频点对应的函数值大于或等于相邻的低频点对应的函数值时,将所述高频点对应的频率作为所述低频起始点;
在所述自然对数曲线中以所述第一最大值或所述第二最大值作为起始点,向高频端进行搜索,当高频点对应的函数值大于或等于相邻低频点对应的函数值时,将所述低频点对应的频率作为所述高频终止点。
进一步地,所述装置的另一个实施例中,所述有效范围确定子单元具体用于:
采用如下公式调整所述低频起始点和高频终止点:
式中,F'1表示调整后的低频起始点,F'2表示调整后的高频终止点,FMAX表示所述第一频率或所述第二频率,F1表示所述低频起始点,F2表示所述高频终止点,X1、X2表示调整系数,X1、X2是(0,1)之间的实数;
将调整后的低频起始点和调整后的高频终止点之间的范围作为所述有效范围。
进一步地,所述装置的另一个实施例中,所述装置还包括品质因子模型构建模块,所述品质因子模型构建模块具体用于:
根据不同时间段、不同位置的地震数据计算所述地层品质因子;
根据不同时间段、不同位置对应的地层品质因子,构建品质因子模型。
再一方面,本申请还提供了一种地层品质因子的计算装置,包括处理器及用于存储处理器可执行指令的存储器,所述指令被所述处理器执行时实现包括以下步骤:
根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱;
对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值;
根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子。
又一方面,本申请还提供了一种地层品质因子的自动估算系统,包括:处理器以及用于存储处理器可执行指令的存储器,所述处理器执行所述指令时实现上述地层品质因子的计算方法。
本申请提供的地层品质因子的计算方法及装置,利用地震数据,获得地震数据的振幅谱,对获得振幅谱进行低通滤波处理,估算出地震子波振幅谱,为后续地层品质因子的计算提供了准确的数据基础,不需要提取地震子波,减少了计算量。在估算出的地震子波振幅谱即第一估计值和第二估计值的基础上,利用谱比法,计算出地层品质因子,减少了现有技术中谱比法估算品质因子时,繁琐的参数设置以及复杂的计算过程,提高了地层品质因子计算结果的准确性。同时,利用低通滤波器对地震数据自动估算地层品质因子Q值,自动化程度高,节省人力,减小人为因素对估算结果的干扰,高效实用。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请提供的一种地层品质因子的计算方法一个实施例的方法流程示意图;
图2是本申请实施例中地震记录振幅谱以及估算出的地震子波振幅谱示意图;
图3是本申请实施例中根据地震记录不同时窗估算获得的地震子波振幅谱;
图4是本申请一个实施例中的振幅谱比的自然对数曲线;
图5是本申请提供的地层品质因子的计算装置一个实施例的模块结构示意图;
图6是本申请一个实施中子波振幅谱估算模块的结构示意图;
图7是本申请一个实施中频宽计算模块的结构示意图;
图8是本申请一个实施例中品质因子计算模块的结构示意图;
图9是本申请提供的一种地层品质因子的计算系统的实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
地层品质因子(以下也可以简称为Q)可以反映地震波衰减,因此地层品质因子Q可以有效的作用于描述地下构造分布和预测油气储层。获取地震波传播介质的品质因子即Q值在地震数据处理中也具有重要意义,如品质因子Q补偿是高分辨率处理的最常用也是最有效的手段,估算地震数据较准确的Q值是该高分辨率处理的前提条件。
通常可以利用谱比法计算地层品质因子的计算方法,本申请实施例采用低通滤波器从地震数据中估算较准确的地震子波振幅谱,而不必先提取地震子波,在谱比法的基础上,计算出地层品质因子。整个过程无繁琐参数设置、可全自动进行,无需处理员人工干预,使谱比法计算品质因子成为一种在生产中可广泛推广应用的实用技术,并且估算的品质因子可直接用于品质因子Q补偿等高分辨率处理。
具体地,图1是本申请提供的一种地层品质因子的计算方法一个实施例的方法流程示意图,本申请提供的地层品质因子的计算方法包括:
S1、根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱。
具体地,地震数据可以通过地震资料、测井资料等获得。可以以第一时刻为中心,设置一定的时窗长度,选取以第一时刻为中心设置的时窗长度内的地震数据,将选取到的地震数据作傅里叶变换,获得第一时刻对应的第一振幅谱。同样的,以第二时刻为中心,设置一定的时窗长度,选取以第二时刻为中心设置的时窗长度内的地震数据,将选取到的地震数据作傅里叶变换,获得第二时刻对应的第二振幅谱。在选取地震数据时,时窗长度可以根据实际需要设置,在选取第一时刻对应的地震数据和第二时刻对应的地震数据时设置的时窗长度可以相同,也可以不同。另外在计算第一振幅谱和第二振幅谱时,可以利用单道地震数据,也可利用多道地震数据分别计算然后累加以提高数据的信噪比。
S2、对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值。
可以对获得的第一时刻对应的第一振幅谱,和第二时刻对应的第二振幅谱分别进行低通滤波处理。具体可以使用低通滤波器对第一振幅谱进行低通滤波处理,获得地震子波的振幅谱的第一估计值,对第二振幅谱进行低通滤波处理,获得地震子波的振幅谱的第二估计值。不需要提取地震子波,直接根据地震数据的振幅谱,估计出地震子波的振幅谱,减少了计算量。
图2是本申请实施例中地震记录振幅谱以及估算出的地震子波振幅谱示意图,如图2所示,地震子波的振幅谱的频率范围是从0Hz开始至奈奎斯特频率,如果采用现有技术中的多项式拟合法估算一般是不能取得这么宽的频带(需要设置拟合频率范围,范围大则拟合结果误差大)。
S3、根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子。
获得第一时刻和第二时刻分别对应的地震子波的振幅谱的第一估计值和第二估计值后,利用第一估计值和第二估计值,可以利用谱比法计算出地层品质因子。谱比法计算地层品质因子的原理可以包括:
设地震波在传播过程中、两个不同时间t1和t2的地震子波振幅谱分别为A1(f)、A2(f),在t1到t2之间介质的地层品质因子为Q,根据频谱比法原理有:
D(f)=ln A2(f)-ln A1(f)=-πfQ-1(t2-t1) (1)
从上式可以看出D(f)与f的关系是一条直线,其斜率为-πQ-1(t2-t1)。
由此可得:
上式中,ΔD(f)/Δf可以表示D(f)对频率f的斜率,D(f)可以表示地震子波的振幅谱比的对数,f可以表示频率,π可以表示圆周率可以取3.14。
第一估计值和第二估计值可以表示上式中的A1(f)、A2(f),因此,获得第一时刻和第二时刻对应的地震子波的第一估计值和第二估计值后,可以用上述公式(1)、(2)计算获得地层品质因子Q。
获得地层品质因子后,可以利用地层品质因子对地震数据进行处理,可以用地层品质因子描述地下构造分布和预测油气储层等。
本申请一个实施例中,还可以利用上述实施例的方法根据不同时刻、不同位置的地震数据,计算获得不同时刻、不同位置对应的地层品质因子。根据不同时刻、不同位置对应的地层品质因子,可以获得地层品质因子随时间和位置的变化关系,构建出品质因子模型。利用品质因子模型,可以准确的描述地下构建发布、预测油气储层等。
图3是本申请实施例中根据地震记录不同时窗估算获得的地震子波振幅谱,图3表示根据实际地震数据,从浅层至深层取不同的时间点,本申请实施例自动估算所得的地震子波的振幅谱。如图3所示,可以看出从浅至深,地震子波振幅谱的频宽逐渐变窄,能量减弱,且频率向低频端移动,这与地震波传播理论是相符的,也与品质因子的表征相符。利用不同时刻、不同位置对应的地震子波的振幅谱,可以采用谱比法计算出对应的地层品质因子,构建出品质因子模型。
本申请实施例提供的地层品质因子的计算方法,利用地震数据,获得地震数据的振幅谱,对获得振幅谱进行低通滤波处理,估算出地震子波的振幅谱,为后续地层品质因子的计算提供了准确的数据基础,不需要提取地震子波,减少了计算量。在估算出的地震子波的振幅谱即第一估计值和第二估计值的基础上,利用谱比法,计算出地层品质因子,减少了现有技术中谱比法估算品质因子时,繁琐的参数设置以及复杂的计算过程,提高了地层品质因子计算结果的准确性。同时,利用低通滤波器对地震数据自动估算地层品质因子Q值,自动化程度高,节省人力,减小人为因素对估算结果的干扰,高效实用。
在上述实施例的基础上,所述对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,可以包括:
计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽;
根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率。
具体地,在对第一振幅谱和第二振幅谱进行低通滤波处理时,可以先计算出第一振幅谱对应的第一频宽和第二振幅谱对应的第二频宽,根据频宽分别设置第一振幅谱和第二振幅谱对应的低通滤波器的截止频率,使用设置好截止频率的低通滤波器分别对第一振幅谱和第二振幅谱进行低通滤波处理。
第一频宽可以是第一振幅谱中振幅值大于一定阈值的频率范围的宽度,第二频宽可以是第二振幅谱中振幅值大于一定阈值的频率范围的宽度。阈值的具体选取可以根据实际需要设置,一般可以选择在-10dB~-20dB之间。
在上述实施例的基础上,本申请一个实施例中,所述计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽,可以包括:
利用第三低通滤波器对所述第一振幅谱进行平滑处理,获得第一光滑振幅谱;
根据所述第一光滑振幅谱计算所述第一频宽;
利用第四低通滤波器对所述第二振幅谱进行平滑处理,获得第二光滑振幅谱;
根据所述第二光滑振幅谱计算所述第二频宽。
具体地,在根据地震数据获得第一时刻和第二时刻分别对应的第一振幅谱和第二振幅谱后,本申请一个实施例可以使用第三低通滤波器对第一振幅谱进行滤波平滑处理,使用第四低通滤波器对第二振幅谱进行滤波平滑处理。去除第一振幅谱和第二振幅谱中的噪声,提高振幅谱的信噪比,分别获得比较光滑的曲线第一光滑振幅谱和第二光滑振幅谱。第三低通滤波器和第四低通滤波器的截止频率可以根据实际需要进行设置,例如:可以以第一振幅谱和第二振幅谱分别对应的地震子波的时间域长度的作为参考,分别设置第三低通滤波器和第四低通滤波器的截止频率。一般第三低通滤波器和第四低通滤波器的截止频率可以取较大值,如可以取160ms,当然根据实际使用的需要,还可以设置为其他数值。
对第一振幅谱和第二振幅谱进行低通滤波后,获得第一光滑振幅谱和第二光滑振幅谱后,利用第一光滑振幅谱和第二光滑振幅谱分别获得第一频宽和第二频宽。关于频宽的计算,可以取第一光滑振幅谱中振幅值高于-14dB(约为振幅最大值的0.2倍,可在-10dB~-20dB之间选取)的频率范围的宽度,作为第一频宽。同样的,可以取第二光滑振幅谱中振幅值高于-14dB(约为振幅最大值的0.2倍,可在-10dB~-20dB之间选取)的频率范围的宽度,作为第二频宽。在计算频宽之前,对获得第一振幅谱和第二振幅谱进行低通滤波处理,可以提高第一频宽和第二频宽计算的准确性,为地震子波振幅谱的低通滤波器的设计做准备,为后续地层品质因子的计算提供了准确的数据基础。
计算出第一振幅谱对应的第一频宽和第二振幅谱计算的第二频宽后,可以根据第一频宽设置第一低通滤波器的截止频率,根据第二频宽设置第二低通滤波器的截止频率。第一低通滤波器和第二低通滤波器的设置原则可以是:频宽大,则截止频率小,频宽小,则截止频率大,一般截止频率的大小可以设置在30~160ms。具体可以根据实际需要,可以通过历史数据进行模拟实验,预先获得不同的频宽对应的截止频率的大小。获得第一频宽和第二频宽后,对应的设置第一低通滤波器的截止频率和第二低通滤波器的截止频率。也可以,预先通过历史数据进行模拟实验,获得频宽和截止频率之间的函数关系,根据获得的函数关系,设置第一低通滤波器和第二低通滤波器对应的截止频率。
设置好第一低通滤波器和第二低通滤波器的截止频率后,可以利用第一低通滤波器对第一时刻对应的第一振幅谱进行滤波处理,获得第一时刻对应的地震子波的振幅谱的第一估计值。利用第二低通滤波器对第二时刻对应的第二振幅谱进行滤波处理,获得第二时刻对应的地震子波的振幅谱的第二估计值。
需要说明的是,本申请实施例中的第一低通滤波器、第二低通滤波器、第三低通滤波器、第四低通滤波器可以是一个低通滤波器,在进行不同步骤的滤波处理时,设置不同的截止频率。也可以使用不同的低通滤波器,对应于不同振幅谱的滤波处理,可以根据实际需要进行选择。
本申请实施例提供的地层品质因子的计算方法,根据第一振幅谱和第二振幅谱的频宽设计低通滤波器的截止频率,利用低通滤波器估算出地震子波振幅谱,再利用谱比法计算出地层品质因子。不需要提取地震子波,降低了地层品质因子计算的复杂程度,提高了地层品质因子计算的准确性和适用性。
在上述实施例的基础上,本申请一个实施例中,所述根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率,可以包括:
分别设置频宽阈值一对应的第一截止频率,和频宽阈值二对应的第二截止频率,所述频宽阈值一小于所述频宽阈值二,利用如下公式计算所述第一低通滤波器的截止频率,和所述第二低通滤波器的截止频率:
式中,FEX1可以表示第一低通滤波器的截止频率,FEX2可以表示第二低通滤波器的截止频率,f1可以表示第一截止频率,f2可以表示第二截止频率,BW1可以表示第一频宽,BW2可以表示第二频宽,B1可以表示频宽阈值一,B2可以表示频宽阈值二。
具体地,可以根据实际需要,分别设置频宽阈值一对应的第一截止频率,和频宽阈值二对应的第二截止频率。频宽阈值一可以表示较低的频宽,频宽阈值二可以表示较高的频宽,频宽阈值一、频宽阈值二以及对应的第一截止频率、第二截止频率可以通过实验或根据历史数据获得,具体数值本申请实施例不作具体限定。本申请一个实施例中,可以将频宽阈值一可以设置为25Hz,对应的第一截止频率可以设置为0.160s,频宽阈值二可以设置为135Hz,对应的第二截止频率可以设置为0.03s。在获得第一振幅谱或第一光滑振幅谱对应的第一频宽,以及第二振幅谱或第二光滑振幅谱对应的第二频宽后,可以利用上述公式(3),分别计算出第一低通滤波器对应的截止频率和第二低通滤波器对应的截止频率。
在上述实施例的基础上,本申请一个实施例中,所述根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率,可以包括:
若所述第一频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第二低通滤波器的截止频率;
若所述第一频宽大于所述频宽阈值二,则将所述第二截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽大于所述频宽阈值一,则将所述第二截止频率作为所述第一低通滤波器的截止频率。
具体地,若获得的第一频宽或第二频宽不在频宽阈值一和频宽阈值二的范围内时,如:若第一频宽小于频宽阈值一,或第二频宽小于频宽阈值一,则可以将频宽阈值一作为第一频宽或第二频宽,对应的可以将频宽阈值一对应的第一截止频率作为第一低通滤波器的截止频率,或可以将频宽阈值一对应的第一截止频率作为第二低通滤波器的截止频率。若第一频宽大于频宽阈值二,或第二频宽大于频宽阈值二,则可以将频宽阈值二作为第一频宽或第二频宽,对应的可以将频宽阈值二对应的第二截止频率作为第一低通滤波器的截止频率,或可以将频宽阈值二对应的第二截止频率作为第二低通滤波器的截止频率。
本申请实施例,根据第一振幅谱(或第一光滑振幅谱)、第二振幅谱(或第二光滑振幅谱)分别对应的第一频宽、第二频宽,合理设计第一低通滤波器和第二低通滤波器的截止频率,通过第一低通滤波器和第二低通滤波器估算出地震子波的振幅谱,为后续地层品质因子的估算提供了准确的数据基础。
在上述实施例的基础上,本申请一个实施例中,所述根据所述第一估计值和所述第二估计值,采用谱比法计算地层的品质因子,可以包括:
根据所述第一估计值和所述第二估计值,计算所述第一时刻和所述第二时刻对应的地震子波的振幅谱比;
根据所述振幅谱比获得地震子波的振幅谱比的自然对数曲线;
根据所述第一估计值或所述第二估计值,确定所述自然对数曲线中的有效范围;
在所述有效范围内对所述自然对数曲线进行最小平方直线拟合,获得拟合后直线的斜率;
根据所述斜率,计算获得所述地层品质因子。
具体地,本申请实施例在估算出地震子波的振幅谱后,利用谱比法进行地层品质因子的计算。可以先根据第一估计值和第二估计值计算出地震子波的振幅谱比,具体可以将第一估计值和第二估计值进行相除处理,获得地震子波的振幅谱比。将获得的振幅谱比进行对数处理,获得地震子波的振幅谱比的自然对数曲线。根据第一估计值或第二估计值,确定自然对数曲线中的有效范围,有效范围可以表示自然对数曲线中能够反映振幅谱比和频率之间的呈直线关系的范围。从公式(1)、(2)可以看出,要计算地层品质因子,需要获得地震子的振幅谱比的对数值与频率之间直线的斜率。本申请实施例利用第一估计值和第二估计值在振幅谱比的自然对数曲线中获取相对呈直线的有效范围,在有效范围内对自然对数曲线进行最小平方直线拟合,获得拟合后直线的斜率,再利用获得的斜率,利用上述公式(2)计算出地层品质因子。
本申请实施例,在估算出地震子波的振幅谱后,利用谱比法,拟合获得地震子波的振幅谱比与频率之间的直线,根据获得的该直线的斜率,计算获得地层品质因子。不需要提取地震子波,并且整个过程可以全自动进行,不需要处理员的人工干预,提高了地层品质因子计算结果的准确性,实现了地层品质因子计算的自动化处理。
在上述实施例的基础上,本申请一个实施例中,所述根据所述第一估计值或所述第二估计值,确定所述自然对数曲线中的有效范围,可以包括:
获取所述第一估计值对应的第一最大值,或者获取所述第二估计值对应的第二最大值;
在所述自然对数曲线中以所述第一最大值对应的第一频率或以所述第二最大值对应的第二频率作为起始点,分别向低频端和高频端搜索,确定低频起始点和高频终止点;
将所述低频起始点和所述高频终止点之间的范围作为所述有效范围。
具体地,在确定自然对数曲线的有效范围时,可以选取第一时刻对应的地震子波的振幅谱(即第一估计值)的最大值作为第一最大值,并获取第一最大值对应的第一频率。或者,也可以选取第二时刻对应的地震子波的振幅谱(即第二估计值)的最大值作为第二最大值,并获取第二最大值对应的第二频率。在获得的地震子波的自然对数曲线中,以第一频率或第二频率为起点,分别向低频端和高频端搜索,获得有效范围的低频起始点和高频终止点,确定出有效范围。低频端可以表示小于第一频率(或第二频率)的自然对数曲线的一端,高频端可以表示大于第一频率(或第二频率)的自然对数曲线的一端。低频起始点可以表示低频端对应的有效范围的起始点,高频终止点可以表示高频端对应的有效范围的终止点。低频起始点和高频终止点的相邻点的导数以及相邻点的自然对数曲线的值的大小会发生变化,可以通过依次搜索,当搜索到相邻点的导数或自然对数曲线的取值的大小关系发生变化时,确定出低频起始点和高频终止点。
本申请一个实施例中,所述分别向低频端和高频端搜索,确定低频起始点和高频终止点,可以包括:
在所述自然对数曲线中以所述第一最大值或所述第二最大值作为起始点,向低频端进行搜索,当高频点对应的函数值大于或等于相邻的低频点对应的函数值时,将所述高频点对应的频率作为所述低频起始点;
在所述自然对数曲线中以所述第一频率或所述第二频率作为起始点,向高频端进行搜索,当高频点对应的函数值大于或等于相邻低频点对应的函数值时,将所述低频点对应的频率作为所述高频终止点。
具体地,可以以第一频率或第二频率作为起点向低频端搜索,当Di≥Di-1(i表示频率,Di可以表示频率i对应的振幅谱比值)时停止,将Di对应的频率作为低频起始点;以第一频率或第二频率作为起点向高频端搜索,当Di≤Di+1时停止,可以将Di对应频率作为高频终止点。
采用自动搜索的方式,确定自然对数曲线中的有效范围,进一步可以获得自然对数曲线中的直线段,为后续的地层品质因子的计算提供了准确的数据基础。
在上述实施例的基础上,本申请一个实施例中,所述将所述低频起始点和所述高频终止点之间的范围作为所述有效范围,可以包括:
采用如下公式调整所述低频起始点和高频终止点:
式中,F'1可以表示调整后的低频起始点,F'2可以表示调整后的高频终止点,FMAX可以表示第一频率或第二频率,F1可以表示低频起始点,F2可以表示高频终止点,X1、X2可以表示调整系数,X1、X2是(0,1)之间的实数;
将调整后的低频起始点和调整后的高频终止点之间的范围作为所述有效范围。
在获得低频起始点和高频终止点,确定出有效范围后,可以通过上述公式(3)调整低频起始点和高频终止点的大小,适当的缩小有效范围的大小,减小了有效范围的两端的低信噪比带来的误差,提高了后续直线拟合的准确性,进一步提高了地层品质因子计算的准确性。
图4是本申请一个实施例中的振幅谱比的自然对数曲线,本申请实施例中的有效范围可以在地震记录主频带范围附近,因此可以选取地震记录主频带范围附近对自然对数曲线D(f)做直线拟合。如图3所示,本申请一个实施例中地震记录的主频为20Hz左右,可以对图4中自然对数曲线在15Hz到25Hz的范围做直线拟合,求得直线斜率后即可根据公式(2)计算地层品质因子Q值。当然,具体拟合的有效范围,也可以利用上述实施例的方法计算获得。
需要说明的是,本申请中的各个公式只是一种实施例,还可以根据需要对公式进行变形或变换,本申请不作具体限定。
下面结合具体的示例,介绍本申请实施例的方案:
(1)输入地震数据,设置地层品质因子Q值估算的第一时刻t1、第二时刻t2,以及时窗长度L(200ms~1000ms)。
(2)计算以第一时刻t1为中心,时窗长度L的地震数据的第一振幅谱A1,具体计算方法可以参考上述实施例的介绍,此处不再赘述。
(3)构建以第一振幅谱A1的频率间隔DF1为采样间隔,截止频率为FE1(取值0.120s,0.08s~0.160s可选)的低通滤波器,对第一振幅谱A1做滤波处理,得到第一光滑振幅谱
(4)计算第一光滑振幅谱的频宽:可以以第一光滑振幅谱的振幅最大值的频率点为起始,分别向低频端和高频端搜索,幅值在MAX1(取值-14dB,-20dB~-10dB可选)之上的频率宽度,即为第一频宽BW1。
(5)根据第一频宽BW1设计第一低通滤波器的截止频率,基本原则是:频宽大则截频小,频宽小则截频大。具体实施方式可以如下:
设置频宽-截止频率的基准值,可由用户根据地震数据情况预设,也可设置缺省值如下:频宽25Hz对应截止频率0.160s,频宽135Hz对应截止频率0.03s,频宽超过25Hz~135Hz的,按25Hz或135Hz处理,第一频宽BW1对应截止频率如下:
FEX1=0.160+(0.03-0.160)×(BW1-25)/(135-25)
(6)对第一振幅谱A1用步骤(5)设计好的第一低通滤波器处理,得到第一时刻t1的地震子波的振幅谱的第一估计值AW1。步骤(3)和步骤(5)中两次低通滤波器设计时可以加斜坡处理,以减弱滤波处理的Gibbs(吉布斯效应)效应。
Gibbs效应可以包括:将具有不连续点的周期函数(如矩形脉冲)进行傅立叶级数展开后,选取有限项进行合成。当选取的项数越多,在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时,该峰起值趋于一个常数,大约等于总跳变值的9%。
(7)对第二时刻t2重复步骤(2)~(6),得第二时刻t2的地震子波的振幅谱的第二估计值AW2。
(8)计算第一时刻t1、第二时刻t2的地震子波的振幅谱比、以及它的自然对数曲线D(f):
D(f)=ln(AW2/AW1)
(9)搜索第一时刻t1地震子的波振幅谱的第一估计值AW1的最大值,其频率为FMAX1;在自然对数曲线D(f)中,以频率FMAX1为起始点,前后搜索确定自然对数曲线的有效范围(F1~F2):
从频率点FMAX1向低频端搜索,当Di≥Di-1(i为频率下标)时停止,Di对应频率为F1;从频率点FMAX1向高频端搜索,当Di≤Di+1时停止,Di对应频率为F2。
当然,根据需要也可以选取第二时刻的地震子波的振幅谱的第二估计值AW2的最大值,以及最大值对应的频率,进行搜索,确定有效范围。
为减小频率两端由于低信噪比带来的误差,可适当缩小(F1~F2)范围,具体可以参考上述公式(4),式中X1、X2为(0,1)之间的实数,通常可取X1=X2=0.8。
(10)在有效范围内(F1~F2)对自然对数曲线做最小平方直线拟合,计算直线斜率a。
(11)根据谱比法,利用上述公式(2)计算地层品质因子Q值:Q=-π(t2-t1)a-1
(12)重复以上步骤(1)~(11),计算其它时刻、其他位置的地层品质因子Q值,构建品质因子模型。
本申请提供的地层品质因子的计算方法,在满足地震数据处理技术的常规假设条件下即可实现,因而方法原理是正确的,无须验证;利用地震记录振幅谱、采用低通滤波器估算子波振幅谱在实际应用中可以验证是可靠的,从而确保地层品质因子Q值估算的正确性。利用低通滤波器对地震数据自动估算地层品质因子Q值,自动化程度高,节省人力,减小人为因素对估算结果的干扰,高效实用。
基于上述所述的地层品质因子的计算方法,本说明书一个或多个实施例还提供一种地层品质因子的计算装置。所述的装置可以包括使用了本说明书实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思,本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似,因此本说明书实施例具体的装置的实施可以参见前述方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
具体地,图5是本申请提供的地层品质因子的计算装置一个实施例的模块结构示意图,如图5所示,本申请中提供的地层品质因子的计算装置包括:振幅谱获取模块51、子波振幅谱估算模块52、品质因子计算模块53。
振幅谱获取模块51,可以用于根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱;
子波振幅谱估算模块52,可以用于对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值;品质因子计算模块53,可以用于根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子。
本申请提供的地层品质因子的计算装置,利用地震数据,获得地震数据的振幅谱,对获得振幅谱进行低通滤波处理,估算出地震子波振幅谱,为后续地层品质因子的计算提供了准确的数据基础,不需要提取地震子波,减少了计算量。在估算出的地震子波振幅谱即第一估计值和第二估计值的基础上,利用谱比法,计算出地层品质因子,减少了现有技术中谱比法估算品质因子时,繁琐的参数设置以及复杂的计算过程,提高了地层品质因子计算结果的准确性。同时,利用低通滤波器对地震数据自动估算地层品质因子Q值,自动化程度高,节省人力,减小人为因素对估算结果的干扰,高效实用。
图6是本申请一个实施中子波振幅谱估算模块的结构示意图,如图6所示,在上述实施例的基础上,所述子波振幅谱估算模块52包括:
频宽计算模块61,可以用于计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽;
滤波器设置模块62,可以用于根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率。
本申请提供的地层品质因子的计算装置,利用第一振幅谱、第二振幅谱的频宽,设计低通滤波器,进一步使用低通滤波器对第一振幅谱、第二振幅谱的频宽进行滤波处理,估算出地震子波的振幅谱。为后续地层品质因子的计算提供了准确的数据基础,同时,不需要提取地震子波,降低了数据处理的复杂程度。
图7是本申请一个实施中频宽计算模块的结构示意图,如图7所示,在上述实施例的基础上,所述频宽计算模块61包括:
第一光滑振幅谱获取单元71,可以用于利用第三低通滤波器对所述第一振幅谱进行平滑处理,获得第一光滑振幅谱;
第一频宽获取单元72,可以用于根据所述第一光滑振幅谱计算所述第一频宽;
第二光滑振幅谱获取单元73,可以用于利用第四低通滤波器对所述第二振幅谱进行平滑处理,获得第二光滑振幅谱;
第二频宽获取单元74,可以用于根据所述第二光滑振幅谱计算所述第二频宽。
本申请提供的地层品质因子的计算装置,在计算频宽之前,对获得第一振幅谱和第二振幅谱进行低通滤波处理,可以提频宽阈值二计算的准确性,为地震子波振幅谱的低通滤波器的设计做准备,为后续地层品质因子的计算提供了准确的数据基础。
在上述实施例的基础上,所述滤波器设置模块具体用于:
分别设置频宽阈值一对应的第一截止频率,和频宽阈值二对应的第二截止频率,所述频宽阈值一小于所述频宽阈值二,利用如下公式计算所述第一低通滤波器的截止频率,和所述第二低通滤波器的截止频率:
式中,FEX1表示所述第一低通滤波器的截止频率,FEX2表示所述第二低通滤波器的截止频率,f1表示所述第一截止频率,f2表示所述第二截止频率,BW1表示所述第一频宽,BW2表示所述第二频宽,B1表示所述频宽阈值一,B2表示所述频宽阈值二。
本申请提供的地层品质因子的计算装置,给出了低通滤波器的截止频率的具体设计方法,通过合理的设计低通滤波器,利用低通滤波器对地震数据对应的第一振幅谱和第二振幅谱进行滤波处理,准确的估算出地震子波的振幅谱,为地层品质因子的计算提供了准确的数据基础。
在上述实施例的基础上,所述滤波器设置模块还用于:
若所述第一频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第二低通滤波器的截止频率;
若所述第一频宽大于所述频宽阈值二,则将所述第二截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽大于所述频宽阈值一,则将所述第二截止频率作为所述第一低通滤波器的截止频率。
本申请提供的地层品质因子的计算装置,给出了频宽不在频宽阈值一和频宽阈值二范围内的情况下,低通滤波器的截止频率的设计方法,提高了地层品质因子计算的准确性和适用范围。
图8是本申请一个实施例中品质因子计算模块的结构示意图,如图8所示,在上述实施例的基础上,所述品质因子计算模块53包括:
振幅谱比计算单元81,用于根据所述第一估计值和所述第二估计值,计算所述第一时刻和所述第二时刻对应的地震子波的振幅谱比;
对数处理单元82,可以用于根据所述振幅谱比获得地震子波的振幅谱比的自然对数曲线;
有效范围获取单元83,可以用于根据所述第一估计值或所述第二估计值,确定所述自然对数曲线中的有效范围;
直线拟合单元84,可以用于在所述有效范围内对所述自然对数曲线进行最小平方直线拟合,获得拟合后直线的斜率;
品质因子计算单元85,可以用于根据所述斜率,计算获得所述地层品质因子。
本申请提供的地层品质因子的计算装置,在估算出地震子波的振幅谱后,利用谱比法,拟合获得地震子波的振幅谱比与频率之间的直线,根据获得该直线的斜率,计算获得地层品质因子。不需要提取地震子波,并且整个过程可以全自动进行,不需要处理员的人工干预,提高了地层品质因子计算结果的准确性,实现了地层品质因子计算的自动化处理。
在上述实施例的基础上,所述有效范围获取单元包括:
最大值获取子单元,用于获取所述第一估计值对应的第一最大值,或者获取所述第二估计值对应的第二最大值;
有效范围搜索子单元,用于在所述自然对数曲线中以所述第一最大值对应的第一频率或以所述第二最大值对应的第二频率作为起始点,分别向低频端和高频端搜索,确定低频起始点和高频终止点;
有效范围确定子单元,用于将所述低频起始点和所述高频终止点之间的范围作为所述有效范围。
本申请提供的地层品质因子的计算装置,给出了自然对数曲线拟合时的有效范围的具体获取方法,提高了自然对数曲线拟合的准确性,为地层品质因子的计算提供了准确的数据基础。
在上述实施例的基础上,所述有效范围搜索子单元具体用于:
在所述自然对数曲线中以所述第一频率或所述第二频率作为起始点,向低频端进行搜索,当高频点对应的函数值大于或等于相邻的低频点对应的函数值时,将所述高频点对应的频率作为所述低频起始点;
在所述自然对数曲线中以所述第一最大值或所述第二最大值作为起始点,向高频端进行搜索,当高频点对应的函数值大于或等于相邻低频点对应的函数值时,将所述低频点对应的频率作为所述高频终止点。
本申请提供的地层品质因子的计算装置,采用自动搜索的方式,确定自然对数曲线中的有效范围,进一步可以获得自然对数曲线中的直线段,为后续的地层品质因子的计算提供了准确的数据基础。
在上述实施例的基础上,所述有效范围确定子单元具体用于:
采用如下公式调整所述低频起始点和高频终止点:
式中,F'1表示调整后的低频起始点,F'2表示调整后的高频终止点,FMAX表示所述第一频率或所述第二频率,F1表示所述低频起始点,F2表示所述高频终止点,X1、X2表示调整系数,X1、X2是(0,1)之间的实数;
将调整后的低频起始点和调整后的高频终止点之间的范围作为所述有效范围。
本申请提供的地层品质因子的计算装置,通过调整低频起始点和高频终止点的大小,适当的缩小有效范围的大小,减小有效范围的两端的低信噪比带来的误差,提高后续直线拟合的准确性,进一步提高地层品质因子计算的准确性。
在上述实施例的基础上,所述装置还包括品质因子模型构建模块,所述品质因子模型构建模块具体用于:
根据不同时间段、不同位置的地震数据计算所述地层品质因子;
根据不同时间段、不同位置对应的地层品质因子,构建品质因子模型。
本申请提供的地层品质因子的计算装置,通过获取不同时间不同位置对应的地层品质因子,构建出品质因子模型。利用品质因子模型,可以准确的描述地下构建发布、预测油气储层等。
需要说明的,上述所述的装置根据方法实施例的描述还可以包括其他的实施方式,具体的实现方式可以参照相关方法实施例的描述,在此不作一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上,所述的存储介质可以计算机读取并执行,实现本说明书实施例所描述方案的效果。因此,本说明书还提供一种地层品质因子的计算装置,包括处理器及存储处理器可执行指令的存储器,所述指令被所述处理器执行时实现包括以下步骤:
根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱;
对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值;
根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子。
所述存储介质可以包括用于存储信息的物理装置,通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括:利用电能方式存储信息的装置如,各式存储器,如RAM、ROM等;利用磁能方式存储信息的装置如,硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘;利用光学方式存储信息的装置如,CD或DVD。当然,还有其他方式的可读存储介质,例如量子存储器、石墨烯存储器等等。
需要说明的,上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述,在此不作一一赘述。
上述实施例所述的一种地层品质因子的计算装置,利用低通滤波器对地震数据自动估算地层品质因子Q值,自动化程度高,节省人力,减小人为因素对估算结果的干扰,高效实用。
本说明书实施例提供的上述地层品质因子的计算方法或装置可以在计算机中由处理器执行相应的程序指令来实现,如使用windows操作系统的c++语言在PC端实现、linux系统实现,或其他例如使用android、iOS系统程序设计语言在智能终端实现,以及基于量子计算机的处理逻辑实现等。本说明书提供的一种地层品质因子的计算系统的一个实施例中,图9是本申请提供的一种地层品质因子的计算系统的实施例的模块结构示意图,如图9所示,本申请另一实施例提供的地层品质因子的自动估算系统可以包括处理器91以及用于存储处理器可执行指令的存储器92,
处理器91和存储器92通过总线93完成相互间的通信;
所述处理器91用于调用所述存储器92中的程序指令,以执行上述各地层品质因子的计算方法实施例所提供的方法,例如包括:根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值;根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子。
需要说明的是说明书上述所述的系统根据相关方法实施例的描述还可以包括其他的实施方式,具体的实现方式可以参照方法实施例的描述,在此不作一一赘述。本申请中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于硬件+程序类实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例,仍然可以属于本说明书实施例的可选实施方案范围之内。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的,计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
虽然本说明书一个或多个实施例提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或终端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。第一,第二等词语用来表示名称,而并不表示任何特定的顺序。
为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
本发明是参照根据本发明实施例的方法、装置(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储、石墨烯存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
本领域技术人员应明白,本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此,本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书一个或多个实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本本说明书一个或多个实施例,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书一个或多个实施例的实施例而已,并不用于限制本本说明书一个或多个实施例。对于本领域技术人员来说,本说明书一个或多个实施例可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在权利要求范围之内。
Claims (18)
1.一种地层品质因子的计算方法,其特征在于,所述方法包括:
根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱;
对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值;
根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子;
其中,所述对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,包括:
计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽;
根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率;
所述根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率,包括:
分别设置频宽阈值一对应的第一截止频率,和频宽阈值二对应的第二截止频率,所述频宽阈值一小于所述频宽阈值二,利用如下公式计算所述第一低通滤波器的截止频率,和所述第二低通滤波器的截止频率:
式中,FEX1表示所述第一低通滤波器的截止频率,FEX2表示所述第二低通滤波器的截止频率,f1表示所述第一截止频率,f2表示所述第二截止频率,BW1表示所述第一频宽,BW2表示所述第二频宽,B1表示所述频宽阈值一,B2表示所述频宽阈值二。
2.如权利要求1所述的一种地层品质因子的计算方法,其特征在于,所述计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽,包括:
利用第三低通滤波器对所述第一振幅谱进行平滑处理,获得第一光滑振幅谱;
根据所述第一光滑振幅谱计算所述第一频宽;
利用第四低通滤波器对所述第二振幅谱进行平滑处理,获得第二光滑振幅谱;
根据所述第二光滑振幅谱计算所述第二频宽。
3.如权利要求1所述的一种地层品质因子的计算方法,其特征在于,所述根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率,包括:
若所述第一频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第二低通滤波器的截止频率;
若所述第一频宽大于所述频宽阈值二,则将所述第二截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽大于所述频宽阈值二,则将所述第二截止频率作为所述第一低通滤波器的截止频率。
4.如权利要求1所述的一种地层品质因子的计算方法,其特征在于,所述根据所述第一估计值和所述第二估计值,采用谱比法计算地层的品质因子,包括:
根据所述第一估计值和所述第二估计值,计算所述第一时刻和所述第二时刻对应的地震子波的振幅谱比;
根据所述振幅谱比获得地震子波的振幅谱比的自然对数曲线;
根据所述第一估计值或所述第二估计值,确定所述自然对数曲线中的有效范围;
在所述有效范围内对所述自然对数曲线进行最小平方直线拟合,获得拟合后直线的斜率;
根据所述斜率,计算获得所述地层品质因子。
5.如权利要求4所述的一种地层品质因子的计算方法,其特征在于,所述根据所述第一估计值或所述第二估计值,确定所述自然对数曲线中的有效范围,包括:
获取所述第一估计值对应的最大振幅值作为第一最大值,或者获取所述第二估计值对应的最大振幅值作为第二最大值;
在所述自然对数曲线中以所述第一最大值对应的第一频率或以所述第二最大值对应的第二频率作为起始点,分别向低频端和高频端搜索,确定低频起始点和高频终止点;
将所述低频起始点和所述高频终止点之间的范围作为所述有效范围。
6.如权利要求5所述的一种地层品质因子的计算方法,其特征在于,所述分别向低频端和高频端搜索,确定低频起始点和高频终止点,包括:
在所述自然对数曲线中以所述第一频率或所述第二频率作为起始点,向低频端进行搜索,当高频点对应的函数值大于或等于相邻的低频点对应的函数值时,将所述高频点对应的频率作为所述低频起始点;
在所述自然对数曲线中以所述第一最大值或所述第二最大值作为起始点,向高频端进行搜索,当高频点对应的函数值大于或等于相邻低频点对应的函数值时,将所述低频点对应的频率作为所述高频终止点。
7.如权利要求5所述的一种地层品质因子的计算方法,其特征在于,所述将所述低频起始点和所述高频终止点之间的范围作为所述有效范围,包括:
采用如下公式调整所述低频起始点和高频终止点:
式中,F'1表示调整后的低频起始点,F'2表示调整后的高频终止点,FMAX表示所述第一频率或所述第二频率,F1表示所述低频起始点,F2表示所述高频终止点,X1、X2表示调整系数,X1、X2是(0,1)之间的实数;
将调整后的低频起始点和调整后的高频终止点之间的范围作为所述有效范围。
8.如权利要求1-7任一项所述的一种地层品质因子的计算方法,其特征在于,所述方法还包括:
根据不同时间段、不同位置的地震数据计算所述地层品质因子;
根据不同时间段、不同位置对应的地层品质因子,构建品质因子模型。
9.一种地层品质因子的计算装置,其特征在于,包括:
振幅谱获取模块,用于根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱;
子波振幅谱估算模块,用于对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值;
品质因子计算模块,用于根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子;
其中,所述子波振幅谱估算模块包括:
频宽计算模块,用于计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽;
滤波器设置模块,用于根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率;
所述滤波器设置模块具体用于:
分别设置频宽阈值一对应的第一截止频率,和频宽阈值二对应的第二截止频率,所述频宽阈值一小于所述频宽阈值二,利用如下公式计算所述第一低通滤波器的截止频率,和所述第二低通滤波器的截止频率:
式中,FEX1表示所述第一低通滤波器的截止频率,FEX2表示所述第二低通滤波器的截止频率,f1表示所述第一截止频率,f2表示所述第二截止频率,BW1表示所述第一频宽,BW2表示所述第二频宽,B1表示所述频宽阈值一,B2表示所述频宽阈值二。
10.如权利要求9所述的一种地层品质因子的计算装置,其特征在于,所述频宽计算模块包括:
第一光滑振幅谱获取单元,用于利用第三低通滤波器对所述第一振幅谱进行平滑处理,获得第一光滑振幅谱;
第一频宽获取单元,用于根据所述第一光滑振幅谱计算所述第一频宽;
第二光滑振幅谱获取单元,用于利用第四低通滤波器对所述第二振幅谱进行平滑处理,获得第二光滑振幅谱;
第二频宽获取单元,用于根据所述第二光滑振幅谱计算所述第二频宽。
11.如权利要求9所述的一种地层品质因子的计算装置,其特征在于,所述滤波器设置模块还用于:
若所述第一频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽小于所述频宽阈值一,则将所述第一截止频率作为所述第二低通滤波器的截止频率;
若所述第一频宽大于所述频宽阈值二,则将所述第二截止频率作为所述第一低通滤波器的截止频率;
若所述第二频宽大于所述频宽阈值二,则将所述第二截止频率作为所述第一低通滤波器的截止频率。
12.如权利要求9所述的一种地层品质因子的计算装置,其特征在于,所述品质因子计算模块包括:
振幅谱比计算单元,用于根据所述第一估计值和所述第二估计值,计算所述第一时刻和所述第二时刻对应的地震子波的振幅谱比;
对数处理单元,用于根据所述振幅谱比获得地震子波的振幅谱比的自然对数曲线;
有效范围获取单元,用于根据所述第一估计值或所述第二估计值,确定所述自然对数曲线中的有效范围;
直线拟合单元,用于在所述有效范围内对所述自然对数曲线进行最小平方直线拟合,获得拟合后直线的斜率;
品质因子计算单元,用于根据所述斜率,计算获得所述地层品质因子。
13.如权利要求12所述的一种地层品质因子的计算装置,其特征在于,所述有效范围获取单元包括:
最大值获取子单元,用于获取所述第一估计值对应的最大振幅值作为第一最大值,或者获取所述第二估计值对应的最大振幅值作为第二最大值;
有效范围搜索子单元,用于在所述自然对数曲线中以所述第一最大值对应的第一频率或以所述第二最大值对应的第二频率作为起始点,分别向低频端和高频端搜索,确定低频起始点和高频终止点;
有效范围确定子单元,用于将所述低频起始点和所述高频终止点之间的范围作为所述有效范围。
14.如权利要求13所述的一种地层品质因子的计算装置,其特征在于,所述有效范围搜索子单元具体用于:
在所述自然对数曲线中以所述第一频率或所述第二频率作为起始点,向低频端进行搜索,当高频点对应的函数值大于或等于相邻的低频点对应的函数值时,将所述高频点对应的频率作为所述低频起始点;
在所述自然对数曲线中以所述第一最大值或所述第二最大值作为起始点,向高频端进行搜索,当高频点对应的函数值大于或等于相邻低频点对应的函数值时,将所述低频点对应的频率作为所述高频终止点。
15.如权利要求13所述的一种地层品质因子的计算装置,其特征在于,所述有效范围确定子单元具体用于:
采用如下公式调整所述低频起始点和高频终止点:
式中,F'1表示调整后的低频起始点,F'2表示调整后的高频终止点,FMAX表示所述第一频率或所述第二频率,F1表示所述低频起始点,F2表示所述高频终止点,X1、X2表示调整系数,X1、X2是(0,1)之间的实数;
将调整后的低频起始点和调整后的高频终止点之间的范围作为所述有效范围。
16.如权利要求9-15任一项所述的一种地层品质因子的计算装置,其特征在于,所述装置还包括品质因子模型构建模块,所述品质因子模型构建模块具体用于:
根据不同时间段、不同位置的地震数据计算所述地层品质因子;
根据不同时间段、不同位置对应的地层品质因子,构建品质因子模型。
17.一种地层品质因子的计算装置,其特征在于,包括处理器及用于存储处理器可执行指令的存储器,所述指令被所述处理器执行时实现包括以下步骤:
根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱;
对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,分别获得对应的地震子波振幅谱的第一估计值、第二估计值;
根据所述第一估计值和所述第二估计值,采用谱比法计算地层品质因子;
其中,所述对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理,包括:
计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽;
根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率;
所述根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率,根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率,包括:
分别设置频宽阈值一对应的第一截止频率,和频宽阈值二对应的第二截止频率,所述频宽阈值一小于所述频宽阈值二,利用如下公式计算所述第一低通滤波器的截止频率,和所述第二低通滤波器的截止频率:
式中,FEX1表示所述第一低通滤波器的截止频率,FEX2表示所述第二低通滤波器的截止频率,f1表示所述第一截止频率,f2表示所述第二截止频率,BW1表示所述第一频宽,BW2表示所述第二频宽,B1表示所述频宽阈值一,B2表示所述频宽阈值二。
18.一种地层品质因子的计算系统,其特征在于,包括处理器以及用于存储处理器可执行指令的存储器,所述处理器执行所述指令时实现如权利要求1至8中任意一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810123034.8A CN108469633B (zh) | 2018-02-07 | 2018-02-07 | 一种地层品质因子的计算方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810123034.8A CN108469633B (zh) | 2018-02-07 | 2018-02-07 | 一种地层品质因子的计算方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108469633A CN108469633A (zh) | 2018-08-31 |
CN108469633B true CN108469633B (zh) | 2019-10-11 |
Family
ID=63266279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810123034.8A Active CN108469633B (zh) | 2018-02-07 | 2018-02-07 | 一种地层品质因子的计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108469633B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110515127B (zh) * | 2019-09-26 | 2020-09-18 | 中国石油大学(北京) | 一种地震品质因子确定方法、装置、设备、介质 |
CN114252910B (zh) * | 2020-09-25 | 2024-07-26 | 中国石油天然气股份有限公司 | 地震品质因子确定方法及装置 |
CN112578448B (zh) * | 2020-12-03 | 2021-10-08 | 成都理工大学 | 一种高精度的地层品质因子的提取方法 |
US11762116B2 (en) | 2021-03-18 | 2023-09-19 | Saudi Arabian Oil Company | System and method of hydrocarbon detection using nonlinear model frequency slope |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7376517B2 (en) * | 2005-05-13 | 2008-05-20 | Chevron U.S.A. Inc. | Method for estimation of interval seismic quality factor |
CN103376464B (zh) * | 2012-04-13 | 2016-04-06 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
CN106547019A (zh) * | 2015-09-17 | 2017-03-29 | 中国石油化工股份有限公司 | 一种确定地层品质因子的方法 |
-
2018
- 2018-02-07 CN CN201810123034.8A patent/CN108469633B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108469633A (zh) | 2018-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108469633B (zh) | 一种地层品质因子的计算方法及装置 | |
CN101860344A (zh) | 一种选频滤波器的构建方法及采用该方法实现fir型、iir型滤波器的构建方法 | |
CN107194180B (zh) | 油藏水驱体积波及系数的确定方法、装置及系统 | |
CN109359120A (zh) | 一种模型训练中的数据更新方法、装置及设备 | |
CN106226812A (zh) | 基于时频二次谱提高地震数据分辨率的方法及装置 | |
CN110532295A (zh) | 一种计算机实施的信息处理、信息查询的方法及装置 | |
CN109460613A (zh) | 模型裁剪方法及装置 | |
CN108897036A (zh) | 一种地震数据处理方法及装置 | |
CN109938740A (zh) | 一种步态周期检测方法、装置及计算机可读存储介质 | |
CN107300519B (zh) | 致密砂岩导电分析方法及系统 | |
CN108107471B (zh) | 一种分方位初至数据体的获取方法及装置 | |
CN108181649A (zh) | 一种介质分界面振幅补偿方法及装置 | |
CN109725351A (zh) | 一种3d弹性波混合吸收边界条件的确定方法、装置及系统 | |
CN104091591A (zh) | 一种音频处理方法及装置 | |
CN109444957B (zh) | 基于决策树评价的水合物储层地震预测方法及处理终端 | |
CN110068860A (zh) | 一种时移地震数据处理方法、装置及系统 | |
CN110532755A (zh) | 一种计算机实现的风险识别的方法及装置 | |
CN104091595A (zh) | 一种音频处理方法及装置 | |
CN108710149A (zh) | 一种薄储层的预测方法及装置 | |
CN109856672B (zh) | 基于深度波数谱的瞬变波包提取方法、存储介质与终端 | |
CN102589551A (zh) | 一种基于小波变换的船用光纤陀螺信号实时滤波方法 | |
CN116224439B (zh) | 强各向异性层状vti介质高效射线追踪方法、设备和介质 | |
CN104459774A (zh) | 一种地质岩性差异识别方法及系统 | |
CN107688197B (zh) | 薄层预测方法和装置 | |
CN106896419B (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 |