CN112904415B - 基于重构源谱一致性的q层析反演方法、系统及处理设备 - Google Patents
基于重构源谱一致性的q层析反演方法、系统及处理设备 Download PDFInfo
- Publication number
- CN112904415B CN112904415B CN202110078478.6A CN202110078478A CN112904415B CN 112904415 B CN112904415 B CN 112904415B CN 202110078478 A CN202110078478 A CN 202110078478A CN 112904415 B CN112904415 B CN 112904415B
- Authority
- CN
- China
- Prior art keywords
- seismic
- reconstructed
- layer
- consistency
- detection point
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000004587 chromatography analysis Methods 0.000 title claims abstract description 16
- 238000001514 detection method Methods 0.000 claims abstract description 37
- 238000010521 absorption reaction Methods 0.000 claims abstract description 29
- 238000012512 characterization method Methods 0.000 claims abstract description 17
- 230000001902 propagating effect Effects 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 16
- 238000003325 tomography Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 239000000523 sample Substances 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000036039 immunity Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation 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/282—Application of seismic models, synthetic seismograms
-
- 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/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
-
- 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
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种基于重构源谱一致性的Q层析反演方法、系统及处理设备,包括以下步骤:S1:获取地震波传播的射线路径,得到地震波穿过每层介质的旅行时Δtij;S2:根据地震波穿过每层介质的旅行时Δtij和给定的层Q模型,计算地震波在地下介质中传播的等效吸收因子ai(f,Q);S3:根据地震波穿过地下地层的等效吸收衰减因子ai(f,Q),确定每个检波点接收的信号的衰减量exp[‑ai(f,Q)],对每个检波点接收的信号进行反向吸收补偿,重构出震源振幅谱Si(f,Q);S4:根据每个检波点处重构的震源振幅谱Si(f,Q),建立重构源谱一致性表征函数J(Q);S5:重复步骤S2至S4,遍历所有的Q模型,得到所有Q模型对应的重构源谱一致性表征函数值,通过比较确定所求的Q模型。
Description
技术领域
本发明涉及一种基于重构源谱一致性的Q层析反演方法、系统、处理设备及存储介质,属于地球物理勘探中的地震信息处理技术领域。
背景技术
地震波在粘弹性介质中传播时会遭受地层的吸收衰减作用,降低了地震资料的分辨率,畸变了地震信号在油气储层上的反射特征,地层的粘弹性效应已经日益成为制约地震资料处理质量的主要障碍。因此,利用地震数据对地下地层的吸收结构进行反演估算,并建立可靠地层吸收结构模型,据此对地震波的吸收进行衰减补偿,已经成为提高地震资料处理质量的关键技术。
目前,品质因子Q的估算方法有十余种,按作用域可分为时间域、频率域和时频域三类方法。1974年,Stacey与Gladwin分析了地震波在传播过程中子波延续时间变长的现象及脉冲增宽现象,提出上升时间法,该方法只需要测量脉冲宽度和两个观测点之间的走时就能求取品质因子。1974年,Bath首次提出利用参考子波衰减前后的谱比来估计Q值,开创了频率域估算Q值的历史。1991年,Tonn利用VSP数据比较了十几种Q值估算方法,发现没有哪一种方法能得到比较满意的结果,其中谱比法的优势是理论精度高,缺点是对噪音过于敏感。1997年,Quan和Harris根据子波衰减前后质心频率的变化估算Q值,该方法利用了频谱的统计学特性,具有良好的抗噪性,且假设震源子波满足高斯分布。2008年,高静怀推导了瞬时频率与品质因子的解析关系,在小波域中发展了一种计算瞬时频率的方法,并以此为基础提出了估算VSP资料的小波域包络峰值处瞬时频率法。2016年,李国发等通过对微测井采集过程中的影响因素分析,提出了新的微测井观测系统,该观测系统可以有效的避免激发因素、检波点耦合等对Q值估算的影响,李国发等应用该观测系统进行了谱比法衰减层析反演。
在众多方法中,谱比法和质心频移法是目前最为流行的Q值反演方法,这两种方法各有优缺点。谱比法的优点在于对震源子波没有前提假设,且理论精度较高,缺点是抗噪性较差;质心频移法的优点是算法稳定、抗噪性较强,但是该方法假设震源子波振幅谱满足高斯分布,这个假设在实际中一般是不成立的,这也在一定程度上降低了该方法的精度。
发明内容
针对上述问题,本发明的目的是提供一种基于重构源谱一致性的Q层析反演方法、系统、处理设备及存储介质,该方法无需已知震源子波,无需对子波进行任何假设,算法理论精度高,并且利用统计学参数来比较重构的震源振幅谱的一致性,具有较好的抗噪性。
为实现上述目的,本发明采用以下技术方案:一种基于重构源谱一致性的Q层析反演方法,包括以下步骤:
S1:获取地震波传播的射线路径和速度模型,计算得到地震波穿过每层介质的旅行时Δtij;
S2:根据地震波穿过每层介质的旅行时Δtij和给定的层Q模型,计算地震波在地下介质中传播的等效吸收因子ai(f,Q),
式中,n表示地层的分层数;Qj表示每层介质的品质因子;i表示第i个检波点,f表示频率;
S3:根据地震波穿过地下地层的等效吸收衰减因子ai(f,Q),确定每个检波点接收的信号的衰减量exp[-ai(f,Q)],进而对每个检波点接收的信号进行反向吸收补偿,重构出震源振幅谱Si(f,Q),
Si(f,Q)=Ri(f)exp(ai(f,Q))
式中,Ri(f)表示第i个检波点接收信号的振幅谱;
S4:根据每个检波点处重构的震源振幅谱Si(f,Q),建立重构源谱一致性表征函数J(Q),
式中,N表示检波点个数,Sa(f,Q)表示每个检测点的平均振幅谱;
S5:重复步骤S2至S4,遍历所有可能的Q模型,得到所有Q模型对应的重构源谱一致性表征函数值,通过比较确定震源振幅谱的一致性最好时对应的Q模型,即为所求的Q模型。
在一些实施例中,在上述步骤S1中,根据微测井地震数据获取初至波旅行时曲线,对初至波旅行时曲线进行解释,获得速度模型;根据实际采集观测系统和速度模型,利用试射法进行射线追踪正演模拟,得到地震波传播的射线路径,进而确定地震波从炮点到检波点穿过每层介质的旅行时Δtij。
在一些实施例中,在上述步骤S4中,平均振幅谱Sa(f,Q)计算公式如下:
本发明还提供一种基于重构源谱一致性的Q层析反演系统,包括:
旅行时计算模块,被配置为根据获取的速度模型和地震波传播的射线路径,计算地震波穿过每层介质的旅行时;
等效吸收衰减因子计算模块,被配置为根据地震波穿过每层介质的旅行时和给定的Q模型,计算地震波在地下介质中传播的等效吸收衰减因子;
震源振幅谱重构模块,被配置为根据地震波穿过地下地层的等效吸收衰减因子,确定每个检波点接收的信号的衰减量,再对每个检波点接收的信号进行反向吸收补偿,重构出震源振幅谱;
重构源谱一致性表征函数建立模块,被配置为根据每个检波点处重构的震源振幅谱建立重构源谱一致性表征函数;
Q模型确定模块,被配置为通过比较所有给定Q模型对应的重构源谱一致性表征函数值,确定所求的Q模型。
一种处理设备,包括存储器、处理器以及存储在所述存储器中计算机程序,所述处理器执行所述计算机程序时实现上述实施例任一项所述的Q层析反演方法。
一种存储介质,存储介质上存储计算机程序,该程序被处理器执行时实现上述实施例中任一项所述的Q层析反演方法。
本发明采用以上技术方案,其具有如下优点:本发明提供的基于重构源谱一致性的Q层析反演方法,根据地震波穿过每层介质的旅行时Δtij以及给定的Q模型,计算地震波在地下介质中传播的等效吸收衰减因子,进而确定每个检波点接收的信号的衰减量,根据衰减量对每个检波点接收的信号进行反向吸收补偿,重构出震源振幅谱,依据震源振幅谱构建重构源谱一致性表征函数J(Q),采用遍历法遍历所有可能的Q模型,通过比较得到的J(Q)确定所求的Q模型,无需已知震源子波,也不需要对子波进行任何假设,算法理论精度很高;此外,该方法利用统计学参数来比较重构的震源振幅谱的一致性,具有较好的抗噪性;利用该方法可以建立可靠地层吸收结构模型,为后期的地层吸收补偿和油气识别提供高质量的基础数据,是提高地震资料处理质量的关键技术。
附图说明
图1为本发明微测井单炮记录的示意图;
图2为本发明初至波旅行时的散点图;
图3为本发明获取的速度模型示意图;
图4为本发明的射线追踪正演模拟的示意图;
图5为本发明给定的Q模型的示意图;
图6为本发明中合成的衰减地震记录示意图;
图7为本发明每个检波点重构出的震源振幅谱;
图8为本发明中重构的震源谱一致性特征示意图;
图9为本发明利用遍历法得到的Q模型示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“上”、“下”、“前”、“后”“内”、“外”、“横”、“竖”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的系统或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,使用术语“第一”、“第二”等词语来限定零部件,仅仅是为了便于对上述零部件进行区别,如没有另行声明,上述词语并没有特殊含义,不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“设置”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
本发明提供一种基于重构源谱一致性的Q层析反演方法,其包括以下步骤:
S1:获取速度模型和地震波传播的射线路径,计算地震波穿过每层介质的旅行时Δtij;
S2:根据地震波穿过每层介质的旅行时Δtij和给定的Q模型,计算地震波在地下介质中传播的等效吸收衰减因子ai(f,Q),
式中,n表示地层的分层数;Qj表示每层介质的品质因子;i表示第i个检波点,f表示频率;
S3:根据地震波穿过地下地层的等效吸收衰减因子ai(f,Q),确定每个检波点接收的信号的衰减量exp[-ai(f,Q)],进而对每个检波点接收的信号进行反向吸收补偿,重构出震源振幅谱Si(f,Q),
Si(f,Q)=Ri(f)exp(ai(f,Q))
式中,Ri(f)表示第i个检波点接收信号的振幅谱;
S4:根据每个检波点处重构的震源振幅谱Si(f,Q),建立重构源谱一致性表征函数J(Q),
式中,N表示检波点个数,Sa(f,Q)表示每个检测点的平均振幅谱;
S5:重复步骤S2至S4,遍历所有可能的Q模型,得到所有Q模型对应的重构源谱一致性表征函数值,通过比较确定震源振幅谱的一致性最好时对应的Q模型,即为所求的Q模型。
进一步地,在上述步骤S1中,根据微测井地震数据获取初至波旅行时曲线,并对初至波旅行时曲线进行解释,获得速度模型;根据实际采集观测系统和速度模型,利用试射法进行射线追踪正演模拟,得到地震波传播的射线路径,进而确定地震波从炮点到检波点穿过每层介质的旅行时Δtij。
进一步地,在上述步骤S4中,平均振幅谱Sa(f,Q)计算公式如下:
下面以具体实施例对本发明作以说明:
步骤1:获取速度模型。
根据微测井单炮记录(如图1所示)拾取初至波旅行时,拟合形成初至波旅行时分段折线(如图2所示),根据每段折线的斜率得到每一层的速度,折线转折处即为分层界面,从而确定速度模型(如图3所示);
步骤2:如图4所示,根据获取的速度模型和实际采集观测系统,利用射线追踪正演模拟可以得到地震波传播的射线路径,进而计算得到地震波穿过每层介质的旅行时Δtij;
步骤4:如图7所示,根据每个检波点的衰减量因子,对每个检波点接收的信号进行反向吸收补偿,重构出震源振幅谱Si(f,Q)=Ri(f)exp(ai(f,Q));
理论上来讲,如果给定的Q模型是准确的,则每个检波点经历的衰减量因子和补偿量是准确的,进而反向补偿后重构出的震源振幅谱肯定的一致的(因为检波点接收到的信号是从同一个震源激发的,反过来补偿回去肯定一样)如图8所示。如果重构出的震源振幅谱一致较差,那肯定是由于给定的Q模型不对。
步骤5:建立重构源谱一致性表征函数J(Q),这个函数的作用就是定量描述重构源谱的一致性好坏。
J(Q)越小,说明重构震源振幅谱的一致性越好;
步骤6:遍历所有可能的Q模型,重复步骤3至步骤5,得到所有Q模型对应的重构源谱一致性表征函数值,通过比较确定源谱的一致性最好时对应的Q模型,即为所求的Q模型,如图9所示。
本公开一实施例还提供一种基于重构源谱一致性的Q层析反演系统,包括:
旅行时计算模块,被配置为根据获取的速度模型和地震波传播的射线路径,计算地震波穿过每层介质的旅行时;
等效吸收衰减因子计算模块,被配置为根据地震波穿过每层介质的旅行时和给定的Q模型,计算地震波在地下介质中传播的等效吸收衰减因子;
震源振幅谱重构模块,被配置为根据地震波穿过地下地层的等效吸收衰减因子,确定每个检波点接收的信号的衰减量,再对每个检波点接收的信号进行反向吸收补偿,重构出震源振幅谱;
重构源谱一致性表征函数建立模块,被配置为根据每个检波点处重构的震源振幅谱建立重构源谱一致性表征函数;
Q模型确定模块,被配置为通过比较所有给定Q模型对应的重构源谱一致性表征函数值,确定所求的Q模型。
本公开一实施例还提供一种处理设备,包括存储器、处理器以及存储在存储器中计算机程序,处理器执行计算机程序时实现上述任一实施例中的Q层析反演方法。
另外地,本公开一实施例还提供一种存储介质,其上存储计算机程序,该程序被处理器执行时实现上述任一实施例中的Q层析反演方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.一种基于重构源谱一致性的Q层析反演方法,其特征在于,包括以下步骤:
S1:获取地震波传播的射线路径和速度模型,计算得到地震波穿过每层介质的旅行时Δtij;
S2:根据地震波穿过每层介质的旅行时Δtij和给定的层Q模型,计算地震波在地下介质中传播的等效吸收因子ai(f,Q),
式中,n表示地层的分层数;Qj表示每层介质的品质因子;i表示第i个检波点,f表示频率;
S3:根据地震波穿过地下地层的等效吸收衰减因子ai(f,Q),确定每个检波点接收的信号的衰减量exp[-ai(f,Q)],进而对每个检波点接收的信号进行反向吸收补偿,重构出震源振幅谱Si(f,Q),
Si(f,Q)=Ri(f)exp(ai(f,Q))
式中,Ri(f)表示第i个检波点接收信号的振幅谱;
S4:根据每个检波点处重构的震源振幅谱Si(f,Q),建立重构源谱一致性表征函数J(Q),
式中,N表示检波点个数,Sa(f,Q)表示每个检测点的平均振幅谱;
S5:重复步骤S2至S4,遍历所有可能的Q模型,得到所有Q模型对应的重构源谱一致性表征函数值,通过比较确定震源振幅谱的一致性最好时对应的Q模型,即为所求的Q模型。
2.如权利要求1所述的基于重构源谱一致性的Q层析反演方法,其特征在于:在上述步骤S1中,根据微测井地震数据获取初至波旅行时曲线,对初至波旅行时曲线进行解释,获得速度模型;根据实际采集观测系统和速度模型,利用试射法进行射线追踪正演模拟,得到地震波传播的射线路径,进而确定地震波从炮点到检波点穿过每层介质的旅行时Δtij。
4.一种处理设备,其特征在于:包括存储器、处理器以及存储在所述存储器中计算机程序,所述处理器执行所述计算机程序时实现如权利要求1至3任一项所述的Q层析反演方法。
5.一种存储介质,其特征在于:存储介质上存储计算机程序,该程序被处理器执行时实现如权利要求1至3任一项所述的Q层析反演方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110078478.6A CN112904415B (zh) | 2021-01-20 | 2021-01-20 | 基于重构源谱一致性的q层析反演方法、系统及处理设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110078478.6A CN112904415B (zh) | 2021-01-20 | 2021-01-20 | 基于重构源谱一致性的q层析反演方法、系统及处理设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112904415A CN112904415A (zh) | 2021-06-04 |
CN112904415B true CN112904415B (zh) | 2022-03-15 |
Family
ID=76117497
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110078478.6A Active CN112904415B (zh) | 2021-01-20 | 2021-01-20 | 基于重构源谱一致性的q层析反演方法、系统及处理设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112904415B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115793058B (zh) * | 2023-02-08 | 2023-04-18 | 成都理工大学 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
CN117724166A (zh) * | 2024-02-07 | 2024-03-19 | 中国石油大学(华东) | 基于大炮初至的近地表三维速度建模方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1467509A (zh) * | 2002-07-12 | 2004-01-14 | 中国石油集团东方地球物理勘探有限责 | 一种时频域大地吸收衰减补偿方法 |
CN102053273A (zh) * | 2009-10-29 | 2011-05-11 | 中国石油化工股份有限公司 | 一种对地震波信号进行反q滤波的方法 |
CN104570108A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 估算等效品质因子方法及用其估算地层品质因子的方法 |
CN106547019A (zh) * | 2015-09-17 | 2017-03-29 | 中国石油化工股份有限公司 | 一种确定地层品质因子的方法 |
CN108845357A (zh) * | 2018-06-13 | 2018-11-20 | 成都信息工程大学 | 一种基于同步挤压小波变换估计地层等效品质因子的方法 |
-
2021
- 2021-01-20 CN CN202110078478.6A patent/CN112904415B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1467509A (zh) * | 2002-07-12 | 2004-01-14 | 中国石油集团东方地球物理勘探有限责 | 一种时频域大地吸收衰减补偿方法 |
CN102053273A (zh) * | 2009-10-29 | 2011-05-11 | 中国石油化工股份有限公司 | 一种对地震波信号进行反q滤波的方法 |
CN104570108A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 估算等效品质因子方法及用其估算地层品质因子的方法 |
CN106547019A (zh) * | 2015-09-17 | 2017-03-29 | 中国石油化工股份有限公司 | 一种确定地层品质因子的方法 |
CN108845357A (zh) * | 2018-06-13 | 2018-11-20 | 成都信息工程大学 | 一种基于同步挤压小波变换估计地层等效品质因子的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112904415A (zh) | 2021-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111208561B (zh) | 基于时变子波与曲波变换约束的地震声波阻抗反演方法 | |
US10705238B2 (en) | Method and apparatus for processing seismic data | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
CN112904415B (zh) | 基于重构源谱一致性的q层析反演方法、系统及处理设备 | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
CN107219554B (zh) | 陆地地震资料的剩余静校正量的自动获取方法 | |
FR2965066A1 (fr) | Procede d'estimation de parametres elastiques par inversion de mesures sismiques 4d | |
WO2015100544A1 (zh) | 基于零偏垂直地震剖面数据估计品质因子的方法和装置 | |
US8929175B2 (en) | Method of marine time-lapse seismic surveying | |
WO2012010790A2 (fr) | Procede d'estimation de parametres elastiques par inversion de mesures sismiques 4d | |
CN111060961B (zh) | 基于多信息约束反演的品质因子确定方法、装置及系统 | |
CN110967749A (zh) | 一种vsp地震资料频变q值估计与反q滤波方法 | |
Liu et al. | Seismic quality factor estimation using frequency-dependent linear fitting | |
KR101693310B1 (ko) | 평면파를 이용한 이방성 매질에서의 전파형 역산방법 | |
CN109143345B (zh) | 基于模拟退火的品质因子q非线性反演方法及系统 | |
CN110244383B (zh) | 基于近地表数据的地质岩性综合模型创建方法 | |
Chen et al. | Improving the Precision of Surface Seismic Data Processing by Walkaway VSP | |
CN112213774B (zh) | 一种浅层q模型估计方法及装置 | |
CN114740528A (zh) | 一种超微分拉普拉斯块约束的叠前多波联合反演方法 | |
CN112213784B (zh) | 复杂地表地震数据一次处理快速静校正方法 | |
CN112415601A (zh) | 表层品质因子q值的确定方法及装置 | |
CN114428343A (zh) | 基于归一化互相关的Marchenko成像方法及系统 | |
Yang et al. | Seismic attenuation estimation from instantaneous frequency | |
CN114114406B (zh) | 一种储层渗透率估计方法及装置 | |
CN110850504B (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 |