CN112379434B - 适用于沙漠区的层析反演静校正的方法及装置 - Google Patents
适用于沙漠区的层析反演静校正的方法及装置 Download PDFInfo
- Publication number
- CN112379434B CN112379434B CN202011193686.2A CN202011193686A CN112379434B CN 112379434 B CN112379434 B CN 112379434B CN 202011193686 A CN202011193686 A CN 202011193686A CN 112379434 B CN112379434 B CN 112379434B
- Authority
- CN
- China
- Prior art keywords
- speed
- desert area
- thickness
- tomographic inversion
- 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
- 230000003068 static effect Effects 0.000 title claims abstract description 86
- 238000012937 correction Methods 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 64
- 239000004576 sand Substances 0.000 claims abstract description 107
- 238000001514 detection method Methods 0.000 claims abstract description 46
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 44
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 21
- 239000010410 layer Substances 0.000 description 33
- 238000010586 diagram Methods 0.000 description 18
- 238000003384 imaging method Methods 0.000 description 10
- 238000011835 investigation Methods 0.000 description 7
- 239000002344 surface layer Substances 0.000 description 7
- 230000009189 diving Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 3
- 238000005457 optimization Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000009933 burial Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process 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
- 230000000704 physical effect 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/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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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)
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种适用于沙漠区的层析反演静校正的方法及装置,其中该方法包括:根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的相应数据,得到沙漠区替换速度模型;以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量。本发明削弱了层析反演的边界效应,同时提高了静校正的精度和效率。
Description
技术领域
本发明涉及地球物理勘探技术领域,尤其涉及适用于沙漠区的层析反演静校正的方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
静校正工作是地震资料处理过程中的一个关键步骤,尤其是在近地表速度结构或者高程剧烈变化的复杂近地表地区,如果不能保证静校正精度,会导致原始单炮记录中有效地震反射同相轴扭曲、错断特征明显,直接影响最终的剖面成像效果。
静校正最常用的方法主要有表层模型(包括时深曲线)静校正方法、折射波静校正方法和层析反演静校正方法,相对于依赖微测井等资料的表层模型静校正方法而言,地震初至波具有丰富的高频成分,理论上层析静校正方法可同时求得长、短波长静校正量,因此是目前主流的是层析反演静校正技术。
但是传统的层析反演静校正方法也存在一定的问题:一是存在边界效应;二是其反演的近地表速度场与实际结果存在着一定的误差,需要准确的表层信息加以约束和标定。因此,为了提高层析反演精度,需要额外开展一定数量的微测井表层调查工作,这将增加野外实施的工作量和勘探投资。
在地震勘探过程中,尤其是在富含饱和水的沙漠探区,如何根据施工区的特点,在减少不必要的工作量、降低勘探成本及投入的基础上,提高层析反演精度和效率,有效改善剖面的成像效果,是一个现有技术无法解决的问题。
发明内容
本发明实施例提供一种适用于沙漠区的层析反演静校正的方法,用以提高层析反演的精度和效率,该方法包括:
根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;
以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的相应数据,得到沙漠区替换速度模型;
以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;
从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;
根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量。
本发明实施例还提供一种适用于沙漠区的层析反演静校正的装置,用以提高层析反演的精度和效率,该装置包括:
沙丘厚度和速度获取模块,用于根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;
模型替换模块,用于以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的相应数据,得到沙漠区替换速度模型;
速度场反演模块,用于以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;
风化层厚度和速度提取模块,用于从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;
层析反演静校正量计算模块,用于根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述适用于沙漠区的层析反演静校正的方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述适用于沙漠区的层析反演静校正的方法的计算机程序。
本发明实施例中,根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的数据,得到沙漠区替换速度模型;以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量,从而可根据少量的静水面调查结果及沙丘曲线就能得到沙丘厚度和速度,可准确描述沙丘的近地表结构,并以其替换沙漠区初始速度模型,改变了层析反演的初始约束条件,从而削弱了层析反演的边界效应,同时反演获得的近地表速度场与实际的浅层地质条件吻合度较高,大大地提高了静校正的精度和效率,有效地改善了剖面的成像效果;与现有技术对比,无需额外开展微测井表层调查工作,节省了人力物力消耗。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中一种适用于沙漠区的层析反演静校正的方法的流程示意图;
图2为本发明实施例中一种适用于沙漠区的层析反演静校正的方法的具体示例图;
图3为本发明实施例中一种适用于沙漠区的层析反演静校正的方法的具体示例图;
图4为本发明实施例中一种适用于沙漠区的层析反演静校正的方法的具体示例图;
图5为本发明实施例中一种适用于沙漠区的层析反演静校正的方法的具体示例图;
图6为本发明实施例中应用沙丘曲线静校正法对沙漠区层析反演静校正的方法的具体示例图;
图7为本发明实施例中应用无约束层析反演静校正法对沙漠区层析反演静校正的方法的具体示例图;
图8为本发明实施例中一种适用于沙漠区的层析反演静校正的方法的具体示例图;
图9为本发明实施例中一种适用于沙漠区的层析反演静校正的装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
图1为本发明实施例中一种适用于沙漠区的层析反演静校正的方法的流程示意图,如图1所示,本发明实施例提供的适用于沙漠区的层析反演静校正的方法,可以包括:
步骤101:根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;
步骤102:以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的相应数据,得到沙漠区替换速度模型;
步骤103:以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;
步骤104:从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;
步骤105:根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量。
本发明实施例中,根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的数据,得到沙漠区替换速度模型;以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量,从而可根据少量的静水面调查结果及沙丘曲线就能得到沙丘厚度和速度,可准确描述沙丘的近地表结构,并以其替换沙漠区初始速度模型,改变了层析反演的初始约束条件,从而削弱了层析反演的边界效应,同时反演获得的近地表速度场与实际的浅层地质条件吻合度较高,大大地提高了静校正的精度和效率,有效地改善了剖面的成像效果;与现有技术对比,无需额外开展微测井表层调查工作,节省了人力物力消耗。
具体实施时,首先根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度。
实施例中,根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度,可以包括:根据沙漠区的静水面调查数据,获得沙漠区的静水面高程;根据沙漠区的静水面高程,获取每个炮点和检波点的沙丘厚度;根据沙漠区的沙丘曲线量板,获取每个炮点和检波点的沙丘速度。
在上述实施例中,利用工区范围内以往的静水面调查数据获得每个炮点和检波点的沙丘厚度,根据沙丘厚度并结合工区的沙丘曲线量板,提取其对应的每个炮点和检波点的沙丘速度。沙丘曲线量板是指采用较短排列横跨一个或者两个典型的沙丘,两端激发、固定排列接收,检波器间距可适当小一些,可实测每个物理点高程,通过微测井调查等方法获得潜水面的高程,地表高程与潜水面高程的差值为沙丘的厚度。根据延迟时的公式通过推导后可以得到沙丘厚度与垂直时(平均速度)关系式,由该关系式可得到不同沙丘厚度对应的垂直时(平均速度),将它们放在同一直角坐标系下,从而生成沙丘曲线量板。
多年勘探经验表明,在富含饱和水的沙漠区,潜水面即为高速顶界面,其形态为稳定的平面或者单斜面。潜水面上覆沙丘为连续介质,由于压实作用和含水度的不同,表层介质的物性表现为在垂向上表层速度随着深度增加而逐渐增大的特点。根据这种特性,可以将沙丘的厚度与平均速度或垂直传播时间的关系统计为沙丘曲线量板或拟合为经验公式,如果通过某种方法能够获得沙丘厚度,就可以根据沙丘曲线量板或拟合的经验公式得到其对应的平均速度。
静水面调查是指在地势低洼处,通过推水坑或者打裸眼(或套管)井后量取静水面高程的方法来求取沙丘厚度。再根据沙漠高速顶界面(潜水面)形态是平面或者单斜面的特点,在二维测线上或者三维工区平面上内插多个静水面调查的结果即可获得所有物理点的高速顶界面高程,每个炮点和检波点的地表高程与高速顶界面高程之差就是对应的沙丘厚度。
具体实施时,本发明实施例提供的适用于沙漠区的层析反演静校正的方法,还可以包括:根据预设置的地表初始速度和梯度因子,预建立沙漠区初始速度模型。
实施例中,预设置的地表初始速度可为500~1500m/s,梯度因子可为3~5s-1。预建立沙漠区初始速度模型,可以包括:根据给定的地表初始速度500~1500m/s,梯度因子3~5s-1,预建立一个沙漠区初始速度模型。
在上述实施例中,通过预建立沙漠区初始速度模型,有助于在后续步骤中对沙漠区初始速度模型中相关数据进行替换。
具体实施时,在根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度后,以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的相应数据,得到沙漠区替换速度模型。
在上述实施例中,将地表至高速顶界面(潜水面)之间的相应数据,以每个炮点和检波点的沙丘厚度和速度进行替代,并通过迭代优化,可使建立的层析反演的沙漠区替换速度模型更符合工区地表实际情况,提高层析反演的精度和效率。
具体实施时,在以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的相应数据,得到沙漠区替换速度模型后,以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场。
实施例中,所述层析反演参数包括沿接收线方向的网格大小、沿垂直接收线方向的网格大小、垂向方向的网格大小、层析反演迭代次数和初至的偏移距范围;以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场,包括:以预设置的沿接收线方向的网格大小、沿垂直接收线方向的网格大小、垂向方向的网格大小、层析反演迭代次数和初至的偏移距范围,对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场。
在上述实施例中,以预设置的层析反演参数开展模型约束层析反演,获得沙漠区近地表速度场,可以包括:对二维测线和/或3D线束,进行模型约束层析反演。层析反演的网格参数可按如下方式预设置:在沿接收线方向(InLine)和沿垂直接收线方向(CrossLine)上,网格大小设置为0.5~8倍道距;纵向上网格大小设置为2~20m;初至的偏移距范围设置为1000~6000m;层析反演迭代次数设置在8~10次之间。然后对沙漠区替换速度模型进行层析反演,反演后获得包含二维测线和/或3D线束的近地表速度场。
具体实施时,在以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场后,从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度。
实施例中,从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度,可以包括:以预设置的沙漠区顶界面速度范围,从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度。
在上述实施例中,预设置的沙漠区顶界面速度范围可为预设置的沙漠区顶界面速度范围;提取每个炮点和检波点的沙丘厚度和速度可以包括:从沙漠区近地表速度场,以给定的1400~2000m/s的速度范围作为高速顶界面的速度,来提取高速顶界面;之后在地表至高速顶界面之间,可根据沙漠区近地表速度场获得每个炮点和检波点的风化层厚度和速度。
具体实施时,在沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度后,根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量。
实施例中,计算沙漠区的层析反演静校正量,可以包括:根据按照要求给定的工区的统一基准面高程和替换速度,结合根据每个炮点和检波点的风化层厚度和速度,计算层析反演静校正量。
在上述实施例中,实现了对层析反演速度场的解释以及最终静校正的计算。
下面给出一个具体实施例,来说明本发明的方法的具体应用。该实施例依托塔里木盆地某富含饱和水的沙漠区二维地震勘探采集项目,针对了其中某一条二维测线开展了如下工作,包括如下步骤:
1、利用沙丘曲线和静水面调查数据获取每个炮点和检波点的沙丘厚度和速度
在二维测线上内插多个静水面调查的数据即可获得所有物理点的高速顶界面高程,每个炮点和检波点的地表高程与高速顶界面高程之差就是对应的沙丘厚度,如图2所示,图2展示了沙漠区的静水面调查位置、地表高程、高速顶界面高程及沙丘厚度的分布;用厚度根据工区的沙丘曲线量板提取其对应的每个炮点和检波点的沙丘速度,如图3所示,图3展示了沙漠区的沙丘速度分布;
2、优化层析反演的沙漠区初始速度模型
首先给定500~1500m/s为地表初始速度、以3~5s-1的梯度因子建立一个沙漠区初始速度模型;然后将地表至高速顶界面(潜水面)之间的模型数据用步骤1获得的沙丘厚度、速度替代,通过优化迭代,使建立的层析反演的初始速度模型更符合工区地表实际情况,如图4所示;
3、选定合适的层析反演参数开展模型约束层析反演获得沙漠区近地表速度场
对二维测线进行模型约束层析反演;层析反演的网格参数可在InLine方向上设置为1倍道距(30m),垂向上网格大小设置为10m;初至的偏移距范围设置为2000m;层析反演迭代次数设置为10次。然后运用步骤2获得的沙漠区初始速度模型进行层析反演,反演后获得二维测线的沙漠区近地表速度场,如图5所示;
4、提取每个炮点和检波点的风化层厚度和速度
从步骤5获得的近地表速度场中,以给定1500m/s的速度范围作为高速顶界面的速度来提取高速顶界面,在地表至高速顶界面之间根据速度场就能够获得每个炮点和检波点的风化层厚度和速度;
5、计算沙漠区的层析反演静校正量
根据按照要求所给定的工区的统一基准面高程和替换速度,再通过步骤4获得的每一个炮点和检波点的沙丘厚度和速度,按如下公式计算层析反演静校正量:
其中,T为沙漠区的每一炮点或检波点的层析反演静校正量,单位为ms;
hi为沙漠区第i层介质的风化层厚度,单位为m;
Vi为沙漠区第i层介质的风化层速度,单位为m/s;
τ为沙漠区的井深或检波器埋深的校正时间,单位为s;
Hd为沙漠区的基准面高程,单位为m;
Hg为沙漠区的高速层的顶界面高程,单位为m;
Vs为沙漠区的基准面的校正速度,单位为m/s。
图6为本发明实施例中应用沙丘曲线静校正法对沙漠区层析反演静校正的方法的具体示例图,如图6所示,在应用沙丘曲线静校正法对沙漠区层析反演静校正时,应用沙丘曲线静校正方法得到的时间水平叠加剖面,整体成像效果相对较差,但无边界效应。图7为本发明实施例中应用无约束层析反演静校正法对沙漠区层析反演静校正的方法的具体示例图,如图7所示,应用无约束层析反演静校正方法得到地时间水平叠加剖面除边界效应外,其成像效果好于沙丘曲线静校正方法。图8为本发明实施例中一种适用于沙漠区的层析反演静校正的方法的具体示例图,如图8所示,本发明实施例中适用于沙漠区的层析反演静校正的方法的整体成像效果最好,且削弱了边界效应。
本发明实施例中,根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的数据,得到沙漠区替换速度模型;以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量,从而可根据少量的静水面调查结果及沙丘曲线就能得到沙丘厚度和速度,可准确描述沙丘的近地表结构,并以其替换沙漠区初始速度模型,改变了层析反演的初始约束条件,从而削弱了层析反演的边界效应,同时反演获得的近地表速度场与实际的浅层地质条件吻合度较高,大大地提高了静校正的精度和效率,有效地改善了剖面的成像效果;与现有技术对比,无需额外开展微测井表层调查工作,节省了人力物力消耗。
如上所述,本发明实施例充分地利用了沙漠区近地表结构的特点,可只根据少量静水面调查数据及沙丘曲线量板,就能得到沙丘的近地表结构;并将其作为层析反演的初始约束条件,从而削弱了层析反演的边界效应,同时提高了层析反演的精度和效率,反演获得的沙漠区近地表速度场与实际的浅层地质条件吻合度较高,大大地提高了静校正的精度,有效地改善了剖面的成像效果。
本发明实施例中还提供了一种适用于沙漠区的层析反演静校正的装置,如下面的实施例所述。由于该装置解决问题的原理与适用于沙漠区的层析反演静校正的方法相似,因此该装置的实施可以参见适用于沙漠区的层析反演静校正的方法的实施,重复之处不再赘述。
图9为本发明实施例中一种适用于沙漠区的层析反演静校正的装置的结构示意图,如图9所示,本发明实施例提供的适用于沙漠区的层析反演静校正的装置,可以包括:
沙丘厚度和速度获取模块01,用于根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;
模型替换模块02,用于以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的相应数据,得到沙漠区替换速度模型;
速度场反演模块03,用于以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;
风化层厚度和速度提取模块04,用于从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;
层析反演静校正量计算模块05,用于根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量。
在一个实施例中,沙丘厚度和速度获取模块,具体用于:
根据沙漠区的静水面调查数据,获得沙漠区的静水面高程;
根据沙漠区的静水面高程,获取每个炮点和检波点的沙丘厚度;
根据沙漠区的沙丘曲线量板,获取每个炮点和检波点的沙丘速度。
在一个实施例中,本发明实施例提供的适用于沙漠区的层析反演静校正的装置,还可以包括:建模模块,用于:根据预设置的地表初始速度和梯度因子,预建立沙漠区初始速度模型。
在一个实施例中,所述层析反演参数包括网格参数及方向、层析反演迭代次数和初至的偏移距范围;速度场反演模块,具体用于:以预设置的网格参数及方向、层析反演迭代次数和初至的偏移距范围,对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场。
在一个实施例中,风化层厚度和速度提取模块,具体用于:以预设置的沙漠区顶界面速度范围,从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述适用于沙漠区的层析反演静校正的方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述适用于沙漠区的层析反演静校正的方法的计算机程序。
本发明实施例中,根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的数据,得到沙漠区替换速度模型;以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量,从而可根据少量的静水面调查结果及沙丘曲线就能得到沙丘厚度和速度,可准确描述沙丘的近地表结构,并以其替换沙漠区初始速度模型,改变了层析反演的初始约束条件,从而削弱了层析反演的边界效应,同时反演获得的近地表速度场与实际的浅层地质条件吻合度较高,大大地提高了静校正的精度和效率,有效地改善了剖面的成像效果;与现有技术对比,无需额外开展微测井表层调查工作,节省了人力物力消耗。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (12)
1.一种适用于沙漠区的层析反演静校正的方法,其特征在于,包括:
根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;
以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的相应数据,得到沙漠区替换速度模型;
以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;
从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;
根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量。
2.如权利要求1所述的方法,其特征在于,根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度,包括:
根据沙漠区的静水面调查数据,获得沙漠区的静水面高程;
根据沙漠区的静水面高程,获取每个炮点和检波点的沙丘厚度;
根据沙漠区的沙丘曲线量板,获取每个炮点和检波点的沙丘速度。
3.如权利要求1所述的方法,其特征在于,还包括:根据预设置的地表初始速度和梯度因子,预建立沙漠区初始速度模型。
4.如权利要求1所述的方法,其特征在于,所述层析反演参数包括沿接收线方向的网格大小、沿垂直接收线方向的网格大小、垂向方向的网格大小、层析反演迭代次数和初至的偏移距范围;
以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场,包括:
以预设置的沿接收线方向的网格大小、沿垂直接收线方向的网格大小、垂向方向的网格大小、层析反演迭代次数和初至的偏移距范围,对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场。
5.如权利要求1所述的方法,其特征在于,从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度,包括:
以预设置的沙漠区顶界面速度范围,从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度。
6.一种适用于沙漠区的层析反演静校正的装置,其特征在于,包括:
沙丘厚度和速度获取模块,用于根据沙漠区的静水面调查数据和沙丘曲线,获取每个炮点和检波点的沙丘厚度和速度;
模型替换模块,用于以每个炮点和检波点的沙丘厚度和速度,替换预建立的沙漠区初始速度模型中的地表至顶界面的相应数据,得到沙漠区替换速度模型;
速度场反演模块,用于以预设置的层析反演参数对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场;
风化层厚度和速度提取模块,用于从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度;
层析反演静校正量计算模块,用于根据每个炮点和检波点的风化层厚度和速度,结合预设的沙漠区基准面高程和替换速度,计算沙漠区的层析反演静校正量。
7.如权利要求6所述的装置,其特征在于,沙丘厚度和速度获取模块,具体用于:
根据沙漠区的静水面调查数据,获得沙漠区的静水面高程;
根据沙漠区的静水面高程,获取每个炮点和检波点的沙丘厚度;
根据沙漠区的沙丘曲线量板,获取每个炮点和检波点的沙丘速度。
8.如权利要求6所述的装置,其特征在于,还包括建模模块,用于:根据预设置的地表初始速度和梯度因子,预建立沙漠区初始速度模型。
9.如权利要求6所述的装置,其特征在于,所述层析反演参数包括沿接收线方向的网格大小、沿垂直接收线方向的网格大小、垂向方向的网格大小、层析反演迭代次数和初至的偏移距范围;
速度场反演模块,具体用于:以预设置的沿接收线方向的网格大小、沿垂直接收线方向的网格大小、垂向方向的网格大小、层析反演迭代次数和初至的偏移距范围,对沙漠区替换速度模型进行模型约束层析反演,得到沙漠区近地表速度场。
10.如权利要求6所述的装置,其特征在于,风化层厚度和速度提取模块,具体用于:以预设置的沙漠区顶界面速度范围,从沙漠区近地表速度场,提取每个炮点和检波点的风化层厚度和速度。
11.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至5任一所述方法。
12.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至5任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011193686.2A CN112379434B (zh) | 2020-10-30 | 2020-10-30 | 适用于沙漠区的层析反演静校正的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011193686.2A CN112379434B (zh) | 2020-10-30 | 2020-10-30 | 适用于沙漠区的层析反演静校正的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112379434A CN112379434A (zh) | 2021-02-19 |
CN112379434B true CN112379434B (zh) | 2024-04-30 |
Family
ID=74576053
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011193686.2A Active CN112379434B (zh) | 2020-10-30 | 2020-10-30 | 适用于沙漠区的层析反演静校正的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112379434B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104133245A (zh) * | 2014-07-22 | 2014-11-05 | 中国石油天然气集团公司 | 一种地震资料的静校正方法及系统 |
CN105467452A (zh) * | 2015-12-18 | 2016-04-06 | 中国石油天然气集团公司 | 一种对过渡区域进行静校正的方法及装置 |
WO2017035104A1 (en) * | 2015-08-24 | 2017-03-02 | Schlumberger Technology Corporation | Velocity model seismic static correction |
CN106990438A (zh) * | 2017-03-28 | 2017-07-28 | 中国石油天然气集团公司 | 复杂地表区域微测井点位及井深确定方法及装置 |
CN109581503A (zh) * | 2017-09-28 | 2019-04-05 | 中国石油化工股份有限公司 | 一种解决沙漠探区层析静校正边界问题的方法及系统 |
CN111399037A (zh) * | 2019-01-02 | 2020-07-10 | 中国石油天然气集团有限公司 | 高速顶界面提取的方法和装置 |
CN111435172A (zh) * | 2019-01-15 | 2020-07-21 | 中国石油天然气集团有限公司 | 层析静校正的方法和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10481287B2 (en) * | 2016-08-05 | 2019-11-19 | Saudi Arabian Oil Company | Surface consistent statics solution and amplification correction |
US10920585B2 (en) * | 2017-12-26 | 2021-02-16 | Saudi Arabian Oil Company | Determining sand-dune velocity variations |
-
2020
- 2020-10-30 CN CN202011193686.2A patent/CN112379434B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104133245A (zh) * | 2014-07-22 | 2014-11-05 | 中国石油天然气集团公司 | 一种地震资料的静校正方法及系统 |
WO2017035104A1 (en) * | 2015-08-24 | 2017-03-02 | Schlumberger Technology Corporation | Velocity model seismic static correction |
CN105467452A (zh) * | 2015-12-18 | 2016-04-06 | 中国石油天然气集团公司 | 一种对过渡区域进行静校正的方法及装置 |
CN106990438A (zh) * | 2017-03-28 | 2017-07-28 | 中国石油天然气集团公司 | 复杂地表区域微测井点位及井深确定方法及装置 |
CN109581503A (zh) * | 2017-09-28 | 2019-04-05 | 中国石油化工股份有限公司 | 一种解决沙漠探区层析静校正边界问题的方法及系统 |
CN111399037A (zh) * | 2019-01-02 | 2020-07-10 | 中国石油天然气集团有限公司 | 高速顶界面提取的方法和装置 |
CN111435172A (zh) * | 2019-01-15 | 2020-07-21 | 中国石油天然气集团有限公司 | 层析静校正的方法和装置 |
Non-Patent Citations (8)
Title |
---|
一种新的沙漠区资料静校正方法;孙天真;尚新民;赵爱国;滕厚华;刘立平;《油气地球物理》;20100131;第8卷(第01期);第17-22页 * |
含压实效应的沙丘曲线静校正方法及其应用;陈可洋;吴清岭;林春华;李来林;裴江云;范兴才;关昕;《矿业工程研究》;20111231;第26卷(第04期);第53-57页 * |
地下饱和含水沙漠区静校正技术应用及效果;周翼;杨洪波;张新东;吕景峰;张立超;邸江伟;《非常规油气》;20170630;第4卷(第03期);第1-7页 * |
塔里木盆地二维地震老资料拼接处理技术及应用;杜娟;;《内蒙古石油化工》;20180228;第44卷(第02期);第92-98页 * |
复杂近地表区走时射线层析的应用及效果;周翼;杨洪波;黄有晖;吕景峰;张立超;陈建国;《天然气勘探与开发》;20171231;第40卷(第04期);第52-58页 * |
折射静校正在苏里格气田三维处理中的应用;王克斌, 赵灵芝, 张旭民;《石油物探》;20030731;第42卷(第02期);第248-251页 * |
约束层析反演及其在地震速度计算中的应用;袁刚;冯心远;蒋波;张涛;《石油物探》;20130131;第52卷(第01期);第55-59页 * |
董1井区静校正方法研究与应用;王增明, 刘超敏, 刘斌, 潘家智, 夏先贵;《石油物探》;20050831;第44卷(第04期);第379-382页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112379434A (zh) | 2021-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105277978B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN102027390B (zh) | 对地震数据的联合内插及消除虚反射 | |
CN102967882B (zh) | 地层的层速度模型的建模方法 | |
CN106526674A (zh) | 一种三维全波形反演能量加权梯度预处理方法 | |
CN104570102B (zh) | 一种近地表速度模型与中深层速度模型的融合方法 | |
US10739481B2 (en) | 2D multiline seismic reflection tomography with seismic-tie constraint | |
CN113552625B (zh) | 一种用于常规陆域地震数据的多尺度全波形反演方法 | |
CN109917454B (zh) | 基于双基准面的真地表叠前深度偏移成像方法及装置 | |
CN107817516B (zh) | 基于初至波信息的近地表建模方法及系统 | |
CN108845350B (zh) | 反演二维速度模型的方法及装置 | |
CN103149588A (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN109946742A (zh) | 一种TTI介质中纯qP波地震数据模拟方法 | |
CN109696705A (zh) | 深海obs上下行波场成像数据校正方法及装置 | |
CN105301638B (zh) | 一种提取风化层底界面的方法和装置 | |
CN106291698B (zh) | 地震相沉积相确定方法和装置 | |
Hole et al. | Interface inversion using broadside seismic refraction data and three‐dimensional travel time calculations | |
CN111399037B (zh) | 高速顶界面提取的方法和装置 | |
CN112379434B (zh) | 适用于沙漠区的层析反演静校正的方法及装置 | |
CN109725354B (zh) | 各向异性速度建模方法及系统 | |
CN113777654B (zh) | 一种基于伴随状态法初至波走时层析的海水速度建模方法 | |
CN115826038A (zh) | 一种vti介质伴随状态法走时多参数层析成像方法及系统 | |
CN115586575A (zh) | 一种3d被动源拟炮集计算方法、系统、设备及存储介质 | |
CN111435172A (zh) | 层析静校正的方法和装置 | |
CN112285782B (zh) | 近地表地震波吸收衰减调查方法及装置 | |
US10401515B2 (en) | Estimation of water properties from seismic data |
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 |