CN112859161A - 基于压缩感知的地震勘探方法及装置 - Google Patents
基于压缩感知的地震勘探方法及装置 Download PDFInfo
- Publication number
- CN112859161A CN112859161A CN201911187300.4A CN201911187300A CN112859161A CN 112859161 A CN112859161 A CN 112859161A CN 201911187300 A CN201911187300 A CN 201911187300A CN 112859161 A CN112859161 A CN 112859161A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- excitation
- seismic
- receiving
- points
- 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 59
- 230000005284 excitation Effects 0.000 claims abstract description 71
- 230000001788 irregular Effects 0.000 claims description 35
- 238000012545 processing Methods 0.000 claims description 24
- 238000005070 sampling Methods 0.000 claims description 23
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 16
- 238000005457 optimization Methods 0.000 claims description 10
- 238000007781 pre-processing Methods 0.000 claims description 10
- 238000013508 migration Methods 0.000 claims description 8
- 230000005012 migration Effects 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 7
- 230000003068 static effect Effects 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 239000002360 explosive Substances 0.000 claims description 4
- 239000004215 Carbon black (E152) Substances 0.000 claims description 3
- 229930195733 hydrocarbon Natural products 0.000 claims description 3
- 150000002430 hydrocarbons Chemical class 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 13
- 238000005516 engineering process Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 6
- 238000010304 firing Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000017105 transposition Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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/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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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
- G01V2210/624—Reservoir parameters
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
技术领域
本发明涉及地球物理勘探技术领域,在地球物理勘探过程中,尤其涉及一种基于压缩感知的地震勘探方法及装置。
背景技术
Shannon/Nyquist采样定理指出,当采样频率大于模拟信号中最高频率的2倍时,采样之后的数字信号能够不失真的恢复出模拟信号。Shannon/Nyquist采样定理作为模拟信号和数字信号之间的桥梁,几十年来一直支撑并引导着现代信号处理各个领域的技术发展,地震勘探技术同样是基于Shannon/Nyquist采样定理一种信号处理技术,并不断向前发展。
近年来,在地震勘探中,为了得到更高成像质量的剖面,广泛应用的是宽频、宽方位和高密度的地震勘探技术,即“两宽一高”地震勘探技术,“两宽一高”技术带来的问题是海量地震数据的采集、存储成本的大幅上升。因此,如何在不降低成像质量的情况下,降低勘探成本成为业内亟需解决的问题。
发明内容
本发明目的在于提供一种完整的基于压缩感知地震勘探方法及装置,以有效降低地震勘探成本。
为达上述目的,本发明所提供的基于压缩感知的地震勘探方法,具体包含:获取待勘探项目的项目需求,根据项目需求获得预存的对应观测系统;对所述观测系统进行稀疏性约束获得不规则化后的激发点和接收点,以所述激发点和所述接收点作为所述观测系统的输出数据;通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据;对所述地震数据进行数据重构恢复为规则的地震数据,并根据规则的地震数据计算获得油气储层信息。
在上述基于压缩感知的地震勘探方法中,优选的,通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据还包含:对所述地震数据进行静校正和/或强噪声衰减处理。
在上述基于压缩感知的地震勘探方法中,优选的,根据项目需求获得预存的对应观测系统中所述观测系统包含:接收线数、接收线距、接收点数、接收点距、激发线数、激发线距、激发点数、激发点距以及接收和激发的起始位置。
在上述基于压缩感知的地震勘探方法中,优选的,对所述观测系统进行稀疏性约束获得不规则化后的激发点和接收点包含:根据所述项目需求获得预设的接收点和激发点的欠采样比例因子;通过预定稀疏性约束长度将接收线和激发线分别为多个长度为预定稀疏性约束长度的子集;根据所述接收点和激发点的欠采样比例因子与所述子集获得不规则化后的接收点和激发点。
在上述基于压缩感知的地震勘探方法中,优选的,通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据包含:将检波器布设于接收点,通过可控震源或爆炸物于激发点激发,并利用记录仪采集记录激发的地震数据。
在上述基于压缩感知的地震勘探方法中,优选的,对所述地震数据进行数据重构恢复为规则的地震数据包含:通过最优化法求解以下公式,将所述地震数据进行数据重构恢复为规则的地震数据:
min||x||1s.t.||ΦCTx-b||2<τ;
在上式中s.t.表示约束条件,b不规则的地震数据,x为稀疏向量,Φ为观测矩阵,CT为稀疏变换矩阵的转置,τ为阈值。
在上述基于压缩感知的地震勘探方法中,优选的,根据规则的地震数据计算获得油气储层信息包含:根据规则的地震数据计算获得叠加或偏移剖面;根据所述叠加或偏移剖面计算获得油气储层信息。
本发明还提供一种基于压缩感知的地震勘探装置,所述装置包含分析模块、处理模块、采集模块和计算模块;所述分析模块用于获取待勘探项目的项目需求,根据项目需求获得预存的对应观测系统;所述处理模块用于对所述观测系统进行稀疏性约束获得不规则化后的激发点和接收点,以所述激发点和所述接收点作为所述观测系统的输出数据;所述采集模块用于通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据;所述计算模块用于对所述地震数据进行数据重构恢复为规则的地震数据,并根据规则的地震数据计算获得油气储层信息。
在上述基于压缩感知的地震勘探装置中,优选的,所述采集模块还包含预处理模块,所述预处理模块用于对所述地震数据进行静校正和/或强噪声衰减处理。
在上述基于压缩感知的地震勘探装置中,优选的,所述处理模块还包含:根据所述项目需求获得预设的接收点和激发点的欠采样比例因子;通过预定稀疏性约束长度将接收线和激发线分别为多个长度为预定稀疏性约束长度的子集;根据所述接收点和激发点的欠采样比例因子与所述子集获得不规则化后的接收点和激发点。
在上述基于压缩感知的地震勘探装置中,优选的,所述计算模块还包含:通过最优化法求解以下公式,将所述地震数据进行数据重构恢复为规则的地震数据:
min||x||1s.t.||ΦCTx-b||2<τ;
在上式中s.t.表示约束条件,b不规则的地震数据,x为稀疏向量,Φ为观测矩阵,CT为稀疏变换矩阵的转置,τ为阈值。
在上述基于压缩感知的地震勘探装置中,优选的,所述计算模块还包含:根据规则的地震数据计算获得叠加或偏移剖面;根据所述叠加或偏移剖面计算获得油气储层信息。
本发明还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述方法。
本发明还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述方法的计算机程序。
本发明的有益技术效果在于:(1)降低勘探成本。在野外采集资料的过程中,应用较少的野外勘探设备进行稀疏采样,完成与常规地震勘探方法相当的勘探效果,从而降低勘探成本;(2)扩大勘探面积。在设备投入和勘探网格不变的情况下,通过扩大勘探区域,利用稀疏采样完成野外资料采集,从而达到在不增加勘探成本的情况下,扩大勘探面积的目的。(3)提高地震勘探的空间分辨率。在不增加勘探成本的情况下,通过缩小炮检点网格,利用稀疏采样完成野外资料采集,从而达到在不增加勘探成本的情况下,提高地震勘探空间分辨率。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1A为本发明一实施例所提供的基于压缩感知地震勘探方法的流程示意图;
图1B为本发明一实施例所提供的基于压缩感知地震勘探方法的原理流程图;
图2为本发明一实施例所提供的常规观测系统检测示意图;
图3为本发明一实施例所提供的稀疏性约束的非规则观测系统检测示意图;
图4为本发明一实施例所提供的非规则观测系统的单炮记录示意图;
图5为本发明一实施例所提供的重构后的单炮记录示意图;
图6为本发明一实施例所提供的处理后的剖面示意图;
图7为本发明一实施例所提供的基于压缩感知地震勘探装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
请参考图1A所示,本发明所提供的基于压缩感知的地震勘探方法,具体包含:S101获取待勘探项目的项目需求,根据项目需求获得预存的对应观测系统;S102对所述观测系统进行稀疏性约束获得不规则化后的激发点和接收点,以所述激发点和所述接收点作为所述观测系统的输出数据;S103通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据;S104对所述地震数据进行数据重构恢复为规则的地震数据,并根据规则的地震数据计算获得油气储层信息。以此,将压缩感知理论引入到地震勘探技术中,使在不降低成像质量的情况下,有效降低勘探成本成为可能。值得说明的是,压缩感知理论突破了Shannon/Nyquist采样理论,它指出对稀疏的或可压缩的信号,采用线性随机观测的方法,获取少量的观测数据,再利用非线性的优化计算方法可以精确地重构信号。本发明利用压缩感知理论,进行特殊方式的非规则采样,与常规采集相比,得到的是更少的非规则地震数据,即实现压缩采样,降低勘探成本,然后通过非线性的优化计算方法精确地重构规则的地震数据,达到与常规勘探相当的成像质量。
在上述实施例中,根据项目需求获得预存的对应观测系统中所述观测系统包含:接收线数、接收线距、接收点数、接收点距、激发线数、激发线距、激发点数、激发点距以及接收和激发的起始位置。当然,实际工作中,也可根据勘探项目需求,确定采用压缩感知地震勘探技术的需求是降低成本、扩大勘探面积或提高空间分辨率;例如:根据项目的需求,设计常规规则网格的观测系统,包括确定接收线数、接收线距、接收点数、接收点距、激发线数、激发线距、激发点数、激发点距以及接收和激发的起始位置等;如果项目需求是降低成本,则常规地震勘探观测系统不需改变;如果是扩大勘探面积,需要根据勘探区域,需要在原观测系统的基础上,增加布设激发线和接收线;如果是提高空间分辨率,则需要缩小激发点距、激发线距,接收点距、接收线距,重新布设炮检点;其后,通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据可包含:将检波器布设于接收点,通过可控震源或爆炸物于激发点激发,并利用记录仪采集记录激发的地震数据。
在本发明一实施例中,通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据之后还需对所述地震数据做预处理,具体可包含:对所述地震数据进行静校正和/或强噪声衰减处理等。当然,实际工作中因数据要求不同,本领域相关技术人员也可根据实际需要对所述地震数据做适应性预处理,本发明在此不做过多限制。
在上述步骤S102中对所述观测系统进行稀疏性约束获得不规则化后的激发点和接收点可包含:根据所述项目需求获得预设的接收点和激发点的欠采样比例因子;通过预定稀疏性约束长度将接收线和激发线分别为多个长度为预定稀疏性约束长度的子集;根据所述接收点和激发点的欠采样比例因子与所述子集获得不规则化后的接收点和激发点。在该实施例中,其主体设计原理如下:
(1)根据勘探需求确定接收点的欠采样比例因子P;
(2)假设某条具有N个接收点的接收线,可以表示为:
X={x(1),x(2),...,x(N)};
(3)给出稀疏性约束长度L,将X分成若干个长度为L的子集S,对每个子集Si,根据欠采样比例因子,确定采样位置,则该子集的实际接收点为:
其中,h(m)随机保留M个接收点位置,0表示空掉的接收点;
(4)激发点的不规则方法与上述方法相同;
(5)将不规则化后的激发点和接收点输出为观测系统文件。
亦即,根据所述接收点和激发点的欠采样比例因子与所述子集,并通过下述公式确定采样位置,获得不规则化后的接收点和激发点;
其中,y等于h[m]为采样函数,f表示接收点和激发点,m表示第m个子区,j接收点或激发点的位置序号;r表示一个随机变量,mod表示求余数,r0表示初始值,p表示归一化的随机变量,1表示子区下限,u表示子区上限。
在上述步骤S104中对所述地震数据进行数据重构恢复为规则的地震数据包含:通过最优化法求解以下公式,将所述地震数据进行数据重构恢复为规则的地震数据:
min||x||1s.t.||ΦCTx-b||2<τ;
在上式中s.t.表示约束条件,b不规则的地震数据,x为稀疏向量,Φ为观测矩阵,CT为稀疏变换矩阵的转置,τ为阈值。
实际工作中,对预处理后的不规则地震数据进行重构,将其恢复为规则的地震数据,重构方法为:
(1)将重构过程看作一个线性问题,
b=Ax
其中b不规则的地震数据,A为欠定矩阵或感知矩阵,x为稀疏向量,此公式也可以写为:
b=ΦCTx
其中Φ为观测矩阵,CT为稀疏变换矩阵的转置,稀疏变换可以是曲波变换、傅里叶变换等。
(2)采用最优化方法求解此方程,得到x
min||x||1s.t.||ΦCTx-b||2<τ
(3)最终得到规则化的地震数据
d=CTx
在上述实施例中,根据规则的地震数据计算获得油气储层信息包含:根据规则的地震数据计算获得叠加或偏移剖面;根据所述叠加或偏移剖面计算获得油气储层信息;亦即,对重构后的规则地震数据进行处理,得到偏移剖面,然后对偏移剖面进行综合解释,以获得油气储层信息。
为更清楚的说明本发明所提供的上述基于压缩感知的地震勘探方法,以下通过整体举例,对上述各实施例做完整说明:
本发明所提供的基于压缩感知的地震勘探方法具体可包含:确定勘探需求;常规地震勘探观测系统设计;稀疏性约束的不规则观测系统优化设计;野外地震资料采集与实施;地震资料预处理;基于研所感知的不规则地震数据规则化重构;地震资料成像处理及地震资料综合解释这七步;其中,所述的确定勘探需求,是根据勘探项目情况,确定采用压缩感知地震勘探技术的需求是降低成本、扩大勘探面积或提高空间分辨率;上述各步骤的具体含义及实施方式如下:
1、确定项目需求。根据勘探项目情况,确定采用压缩感知地震勘探技术的需求是降低成本、扩大勘探面积或提高空间分辨率。
2、布设常规观测系统。图2是根据项目需求,设计的常规地震勘探观测系统,包括确定接收线数、接收线距、接收点数、接收点距、激发线数、激发线距、激发点数、激发点距以及接收和激发的起始位置等。
如果是降低成本,只需按图1B的流程继续实施即可;如果是扩大勘探面积,需要根据勘探区域,需要在原观测系统的基础上,增加布设激发线和接收线;如果是提高空间分辨率,则需要缩小激发点距、激发线距,接收点距、接收线距,重新布设炮检点。
3、设计稀疏性约束的不规则观测观测系统。图3是稀疏性约束的不规则观测观测系统,具体不规则观测系统优化设计方法为:
1)根据勘探需求确定接收点的欠采样比例因子P=50%,
2)假设某条具有N个接收点的接收线,可以表示为X={x(1),x(2),...,x(N)}
3)给出稀疏性约束长度L=6,表示约束长度为6个激发点或接收点,将X分成若干个长度为L的子集S,对每个子集Si,根据欠采样比例因子,确定采样位置,则该子集的实际接收点为:
其中,h(m)随机保留M个接收点位置,0表示空掉的接收点。
4)激发点的不规则方法与上述方法相同。
5)将不规则化后的激发点和接收点输出为观测系统文件。
4、野外地震资料采集与实施,是根据设计的不规则观测系统进行野外资料采集,将检波器布设于接收点,采用可控震源或者炸药,在激发点位置进行激发,同时有记录仪器进行记录激发的地震数据。图4是采样不规则观测系统接收到单炮地震记录。
5、地震资料预处理,是对记录的不规则地震数据进行预处理,预处理包括但不限于静校正、强噪声衰减等步骤。
6、不规则地震数据重构,是对预处理后的不规则地震数据进行重构,将其恢复为规则的地震数据,图5是对不规则地震数据数据重构后的规则单炮记录,具体重构方法为:
(1)将重构过程看作一个线性问题,
b=Ax
其中b不规则的地震数据,A为欠定矩阵或感知矩阵,x为稀疏向量,此公式也可以写为:
b=ΦCTx
其中Φ为观测矩阵,CT为稀疏变换矩阵的转置,稀疏变换可以是曲波变换、傅里叶变换等。
(2)采用最优化方法求解此方程,得到x
min||x||1s.t.||ΦCTx-b||2<τ
(3)最终得到规则化的地震数据
d=CTx
7、地震资料成像处理及地震资料综合解释,是对重构后的规则地震数据进行处理,得到叠加或偏移剖面,然后对偏移剖面进行综合解释,以获得油气储层信息;图6是对重构数据处理后的叠加剖面,可以对其进行综合解释,
请参考图7所示,本发明还提供一种基于压缩感知的地震勘探装置,所述装置包含分析模块、处理模块、采集模块和计算模块;所述分析模块用于获取待勘探项目的项目需求,根据项目需求获得预存的对应观测系统;所述处理模块用于对所述观测系统进行稀疏性约束获得不规则化后的激发点和接收点,以所述激发点和所述接收点作为所述观测系统的输出数据;所述采集模块用于通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据;所述计算模块用于对所述地震数据进行数据重构恢复为规则的地震数据,并根据规则的地震数据计算获得油气储层信息。
在上述本发明一实施例中,所述采集模块还包含预处理模块,所述预处理模块用于对所述地震数据进行静校正和/或强噪声衰减处理。所述处理模块还包含:根据所述项目需求获得预设的接收点和激发点的欠采样比例因子;通过预定稀疏性约束长度将接收线和激发线分别为多个长度为预定稀疏性约束长度的子集;根据所述接收点和激发点的欠采样比例因子与所述子集获得不规则化后的接收点和激发点。所述计算模块还包含:通过最优化法求解以下公式,将所述地震数据进行数据重构恢复为规则的地震数据:
min||x||1s.t.||ΦCTx-b||2<τ;
在上式中s.t.表示约束条件,b不规则的地震数据,x为稀疏向量,Φ为观测矩阵,CT为稀疏变换矩阵的转置,τ为阈值。
在上述实施例中,所述计算模块还包含:根据规则的地震数据计算获得叠加或偏移剖面;根据所述叠加或偏移剖面计算获得油气储层信息。
本发明还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述方法。
本发明还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述方法的计算机程序。
本发明的有益技术效果在于:(1)降低勘探成本。在野外采集资料的过程中,应用较少的野外勘探设备进行稀疏采样,完成与常规地震勘探方法相当的勘探效果,从而降低勘探成本;(2)扩大勘探面积。在设备投入和勘探网格不变的情况下,通过扩大勘探区域,利用稀疏采样完成野外资料采集,从而达到在不增加勘探成本的情况下,扩大勘探面积的目的。(3)提高地震勘探的空间分辨率。在不增加勘探成本的情况下,通过缩小炮检点网格,利用稀疏采样完成野外资料采集,从而达到在不增加勘探成本的情况下,提高地震勘探空间分辨率。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (14)
1.一种基于压缩感知的地震勘探方法,其特征在于,所述方法包含:
获取待勘探项目的项目需求,根据项目需求获得预存的对应观测系统;
对所述观测系统进行稀疏性约束获得不规则化后的激发点和接收点,以所述激发点和所述接收点作为所述观测系统的输出数据;
通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据;
对所述地震数据进行数据重构恢复为规则的地震数据,并根据规则的地震数据计算获得油气储层信息。
2.根据权利要求1所述的基于压缩感知的地震勘探方法,其特征在于,通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据还包含:对所述地震数据进行静校正和/或强噪声衰减处理。
3.根据权利要求1所述的基于压缩感知的地震勘探方法,其特征在于,根据项目需求获得预存的对应观测系统中所述观测系统包含:接收线数、接收线距、接收点数、接收点距、激发线数、激发线距、激发点数、激发点距以及接收和激发的起始位置。
5.根据权利要求1所述的基于压缩感知的地震勘探方法,其特征在于,通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据包含:将检波器布设于接收点,通过可控震源或爆炸物于激发点激发,并利用记录仪采集记录激发的地震数据。
6.根据权利要求1所述的基于压缩感知的地震勘探方法,其特征在于,对所述地震数据进行数据重构恢复为规则的地震数据包含:通过最优化法求解以下公式,将所述地震数据进行数据重构恢复为规则的地震数据:
min||x||1s.t.||ΦCTx-b||2<τ;
在上式中s.t.表示约束条件,b不规则的地震数据,x为稀疏向量,Φ为观测矩阵,CT为稀疏变换矩阵的转置,τ为阈值。
7.根据权利要求1所述的基于压缩感知的地震勘探方法,其特征在于,根据规则的地震数据计算获得油气储层信息包含:根据规则的地震数据计算获得叠加或偏移剖面;根据所述叠加或偏移剖面计算获得油气储层信息。
8.一种基于压缩感知的地震勘探装置,其特征在于,所述装置包含分析模块、处理模块、采集模块和计算模块;
所述分析模块用于获取待勘探项目的项目需求,根据项目需求获得预存的对应观测系统;
所述处理模块用于对所述观测系统进行稀疏性约束获得不规则化后的激发点和接收点,以所述激发点和所述接收点作为所述观测系统的输出数据;
所述采集模块用于通过所述观测系统于所述待勘探项目对应区域采集激发产生的地震数据;
所述计算模块用于对所述地震数据进行数据重构恢复为规则的地震数据,并根据规则的地震数据计算获得油气储层信息。
9.根据权利要求8所述的基于压缩感知的地震勘探装置,其特征在于,所述采集模块还包含预处理模块,所述预处理模块用于对所述地震数据进行静校正和/或强噪声衰减处理。
11.根据权利要求8所述的基于压缩感知的地震勘探装置,其特征在于,所述计算模块还包含:通过最优化法求解以下公式,将所述地震数据进行数据重构恢复为规则的地震数据:
min||x||1s.t.||ΦCTx-b||2<τ;
在上式中s.t.表示约束条件,b不规则的地震数据,x为稀疏向量,Φ为观测矩阵,CT为稀疏变换矩阵的转置,τ为阈值。
12.根据权利要求8所述的基于压缩感知的地震勘探装置,其特征在于,所述计算模块还包含:根据规则的地震数据计算获得叠加或偏移剖面;根据所述叠加或偏移剖面计算获得油气储层信息。
13.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7任一所述方法。
14.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至7任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911187300.4A CN112859161B (zh) | 2019-11-28 | 基于压缩感知的地震勘探方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911187300.4A CN112859161B (zh) | 2019-11-28 | 基于压缩感知的地震勘探方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112859161A true CN112859161A (zh) | 2021-05-28 |
CN112859161B CN112859161B (zh) | 2024-09-27 |
Family
ID=
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140372044A1 (en) * | 2013-06-17 | 2014-12-18 | Westerngeco L.L.C. | Seismic data processing |
CN104422954A (zh) * | 2013-08-20 | 2015-03-18 | 中国石油化工股份有限公司 | 一种地震数据规则化中稀疏采样方法 |
CN106291675A (zh) * | 2015-05-22 | 2017-01-04 | 中国石油化工股份有限公司 | 一种基于基追踪技术的地震数据重构方法 |
US20170090053A1 (en) * | 2015-09-28 | 2017-03-30 | Conocophillips Company | 3d seismic acquisition |
CN107817527A (zh) * | 2017-10-31 | 2018-03-20 | 吉林大学 | 基于块稀疏压缩感知的沙漠地震勘探随机噪声压制方法 |
GB201802978D0 (en) * | 2018-02-23 | 2018-04-11 | Statoil Petroleum As | Improved seismic source firing sequence and reciever arrangement |
CN109407143A (zh) * | 2017-08-16 | 2019-03-01 | 中石化石油工程技术服务有限公司 | 基于压缩感知的地震勘探非规则观测系统设计方法 |
US20190129050A1 (en) * | 2016-11-02 | 2019-05-02 | Conocophillips Company | Use nuos technology to acquire optimized 2d data |
RU2017146989A (ru) * | 2017-12-29 | 2019-07-03 | федеральное государственное автономное образовательное учреждение высшего образования "Казанский (Приволжский) федеральный университет" (ФГАОУ ВО КФУ) | Способ сейсмического мониторинга разработки мелкозалегающих залежей сверхвязкой нефти |
RU2018105637A (ru) * | 2018-02-14 | 2019-08-14 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Уральский государственный горный университет" (ФГБОУ ВО "УГГУ") | Способ получения мигрированных сейсмических изображений геологических сред по данным сейсморазведки 2d |
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140372044A1 (en) * | 2013-06-17 | 2014-12-18 | Westerngeco L.L.C. | Seismic data processing |
CN104422954A (zh) * | 2013-08-20 | 2015-03-18 | 中国石油化工股份有限公司 | 一种地震数据规则化中稀疏采样方法 |
CN106291675A (zh) * | 2015-05-22 | 2017-01-04 | 中国石油化工股份有限公司 | 一种基于基追踪技术的地震数据重构方法 |
US20170090053A1 (en) * | 2015-09-28 | 2017-03-30 | Conocophillips Company | 3d seismic acquisition |
US20190129050A1 (en) * | 2016-11-02 | 2019-05-02 | Conocophillips Company | Use nuos technology to acquire optimized 2d data |
CN109407143A (zh) * | 2017-08-16 | 2019-03-01 | 中石化石油工程技术服务有限公司 | 基于压缩感知的地震勘探非规则观测系统设计方法 |
CN107817527A (zh) * | 2017-10-31 | 2018-03-20 | 吉林大学 | 基于块稀疏压缩感知的沙漠地震勘探随机噪声压制方法 |
RU2017146989A (ru) * | 2017-12-29 | 2019-07-03 | федеральное государственное автономное образовательное учреждение высшего образования "Казанский (Приволжский) федеральный университет" (ФГАОУ ВО КФУ) | Способ сейсмического мониторинга разработки мелкозалегающих залежей сверхвязкой нефти |
RU2018105637A (ru) * | 2018-02-14 | 2019-08-14 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Уральский государственный горный университет" (ФГБОУ ВО "УГГУ") | Способ получения мигрированных сейсмических изображений геологических сред по данным сейсморазведки 2d |
GB201802978D0 (en) * | 2018-02-23 | 2018-04-11 | Statoil Petroleum As | Improved seismic source firing sequence and reciever arrangement |
Non-Patent Citations (8)
Title |
---|
GILLES HENNENFENT等: "Simply denoise: Wavefield reconstruction via jittered undersampling", GEOPHYSICS, vol. 73, no. 03, 30 June 2008 (2008-06-30), pages 19 - 28 * |
MEN ZHE等: "A method and application of irregular geometry design based on compressive sensing", SEG TECHNICAL PROGRAM EXPANDED ABSTRACTS, 10 August 2019 (2019-08-10), pages 268 - 272 * |
宁宏晓等: "2015年我国石油地球物理勘探成果述评", 地球物理学进展, vol. 33, no. 01, 28 February 2018 (2018-02-28), pages 385 - 405 * |
张华;陈小宏;: "基于jitter采样和曲波变换的三维地震数据重建", 地球物理学报, vol. 56, no. 05, 31 May 2013 (2013-05-31), pages 1637 - 1649 * |
张良;韩立国;刘争光;李宇;: "基于压缩感知和Contourlet变换的地震数据重建方法", 石油物探, vol. 56, no. 06, pages 805 * |
李鹏等: "基于jittered采样的浅层三维地震数据处理及应用", 地球物理学进展, vol. 32, no. 04, 31 August 2017 (2017-08-31), pages 1784 - 1790 * |
舒国旭;吕公河;吕尧;石太昆;邸志欣;霍守东;: "基于压缩感知的地震数据重建", 石油物探, no. 04, pages 549 - 554 * |
门哲等: "基于压缩感知的随机采样与稀疏性约束采样重构效果对比", 中国石油学会2019年物探技术研讨会文集, pages 1130 - 1132 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11740375B2 (en) | Methods for simultaneous source separation | |
US7599798B2 (en) | Migrating composite seismic response data to produce a representation of a seismic volume | |
US8737165B2 (en) | Interferometric seismic data processing for a towed marine survey | |
US10520623B2 (en) | Methods and systems for marine survey acquisition | |
US20120271550A1 (en) | Seismic Imaging Systems and Methods Employing a 3D Reverse Time Migration with Tilted Transverse Isotropy | |
EP2144090A2 (en) | Flexural wave attenuation | |
CN104062683A (zh) | 一种基于曲波变换和全变差的联合衰减随机噪声处理方法 | |
WO2014165129A1 (en) | Time-lapse monitoring | |
CN114839673B (zh) | 多震源高效采集波场分离方法、分离系统及计算机设备 | |
WO2013150464A1 (en) | Seismic data processing with frequency diverse de-aliasing filtering | |
US20160320507A1 (en) | Time lapse seismic data processing | |
US10330809B2 (en) | Device and method for optimization of 4D and 3D seismic data | |
CN110907989A (zh) | 重建拟地面地震反射波成像方法及系统 | |
CN112859161A (zh) | 基于压缩感知的地震勘探方法及装置 | |
CN105319594B (zh) | 一种基于最小二乘参数反演的傅里叶域地震数据重构方法 | |
US20230051004A1 (en) | Method and apparatus for performing efficient modeling of extended-duration moving seismic sources | |
US20180095186A1 (en) | Noise models by selection of transform coefficients | |
CN113743193B (zh) | 一种叠前地震资料线性干扰压制方法及系统 | |
US20080232195A1 (en) | Apparatus and Method for Processing Geophysical Information | |
CN104166160B (zh) | 一种地震数据组合方法、装置及系统 | |
CN116413801B (zh) | 一种各向异性介质弹性波高精度成像方法和系统 | |
US20220413175A1 (en) | Enhanced projection on convex sets for interpolation and deblending | |
Falowo et al. | The Application of Wavelet Deconvolution for Noise Reduction in Seismic Data | |
Zizi et al. | Low‐frequency seismic deghosting in a compressed domain using parabolic dictionary learning | |
CN116070055A (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 |