CN101737038A - 一种测量地层岩性和密度的方法 - Google Patents
一种测量地层岩性和密度的方法 Download PDFInfo
- Publication number
- CN101737038A CN101737038A CN200910215889A CN200910215889A CN101737038A CN 101737038 A CN101737038 A CN 101737038A CN 200910215889 A CN200910215889 A CN 200910215889A CN 200910215889 A CN200910215889 A CN 200910215889A CN 101737038 A CN101737038 A CN 101737038A
- Authority
- CN
- China
- Prior art keywords
- density
- stratum
- equation
- formula
- rho
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种测量地层岩性和密度的方法,属于石油勘测技术领域。根据康普顿散射理论,计算每个能窗的计数率响应:ln(C)=a1+a2ρe+a3(ρe)2+a4(ρe)3+a5L+a6x3+a7x4,建立包含4个未知数(ρe、L、x3、x4)的若干个非线性方程;每个测量点运用包含4个未知数的非线性方程组进行计算;用最小二乘法求解方程组,得到ρe、L以及x3、x4;地层密度ρb的公式ρb=1.0704-0.1883,x3、x4与密度补偿值Δρ之间的关系为Δρ=d6x3+d7x4,在全谱岩性密度测井中,用地层和泥饼之间的密度差ρp=d6x3和地层和密度之间的岩性差别Δρm=d7x4作为曲线输出,地层光电指数Pe=4.16L+0.464。采用本发明得到的固定肋角在不同井况条件下测量的准确性高。
Description
技术领域
本发明涉及一种测量地层岩性和密度的方法,属于石油勘测技术领域。
技术背景
目前,用于石油测井的岩性密度测井算法普遍基于康普顿散射理论采用脊肋图法。该算法认为在没有泥饼的条件下,不同密度地层,长、短源距探测器测得的数据,在脊肋图上形成一条直线,即脊线。存在泥饼的时候,短源距探测器源距小,受泥饼影响大,长源距探测器源距大,受泥饼影响小。长、短源距探测器测得的数据,在脊肋图上偏离脊线。一定密度地层,一定密度泥饼,不同泥饼厚度,长、短源距探测器测得的数据,在脊肋图上形成一条曲线,即肋线。通过标定脊线上的点与密度值一一对应。通过刻度,把仪器规一化。测井时,如果测得的数据正好落在脊线上,即可立即读取密度值;如果测得的数据没有落在脊线上,则沿肋线回归到脊线上,读取密度值。
脊肋图的一个演变是双脊双肋法,该算法认为轻泥饼肋线实际上是两个线性段,靠近脊线的线段角度与重泥饼肋线相同,另一线段的角度不同,两条肋线的交点构成第二条脊线。
脊肋图算法或双脊双肋图算法都是在算法里固定了肋线的角度,即便是肋线可以用折线表示,也是对厚泥饼改变了一个固定的角度。实际情况是地层密度在不断变化,泥饼厚度和泥饼密度也是各处不同的,也就是说实际的肋线的肋角是不固定的。
传统算法当泥浆中含有重晶石成分时,岩性不可测。
发明内容
本发明的目的是为了解决固定肋角在不同井况条件下带来测量误差的问题,提出一种测量地层岩性和密度的方法。
本发明的目的是通过以下技术方案实现的。
本发明的一种测量地层岩性和密度的方法,其具体测量步骤为:
根据康普顿散射理论,计算每个能窗的计数率响应如公式(1)所示:
ln(C)=a1+a2ρe+a3(ρe)2+a4(ρe)3+a5L+a6x3+a7x4 (1)
式中C为探测器能窗计数率;ρe为地层电子密度;L为地层岩性系数;a1、a2、a3、a4、a5、a6、a7为系数;x3和x4泥饼影响项;
1)地层密度的计算
根据公式(1),建立包含4个未知数:ρe为地层电子密度、L为地层岩性系数、x3泥饼影响项和x4泥饼影响项的若干个非线性方程,方程数量取决于能谱分能窗的数量;每个测量点运用包含以上4个未知数的非线性方程组进行计算;用最小二乘法求解方程组,得到ρe、L、x3和x4;方程组中包含4个未知数,故称为“四维算法”;其中地层电子密度ρe转换成地层视体积密度即地层密度ρb的公式如公式(2)所示:
ρb=1.0704-0.1883 (2)
上述方法得到的x3、x4与密度补偿值Δρ之间的关系如公式(3)所示:
Δρ=d6x3+d7x4 (3)
式中d6、d7为常数;
在全谱岩性密度测井中,用地层和泥饼之间的密度差ρp=d6x3和地层和密度之间的岩性差别Δρm=d7x4分别作为曲线输出;
2)此外地层岩性计算
利用下列公式计算地层光电指数Pe,如公式(4)所示:其中L为地层岩性系数;
Pe=4.16L+0.464 (4)
由此可得到地层岩性和密度。
有益效果
采用本发明得到的固定肋角在不同井况条件下测量的准确性高。
具体实施方式
一种测量地层岩性和密度的测井仪,其具体测量步骤为:
根据康普顿散射理论,计算每个能窗的计数率响应如公式(1)所示:
ln(C)=a1+a2ρe+a3(ρe)2+a4(ρe)3+a5L+a6x3+a7x4 (1)
式中C为探测器能窗计数率;ρe为地层电子密度;L为地层岩性系数;a1、a2、a3、a4、a5、a6、a7为系数;x3和x4泥饼影响项;
1)地层密度的计算
根据公式(1),建立包含4个未知数(ρe、L、x3、x4)的4个非线性方程,每个测量点运用包含4个未知数的非线性方程组进行计算;用最小二乘法求解方程组,得到ρe、L以及x3、x4;方程组中包含4个未知数,故称为“四维算法”;得到地层电子密度ρe,地层电子密度转换成地层视体积密度即地层密度ρb的公式如公式(2)所示:
ρb=1.0704-0.1883 (2)
上述方法得到的x3、x4与密度补偿值Δρ之间的关系如公式(3)所示:
Δρ=d6x3+d7x4 (3)
式中d6、d7为常数;
在全谱岩性密度测井中,用地层和泥饼之间的密度差ρp=d6x3和地层和密度之间的岩性差别Δρm=d7x4分别作为曲线输出;
2)地层岩性计算
利用下列公式计算地层光电指数Pe,如公式(4)所示:
Pe=4.16L+0.464 (4)
Claims (1)
1.一种测量地层岩性和密度的方法,其特征在于其具体测量步骤为:
根据康普顿散射理论,计算每个能窗的计数率响应如公式(1)所示:
1n(C)=a1+a2ρe+a3(ρe)2+a4(ρe)3+a5L+a6x3+a7x4 (1)
式中C为探测器能窗计数率;ρe为地层电子密度;L为地层岩性系数;a1、a2、a3、a4、a5、a6、a7为系数;x3和x4泥饼影响项;
1)地层密度的计算
根据公式(1),建立包含4个未知数:ρe为地层电子密度、L为地层岩性系数、x3泥饼影响项和x4泥饼影响项的若干个非线性方程,方程数量取决于能谱分能窗的数量;每个测量点运用包含4个未知数的非线性方程组进行计算;用最小二乘法求解方程组,得到ρe、L以及x3、x4;得到地层电子密度ρe,地层电子密度转换成地层视体积密度即地层密度ρb的公式如公式(2)所示:
ρb=1.0704-0.1883 (2)
上述方法得到的x3、x4与密度补偿值Δρ之间的关系如公式(3)所示:
Δρ=d6x3+d7x4 (3)
式中d6、d7为常数;
在全谱岩性密度测井中,用地层和泥饼之间的密度差ρp=d6x3和地层和密度之间的岩性差别Δρm=d7x4分别作为曲线输出;
2)此外地层岩性计算
利用下列公式计算地层光电指数Pe,如公式(4)所示:其中L为地层岩性系数
Pe=4464 (4)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910215889A CN101737038A (zh) | 2009-12-31 | 2009-12-31 | 一种测量地层岩性和密度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910215889A CN101737038A (zh) | 2009-12-31 | 2009-12-31 | 一种测量地层岩性和密度的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101737038A true CN101737038A (zh) | 2010-06-16 |
Family
ID=42461005
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910215889A Pending CN101737038A (zh) | 2009-12-31 | 2009-12-31 | 一种测量地层岩性和密度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101737038A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111894564A (zh) * | 2020-07-07 | 2020-11-06 | 中国石油大学(华东) | 一种基于x射线岩性密度测井计算地层光电吸收截面指数的方法 |
-
2009
- 2009-12-31 CN CN200910215889A patent/CN101737038A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111894564A (zh) * | 2020-07-07 | 2020-11-06 | 中国石油大学(华东) | 一种基于x射线岩性密度测井计算地层光电吸收截面指数的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102562048B (zh) | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 | |
CN104536056B (zh) | 小口径伽马能谱测井装置及数据采集传输和自稳方法 | |
CN106597544B (zh) | 致密油气藏脆性的预测方法及装置 | |
CN110517794A (zh) | 一种建立页岩气藏的埋藏-热演化史图的方法 | |
CN105545284B (zh) | 一种随钻伽马成像数据处理方法 | |
CA2349763A1 (en) | Standoff compensation for nuclear measurements | |
CN115311821A (zh) | 一种基于数字孪生技术的地质灾害智能监测云平台 | |
CN115267928B (zh) | 一种用于随钻元素测井的能谱智能处理方法 | |
CN109669424B (zh) | 一种复杂条件下岩石力学性能试验控制系统及方法 | |
CN110160435A (zh) | 滑坡含水率增载响应比预测参数与方法 | |
CN103674065B (zh) | 一种滑动式测斜仪检测装置 | |
US20060229815A1 (en) | Method and apparatus for shale bed detection in deviated and horizontal wellbores | |
CN109141271A (zh) | 多点式光纤光栅孔底应变计 | |
NO844239L (no) | Fremgangsmaate for bestemmelse av restoljemetning i en undergrunnsformasjon | |
Massey et al. | Relative sea‐level change and postglacial isostatic adjustment along the coast of south Devon, United Kingdom | |
CN103590828A (zh) | 一种录井dck指数法评价地层压力的方法 | |
CN111208564A (zh) | 一种深度域层位标定方法及装置 | |
Hawley et al. | Borehole optical stratigraphy and neutron-scattering density measurements at Summit, Greenland | |
CN111812746B (zh) | 基于地震表面波的钻孔应变仪原位标定方法 | |
CN109444972A (zh) | 一种双能x射线测量地层密度测井装置和方法 | |
CN101737038A (zh) | 一种测量地层岩性和密度的方法 | |
CN107290502B (zh) | 矿山边坡岩体潜在滑移面野外快速判别方法 | |
US8141259B2 (en) | Method of determining the dip of a formation | |
CN108915677A (zh) | 一种岩石切削强度及抗切削强度指标快速确定的方法 | |
US20060208184A1 (en) | Method for determining shale bed boundaries and gamma ray activity with gamma ray instrument |
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 |
Application publication date: 20100616 |