CN101625419A - 利用井口时间提高静校正量精度的方法 - Google Patents
利用井口时间提高静校正量精度的方法 Download PDFInfo
- Publication number
- CN101625419A CN101625419A CN200810116290A CN200810116290A CN101625419A CN 101625419 A CN101625419 A CN 101625419A CN 200810116290 A CN200810116290 A CN 200810116290A CN 200810116290 A CN200810116290 A CN 200810116290A CN 101625419 A CN101625419 A CN 101625419A
- Authority
- CN
- China
- Prior art keywords
- static correction
- speed
- uphole time
- uphole
- time
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是地震勘探中利用井口时间提高静校正量精度的方法,收集记录井口时间,剔除其中的异常值,用已知的表层厚度和速度模型,直接用井口时间计算炮点静校正量,当野外表层调查数据很少或没有表层调查数据,没有已知的表层厚度和速度模型时,提取每一激发井的深度和井口时间,计算等效平均速度,建立高密度的表层厚度模型和速度模型,计算静校正量。本发明直接利用井口时间提高了静校正量精度,优化了炮点静校正量计算结果,也可以通过井口时间在表层调查控制点稀少或没有的工区建立高精度的静校正量计算模型,得到高精度静校正量。同时得到的表层速度数据模型能够在初至折射波、层析成像、波动方程等静校正方法中应用,经济效益可观。
Description
技术领域
本发明涉及地震勘探技术,具体是利用井口时间提高静校正量精度的方法。
背景技术
石油地震勘探是通过人工方法产生地震波,通过研究地震波在地下岩层中的传播情况,勘察地下地质构造和岩层变化,寻找油气田的一种地球物理勘探方法。地震勘探可以分为地震资料采集、地震资料处理、地震资料解释三个大的环节。
地震资料采集是通过人工的方法产生地震波,在陆地上通常使用炸药震源,在选定的炮点打一口浅井,井深一般为6-30米。炸药的药量一般为1-25公斤,装入电雷管放到井底,引爆炸药产生地震波。然后用地面地震接收仪器记录由地下反射回来的地震波所引起的地面振动信号,并记录在磁带上,提供给室内计算中心进行后续处理。
但在地震资料采集过程中,在近地表存在一种速度很低的地层结构,造成地震波在通过它时花费更多的时间。这种结构在纵向和横向上变化很大,因此不同位置地震波消耗时间也不一样,这种时间上的差异会给后续地震资料处理带来很大的麻烦,不能同相叠加,地下成像困难。因此,要得到能反映地下真实情况的地震剖面,就必须消除掉表层结构造成的横向时差影响。地震勘探专业中将这个消除时差过程称之为静校正。
为了实现静校正,即得到时差并加以消除,就必须要采用一些技术和方法。多年的发展推出了许多方法,有大炮初至折射,层析成像,小折射,微测井,地质露头,陆地声纳,面波勘探,电磁法等十多种地震和非地震方法,这些方法主要是通过建立表层结构模型(包括厚度模型和速度模型)来计算得到地震波穿过表层的时间,属于间接方法。这些方法在应用时都需要计算一些相关参数,有些参数能够运用数学中的计算方法准确求得,有些参数则不能。这其中有一个很重要的参数就是表层速度,是静校正量计算不可或缺的。目前应用最广泛的几种静校正方法中,野外表层调查可以提取若干表层速度信息,但通常是点控制,密度不会太大,精度受到限制。大炮初至折射不能直接得到表层速度,层析成像虽然可以反演速度模型,但往往出现多解性,结果很不稳定。因此,表层速度的精度对静校正效果产生了很大影响。
井口时间是野外井炮生产时产生的一种信号,是地震勘探中的重要参数,它所包含的信息相当丰富,其具有炮点数量多、密度高、信息量大、覆盖完整的特点,每一井位的井口时间、井深、井内岩性均直接地反映了井深段这一距离内的信息,所表现的表层信息具有较好的客观对应性和直观性。
目前,野外生产中井口时间一直都在记录并没有废弃,但也没有利用起来。近年来,随着激发井深的不断加大,井口时间的应用价值也有了极大提高。
发明内容
本发明目的在于提供一种优化炮点静校正量计算结果的利用井口时间提高静校正量精度的方法。
本发明具体步骤如下:
1)收集野外地震采集施工时记录井口时间,剔除其中的异常值;
步骤1)所述的异常值是为零值或异常大于正常井口时间值的数值。
2)利用已知的表层厚度和速度模型,直接用井口时间按以下公式计算炮点静效正量:
式中:
T-炮点静校正量(ms),
hi-第i层介质的厚度(m),
Vi-第i层介质的速度(m/s),
Hd-基准面高程(m),
Hg-高速顶界面高程(m),
Vs-基准面校正速度(m/s),
τ-井口时间(ms),
dj-井深在第j层介质的厚度(m)。
3)当野外表层调查数据很少或没有表层调查数据,没有已知的表层厚度和速度模型时,提取每一激发井的深度(d)和井口时间(τ),计算等效平均速度v=d/τ,每一个激发点是一个空间样点;
4)根据深度d和速度v两个属性参数,建立高密度的表层厚度模型和速度模型;
步骤4)所述的模型建立方法采用常规空间建模方法。
5)在上面建立的模型基础上,用以下公式计算静校正量:
检波点静校正量:
炮点静校正量:
上式中:
T-炮点、检波点静校正量(ms),
h-表层介质厚度(m),
V-表层速度(m/s),
Hd-基准面高程(m),
Hg-高速顶界面高程(m),
Vs-基准面校正速度(m/s)。
本发明直接利用井口时间提高了静校正量精度,优化了炮点静校正量计算结果。还可以通过井口时间在表层调查控制点稀少或没有的工区建立高精度的静校正量计算模型,得到高精度静校正量。同时得到的表层速度数据模型能够在初至折射波、层析成像、波动方程等静校正方法中应用,经济效益可观。
附图说明
图1是炮点、检波点分布图;
图2是本发明检波点静校正量;
图3是本发明炮点静校正量。
具体实施方式
本发明主要通过两个方面实施,一是已有相对准确的表层模型,用井口时间可以优化炮点静校正量计算中的井深校正,提高炮点静校正量精度。炮点静校正量由低降速带校正、井深校正和基准面校正组成。在以往井深校正时,不是直接应用井口时间,而是用表层速度和井深进行时间转换后再应用,由于速度是通过点内插得来或采用常速,精度受到了影响。如果直接采用井口时间这一时间量,炮点静校正量的计算精度就会得到明显提高。
通常静校正量计算公式:
采用井口时间τ后可变为:
式中:
T-炮点静校正量(ms)
hi-第i层介质的厚度(m)
Vi-第i层介质的速度(m/s)
Hd-基准面高程(m)
Hg-高速顶界面高程(m)
Vs-基准面校正速度(m/s)
τ-井口时间(ms)
dj-井深在第j层介质的厚度(m)
公式(1)的物理意义是模型法炮点静校正量计算,公式右边第二项做井深校正时采用层厚度除以层速度累加的方法,为多年常规方法。公式(2)右边第二项直接引入了井口时间对井深加以校正。
二是野外表层调查数据很少或没有表层调查数据时,用井口时间提取深度(厚度)、时间、速度参数,将生产中的炮点作为表层调查控制点,进而建立表层模型。由于生产炮密度大,数量多,因而其建立的表层模型精度也高,静校正量结果也比较精确。
目前,生产中激发井都比较深,通常打穿低降速层到高速层,用高密度的井口时间建立静校正计算模型,提取每一激发井的深度d和井口时间τ,用d和τ计算等效速度,计算公式是v=d/τ,这样,每一个激发点就都是一个空间样点(图1),有了深度(厚度)d和速度v两个属性,用它们就可以建立高密度的表层厚度模型和速度模型。模型建立方法采用常规空间建模方法。
模型建好后,静校正量计算公式就变成下面的两个式子:
检波点静校正量:
炮点静校正量:
式中:
T-炮点、检波点静校正量(ms)
h-表层介质厚度(m)
V-表层速度(m/s)
Hd-基准面高程(m)
Hg-高速顶界面高程(m)
Vs-基准面校正速度(m/s)
如图2、图3为本发明计算的检波点和炮点静校正量。
另外,这种方法同时得到了表层速度数据模型,可以在初至折射波、层析成像、波动方程等静校正方法中加以应用,从而提高静校正精度。
Claims (3)
1、一种利用井口时间提高静校正量精度的方法,其特征在于具体步骤如下:
1)收集野外地震采集施工时记录井口时间,剔除其中的异常值;
2)利用已知的表层厚度和速度模型,直接用井口时间按以下公式计算炮点静效正量:
式中:
T-炮点静校正量(ms),
hi-第i层介质的厚度(m),
Vi-第i层介质的速度(m/s),
Hd-基准面高程(m),
Hg-高速顶界面高程(m),
Vs-基准面校正速度(m/s),
τ-井口时间(ms),
dj-井深在第j层介质的厚度(m)。
3)当野外表层调查数据很少或没有表层调查数据,没有已知的表层厚度和速度模型时,提取每一激发井的深度(d)和井口时间(τ),计算等效平均速度v=d/τ,每一个激发点是一个空间样点;
4)根据深度d和速度v两个属性参数,建立高密度的表层厚度模型和速度模型;
5)在上面建立的模型基础上,用以下公式计算静校正量:
检波点静校正量:
炮点静校正量:
上式中:
T-炮点、检波点静校正量(ms),
h-表层介质厚度(m),
V-表层速度(m/s),
Hd-基准面高程(m),
Hg-高速顶界面高程(m),
Vs-基准面校正速度(m/s)。
2、根据权利要求1所述的利用井口时间提高静校正量精度的方法,其特征在于步骤1)所述的异常值是为零值或异常大于正常井口时间值的数值。
3、根据权利要求1所述的利用井口时间提高静校正量精度的方法,其特征在于步骤4)所述的模型建立方法采用常规空间建模方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200810116290A CN101625419A (zh) | 2008-07-08 | 2008-07-08 | 利用井口时间提高静校正量精度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200810116290A CN101625419A (zh) | 2008-07-08 | 2008-07-08 | 利用井口时间提高静校正量精度的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101625419A true CN101625419A (zh) | 2010-01-13 |
Family
ID=41521328
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200810116290A Pending CN101625419A (zh) | 2008-07-08 | 2008-07-08 | 利用井口时间提高静校正量精度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101625419A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102236104A (zh) * | 2010-05-06 | 2011-11-09 | 中国石油天然气集团公司 | 三维密度变速地震反射空校方法 |
CN102967882A (zh) * | 2012-11-16 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地层的层速度模型的建模方法 |
CN103116184A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种表层调查中计算横波静校正量的方法 |
CN103616720A (zh) * | 2013-10-28 | 2014-03-05 | 中国石油天然气集团公司 | 一种折射速度获取方法及装置 |
CN103869368A (zh) * | 2012-12-10 | 2014-06-18 | 中国石油天然气集团公司 | 一种无表层调查资料约束的大炮初至综合建模静校正方法 |
CN105467452A (zh) * | 2015-12-18 | 2016-04-06 | 中国石油天然气集团公司 | 一种对过渡区域进行静校正的方法及装置 |
CN111323818A (zh) * | 2018-12-13 | 2020-06-23 | 中国石油天然气股份有限公司 | 陆上地震资料的静校正方式的筛选方法及装置 |
-
2008
- 2008-07-08 CN CN200810116290A patent/CN101625419A/zh active Pending
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102236104B (zh) * | 2010-05-06 | 2014-12-24 | 中国石油天然气集团公司 | 三维密度变速地震反射空校方法 |
CN102236104A (zh) * | 2010-05-06 | 2011-11-09 | 中国石油天然气集团公司 | 三维密度变速地震反射空校方法 |
CN103116184B (zh) * | 2011-11-16 | 2016-06-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种表层调查中计算横波静校正量的方法 |
CN103116184A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种表层调查中计算横波静校正量的方法 |
CN102967882B (zh) * | 2012-11-16 | 2015-08-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地层的层速度模型的建模方法 |
CN102967882A (zh) * | 2012-11-16 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地层的层速度模型的建模方法 |
CN103869368A (zh) * | 2012-12-10 | 2014-06-18 | 中国石油天然气集团公司 | 一种无表层调查资料约束的大炮初至综合建模静校正方法 |
CN103616720A (zh) * | 2013-10-28 | 2014-03-05 | 中国石油天然气集团公司 | 一种折射速度获取方法及装置 |
CN103616720B (zh) * | 2013-10-28 | 2016-04-06 | 中国石油天然气集团公司 | 一种折射速度获取方法及装置 |
CN105467452A (zh) * | 2015-12-18 | 2016-04-06 | 中国石油天然气集团公司 | 一种对过渡区域进行静校正的方法及装置 |
CN105467452B (zh) * | 2015-12-18 | 2018-01-05 | 中国石油天然气集团公司 | 一种对过渡区域进行静校正的方法及装置 |
CN111323818A (zh) * | 2018-12-13 | 2020-06-23 | 中国石油天然气股份有限公司 | 陆上地震资料的静校正方式的筛选方法及装置 |
CN111323818B (zh) * | 2018-12-13 | 2022-06-03 | 中国石油天然气股份有限公司 | 陆上地震资料的静校正方式的筛选方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102176054B (zh) | 近地表综合信息处理解释方法 | |
CN103233727B (zh) | 一种反演地层横波速度径向剖面的方法 | |
US9075158B2 (en) | Using a drill bit as a seismic source for SET velocity analysis | |
AU2011296498B2 (en) | Method for detection of subsurface seismic events in vertically transversely isotropic media | |
CN103116184B (zh) | 一种表层调查中计算横波静校正量的方法 | |
CN101625419A (zh) | 利用井口时间提高静校正量精度的方法 | |
CN102053263B (zh) | 调查表层结构的方法 | |
CN106597533A (zh) | 一种用于山前带地震资料处理的深度域速度建模方法 | |
CN106154334A (zh) | 基于网格搜索的井下微地震事件实时反演定位方法 | |
CN103675944B (zh) | 一种沉积盆地铀资源勘探的综合地球物理方法 | |
CN109164487A (zh) | 一种基于模型建立平均速度场的方法及精细成图的方法 | |
EP2812730A2 (en) | Integrated passive and active seismic surveying using multiple arrays | |
WO2008124759A1 (en) | Location of a microseismic event using the combination of a vsp array and a surface array | |
CN103345001A (zh) | 一种测定古湖泊水深的方法 | |
Gajek et al. | Results of the downhole microseismic monitoring at a pilot hydraulic fracturing site in Poland—Part 1: Event location and stimulation performance | |
CN107422384B (zh) | 一种利用测井资料仿垂直地震剖面资料的方法 | |
McCann et al. | Application of cross-hole seismic measurements in site investigation surveys | |
CN109283577A (zh) | 一种地震层位标定方法 | |
CN106908838A (zh) | 构建目标区域地层非弹性衰减品质因子立体模型的方法 | |
CN104181610B (zh) | 一种射线路径弹性反演方法以及系统 | |
CN110989034B (zh) | 一种回归-分形插值法反演测井横波时差方法 | |
CN115586569A (zh) | 基于理论模型约束下数据驱动的地层水平地应力计算方法 | |
CN108897046B (zh) | 远探测声波逆时偏移成像方法及装置 | |
CN107894611B (zh) | 一种水平井数据转直井数据的方法及用于地震反演的方法 | |
CN112031755A (zh) | 一种压裂监测装置、方法、系统、电子设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20100113 |