CN112526611A - 表层地震波品质因子的提取方法及装置 - Google Patents
表层地震波品质因子的提取方法及装置 Download PDFInfo
- Publication number
- CN112526611A CN112526611A CN201910879834.7A CN201910879834A CN112526611A CN 112526611 A CN112526611 A CN 112526611A CN 201910879834 A CN201910879834 A CN 201910879834A CN 112526611 A CN112526611 A CN 112526611A
- Authority
- CN
- China
- Prior art keywords
- seismic wave
- quality factor
- wave quality
- measuring point
- data
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 239000002344 surface layer Substances 0.000 claims abstract description 63
- 230000003595 spectral effect Effects 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims description 45
- 238000001228 spectrum Methods 0.000 claims description 33
- 238000001514 detection method Methods 0.000 claims description 30
- 230000006870 function Effects 0.000 claims description 27
- 238000004590 computer program Methods 0.000 claims description 16
- 238000007781 pre-processing Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 7
- 230000001131 transforming effect Effects 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims 1
- 238000010521 absorption reaction Methods 0.000 abstract description 11
- 230000006872 improvement Effects 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 abstract description 7
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 238000003384 imaging method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 19
- 230000005284 excitation Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种表层地震波品质因子的提取方法及装置，其中该方法包括：获取待测区域内每个测点的Walkaway‑VSP数据；基于谱比法，根据每个测点的Walkaway‑VSP数据，计算每个测点的表层地震波品质因子的值；根据每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场。该方法在求解表层地震波品质因子的值时，应用Walkaway‑VSP技术，解决了因近地表的低、降速带导致的地震波吸收和衰减带来的地震处理补偿问题，从而得到可直接用于地面地震处理的表层地震波品质因子；容易实现，可靠性高，有助于后续的地震波准确成像、提升分辨率等方面的技术改进。
Description
技术领域
本发明涉及地球物理勘探中的地震数据采集、处理及分析技术，属于井中地球物理技术领域，尤其涉及一种表层地震波品质因子的提取方法及装置。
背景技术
表层地震波品质因子(表层Q因子)分析一直以来是地震勘探中的一个重点和难点课题，如何准确的提取表层地震波品质因子，对地震数据进行有效补偿是近年来地球物理学家研究的目标之一。早在上世纪七十年代，学者们已认识到，要提高地震记录的质量，衰减补偿是必须考虑的地震参数。所以除地震波速度外，衰减补偿是另外一个可用于岩性研究的地震参数。Hamilton和Hauge最早进行了衰减因子与岩性间的关系研究(Hamilton，1972；Hauge，l981，1983)。
目前地球物理界公认，近地表低、降速会给地震勘探带来了很多不利影响。相对致密的深部地层来说，近地表疏松的低、降速带地层对地震高频信号的吸收和衰减非常严重，导致了地震子波振幅衰减和相位畸变，降低了地震勘探的可靠性。对于深地层的吸收衰减，目前做法是通过零井源距VSP测井来计算吸收因子，但利用该VSP测井方式无法计算浅地表的吸收因子。
因此需要一种方法来精确计算表层地震波品质因子值，消除近地表的低、降速带导致的地震波衰减，提高地震资料的可靠性。
发明内容
本发明实施例提供一种表层地震波品质因子的提取方法，用以解决因近地表的低、降速带导致的地震波吸收和衰减带来的地震处理补偿问题，得到可直接用于地面地震处理的表层地震波品质因子，该方法包括：
获取待测区域内每个测点的Walkaway-VSP数据；
基于谱比法，根据所述每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；
根据所述每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；所述待测区域的表层地震波品质因子场，用于表征待测区域各处的表层地震波品质因子的值。
本发明实施例还提供一种表层地震波品质因子的提取装置，用以解决因近地表的低、降速带导致的地震波吸收和衰减带来的地震处理补偿问题，得到可直接用于地面地震处理的表层地震波品质因子，该装置包括：
数据获取模块，用于获取待测区域内每个测点的Walkaway-VSP数据；
数据计算模块，用于基于谱比法，根据所述每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；
表层Q场建立模块，用于根据所述每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；所述待测区域的表层地震波品质因子场，用于表征待测区域各处的表层地震波品质因子的值。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述表层地震波品质因子的提取方法。
本发明实施例也提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述表层地震波品质因子的提取方法的计算机程序。
本发明实施例中，通过获取待测区域内每个测点的Walkaway-VSP数据；基于谱比法，根据所述每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；根据所述每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；求解表层地震波品质因子的值时，应用Walkaway-VSP技术，解决了因近地表的低、降速带导致的地震波吸收和衰减带来的地震处理补偿问题，从而得到可直接用于地面地震处理的表层地震波品质因子；该方法容易实现，可靠性高，有助于后续的地震波准确成像、提升分辨率等方面的技术改进。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中表层地震波品质因子的提取方法示意图。
图2为本发明具体实施例中表层地震波品质因子的提取方法示意图。
图3为本发明具体实施例中步骤201的具体实施流程示意图。
图4为本发明具体实施例中的共检波点道集数据曲线示意图。
图5为本发明具体实施例中的表层地震波品质因子Q与表层厚度的关系函数拟合曲线示意图。
图6为本发明具体实施例中的建立的表层地震波品质因子场的示意图。
图7为本发明具体实施例中的观测系统的示意图。
图8为本发明实施例中表层地震波品质因子的提取装置示意图。
图9为本发明具体实施例中表层地震波品质因子的提取装置示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为了解决因近地表的低、降速带导致的地震波吸收和衰减带来的地震处理补偿的问题，本发明实施例提供了一种表层地震波品质因子的提取方法，用以得到可直接用于地面地震处理的表层地震波品质因子，如图1所示，该方法包括：
步骤101：获取待测区域内每个测点的Walkaway-VSP数据；
步骤102：基于谱比法，根据每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；
步骤103：根据每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；待测区域的表层地震波品质因子场，用于表征待测区域各处的表层地震波品质因子的值。
由图1所示的流程可以看出，通过获取待测区域内每个测点的Walkaway-VSP数据；基于谱比法，根据每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；根据每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；求解表层地震波品质因子的值时，应用Walkaway-VSP技术，解决了因近地表的低、降速带导致的地震波吸收和衰减带来的地震处理补偿问题，从而得到可直接用于地面地震处理的表层地震波品质因子；该方法容易实现，可靠性高，有助于后续的地震波准确成像、提升分辨率等方面的技术改进。
具体实施时，首先获取待测区域内每个测点的Walkaway-VSP数据，上述Walkaway-VSP数据是基于Walkaway-VSP(线性变井源距垂直地震剖面法，Walkaway-VerticalSeismic Profiling)技术采集获取的地震资料，例如包括：Walkaway-VSP原始记录数据、观测班报、测量成果等，这些数据整理后得到Walkaway-VSP的共激发点道集数据。
获取每个测点的Walkaway-VSP数据后，在图1所示的方法的基础上还需要进行数据预处理，如图2所示，在图1的基础上还包括：
步骤201：对每个测点的Walkaway-VSP数据预处理。
上述步骤201的具体处理过程例如可以如图3所示，具体例如可以包括：
步骤301：对每个测点的Walkaway-VSP数据中的共激发点道集数据进行去噪处理；
步骤302：根据去噪处理后的数据，以检波点序号重新排序，变换得到共检波点道集数据。
具体实施时，首先对共激发点道集数据进行去噪处理，对每个共激发点道集数据，进行例如去随机噪声、去谐波等噪声压制的处理。进行上述处理后，评价该道集的数据品质，进行数据筛选，例如舍弃品质较差的道集数据，品质较差的道集数据例如可以包括：受面波影响的近地表数据、受套管波影响的数据、受电(光)缆影响的数据以及其他信噪比较低的数据。由于本发明实施例基于大量的数据统计，因此，舍弃较差的部分数据可以使得最终的计算结果更加准确，提高可靠性。
其次，根据去噪处理后的数据，以检波点序号重新排序，变换得到共检波点道集数据。例如图4所示的为一具体实施例中的共检波点道集数据曲线示意图，其中，横坐标代表共检波点道集道号(CDP)，纵坐标代表记录时间t。
进行数据预处理后，基于谱比法，根据每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子，具体步骤包括：
根据每个测点的Walkaway-VSP数据变换得到的共检波点道集数据，拾取初至时间及续至时间，得到计算时间；
根据每个测点的Walkaway-VSP数据变换得到的共检波点道集数据，获得初至频谱及续至频谱，得到计算振幅谱；
根据计算振幅谱的对数比值计算结果，进行直线拟合得到谱比计算斜率；
根据计算时间及谱比计算斜率，按照如下公式计算得到各个续至波处的地震波品质因子的值：
Qn＝-π(tn+1-t1)/kn
其中，Qn代表第n个续至波处的地震波品质因子值；n≥1代表可有效识别的续至波的个数；
t1代表第一计算时间；tn+1代表第n+1计算时间；
kn代表第n个续至波处的谱比计算斜率；
将每个测点对应的各个续至波处的地震波品质因子的值取平均，得到每个测点的表层地震波品质因子的值。
具体实施时，首先根据上述共检波点道集数据，拾取初至时间及续至时间，得到计算时间。在上述共检波点道集数据上，识别并拾取首波时距曲线，得到初至时间，例如可记作第一计算时间t1；识别并拾取第一组续至波时距曲线，得到第一组续至波的续至时间，例如记作第二计算时间t2；如果有更多的续至波可识别，则可用相同方法得到第三计算时间t3、第四计算时间t4……第n计算时间tn。对上述时距曲线的拾取，严格要求拾取波峰位置对应时间。
根据上述共检波点道集数据,获得初至频谱及续至频谱，得到计算振幅谱。在上述共检波点道集数据上，以首波的平均周期为时窗，对时窗内数据进行快速度傅立叶变换，得到初至频谱，记作第一计算振幅谱，例如可记作A1(f)；以第一组续至波的平均周期为时窗，对时窗内数据进行快速度傅立叶变换，得到第一组续至波的续至频谱，记作第二计算振幅谱，例如可记作A2(f)；如果有更多的续至波可识别，利用相同方法可得到第三计算振幅谱A3(f)、第四计算振幅谱A4(f)……第n计算振幅谱An(f)。其中A代表振幅，f代表频率。
本领域技术人员可以理解，得到计算时间与得到计算振幅谱的步骤之间不存在特定的顺序，可以同时进行，也可以先得到计算振幅谱，再得到计算时间。上述步骤流程仅为其中一种具体实施，不用于限制本发明。
接着，利用第一计算振幅谱A1(f)、第二计算振幅谱A2(f)，例如可以取对数后求取比值，即ln(A2(f))/ln(A1(f))。将频率f作为自变量，ln(A2(f))/ln(A1(f))作为因变量，进行拟合得到谱比计算斜率，例如采用直线拟合的方法得到谱比计算斜率k1；利用第一计算振幅谱A1(f)、第三计算振幅谱A3(f)取对数后求取比值，即ln(A3(f))/ln(A1(f))，例如采用直线拟合的方法得到谱比计算斜率k2；利用相同的方法，计算得到k3、k4……kn。
根据计算时间及谱比计算斜率，按照如下公式计算得到各个续至波处的地震波品质因子的值：
Qn＝-π(tn+1-t1)/kn
其中，Qn代表第n个续至波处的地震波品质因子值；n≥1，代表可有效识别的续至波的个数；
t1代表第一计算时间；tn+1代表第n+1计算时间；
kn代表第n个续至波处的谱比计算斜率；
举例说明，利用第一计算时间t1、第二计算时间t2，计算得到第1个续至波处的地震波品质因子Q1＝-π(t2-t1)/k1；利用第一计算时间t1、第三计算时间t3，计算得到品质因子第2个续至波处的地震波Q2＝-π(t3-t1)/k2；利用相同的方法，计算得到Q3、Q4……Qn。
计算得到各个续至波处的地震波品质因子的值后，将每个测点对应的各个续至波处的地震波品质因子的值取平均，得到每个测点的表层地震波品质因子的值Q＝(Q1+Q2+…+Qn)/n。
可以理解的是，上述计算公式仅为一例，实施时可以对该公式进行变形，或采用其它公式或方法，基于谱比法，根据每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值的公式或方法均落入本发明的保护范围，实施例中不再赘述。
计算得到每个测点的表层地震波品质因子的值后，根据每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场，具体步骤包括：
根据每个测点的表层地震波品质因子的值以及每个测点的表层厚度，拟合得到表层地震波品质因子与表层厚度的关系函数；该关系函数用于表征表层地震波品质因子与表层厚度之间的对应关系；
获取待测区域的表层厚度模型，将待测区域的表层厚度模型代入上述关系函数，得到表层地震波品质因子场。
具体实施时，根据上述每个测点的表层地震波品质因子Q的值以及每个测点的表层厚度h，进行拟合，例如进行一次函数拟合得到关系函数Q(h)，该函数建立了表层地震波品质因子Q以及表层厚度h的数值之间的对应关系。如图5所示为一具体实施例中表层地震波品质因子Q与表层厚度h的关系拟合曲线示意图，纵坐标代表Q值，横坐标代表表层厚度h。本领域技术人员可以理解，此处进行一次函数拟合得到关系函数Q(h)仅为一例，只要根据上述表层地震波品质因子Q以及对应此处的表层厚度h进行拟合得到关系函数，例如二次关系函数、数值映射、散点图等形式，都落入本发明的保护范围，此处不再赘述。
得到关系函数Q(h)后，获取待测区域的表层厚度模型，将待测区域的表层厚度模型代入关系函数Q(h)中，得到待测区域的表层地震波品质因子场；该待测区域的表层地震波品质因子场，用于表征待测区域各处的表层地震波品质因子的值。如图6所示为一具体实施例中的建立的表层地震波品质因子场的示意图，纵坐标代表南北向相对地理坐标；横坐标代表东西向相对地理坐标。
下面给出一具体实例说明如何进行表层地震波品质因子的提取：
具体实施时，首先，确定Walkaway-VSP的观测井。在待测区域内，例如可以选择一口井中地震观测井。井中地震观测井的具体选择，例如可以包括以下几方面要求：一是井况条件适宜于井中观测，至少保证1000米井段内可下放检波器；二是井区适宜地面可控震源激发；三是所在井区地表条件对于全区具有代表性。
其次，确定Walkaway-VSP的观测系统。例如包括以下确定要求：一是设计接收井段，综合考虑作业成本和数据需求，以500～1000米深度作为接收井段为较佳实施例；采用10米或20米检波点间距。二是设计激发点测线，该测线的范围应涵盖有代表性的不同表层结构地区，测线长度和激发点间距可根据区域环境适当调整，该观测系统一具体实施示意图如图7所示。
利用上述观测井及观测系统，基于Walkaway-VSP进行数据采集，获取待测区域内每个测点的Walkaway-VSP数据。
本领域技术人员可以理解，上述对观测井以及观测系统的选择及确定仅为举例，可根据实际需要，例如实际工况以及地质条件等，进行调整，此处不再赘述。
获取待测区域内每个测点的Walkaway-VSP数据后，基于谱比法，根据每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；
计算每个测点的表层地震波品质因子的值后，根据每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；该待测区域的表层地震波品质因子场，用于表征待测区域各处的表层地震波品质因子的值。
基于同一发明构思，本发明实施例还提供一种表层地震波品质因子的提取装置，由于表层地震波品质因子的提取装置所解决问题的原理与表层地震波品质因子的提取方法相似，因此表层地震波品质因子的提取装置的实施可以参见表层地震波品质因子的提取方法的实施，重复之处不再赘述，具体结构如图8所示：
数据获取模块801，用于获取待测区域内每个测点的Walkaway-VSP数据；
数据计算模块802，用于基于谱比法，根据每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；
表层Q场建立模块803，用于根据每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；待测区域的表层地震波品质因子场，用于表征待测区域内各处的表层地震波品质因子的值。
本发明一具体实施中的表层地震波品质因子的提取装置，如图9所示，在图8的基础上还包括：
预处理模块901，用于对每个测点的Walkaway-VSP数据预处理；
预处理模块901具体用于：
对每个测点的Walkaway-VSP数据中的共激发点道集数据进行去噪处理；
根据去噪处理后的数据，以检波点序号重新排序，变换得到共检波点道集数据；
具体实施时，数据计算模块802具体用于：
根据每个测点的Walkaway-VSP数据变换得到的共检波点道集数据，拾取初至时间及续至时间，得到计算时间；
根据每个测点的Walkaway-VSP数据变换得到的共检波点道集数据，获得初至频谱及续至频谱，得到计算振幅谱；
根据计算振幅谱的对数比值计算结果，进行直线拟合得到谱比计算斜率；
根据计算时间及谱比计算斜率，按照如下公式计算得到各个续至波处的地震波品质因子的值：
Qn＝-π(tn+1-t1)/kn
其中，Qn代表第n个续至波处的地震波品质因子值；n≥1代表可有效识别的续至波的个数；
t1代表第一计算时间；tn+1代表第n+1计算时间；
kn代表第n个续至波处的谱比计算斜率；
将每个测点对应的各个续至波处的地震波品质因子的值取平均，得到每个测点的表层地震波品质因子的值。
具体实施时，表层Q场建立模块803具体用于：
根据每个测点的表层地震波品质因子的值以及每个测点的表层厚度，拟合得到表层地震波品质因子与表层厚度的关系函数；关系函数用于表征表层地震波品质因子与表层厚度之间的对应关系；
将待测区域的表层厚度模型，代入关系函数，得到待测区域的表层地震波品质因子场。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述表层地震波品质因子的提取方法。
本发明实施例还提供一种计算机可读存储介质，存储有执行上述表层地震波品质因子的提取方法的计算机程序。
综上所述，本发明实施例提供的表层地震波品质因子的提取方法及装置具有如下优点：
通过通过获取待测区域内每个测点的Walkaway-VSP数据；基于谱比法，根据每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；根据每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；求解表层地震波品质因子的值时，应用Walkaway-VSP技术，解决了因近地表的低、降速带导致的地震波吸收和衰减带来的地震处理补偿问题，从而得到可直接用于地面地震处理的表层地震波品质因子；该方法容易实现，可靠性高，有助于后续的地震波准确成像、提升分辨率等方面的技术改进；通过对Walkaway-VSP数据预处理，筛选出较差的共激发点道集数据后舍弃，使得计算结果更加准确，提高了可靠性。
本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、装置和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种表层地震波品质因子的提取方法，其特征在于，包括：
获取待测区域内每个测点的Walkaway-VSP数据；
基于谱比法，根据所述每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；
根据所述每个测点的表层地震波品质因子的值，以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；所述待测区域的表层地震波品质因子场，用于表征待测区域各处的表层地震波品质因子的值。
2.如权利要求1所述的方法，其特征在于，还包括：对所述每个测点的Walkaway-VSP数据预处理，具体处理包括：
对所述每个测点的Walkaway-VSP数据中的共激发点道集数据进行去噪处理；
根据去噪处理后的数据，以检波点序号重新排序，变换得到共检波点道集数据。
3.如权利要求2所述的方法，其特征在于，基于谱比法，根据所述每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值，包括：
根据每个测点的Walkaway-VSP数据变换得到的所述共检波点道集数据，拾取初至时间及续至时间，得到计算时间；
根据每个测点的Walkaway-VSP数据变换得到的所述共检波点道集数据，获得初至频谱及续至频谱，得到计算振幅谱；
根据所述计算振幅谱的对数比值计算结果，进行直线拟合得到谱比计算斜率；
根据所述计算时间及所述谱比计算斜率，按照如下公式计算得到各个续至波处的地震波品质因子的值：
Qn＝-π(tn+1-t1)/kn
其中，Qn代表第n个续至波处的地震波品质因子值；n≥1代表可有效识别的续至波的个数；
t1代表第一计算时间；tn+1代表第n+1计算时间；
kn代表第n个续至波处的谱比计算斜率；
将每个测点对应的各个续至波处的地震波品质因子的值取平均，得到每个测点的表层地震波品质因子的值。
4.如权利要求1所述的方法，其特征在于，根据所述每个测点的表层地震波品质因子的值，以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场，包括：
根据所述每个测点的表层地震波品质因子的值，以及每个测点的表层厚度，拟合得到表层地震波品质因子与表层厚度的关系函数；所述关系函数用于表征表层地震波品质因子与表层厚度之间的对应关系；
将待测区域的表层厚度模型，代入所述关系函数，得到所述待测区域的表层地震波品质因子场。
5.一种表层地震波品质因子的提取装置，其特征在于，包括：
数据获取模块，用于获取待测区域内每个测点的Walkaway-VSP数据；
数据计算模块，用于基于谱比法，根据所述每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；
表层Q场建立模块，用于根据所述每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；所述待测区域的表层地震波品质因子场，用于表征待测区域内各处的表层地震波品质因子的值。
6.如权利要求5所述的装置，其特征在于，还包括：
预处理模块，对所述每个测点的Walkaway-VSP数据预处理；
具体用于：
对所述每个测点的Walkaway-VSP数据中的共激发点道集数据进行去噪处理；
根据去噪处理后的数据，以检波点序号重新排序，变换得到共检波点道集数据。
7.如权利要求6所述的装置，其特征在于，所述数据计算模块具体用于：
根据每个测点的Walkaway-VSP数据变换得到的所述共检波点道集数据，拾取初至时间及续至时间，得到计算时间；
根据每个测点的Walkaway-VSP数据变换得到的所述共检波点道集数据，获得初至频谱及续至频谱，得到计算振幅谱；
根据所述计算振幅谱的对数比值计算结果，进行直线拟合得到谱比计算斜率；
根据所述计算时间及所述谱比计算斜率，按照如下公式计算得到各个续至波处的地震波品质因子的值：
Qn＝-π(tn+1-t1)/kn
其中，Qn代表第n个续至波处的地震波品质因子值；n≥1代表可有效识别的续至波的个数；
t1代表第一计算时间；tn+1代表第n+1计算时间；
kn代表第n个续至波处的谱比计算斜率；
将每个测点对应的各个续至波处的地震波品质因子的值取平均，得到每个测点的表层地震波品质因子的值。
8.如权利要求5所述的装置，其特征在于，所述表层Q场建立模块具体用于：
根据所述每个测点的表层地震波品质因子的值以及每个测点的表层厚度，拟合得到表层地震波品质因子与表层厚度的关系函数；所述关系函数用于表征表层地震波品质因子与表层厚度之间的对应关系；
将待测区域的表层厚度模型，代入所述关系函数，得到所述待测区域的表层地震波品质因子场。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一所述方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910879834.7A CN112526611A (zh) | 2019-09-18 | 2019-09-18 | 表层地震波品质因子的提取方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910879834.7A CN112526611A (zh) | 2019-09-18 | 2019-09-18 | 表层地震波品质因子的提取方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112526611A true CN112526611A (zh) | 2021-03-19 |
Family
ID=74974985
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910879834.7A Pending CN112526611A (zh) | 2019-09-18 | 2019-09-18 | 表层地震波品质因子的提取方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112526611A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114089416A (zh) * | 2021-11-17 | 2022-02-25 | 成都理工大学 | 一种利用薛定谔方程进行地震波衰减梯度估计的方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020185329A1 (en) * | 2001-04-30 | 2002-12-12 | The Regents Of The University Of California | Frequency-dependent processing and interpretation (FDPI) of seismic data for identifiying, imaging and monitoring fluid-saturated underground reservoirs |
| US20060104158A1 (en) * | 2004-11-17 | 2006-05-18 | Joel Walls | Method for determining formation quality factor from well log data and its application to seismic reservoir characterization |
| AU2009243472A1 (en) * | 2009-12-02 | 2011-06-16 | Zuercher, Hannes Dr | Locating oil or gas actively or passively by observing a porous oil and gas saturated system giving off its characteristic resonance response to artificial excitation or ambient background noise, including optional differentiation of oil, gas and water |
| CN102540253A (zh) * | 2011-12-30 | 2012-07-04 | 中国石油天然气股份有限公司 | 一种陡倾角地层及裂缝地震偏移方法和装置 |
| CN102692648A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种基于大炮初至的折射波剩余静校正方法 |
| CN105607124A (zh) * | 2016-03-09 | 2016-05-25 | 蒋立 | 地震波近地表地层品质因子的补偿方法及装置 |
| CN106291693A (zh) * | 2015-05-21 | 2017-01-04 | 中国石油化工股份有限公司 | 一种基于广义s变换的叠前q值反演方法及系统 |
| CN106443786A (zh) * | 2016-11-14 | 2017-02-22 | 中国科学院地质与地球物理研究所 | 基于地面接收的反射地震资料的q值场建模方法 |
| CN106646601A (zh) * | 2016-12-28 | 2017-05-10 | 中国石油化工股份有限公司 | 多信息联合约束的浅中深层三维q体建立方法 |
| CN107272062A (zh) * | 2017-07-05 | 2017-10-20 | 西安交通大学 | 一种数据驱动的地下介质q场估计方法 |
-
2019
- 2019-09-18 CN CN201910879834.7A patent/CN112526611A/zh active Pending
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020185329A1 (en) * | 2001-04-30 | 2002-12-12 | The Regents Of The University Of California | Frequency-dependent processing and interpretation (FDPI) of seismic data for identifiying, imaging and monitoring fluid-saturated underground reservoirs |
| US20060104158A1 (en) * | 2004-11-17 | 2006-05-18 | Joel Walls | Method for determining formation quality factor from well log data and its application to seismic reservoir characterization |
| AU2009243472A1 (en) * | 2009-12-02 | 2011-06-16 | Zuercher, Hannes Dr | Locating oil or gas actively or passively by observing a porous oil and gas saturated system giving off its characteristic resonance response to artificial excitation or ambient background noise, including optional differentiation of oil, gas and water |
| CN102692648A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种基于大炮初至的折射波剩余静校正方法 |
| CN102540253A (zh) * | 2011-12-30 | 2012-07-04 | 中国石油天然气股份有限公司 | 一种陡倾角地层及裂缝地震偏移方法和装置 |
| CN106291693A (zh) * | 2015-05-21 | 2017-01-04 | 中国石油化工股份有限公司 | 一种基于广义s变换的叠前q值反演方法及系统 |
| CN105607124A (zh) * | 2016-03-09 | 2016-05-25 | 蒋立 | 地震波近地表地层品质因子的补偿方法及装置 |
| CN106443786A (zh) * | 2016-11-14 | 2017-02-22 | 中国科学院地质与地球物理研究所 | 基于地面接收的反射地震资料的q值场建模方法 |
| CN106646601A (zh) * | 2016-12-28 | 2017-05-10 | 中国石油化工股份有限公司 | 多信息联合约束的浅中深层三维q体建立方法 |
| CN107272062A (zh) * | 2017-07-05 | 2017-10-20 | 西安交通大学 | 一种数据驱动的地下介质q场估计方法 |
Non-Patent Citations (9)
| Title |
|---|
| S DANGEL ET AL: "《Phenomenology of tremor-like signals observed over hydrocarbon reservoirs》", JOURNAL OF VOLCANOLOGY AND GEOTHERMAL RESEARCH, vol. 128, 15 November 2003 (2003-11-15), pages 135 - 158 * |
| 李生杰等: "《地层衰减数据体的建立》", 新疆地质, vol. 19, no. 2, 30 June 2001 (2001-06-30), pages 146 - 149 * |
| 杨宇山等: "《Walkaway井地联合地震资料Q值波形反演方法》", 石油天然气学报, vol. 31, no. 2, 30 April 2009 (2009-04-30), pages 53 - 58 * |
| 杨智超等: "《近地表地层品质因子在川西北部地区地震勘探中的应用》", 天然气勘探与开发, vol. 41, no. 3, 30 September 2018 (2018-09-30), pages 14 - 19 * |
| 王晓涛等: "《准噶尔盆地厚沙漠区表层实测 Q 采集及计算方法研究》", 中国石油学会2015年物探技术研讨会论文集, 30 April 2015 (2015-04-30), pages 246 * |
| 王晓涛等: "《强衰减条件下能量比法的修正及其应用》", 石油地球物理勘探, vol. 52, no. 1, 28 February 2017 (2017-02-28), pages 42 - 47 * |
| 王静等: "《VSP全井段Q值估计》", 石油地球物理勘探, vol. 53, no. 2, 31 December 2018 (2018-12-31), pages 59 - 62 * |
| 罗勇等: "《利用多项式拟合质心频移法估算近地表Q值》", 石油地球物理勘探, vol. 51, no. 3, 30 June 2016 (2016-06-30), pages 589 - 595 * |
| 陈志德等: "《近地表Q值模型建立方法及其地震叠前补偿应用》", 石油学报, vol. 36, no. 2, 28 February 2015 (2015-02-28), pages 188 - 196 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114089416A (zh) * | 2021-11-17 | 2022-02-25 | 成都理工大学 | 一种利用薛定谔方程进行地震波衰减梯度估计的方法 |
| CN114089416B (zh) * | 2021-11-17 | 2023-02-21 | 成都理工大学 | 一种利用薛定谔方程进行地震波衰减梯度估计的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2579164C1 (ru) | Способ обращения для определения добротности геологической среды | |
| DK177865B1 (da) | Fremgangsmåde til detektering eller monitorering af en subsurface-struktur af carbonhydridreservoirstørrelse | |
| US7917295B2 (en) | Modeling and filtering coherent noise in seismic surveying | |
| CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
| CN101551465B (zh) | 一种自适应识别和消除地震勘探单频干扰的方法 | |
| CN113640881B (zh) | 多偏移距二维横向高分辨率瞬态面波探测方法 | |
| CN110907995B (zh) | 井中vsp地震数据的逆时偏移方法及装置 | |
| CN102176054A (zh) | 近地表综合信息处理解释方法 | |
| CN116520419B (zh) | 一种热流体裂缝通道识别方法 | |
| CN106199698A (zh) | 基于多次波信息的频率域地震数据重构方法 | |
| CN114415234B (zh) | 基于主动源面波频散和h/v确定浅地表横波速度的方法 | |
| CN104597485B (zh) | 一种微小断层检测方法及断层检测装置 | |
| CN109143345B (zh) | 基于模拟退火的品质因子q非线性反演方法及系统 | |
| CN112526611A (zh) | 表层地震波品质因子的提取方法及装置 | |
| EP2249183A1 (en) | Modeling and filtering coherent noise in seismic surveying | |
| CN110244383B (zh) | 基于近地表数据的地质岩性综合模型创建方法 | |
| Gao et al. | Acquisition and processing pitfall with clipped traces in surface-wave analysis | |
| CN110780346A (zh) | 一种隧道超前探测复杂地震波场的分离方法 | |
| CN112415601B (zh) | 表层品质因子q值的确定方法及装置 | |
| CN112379429B (zh) | 地震数据的振幅补偿方法及装置 | |
| CN110568491B (zh) | 一种品质因子q的估算方法 | |
| CN114137606A (zh) | 一种稳健的谱模拟反褶积方法 | |
| CN108761534B (zh) | 陆上地震加速度信号应用新方法 | |
| Yao et al. | Microseismic signal denoising using simple bandpass filtering based on normal time–frequency transform | |
| CN103760599B (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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Cai Zhidong Inventor after: Zhu Yang Inventor after: Zhao Xiaohui Inventor after: Li Fei Inventor after: Wang Jing Inventor after: Wang Yong Inventor before: Cai Zhidong Inventor before: Zhu Yang Inventor before: Zhao Xiaohui Inventor before: Li Fei Inventor before: Wang Jing Inventor before: Wang Yong |