CN110513098B - 一种提高元素录井曲线纵向分辨率的方法 - Google Patents
一种提高元素录井曲线纵向分辨率的方法 Download PDFInfo
- Publication number
- CN110513098B CN110513098B CN201910634199.6A CN201910634199A CN110513098B CN 110513098 B CN110513098 B CN 110513098B CN 201910634199 A CN201910634199 A CN 201910634199A CN 110513098 B CN110513098 B CN 110513098B
- Authority
- CN
- China
- Prior art keywords
- data
- curve
- logging
- target
- calcium
- 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
- 238000000034 method Methods 0.000 title claims abstract description 15
- BHPQYMZQTOCNFJ-UHFFFAOYSA-N Calcium cation Chemical group [Ca+2] BHPQYMZQTOCNFJ-UHFFFAOYSA-N 0.000 claims description 42
- 238000005553 drilling Methods 0.000 claims description 21
- 238000005259 measurement Methods 0.000 claims description 7
- 239000011435 rock Substances 0.000 description 13
- 238000005070 sampling Methods 0.000 description 8
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 5
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 4
- 239000011575 calcium Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 229910052791 calcium Inorganic materials 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种提高元素录井曲线纵向分辨率的方法:收集目标井的录井工程资料、元素录井资料;选取目标元素曲线的标准数据;以标准数据为基础,建立目标元素曲线与录井工程曲线之间的多元回归方程;在相邻实测目标元素数据之间,将录井工程曲线值代入多元回归方程,完成对目标元素曲线的数据加密,从而提高目标元素曲线的纵向分辨率。本发明提高了元素录井曲线的纵向分辨率,操作简单,易于推广应用,适用于现场地质人员对地层信息进行快速评价。
Description
技术领域
本发明涉及石油天然气勘探开发技术领域,更具体的说,是涉及一种提高元素录井曲线纵向分辨率的方法。
背景技术
随着油田勘探领域的扩大和钻井新工艺的应用,钻遇岩性复杂、录井岩屑细小、钻井速度快,给传统的录井岩性识别带来了很大挑战。
元素录井技术通过测量出录井岩屑中的元素含量从而为岩性识别提供了一种新的技术手段。在《特种油气藏》,2017(5):78-84公开的利用元素录井技术进行岩性识别中,其公开了如下主内容:渤海海域介于郯庐断裂带与沧东断裂带之间,目前探井已经揭示的潜山地层和岩石类型主要为中生界火成岩类、古生界碳酸盐岩类和元古-太古界变质岩类,潜山岩性复杂且难以识别。以元素录井数据为基础,结合钻井、岩心、全岩分析等数据,建立了渤海海域潜山岩性-元素数据库,并通过分析对比不同潜山岩性的元素含量差异,提取了潜山各类岩性的敏感元素,运用敏感元素含量交会图的方法,建立了3类主要岩性的元素区分图版和火山岩、碳酸盐岩的元素成分命名图版。
但在钻井新工艺推广应用条件下,钻井速度快,岩屑加密采样难,导致元素录井曲线的采样间距通常为5米/次,曲线的纵向分辨率较低,难以有效识别出薄层的岩性变化。录井工程曲线的采样间距为1米/次、纵向分辨率相对较高,且能一定程度上反应出薄层岩性的变化规律。因此,通过建立录井工程曲线与目标元素曲线之间的数学关系式,即可在录井阶段得到一条采样间距为1米/次的模拟元素曲线,从而有效的提高了元素录井曲线的纵向分辨率。
发明内容
本发明的目的是为了克服现有技术中的不足,提出一种提高元素录井曲线纵向分辨率的方法,通过建立元素曲线与录井工程曲线之间多元回归方程,实现录井阶段对元素曲线的数据加密,从而提高了元素曲线的纵向分辨率,解决了元素录井曲线薄层识别难的问题。
本发明的目的是通过以下技术方案实现的。
本发明提高元素录井曲线纵向分辨率的方法,包括以下步骤:
1)收集目标井的录井工程资料、元素录井资料;
2)选取目标元素曲线的标准数据;
3)以标准数据为基础,建立目标元素曲线与录井工程曲线之间的多元回归方程;
4)在相邻实测目标元素数据之间,将录井工程曲线值代入多元回归方程,完成对目标元素曲线的数据加密,从而提高目标元素曲线的纵向分辨率。
所述录井工程资料具体为:钻时数据、钻压数据、扭矩数据、转盘转速数据。
所述目标元素曲线的标准数据具体为:同一层位目标元素的录井实测数据及其等深度的录井工程数据。
所述多元回归方程具体为:通过EXCEL软件建立的目标元素曲线与录井工程曲线的关系式。
通过纵向分辨率更高的测井元素曲线对加密后的目标元素曲线进行评价。
与现有技术相比,本发明的技术方案所带来的有益效果是:
(1)本发明利用录井工程数据对元素录井数据进行加密处理属于录井资料处理方法创新,是在目前钻井新工艺条件下,钻井速度快,加密采样难的条件下,应用录井工程数据使得元素录井数据的采样间距由5米/次提高到1米/次。
(2)本发明具有实时性,以元素录井、工程录井资料为基础,故能在录井阶段得到一条分辨率较高的元素曲线进行岩性解释。
(3)本发明具有低成本,无需任何实验分析的额外费用。
(4)本发明具有易操作性,利用EXCEL软件即可运行,简单实用,快速高效。
附图说明
图1为加密钙元素(Ca)曲线效果分析图;
图2为加密钙元素数据与测井实测钙元素数据相关性分析图。
具体实施方式
为使本发明的目的、技术方案和有点更加清楚,下面结合附图对本发明实施方式作进一步的详细描述。
本发明提高元素录井曲线纵向分辨率的方法,包括以下步骤:
1)收集目标井的录井工程资料、元素录井资料。其中,录井工程资料具体为:钻时数据、钻压数据、扭矩数据、转盘转速数据。
2)选取目标元素曲线的标准数据:同一层位目标元素的录井实测数据及其等深度的录井工程数据。
3)以标准数据为基础,建立目标元素曲线与录井工程曲线之间的多元回归方程。多元回归方程具体为:通过EXCEL软件建立的目标元素曲线与录井工程曲线的关系式。
4)在相邻实测目标元素数据之间,将录井工程曲线值代入多元回归方程,即可完成对目标元素曲线的数据加密,从而提高目标元素曲线的纵向分辨率。另外还可通过纵向分辨率更高的测井元素曲线对加密后的目标元素曲线进行评价。具体实施例
本发明实施例以渤海B工区B22井古生界碳酸盐岩潜山地层元素录井测量的钙元素(Ca)曲线为例,该井钙元素曲线的采样间距为5米/次,纵向分辨率相对较低,现对该井钙元素曲线进行高分辨率处理,以得到一条采样间距为1米/次、纵向分辨率更高的模拟钙元素曲线,具体实现过程如下:
101:收集目标井的录井工程资料、元素录井资料;
收集B22井的录井工程资料、元素录井资料,其中录井工程资料包括:钻时数据、钻压数据、扭矩数据、转盘转速数据;元素录井资料为:钙元素数据。
102:选取目标元素曲线的标准数据;
B22井钻遇古生界地层深度范围为4346米至4611米,自4350米开始,每隔5米录井测量的钙元素数据以及相同深度点的钻时(ROP)、钻压(WOB)、扭矩(TRQ)、转盘转速(RPM)数据作为标准数据。
103:建立钙元素曲线与录井工程曲线的多元回归方程;
以B22井的标准数据为基础,利用EXCEL软件建立钙元素曲线与录井工程曲线的多元回归方程:
y=-125.936-1.319x1+5.929x2+3.27x3-0.117x4-0.059x2 2-0.11x3 2 (1)
其中式(1)中y为钙元素的质量分数,x1为钻压数据,x2为转盘转速数据,x3为扭矩数据,x4为钻时数据。
104:高分辨率钙元素曲线的模拟。
B22井古生界碳酸盐岩潜山段地层相邻的实测录井钙元素之间,以录井工程参数为基础,利用公式(1)计算出对应深度钙元素的质量分数,将其与实测录井钙元素的质量分数按照深度顺序组合起来,从而实现了钙元素数据的加密,进而得到一条采样间距为1米/次的钙元素曲线。
为保证加密钙元素数据的质量,还需对其计算结果进行可靠性检验。本发明实施例以复相关系数R来评价计算得到的钙元素数据的质量,详见下文描述:
一、回归方程可靠性检验。
检验回归方程的可靠性主要依靠复相关系数R。其中,R表明用所建立的回归方程进行钙元素计算的可靠程度,其值越接近1,越可靠。R可由EXCEL软件在进行多元回归时自动计算得出。本发明实施例中的R值为0.736。
二、计算钙元素数据与元素测井测量的钙元素数据对比分析。
在可靠程度检验之后,若满足检验,则可对建立的回归方程进行实际应用。对于B22井古生界碳酸盐岩潜山段地层,即可对该井段的钙元素录井曲线进行加密处理。提取完钻后加测的元素测井钙元素曲线与加密处理钙元素曲线、录井测量的钙元素曲线进行对比,见图1。三者的整体变化趋势相近,但是加密处理的钙元素曲线相比录井测量的钙元素曲线,曲线的变化形态更接近于高分辨率的元素测井钙元素曲线,表明经过加密处理后,加密处理的钙元素曲线的纵向分辨率有了较明显的提高。
同时,如图2所示的交会图分析计算得到钙元素数据与测井测量钙元素数据的复相关系数可达0.71,说明计算得到的钙元素数据与测井测量的钙元素数据有较好的相关性,表明本发明实施例所计算得到钙元素数据较为准确,能够有效的提高钙元素录井曲线的纵向分辨率。
尽管上面结合附图对本发明的功能及工作过程进行了描述,但本发明并不局限于上述的具体功能和工作过程,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还以做出很多形式,这些均属于本发明的保护之内。
Claims (2)
1.一种提高元素录井曲线纵向分辨率的方法,其特征在于,包括以下步骤:
1)收集目标井的录井工程资料、元素录井资料;其中,所述录井工程资料具体为:钻时数据、钻压数据、扭矩数据、转盘转速数据;
2)选取目标元素曲线的标准数据:同一层位目标元素的录井实测数据及其等深度的录井工程数据;
3)以标准数据为基础,建立目标元素曲线与录井工程曲线之间的多元回归方程;多元回归方程具体为:通过EXCEL软件建立的目标元素曲线与录井工程曲线的关系式;
目标元素为钙元素时,多元回归方程为:
y=-125.936-1.319x1+5.929x2+3.27x3-0.117x4-0.059x2 2-0.11x3 2 (1)
式中,y为钙元素的质量分数,x1为钻压数据,x2为转盘转速数据,x3为扭矩数据,x4为钻时数据;
4)在相邻实测目标元素数据之间,将录井工程曲线值代入多元回归方程,完成对目标元素曲线的数据加密,从而提高目标元素曲线的纵向分辨率;
目标元素为钙元素时,以录井工程参数为基础,利用公式(1)计算出对应深度钙元素的质量分数,将其与实测录井钙元素的质量分数按照深度顺序组合起来,从而实现了钙元素数据的加密。
2.根据权利要求1所述的提高元素录井曲线纵向分辨率的方法,其特征在于,通过纵向分辨率更高的测井元素曲线对加密后的目标元素曲线进行评价。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910634199.6A CN110513098B (zh) | 2019-07-15 | 2019-07-15 | 一种提高元素录井曲线纵向分辨率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910634199.6A CN110513098B (zh) | 2019-07-15 | 2019-07-15 | 一种提高元素录井曲线纵向分辨率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110513098A CN110513098A (zh) | 2019-11-29 |
CN110513098B true CN110513098B (zh) | 2022-07-01 |
Family
ID=68623181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910634199.6A Active CN110513098B (zh) | 2019-07-15 | 2019-07-15 | 一种提高元素录井曲线纵向分辨率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110513098B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101344001A (zh) * | 2008-08-05 | 2009-01-14 | 中国石化集团华北石油局 | 石油钻井中x射线荧光陆源碎屑岩孔隙度的分析方法 |
CN104612675A (zh) * | 2015-02-02 | 2015-05-13 | 中石化西南石油工程有限公司地质录井分公司 | 一种碳酸盐岩地层随钻岩性快速识别方法 |
CN105064986A (zh) * | 2015-07-24 | 2015-11-18 | 中国石油天然气集团公司 | 利用常规测录井资料构建储层四性关系谱的方法 |
CN105675635A (zh) * | 2015-12-31 | 2016-06-15 | 中国石油天然气股份有限公司 | 致密岩石组分相对含量及脆性指数确定方法和装置 |
CN105700034A (zh) * | 2016-02-19 | 2016-06-22 | 中国石油集团川庆钻探工程有限公司 | 利用地层元素迁移指数进行储层识别的方法 |
CN107313770A (zh) * | 2017-07-07 | 2017-11-03 | 中国海洋石油总公司 | 一种利用元素录井特征剖面进行地层划分的新方法 |
CN107489417A (zh) * | 2017-07-25 | 2017-12-19 | 中国海洋石油总公司 | 基于xrf录井的岩矿组成反演法砂泥岩岩性识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7857047B2 (en) * | 2006-11-02 | 2010-12-28 | Exxonmobil Upstream Research Company | Method of drilling and producing hydrocarbons from subsurface formations |
-
2019
- 2019-07-15 CN CN201910634199.6A patent/CN110513098B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101344001A (zh) * | 2008-08-05 | 2009-01-14 | 中国石化集团华北石油局 | 石油钻井中x射线荧光陆源碎屑岩孔隙度的分析方法 |
CN104612675A (zh) * | 2015-02-02 | 2015-05-13 | 中石化西南石油工程有限公司地质录井分公司 | 一种碳酸盐岩地层随钻岩性快速识别方法 |
CN105064986A (zh) * | 2015-07-24 | 2015-11-18 | 中国石油天然气集团公司 | 利用常规测录井资料构建储层四性关系谱的方法 |
CN105675635A (zh) * | 2015-12-31 | 2016-06-15 | 中国石油天然气股份有限公司 | 致密岩石组分相对含量及脆性指数确定方法和装置 |
CN105700034A (zh) * | 2016-02-19 | 2016-06-22 | 中国石油集团川庆钻探工程有限公司 | 利用地层元素迁移指数进行储层识别的方法 |
CN107313770A (zh) * | 2017-07-07 | 2017-11-03 | 中国海洋石油总公司 | 一种利用元素录井特征剖面进行地层划分的新方法 |
CN107489417A (zh) * | 2017-07-25 | 2017-12-19 | 中国海洋石油总公司 | 基于xrf录井的岩矿组成反演法砂泥岩岩性识别方法 |
Non-Patent Citations (1)
Title |
---|
基于神经网络的元素录井岩性识别方法研究及应用;张志虎等;《石油工业技术监督》;20190228;第35卷(第2期);第57-60页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110513098A (zh) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101782294B1 (ko) | 탄층 메탄가스 저류층의 가스 생산 모델링방법 | |
EP3414566B1 (en) | Thermal maturity determination of rock formations using mud gas isotope logging | |
CN105445441B (zh) | 储层含气饱和度确定方法及装置 | |
CN110965991B (zh) | 基于人工智能的随钻岩石矿物组分识别方法及装置 | |
CN107505344A (zh) | 利用“最小二乘积”法的岩性解释方法 | |
CN108374657A (zh) | 井断点自动识别方法 | |
CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
CN110593855B (zh) | 一种测井电阻率的校正方法、油层的识别方法 | |
CN103306672B (zh) | 一种预测页岩地层不同钻进方向的研磨性的方法 | |
CN107727829B (zh) | 一种花岗岩型铀矿剥蚀程度半定量估算方法 | |
CN110513098B (zh) | 一种提高元素录井曲线纵向分辨率的方法 | |
CN110529106B (zh) | 一种利用测井资料确定煤层显微组分含量的方法 | |
CN109738955B (zh) | 一种基于成分-结构分类下的变质岩岩性综合判别方法 | |
Guidry et al. | Devonian shale formation evaluation model based on logs, new core analysis methods, and production tests | |
CN116721223A (zh) | 一种基于倾斜摄影和随钻测量技术的岩体稳定性分析方法 | |
CN113361111B (zh) | 一种低渗透油藏储层特征模型的划分方法 | |
CN112343588B (zh) | 一种基于xrf元素录井随钻中获取岩石泊松比的方法 | |
CN115929285A (zh) | 一种基于拉格朗日支持向量机算法的地温梯度预测方法 | |
CN106285653A (zh) | 确定页岩气储层的有机孔隙度的方法 | |
CN112184033B (zh) | 一种碳酸盐岩岩层断层封闭性评价方法 | |
CN110924942B (zh) | 一种录井作业中随钻评价岩层孔渗性的方法 | |
CN111206923B (zh) | 一种利用钻能确定节理岩体模量比与强度比的测试方法 | |
Zhao et al. | Pore-scale hydraulic properties of virtual sandstone microstructures: spatial variations and voxel scale effects | |
FANG et al. | Intelligent recognition of tunnel stratum based on advanced drilling tests | |
CN113311502A (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 |