CN113835120B - 不同激发能量大小的正演模拟分析方法 - Google Patents
不同激发能量大小的正演模拟分析方法 Download PDFInfo
- Publication number
- CN113835120B CN113835120B CN202010593833.9A CN202010593833A CN113835120B CN 113835120 B CN113835120 B CN 113835120B CN 202010593833 A CN202010593833 A CN 202010593833A CN 113835120 B CN113835120 B CN 113835120B
- Authority
- CN
- China
- Prior art keywords
- energy
- cannon
- work area
- forward modeling
- noise ratio
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 35
- 230000005284 excitation Effects 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 44
- 238000003384 imaging method Methods 0.000 claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims description 22
- 238000011156 evaluation Methods 0.000 claims description 4
- 230000002093 peripheral effect Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 10
- 238000010276 construction Methods 0.000 abstract description 6
- 239000003814 drug Substances 0.000 description 7
- 229940079593 drug Drugs 0.000 description 5
- 230000008859 change Effects 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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. analysis, for interpretation, for correction
- 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
- G01V2210/624—Reservoir parameters
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Abstract
本发明提供一种不同激发能量大小的正演模拟分析方法,该不同激发能量大小的正演模拟分析方法包括:步骤1,收集工区不同药量单炮资料;步骤2,提取不同药量单炮能量信噪比信息;步骤3,根据工区地质构造情况和成像目标需求,建立地质模型进行正演模拟;步骤4,将模拟单炮乘以相应的能量系数;步骤5,采用模拟单炮进行处理分析。该不同激发能量大小的正演模拟分析方法获得的模拟单炮能够更真实的反映实际单炮情况,对于指导地震采集方法具有重要意义,克服了目前正演模拟分析中单炮能量与实际差别大的问题,对于西部复杂地表及东部城区施工具有重要的指导作用。
Description
技术领域
本发明涉及油田开发技术领域,特别是涉及到一种不同激发能量大小的正演模拟分析方法。
背景技术
地震正演模拟分析是分析地震波场特征,在正式地震采集前是指导地震勘探方法重要的方法,可以预知采用该方法到底能够获得什么样的资料,以及预知能不能完成地质任务。但是受目前正演模拟方法的限制,很难做到与实际资料相似的程度,如果与实际情况相差太远,则失去了预知地震勘探效果,指导地震勘探方法的作用,从而得到错误的分析结论。目前地震模拟方法中非常重要的问题是,模拟出来的单炮能量都是一样的,而信噪比也与实际单炮相差很大。而实际施工中,收到地表各种障碍物的影响,不得不采用各种药量进行激发,从而使得单炮的能量和信噪比差异较大。
目前地震正演模拟过程中,模拟的单炮能量、信噪比等信息与实际单炮相差较大,从而无法真正采用正演模拟单炮对地震采集方法进行分析评价。地震勘探正演模拟模拟方法对于指导地震采集方法非常重要,野外地震采集成本较大,而一旦采集完成将无法弥补由于方法的缺陷造成的成像损失,而正演模拟分析就能够发挥较大的作用,从室内就可以预知各种地震采集方法可能获得的成像效果,从而对地震采集方法进行优化。而目前正演模拟方法还存在模拟能量无法满足实际需求的问题。
为此我们发明了一种新的不同激发能量大小的正演模拟分析方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种获得的模拟单炮能够更真实的反映实际单炮情况的不同激发能量大小的正演模拟分析方法。
本发明的目的可通过如下技术措施来实现:不同激发能量大小的正演模拟分析方法,该不同激发能量大小的正演模拟分析方法包括:包括:步骤1,收集工区不同药量单炮资料;步骤2,提取不同药量单炮能量信噪比信息;步骤3,根据工区地质构造情况和成像目标需求,建立地质模型进行正演模拟;步骤4,将模拟单炮乘以相应的能量系数;步骤5,采用模拟单炮进行处理分析。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,收集工区中以往不同药量试验资料,若该工区中没有以往的不同药量试验资料,则收集周边相似工区的不同药量试验资料,若没有周边资料,则在工区中某个位置进行不同药量的试验,从小药量一直测试到大药量,要求等间隔试验。
在步骤2中,在获得不同药量单炮资料之后,分别提取每种药量的近道初至能量,计算出这些单炮能量的比例,同时分别提取每种药量单炮目的层位置的信噪比。
在步骤3中,提取工区单炮资料中的目的层子波,采用该子波作为正演模拟的激发子波进行正演模拟。
在步骤4中,在获得模拟单炮后,利用前面分析得到的不同药量的能量比值对单炮乘以相应的增益因子,改变单炮的能量大小,在调整完增益大小后,根据前面分析的不同药量单炮目的层的信噪比大小,增加相应的噪音能量,使得模拟获得的某个药量单炮的信噪比与实际药量单炮信噪比相当;按照这样的方法增加相应的噪音,使得模拟单炮信噪比与实际单炮对应。
在步骤5,在处理完模拟单炮的信噪比和能量后,对模拟单炮进行处理成像,获得模拟剖面评价预将采用的地震采集方法的优劣,优化地震采集方法。
本发明中的不同激发能量大小的正演模拟分析方法,是一种充分考虑实际地震波场能量特征,根据工区中可能遇到的激发药量而获得的单炮情况,对模拟数据进行调整,从而更加真实的获得实际资料情况,并通过该模拟单炮可以分析各种观测系统在地震勘探施工后可以获得什么样的成像效果,从而指导地震采集方法的调整。经实际资料验证了该方法的合理性,采用该方法获得的模拟单炮能够更真实的反映实际单炮情况,对于指导地震采集方法具有重要意义,克服了目前正演模拟分析中单炮能量与实际差别大的问题,对于西部复杂地表及东部城区施工具有重要的指导作用。
附图说明
图1为本发明的不同激发能量大小的正演模拟分析方法的一具体实施例的流程图;
图2为本发明的一具体实施例中利用本方法获得的观测系统图。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明提供一种不同激发能量大小的正演模拟分析方法。其方法包括根据工区中可能遇到的激发药量而获得的单炮情况,对模拟数据进行调整,从而更加真实的获得实际资料情况,并通过该模拟单炮可以分析各种观测系统在地震勘探施工后可以获得什么样的成像效果,从而指导地震采集方法的调整。经实际资料验证了该方法的合理性,采用该方法获得的模拟单炮能够更真实的反映实际单炮情况,对于指导地震采集方法具有重要意义,克服了目前正演模拟分析中单炮能量与实际差别大的问题,对于西部复杂地表及东部城区施工具有重要的指导作用。
本发明的不同激发能量大小的正演模拟分析方法包括以下步骤:
步骤1,首先收集工区中以往不同药量试验资料,例如药量从1kg,2kg………8kg,10kg等不同药量的单炮资料,如果该工区中没有以往的不同药量试验资料,可以收集周边相似工区的不同药量试验资料,如果周边资料也没有,就需要在工区中某个位置进行不同药量的,从小药量一直测试到大药量,要求等间隔试验,例如1kg,2kg………8kg等等;
步骤2,在获得不同药量之后,分别提取每种药量的近道初至能量,计算出这些单炮能量的比例,同时分别提取每种药量单炮目的层位置的信噪比;
步骤3,根据工区地质构造情况和成像目标需求,建立地质模型进行正演模拟。提取工区单炮资料中的目的层子波,采用该子波作为正演模拟的激发子波进行正演模拟;
步骤4,在获得模拟单炮后,将单炮输出sgy数据,输入到处理系统中,利用前面分析得到的不同药量的能量比值对单炮乘以相应的增益因子,改变单炮的能量大小。例如6kg药量为一个工区正常生产单炮的激发因素,其增益为1,而低于6kg的药量则增益小于1,大于6kg的药量单炮增益大于1,其具体数值大小根据前面分析的结果确定。在调整完增益大小后,根据前面分析的不同药量单炮目的层的信噪比大小,增加相应的噪音能量,使得模拟获得的某个药量单炮的信噪比与实际药量单炮信噪比相当。例如信噪比为1.2的6kg实际单炮,对应的模拟单炮能量为1,应该增加的随机噪音能量为1/1.2=0.833。按照这样的方法增加相应的噪音,使得模拟单炮信噪比与实际单炮对应;
步骤5,在处理完模拟单炮的信噪比和能量后,对模拟单炮进行处理成像,获得模拟剖面评价预将采用的地震采集方法的优劣,优化地震采集方法。
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
为使本发明的目的、特征和优点能更明显易懂,下文配合所附图式,作详细说明如下:
如图1所示,图1为一种波场连续采集的观测系统设计方法的流程图。
在步骤101,首先收集工区中以往不同药量试验资料,例如药量从1kg,2kg………8kg,10kg等不同药量的单炮资料,如果该工区中没有以往的不同药量试验资料,可以收集周边相似工区的不同药量试验资料,如果周边资料也没有,就需要在工区中某个位置进行不同药量的,从小药量一直测试到大药量,要求等间隔试验,例如1kg,2kg………8kg等等。流程进入到步骤102。
在步骤102,在获得不同药量之后,分别提取每种药量的近道初至能量,计算出这些单炮能量的比例,同时分别提取每种药量单炮目的层位置的信噪比。进入到步骤103。
在步骤103,根据工区地质构造情况和成像目标需求,建立地质模型进行正演模拟。提取工区单炮资料中的目的层子波,采用该子波作为正演模拟的激发子波进行正演模拟。流程进入到步骤104。
在步骤104,在获得模拟单炮后,将单炮输出sgy数据,输入到处理系统中,利用前面分析得到的不同药量的能量比值对单炮乘以相应的增益因子,改变单炮的能量大小。例如6kg药量为一个工区正常生产单炮的激发因素,其增益为1,而低于6kg的药量则增益小于1,大于6kg的药量单炮增益大于1,其具体数值大小根据前面分析的结果确定。在调整完增益大小后,根据前面分析的不同药量单炮目的层的信噪比大小,增加相应的噪音能量,使得模拟获得的某个药量单炮的信噪比与实际药量单炮信噪比相当。例如信噪比为1.2的6kg实际单炮,对应的模拟单炮能量为1,应该增加的随机噪音能量为1/1.2=0.833。按照这样的方法增加相应的噪音,使得模拟单炮信噪比与实际单炮对应。流程进入到步骤105。
在步骤105,在处理完模拟单炮的信噪比和能量后,对模拟单炮进行处理成像,获得模拟剖面评价预将采用的地震采集方法的优劣,优化地震采集方法。流程结束。
图2所示是利用该方法计算获得的不同药量大小的模拟单炮。
本发明根据实际资料中不同药量单炮进行分析,分析出不同药量单炮的能量和信噪比,并将该分析结论应用到正演模拟单炮中,从而获得更加真实的模拟单炮,然后采用该单炮进行处理成像,从成像效果来判断地震采集方法的优劣。采用这种方法可以获得更加逼近真实情况的模拟单炮,对于模拟分析实际工区障碍物对成像的影响具有较好的指导意义。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (5)
1.不同激发能量大小的正演模拟分析方法,其特征在于,该不同激发能量大小的正演模拟分析方法包括:
步骤1,收集工区不同药量单炮资料;
步骤2,提取不同药量单炮能量信噪比信息;
步骤3,根据工区地质构造情况和成像目标需求,建立地质模型进行正演模拟;
步骤4,将模拟单炮乘以相应的能量系数;
步骤5,采用模拟单炮进行处理分析;
在步骤4中,在获得模拟单炮后,利用前面分析得到的不同药量的能量比值对单炮乘以相应的增益因子,改变单炮的能量大小,在调整完增益大小后,根据前面分析的不同药量单炮目的层的信噪比大小,增加相应的噪音能量,使得模拟获得的某个药量单炮的信噪比与实际药量单炮信噪比相当;按照这样的方法增加相应的噪音,使得模拟单炮信噪比与实际单炮对应。
2.根据权利要求1所述的不同激发能量大小的正演模拟分析方法,其特征在于,在步骤1中,收集工区中以往不同药量试验资料,若该工区中没有以往的不同药量试验资料,则收集周边相似工区的不同药量试验资料,若没有周边资料,则在工区中某个位置进行不同药量的试验,从小药量一直测试到大药量,要求等间隔试验。
3.根据权利要求1所述的不同激发能量大小的正演模拟分析方法,其特征在于,在步骤2中,在获得不同药量单炮资料之后,分别提取每种药量的近道初至能量,计算出这些单炮能量的比例,同时分别提取每种药量单炮目的层位置的信噪比。
4.根据权利要求1所述的不同激发能量大小的正演模拟分析方法,其特征在于,在步骤3中,提取工区单炮资料中的目的层子波,采用该子波作为正演模拟的激发子波进行正演模拟。
5.根据权利要求1所述的不同激发能量大小的正演模拟分析方法,其特征在于,在步骤5,在处理完模拟单炮的信噪比和能量后,对模拟单炮进行处理成像,获得模拟剖面评价预将采用的地震采集方法的优劣,优化地震采集方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010593833.9A CN113835120B (zh) | 2020-06-24 | 2020-06-24 | 不同激发能量大小的正演模拟分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010593833.9A CN113835120B (zh) | 2020-06-24 | 2020-06-24 | 不同激发能量大小的正演模拟分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113835120A CN113835120A (zh) | 2021-12-24 |
CN113835120B true CN113835120B (zh) | 2024-02-09 |
Family
ID=78965050
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010593833.9A Active CN113835120B (zh) | 2020-06-24 | 2020-06-24 | 不同激发能量大小的正演模拟分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113835120B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012145113A1 (en) * | 2011-04-18 | 2012-10-26 | Siemens Corporation | Methods for deblending of seismic shot gathers |
CN103543465A (zh) * | 2012-07-12 | 2014-01-29 | 中国石油天然气集团公司 | 基于目的层照明能量确定有效激发点的方法 |
CN105093296A (zh) * | 2015-06-30 | 2015-11-25 | 中国石油天然气集团公司 | 一种优化观测系统的方法及装置 |
CN107526101A (zh) * | 2017-07-07 | 2017-12-29 | 中国地质调查局油气资源调查中心 | 一种获取地震反射波的采集和处理方法 |
CN110082820A (zh) * | 2018-01-26 | 2019-08-02 | 中石化石油工程技术服务有限公司 | 炸药震源混合分布式宽频激发的方法 |
CN110501757A (zh) * | 2019-06-27 | 2019-11-26 | 中国石油化工股份有限公司 | 点线面三维逐点井深设计方法 |
-
2020
- 2020-06-24 CN CN202010593833.9A patent/CN113835120B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012145113A1 (en) * | 2011-04-18 | 2012-10-26 | Siemens Corporation | Methods for deblending of seismic shot gathers |
CN103543465A (zh) * | 2012-07-12 | 2014-01-29 | 中国石油天然气集团公司 | 基于目的层照明能量确定有效激发点的方法 |
CN105093296A (zh) * | 2015-06-30 | 2015-11-25 | 中国石油天然气集团公司 | 一种优化观测系统的方法及装置 |
CN107526101A (zh) * | 2017-07-07 | 2017-12-29 | 中国地质调查局油气资源调查中心 | 一种获取地震反射波的采集和处理方法 |
CN110082820A (zh) * | 2018-01-26 | 2019-08-02 | 中石化石油工程技术服务有限公司 | 炸药震源混合分布式宽频激发的方法 |
CN110501757A (zh) * | 2019-06-27 | 2019-11-26 | 中国石油化工股份有限公司 | 点线面三维逐点井深设计方法 |
Non-Patent Citations (2)
Title |
---|
基于目的层照明能量的炮点设计方法;赵虎;尹成;李瑞;徐峰;陈光明;李轩波;;石油物探(第05期);全文 * |
山前褶皱带模型正演模拟及波场分析;王涛;尹川;陈景阳;蔡文杰;李萍;曹自强;段晓梦;;地球物理学进展(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113835120A (zh) | 2021-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Douglas et al. | Recent and future developments in earthquake ground motion estimation | |
CA2795340C (en) | Artifact reduction in iterative inversion of geophysical data | |
Pousse et al. | Nonstationary stochastic simulation of strong ground motion time histories including natural variability: Application to the K-net Japanese database | |
CN105954804B (zh) | 页岩气储层脆性地震预测方法及装置 | |
CN105277978B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN108369289A (zh) | 使用全波场反演点扩展函数分析设计地球物理勘测的方法 | |
CN105319589B (zh) | 一种利用局部同相轴斜率的全自动立体层析反演方法 | |
CN103792246A (zh) | 一种基于lm-bp神经网络的x射线荧光光谱定量分析方法 | |
CN104749624A (zh) | 一种地震岩相识别及其不确定性定量评价同步实现方法 | |
RU2737846C2 (ru) | Система установки наземных сейсмических датчиков с парами смежных многокомпонентных сейсмических датчиков на расстоянии в среднем по меньшей мере двадцать метров | |
RU2631407C1 (ru) | Способ и устройство для обработки сейсмических сигналов | |
CN105353405A (zh) | 一种全波形反演方法和系统 | |
CN111123374A (zh) | 一种基于匹配滤波的探地雷达全波形反演方法 | |
CN103454681A (zh) | 评价三维地震观测系统成像效果的方法和设备 | |
CN102854532A (zh) | 三维叠前炮检域随机噪声压制方法 | |
CA2412995C (en) | Seismic survey system | |
CN113835120B (zh) | 不同激发能量大小的正演模拟分析方法 | |
CN103513279B (zh) | 一种基于地震波波动方程的照明分析计算方法及计算装置 | |
Schnurr et al. | Improved parametric models for explosion pressure signals derived from large datasets | |
CN104166162B (zh) | 基于迭代三参数小波变换的缝洞发育带检测方法 | |
Parker et al. | Ground‐Motion Variability from Kinematic Rupture Models and the Implications for Nonergodic Probabilistic Seismic Hazard Analysis | |
CN112611275B (zh) | 一种爆破盲炮检测方法 | |
Stroujkova et al. | Using machine learning for explosion yield estimation | |
CN115421190A (zh) | 一种地下结构弹性波正演声波反演成像方法及设备 | |
Parasyris et al. | Near surface full waveform inversion via deep learning for subsurface imaging |
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 |