CN113123783B - 一种精细评价地层体积含量的方法 - Google Patents
一种精细评价地层体积含量的方法 Download PDFInfo
- Publication number
- CN113123783B CN113123783B CN201911412325.XA CN201911412325A CN113123783B CN 113123783 B CN113123783 B CN 113123783B CN 201911412325 A CN201911412325 A CN 201911412325A CN 113123783 B CN113123783 B CN 113123783B
- Authority
- CN
- China
- Prior art keywords
- log
- fluid
- mineral
- volume
- objective function
- 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 63
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 82
- 239000011707 mineral Substances 0.000 claims abstract description 82
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 44
- 230000004044 response Effects 0.000 claims abstract description 14
- 239000012530 fluid Substances 0.000 claims description 74
- 238000004364 calculation method Methods 0.000 claims description 22
- 239000002002 slurry Substances 0.000 claims description 17
- 229910052770 Uranium Inorganic materials 0.000 claims description 16
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 claims description 16
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 claims description 15
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 claims description 15
- 229910052776 Thorium Inorganic materials 0.000 claims description 15
- 238000010521 absorption reaction Methods 0.000 claims description 15
- 239000011591 potassium Substances 0.000 claims description 15
- 229910052700 potassium Inorganic materials 0.000 claims description 15
- 238000012360 testing method Methods 0.000 claims description 7
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 239000004698 Polyethylene Substances 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 125000003275 alpha amino acid group Chemical group 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims description 2
- 239000000523 sample Substances 0.000 claims description 2
- 238000005553 drilling Methods 0.000 abstract description 6
- 230000009286 beneficial effect Effects 0.000 abstract description 5
- 238000004422 calculation algorithm Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 5
- 238000011161 development Methods 0.000 abstract description 4
- 238000002474 experimental method Methods 0.000 abstract description 3
- 238000005755 formation reaction Methods 0.000 description 38
- 235000019738 Limestone Nutrition 0.000 description 13
- 239000006028 limestone Substances 0.000 description 13
- 238000005457 optimization Methods 0.000 description 11
- 239000011159 matrix material Substances 0.000 description 7
- 239000013598 vector Substances 0.000 description 6
- 238000013178 mathematical model Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 4
- 239000011148 porous material Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000011028 pyrite Substances 0.000 description 2
- 229910052683 pyrite Inorganic materials 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 229910021532 Calcite Inorganic materials 0.000 description 1
- 230000003416 augmentation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000010433 feldspar Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000002945 steepest descent method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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
- 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 And Detection Of Objects (AREA)
Abstract
本发明提供一种精细评价地层体积含量的方法。利用常规测井资料，通过构造地层体积模型下的常规响应方程，利用线性规划和非线性相结合方法计算页岩地层体积含量。通过构建适当的目标函数，利用测井响应方程的不等式约束，应用单纯性和变尺度算法有利于消除扩径条件下的钻井泥浆的不确定影响，快速、准确地反演地层体积含量，减少了用元素测井等费用昂贵的测井项目，减少了岩心实验确定页岩矿物类型和含量的过程，节省了费用和时间，为加快国内页岩地层勘探、开发起积极作用。
Description
技术领域
本发明涉及石油勘探开发领域，具体涉及一种精细评价地层体积含量的方法。
背景技术
矿物体积含量为每种矿物所占地层中的体积比，在测井解释模型中通常用到的是体积含量模型。矿物体积含量解释方法在Schlumberger、Halliburton和Western-Atlas公司有用于矿物体积含量反演的计算方法。Quirein等(2010年)用优化方法直接建立了页岩地层元素到矿物体积含量的计算；雍世和(1995年)应用优化算法实现了用常规测井资料计算地层矿物体积含量。这些反演通常用在渗透性地层中评价其矿物含量，而对于页岩等非渗透性地层或渗透性异常低的地层，反演效果很差，反演出的矿物类型和体积含量误差较大，因此需要寻找新的方法来准确地进行地层矿物体积含量反演。
线性规划(Linear Programming，缩写为LP)是运筹学的重要分支之一，在实际中应用得较广泛，其方法也较成熟，它是辅助人们进行科学管理的一种数学方法。利用线性规划方法解决问题的目标函数是多个决策变量的线性函数，求目标函数的最大值或最小值；解决问题的约束条件是一组多个决策变量的线性不等式或等式。这些问题与测井过程中所遇到的问题是一致的，运用线性规划方法为解决地层矿物体积含量提供了一种新的途径。
建立数学模型是线性规划的关键一步。线性规划数学模型(Mathematical Modelof Linear Programming)由决策变量(Decision variables)、目标函数(Objectivefunction)和约束条件(Constraints)三个要素构成。一般假设线性规划数学模型中有m个约束、n个决策变量xj，j＝1，2，…，n，目标函数的变量系数用cj表示，cj称为价值系数。约束条件的变量系数用aij表示，aij称为工艺系数。约束条件右端的常数用bi表示，bi称为资源系数。则线性规划数学模型的一般表达式可写成：
max(min) Z＝c1x1+c2x2+...+cnxn
该条件可以写成：
xj≥0(j＝1,2,…,n) (3)
c＝(c1,c2,...,cn,)为价值变量，
x＝(x1,x2,...,xn,)为决策变量。
由于线性规划问题包含等式约束、不等式约束和目标函数，这三个函数组合成增广矩阵。其系数矩阵A的任意一个m×m阶的非奇异的子方阵B(即|B|≠0)，称为LP问题的一个基阵或基：
系数矩阵A中的其他部分称为非基阵：
其中，Pi(i＝1,2,…,m)为基向量；Xi(i＝1,2,…,m)为基变量；Pj(j＝m+1,…,n)为非基向量；Xj(j＝m+1,…,n)为非基变量。
采用目标函数最大值，采用增加变量的方式，把不等式转换成等式，因此线性规划问题变为标准形式为：
xj≥0(i＝1,2,…,n)
其中，Z(0)＝cBB-1b
σj＝cj-cBB-1N
σ′ij＝B-1N
bi′＝B-1bi
单纯形法(Simplex Method)是1947年由G.B.Dantzig提出，是解线性规划问题最有效的算法之一，且已成为整数规划和非线性规划某些算法的基础。在用单纯法求解线性规划问题时，为了讨论问题方便，需将线性规划模型化为统一的标准形式。
单纯形法求解线性规划问题的计算步骤如下：
步骤1：将一般形式转化为标准形式；
步骤2：从标准形式中求出初始基可行解，建立初始单纯形表。对标准形式的LP，在约束条件式的变量的系数矩阵中总会存在一个单位矩阵：
其中：P1,P2,.....,Pm称为基向量，同其对应的变量x1,x2,........,xm称为基变量，模型中其他变量xm+1,xm+2,......xn称为非基变量。若令所有非基变量为0，求出基变量的值，可以得到初始其可行解，将其数据代入单纯形表中，可以得到初始单纯形表；
步骤3：检验各非基变量Xj的检验数σj,如果所有的σj≤0(j＝1,2,…,n)，则已求得最优解，停止计算。否则转入下一步；
步骤4：在所有的σj>0中，如果有某个σk>0，所对应的xk的系数列向量p’k≤0(即a’ik≤0，i＝1,2,…,m)，则此问题解无界，停止计算。否则转入下一步；
步骤5：当有一个以上的检验数大于0时，根据σk＝max{σj|aj＞0,1≤j≤n},其对应的变量xk作为换入的非基变量，称为换入基变量。又根据最小比值法则计算:确定xr是换出的基变量，元素ark决定了从一个基可行解到相邻基可行解的转移去向，称为(取名)主元素。转入下一步；
步骤6：以ark为主元进行换基变换，用初等行变换将Xk所对应的列向量变换成单位列向量，即同时将检验数行中的第k个元素也变换为零，这样从一个基可行解转换到相邻的目标函数值更大的基可行解，可得新的单纯形表。返回步骤2。
单纯性解法当所求未知数较少时，通常能得出较好的结果，但如果未知数较多时，线性规划问题常常出现无解的情况。利用常规测井资料来评价地层矿物含量，除了地层矿物未知外，还包含泥浆性质、井眼环境等，未知数较大(大于10以上)，对线性规划方法带来一定的挑战。
非线性规划的解法常用的有梯度法(又称最速下降法)、牛顿法、抛物线插值法和变尺度法。
发明内容
鉴于上述现有技术中存在的问题，本发明的目的在于提供一种精细评价地层体积含量，尤其是地层中矿物体积含量和/或流体体积含量的方法。本发明利用常规测井资料，考虑在不规则井眼条件下，构造地层体积模型下的常规响应不等式方程，通过构建适当的目标函数，设定线性规划中的决策变量、约束条件和价值系数，首先利用线性规划方法中的单纯性方法计算页岩地层体积含量，把该计算结果作为非线性规划反演的初始值，最后利用非线性规划方法进行无约束极值解。应用线性规划和非线性规划相结合有利于消除扩径条件下的钻井泥浆的不确定影响，快速、准确地反演地层体积含量。
本发明提供一种精细评价地层体积含量的方法，包括：
(1)获得包括实际测井曲线和/或测试数据在内的常规测井资料；
(2)根据地层中实际含有和/或可能含有的矿物和/或流体确定决策变量；
(3)不考虑井眼环境时，将所述常规测井资料和所述决策变量响应为不等式方程，并以所述不等式方程作为约束条件；
(4)考虑井眼环境时，将所述常规测井资料和所述决策变量响应为等式方程，并以所述等式方程作为约束条件；
(5)利用实际测井曲线与理论计算曲线之差的累计和，构建线性目标函数；
(6)将步骤(3)或步骤(4)中所述约束条件和步骤(5)中所述线性目标函数转换为线性规划标准形式；
(7)对所述线性规划标准形式进行求解，从而初步得出地层中矿物体积含量和/或流体体积含量，优选为利用单纯形法对所述线性规划标准形式进行求解；
(8)判断步骤(7)中所得出的地层中矿物体积含量和/或流体体积含量是否为最优求解结果；
(9)当步骤(7)中所得出的地层中矿物体积含量和/或流体体积含量不是最优求解结果时，修改所述决策变量，并重复步骤(3)～(8)；
(10)当步骤(7)中所得出的地层中矿物体积含量和/或流体体积含量是最优求解结果时，输出线性规划的计算结果；
(11)增加实际矿物和/或流体类型，并将步骤(10)中输出的线性规划的计算结果作为非线性规划的初始值；
(12)利用实际测井曲线与理论计算曲线之差的累计平方和，构建非线性目标函数；
(13)将步骤(3)或步骤(4)中所述约束条件和步骤(12)中所述非线性目标函数转换为非线性规划标准形式；
(14)对所述非线性规划标准形式进行求解，从而得出地层中矿物体积含量和/或流体体积含量，优选为利用变尺度法对所述非线性规划标准形式进行求解。
本申请的发明人经研究发现，地层矿物类型、含量复杂，通常要用元素测井资料来准确确定。但由于地层钻开后，存在扩径或缩径现象，尤其是扩径时，测井仪器所探测到测井曲线受到钻井泥浆的影响，同时钻井泥浆侵入地层后使地层的性质发生改变，此时测井曲线不能真实地反映地层信息。通常的做法是对这些曲线进行井眼环境校正，由于泥浆的性质难以确定，所以很难准确地完成井眼环境校正。在此情况下，利用常规测井资料进行最优化反演，实际是进行无约束下的非线性规划反演，得出的地层体积含量存在较大的误差，甚至得出错误的结论。
本发明的技术方案利用常规测井资料，通过构造地层体积模型下的常规响应方程，首先利用线性规划方法中的单纯性方法计算页岩地层体积含量，把该计算结果作为非线性规划反演的初始值，最后利用非线性规划方法进行无约束极值解。通过构建适当的目标函数，利用测井响应方程的不等式约束，应用单纯性算法有利于消除扩径条件下的钻井泥浆的不确定影响，快速、准确地提供地层矿物体积含量初始值，为了获得很多的矿物数量和流体性质，利用非线性规划方法进一步求解。减少了用元素测井等费用昂贵的测井项目，减少了岩心实验确定页岩矿物类型和含量的过程，节省了费用和时间，为加快国内页岩地层勘探、开发起积极作用。
在本发明的一些优选的实施方式中，所述测试数据包括探测深度和井径；所述实际测井曲线包括自然伽马测井曲线、自然电位测井曲线、补偿声波测井曲线、补偿中子测井曲线、补偿密度测井曲线、光电吸收截面指数测井曲线、铀测井曲线、钍测井曲线和钾测井曲线中的至少一种。
在本发明的一些优选的实施方式中，所述不等式方程如式(1)所示：
式(1)中，为理论计算曲线，LOGj表示第j种测井曲线，Vi表示第i种矿物或流体的体积，MCij表示第i种矿物或流体所对应的第j种测井理论数值，其中，1≥Vi≥0，
在本发明的一些优选的实施方式中，
所述不等式方程选自下述式(2)至式(8)中的至少一种：
其中，DEN表示补偿密度测井曲线，AC表示补偿声波测井曲线，CNL表示补偿中子测井曲线，PE表示光电吸收截面测井曲线，U表示铀测井曲线，TH表示钍测井曲线，K表示钾测井曲线，Vi表示第i种矿物或流体的体积，ρi表示第i种矿物或流体所对应的补偿密度，ΔTi表示第i种矿物或流体所对应的补偿声波，CNLi表示第i种矿物或流体所对应的补偿中子，PEi表示第i种矿物或流体所对应的光电吸收截面波，Ui表示第i种矿物或流体所对应的铀值，THi表示第i种矿物或流体所对应的钍值，Ki表示第i种矿物或流体所对应的对应的钾值，且满足Vi≥0。
在本发明的一些优选的实施方式中，
所述等式方程如式(9)所示：
式(9)中，为理论计算曲线，LOGj表示第j种测井曲线，Vi表示第i种矿物或流体的体积，MCij表示第i种矿物或流体所对应的第j种测井理论数值，MCmf表示泥浆对应的测井曲线的响应值，ΔV表示井眼扩径或缩径所形成的体积占测井曲线探测体积的比例，其中，/>1≥Vi,ΔV≥0。
在本发明的一些优选的实施方式中，
所述等式方程选自下述式(10)至式(16)中的至少一种：
其中，DEN表示补偿密度测井曲线，AC表示补偿声波测井曲线，CNL表示补偿中子测井曲线，PE表示光电吸收截面测井曲线，U表示铀测井曲线，TH表示钍测井曲线，K表示钾测井曲线，Vi表示第i种矿物或流体的体积，ΔV表示井眼扩径或缩径所形成的体积占测井曲线探测体积的比例，ρi表示第i种矿物或流体所对应的补偿密度，ΔTi表示第i种矿物或流体所对应的补偿声波，CNLi表示第i种矿物或流体所对应的补偿中子，PEi表示第i种矿物或流体所对应的光电吸收截面波，Ui表示第i种矿物或流体所对应的铀值，THi表示第i种矿物或流体所对应的钍值，Ki表示第i种矿物或流体所对应的对应的钾值，ρmf表示泥浆对应的补偿密度，ΔTmf表示泥浆对应的补偿声波，CNLmf表示泥浆对应的补偿中子，PEmf表示泥浆对应的光电吸收截面波，Umf表示泥浆对应的铀值，THmf表示泥浆对应的钍值，Kmf表示泥浆对应的钾值，且满足Vi≥0。
在本发明的一些优选的实施方式中，所述线性目标函数如式(17)所示：
在本发明的一些优选的实施方式中，所述非线性目标函数如式(18)所示：
在本发明的一些优选的实施方式中，上述方法还包括：
确定所述线性目标函数中的价值系数，并将所述线性目标函数中的价值系数应用于所述线性规划标准形式；和/或
确定所述非线性目标函数中的价值系数，并将所述非线性目标函数中的价值系数应用于所述非线性规划标准形式，其中，
所述线性目标函数中的价值系数如式(19)所示；和/或所述非线性目标函数中的价值系数如式(20)所示：
根据本发明，所述线性规划标准形式有两种方式，
其一为：
1≥Vi≥0
任选地，
其二为：
1≥Vi,ΔV≥0
任选地，根据本发明，所述线性规划标准形式还可以用下述两种方式表示，其一为：
1≥Vi≥0
任选地，
其二为：
1≥Vi,ΔV≥0
任选地，
根据本发明，所述非线性规划标准形式有两种方式，
其一为：
1≥Vi≥0
任选地，
其二为：
1≥Vi,ΔV≥0
任选地，
根据本发明，所述非线性规划标准形式还可以用下述两种方式表示，
其一为：
1≥Vi≥0
任选地，
其二为：
1≥Vi,ΔV≥0
任选地，
在本发明的一些优选的实施方式中，步骤(8)中，若所得出的地层中矿物体积含量和/或流体体积含量均大于等于0，则将所得出的地层中矿物体积含量和/或流体体积含量判断为最优解。
本发明的有益效果至少在于：
本专利利用常规测井资料，通过构建适当的目标函数，利用测井响应方程的不等式约束，应用线性规划和非线性规划相结合有利于消除扩径条件下的钻井泥浆的不确定影响，快速、准确地反演地层体积含量，减少了用元素测井等费用昂贵的测井项目，减少了岩心实验确定页岩矿物类型和含量的过程，节省了费用和时间，为加快国内页岩地层勘探、开发起积极作用。
附图说明
图1是实施例1的评价地层中矿物体积含量的方法流程图。
图2是实施例1中砂岩地层的常规测井曲线。
图3是实施例1中砂岩地层反演体积含量结果图示。
图4是实施例1中利用国外软件对砂岩地层反演体积含量结果图示。
图5是实施例1中页岩地层的常规测井曲线。
图6是实施例1中页岩地层反演体积含量结果图示。
图7是利用国外软件对页岩地层反演体积含量结果图示。
图8是实施例1中灰岩地层的常规测井曲线。
图9是实施例1中灰岩地层反演体积含量结果图示。
图10是利用国外软件对灰岩地层反演体积含量结果图示。
附图中，岩性体积分析一列中涉及到的缩写及其代表的含义为：
PORT-总孔隙度；VDOL-白云石；VCLC-方解石；VSML-其它；VPYR-黄铁矿；VPRD-长石；VQVA-石英；VSH-泥质。
具体实施方式
以下通过实施例对本发明进行详细说明，但本发明的保护范围并不限于下述说明。
下述实施方式中，所涉及到的缩写，例如，DEN、AC和CNL等均具有与上文相同的指代意义，在此不做赘述。
实施例1
(1)获得常规测井资料，包括深度(图2、图5和图8中第一道)、地层分析(图2、图5和图8中第二道，本实施例在计算时未使用该项数据)、自然伽马测井曲线(图2、图5和图8中第三道)、补偿声波测井曲线(图2、图5和图8中第四道)、补偿中子测井曲线(图2、图5和图8中第五道)、补偿密度测井曲线(图2、图5和图8中第六道)、光电吸收截面指数测井曲线(图2、图5和图8中第七道)、铀测井曲线(图2、图5和图8中第八道)、钍测井曲线(图2、图5和图8中第九道)、钾测井曲线(图2、图5和图8中第十道)和电阻率曲线(图2、图5和图8中第十一道，本实施例在计算时未使用该项数据)。
(2)假设某地层包含泥岩、砂岩、灰岩、有机质、黄铁矿和孔隙度等决策变量，即线性规划中基解变量有6个未知数，也即有6个决策变量，具体地，为V1、V2、V3、V4、V5和V6。
(3)不考虑井眼环境时，将所述常规测井资料和所述决策变量响应为不等式方程，并以所述不等式方程作为约束条件，具体地为：
/>
1≥Vi≥0
(5)利用实际测井曲线与理论计算曲线之差的累计和，构建线性目标函数：
其中，MCij值如表1所示。
(6)确定所述线性目标函数中的价值系数为
(7)把测井线性规划问题转换为线性规划标准形式：
/>
1≥Vi≥0
(8)利用单纯形法对所述线性规划标准形式进行求解，从而初步得出地层中矿物体积含量和/或流体体积含量，如图3、图6和图8中第一道所示。
(9)增加实际矿物/和或流体类型，并将步骤(8)中输出的线性规划的计算结果作为非线性规划的初始值；
(10)利用实际测井曲线与理论计算曲线之差的累计平方和，构建非线性目标函数：
(11)确定非线性目标函数中的价值系数
(12)把测井非线性规划问题转换为非线性规划标准形式：
1≥Vi≥0
(13)利用变尺度法对所述线性规划标准形式进行求解，从而得出地层中矿物体积含量和/或流体体积含量，如图3、图6和图8中第二道所示。
表1.地层矿物和流体含量测井曲线响应骨架值
该地层包含砂岩地层、页岩地层和灰岩地层，以下是各自的处理结果。
计算结果显示砂岩地层、页岩地层和灰岩地层的曲线道是一致的。图4、图7和图10为国外软件利用元素测井资料和常规测井资料联合解释的结果。图2、图5和图8中第十一道为优化反演的目标函数。图3、图6和图9中第一道为线性规划的计算结果，第二道为以线性规划结果作为约束条件，运用非线性规划优化解释方法得出的评价结果。详细说明如下：
1.1砂岩地层反演体积含量结果
图3为砂泥岩地层利用线性规划求解初始矿物和流体条件，从图中看出，第一道线性规划方法计算出以砂岩、泥岩和孔隙流体为主，以线性规划结果作为约束条件，运用最优化解释方法得出如第二道所示的矿物成分和流体性质；线性规划计算结果与结合非线性规划优化解释结果相比较，非线性规划优化解释结果增加了束缚水体积含量和少量的灰岩矿物，孔隙度相对减小，而砂岩的体积含量相对增加，与图4中国外软件利用元素测井资料和常规测井资料联合解释的结果一致性较好。
1.2页岩地层反演体积含量结果
图6为页岩地层利用线性规划求解初始矿物和流体条件，从图中看出，第一道线性规划方法计算出以砂岩、泥岩为主，还包含有机质、灰岩和孔隙流体，以线性规划结果作为约束条件，运用最优化解释方法得出如第二道所示的矿物成分和流体性质；线性规划计算结果与结合非线性规划优化解释结果相比较，非线性规划优化解释结果增加了束缚水体积含量和孔隙度相对减小，而砂岩的体积含量相对增加，声波时差、补偿中子、密度、光电吸收截面指数、铀、钍和钾测井曲线的实际测量数值与构造曲线数值匹配良好，并且优化结果与图7中国外软件利用元素测井资料和常规测井资料联合解释的结果一致性较好。
1.3灰岩地层反演体积含量结果
图9为灰岩地层利用线性规划求解初始矿物和流体条件，从图中看出，第十三道线性规划方法计算出以灰岩为主，还包含少量的砂岩、泥岩和孔隙流体，以线性规划结果作为约束条件，运用非线性规划优化解释方法得出如十四道所示的矿物成分和流体性质；线性规划计算结果与结合非线性规划优化解释结果相比较，非线性规划优化解释结果砂岩的体积含量相对增加，灰岩含量和孔隙度相对减小，声波时差、补偿中子、密度、光电吸收截面指数、铀、钍和钾测井曲线的实际测量数值与构造曲线数值匹配良好。该段地层井壁稳定性较差，存在一定的扩径现象，导致实际测井曲线受井眼环境影响，利用线性规划方法计算出地层的主要灰岩含量，以及砂、泥岩含量，为非线性规划优化解释提供了准确的初始解释值，有利于准确优化地层矿物含量和流体性质。
应当注意的是，以上所述的实施例仅用于解释本发明，并不构成对本发明的任何限制。通过参照典型实施例对本发明进行了描述，但应当理解为其中所用的词语为描述性和解释性词汇，而不是限定性词汇。可以按规定在本发明权利要求的范围内对本发明作出修改，以及在不背离本发明的范围和精神内对本发明进行修订。尽管其中描述的本发明涉及特定的方法、材料和实施例，但是并不意味着本发明限于其中公开的特定例，相反，本发明可扩展至其他所有具有相同功能的方法和应用。
Claims (10)
1.一种精细评价地层体积含量的方法，包括：
(1)获得包括实际测井曲线和/或测试数据在内的常规测井资料；
(2)根据地层中实际含有和/或可能含有的矿物和/或流体确定决策变量；
(3)不考虑井眼环境时，将所述常规测井资料和所述决策变量响应为不等式方程，并以所述不等式方程作为约束条件；
(4)考虑井眼环境时，将所述常规测井资料和所述决策变量响应为等式方程，并以所述等式方程作为约束条件；
(5)利用所述实际测井曲线与理论计算曲线之差的累计和，构建线性目标函数；
(6)将步骤(3)或步骤(4)中所述约束条件和步骤(5)中所述线性目标函数转换为线性规划标准形式；
(7)对所述线性规划标准形式进行求解，从而初步得出地层中矿物体积含量和/或流体体积含量，优选为利用单纯形法对所述线性规划标准形式进行求解；
(8)判断步骤(7)中所得出的地层中矿物体积含量和/或流体体积含量是否为最优求解结果；
(9)当步骤(7)中所得出的地层中矿物体积含量和/或流体体积含量不是最优求解结果时，修改所述决策变量，并重复步骤(3)～(8)；
(10)当步骤(7)中所得出的地层中矿物体积含量和/或流体体积含量是最优求解结果时，输出线性规划的计算结果；
(11)增加实际矿物和/或流体类型，并将步骤(10)中输出的线性规划的计算结果作为非线性规划的初始值；
(12)利用实际测井曲线与理论计算曲线之差的累计平方和，构建非线性目标函数；
(13)将步骤(3)或步骤(4)中所述约束条件和步骤(12)中所述非线性目标函数转换为非线性规划标准形式；
(14)对所述非线性规划标准形式进行求解，从而得出地层中矿物体积含量和/或流体体积含量，优选为利用变尺度法对所述非线性规划标准形式进行求解。
2.根据权利要求1所述的方法，其特征在于，所述测试数据包括探测深度和井径中的至少一种；所述实际测井曲线包括自然伽马测井曲线、自然电位测井曲线、补偿声波测井曲线、补偿中子测井曲线、补偿密度测井曲线、光电吸收截面指数测井曲线、铀测井曲线、钍测井曲线和钾测井曲线中的至少一种。
3.根据权利要求1或2所述的方法，其特征在于，所述不等式方程如式(1)所示：
式(1)中，为理论计算曲线，LOGj表示第j种测井曲线，Vi表示第i种矿物或流体的体积，MCij表示第i种矿物或流体所对应的第j种测井理论数值，其中，1≥Vi≥0，/>
4.根据权利要求1-3中任一项所述的方法，其特征在于，所述不等式方程选自式(2)至式(8)中的至少一种：
其中，DEN表示补偿密度测井曲线，AC表示补偿声波测井曲线，CNL表示补偿中子测井曲线，PE表示光电吸收截面测井曲线，U表示铀测井曲线，TH表示钍测井曲线，K表示钾测井曲线，Vi表示第i种矿物或流体的体积，ρi表示第i种矿物或流体所对应的补偿密度，ΔTi表示第i种矿物或流体所对应的补偿声波，CNLi表示第i种矿物或流体所对应的补偿中子，PEi表示第i种矿物或流体所对应的光电吸收截面波，Ui表示第i种矿物或流体所对应的铀值，THi表示第i种矿物或流体所对应的钍值，Ki表示第i种矿物或流体所对应的钾值，且满足Vi≥0。
5.根据权利要求1-4中任一项所述的方法，其特征在于，所述等式方程如式(9)所示：
式(9)中，为理论计算曲线，LOGj表示第j种测井曲线，Vi表示第i种矿物或流体的体积，MCij表示第i种矿物或流体所对应的第j种测井理论数值，MCmf表示泥浆对应的测井曲线的响应值，ΔV表示井眼扩径或缩径所形成的体积占测井曲线探测体积的比例，其中，1≥Vi,ΔV≥0。
6.根据权利要求1-5中任一项所述的方法，其特征在于，所述等式方程选自式(10)至式(16)中的至少一种：
其中，DEN表示补偿密度测井曲线，AC表示补偿声波测井曲线，CNL表示补偿中子测井曲线，PE表示光电吸收截面测井曲线，U表示铀测井曲线，TH表示钍测井曲线，K表示钾测井曲线，Vi表示第i种矿物或流体的体积，ΔV表示井眼扩径或缩径所形成的体积占测井曲线探测体积的比例，ρi表示第i种矿物或流体所对应的补偿密度，ΔTi表示第i种矿物或流体所对应的补偿声波，CNLi表示第i种矿物或流体所对应的补偿中子，PEi表示第i种矿物或流体所对应的光电吸收截面波，Ui表示第i种矿物或流体所对应的铀值，THi表示第i种矿物或流体所对应的钍值，Ki表示第i种矿物或流体所对应的对应的钾值，ρmf表示泥浆对应的补偿密度，ΔTmf表示泥浆对应的补偿声波，CNLmf表示泥浆对应的补偿中子，PEmf表示泥浆对应的光电吸收截面波，Umf表示泥浆对应的铀值，THmf表示泥浆对应的钍值，Kmf表示泥浆对应的钾值，且满足Vi≥0。
7.根据权利要求1-6中任一项所述的方法，其特征在于，所述线性目标函数如式(17)所示：
8.根据权利要求1-7中任一项所述的方法，其特征在于，所述非线性目标函数如式(18)所示：
9.根据权利要求1-8中任一项所述的方法，其特征在于，还包括：
确定所述线性目标函数中的价值系数，并将所述线性目标函数中的价值系数应用于所述线性规划标准形式；和/或
确定所述非线性目标函数中的价值系数，并将所述非线性目标函数中的价值系数应用于所述非线性规划标准形式，其中，
所述线性目标函数中的价值系数如式(19)所示；和/或所述非线性目标函数中的价值系数如式(20)所示：
10.根据权利要求1-9中任一项所述的方法，其特征在于，步骤(8)中，若所得出的地层中矿物体积含量和/或流体体积含量均大于等于0，则将所得出的地层中矿物体积含量和/或流体体积含量判断为最优解。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911412325.XA CN113123783B (zh) | 2019-12-31 | 2019-12-31 | 一种精细评价地层体积含量的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911412325.XA CN113123783B (zh) | 2019-12-31 | 2019-12-31 | 一种精细评价地层体积含量的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113123783A CN113123783A (zh) | 2021-07-16 |
| CN113123783B true CN113123783B (zh) | 2023-10-13 |
Family
ID=76770247
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911412325.XA Active CN113123783B (zh) | 2019-12-31 | 2019-12-31 | 一种精细评价地层体积含量的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113123783B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113123782B (zh) * | 2019-12-31 | 2024-05-14 | 中国石油化工股份有限公司 | 一种利用线性规划评价地层组分的方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105221142A (zh) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物质量含量的方法 |
| CN105298476A (zh) * | 2014-07-04 | 2016-02-03 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物体积含量的方法 |
| WO2017024530A1 (zh) * | 2015-08-11 | 2017-02-16 | 深圳朝伟达科技有限公司 | 一种计算烃源岩中有机碳含量的方法 |
-
2019
- 2019-12-31 CN CN201911412325.XA patent/CN113123783B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105221142A (zh) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物质量含量的方法 |
| CN105298476A (zh) * | 2014-07-04 | 2016-02-03 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物体积含量的方法 |
| WO2017024530A1 (zh) * | 2015-08-11 | 2017-02-16 | 深圳朝伟达科技有限公司 | 一种计算烃源岩中有机碳含量的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 袁超 ; 马萌 ; 周灿灿 ; 冯周 ; 施宇峰 ; .结合地层元素测井和密度测井提高孔隙度评价精度.测井技术.2018,(02),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113123783A (zh) | 2021-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gong et al. | A new approach for characterization and prediction of natural fracture occurrence in tight oil sandstones with intense anisotropy | |
| Santillán et al. | Phase field model of fluid‐driven fracture in elastic media: Immersed‐fracture formulation and validation with analytical solutions | |
| US6662109B2 (en) | Method of constraining by dynamic production data a fine model representative of the distribution in the reservoir of a physical quantity characteristic of the subsoil structure | |
| Wang et al. | Probabilistic identification of underground soil stratification using cone penetration tests | |
| Zhang et al. | Fluid flow through anisotropic and deformable double porosity media with ultra-low matrix permeability: A continuum framework | |
| CN106223942A (zh) | 一种基于测井曲线重构的砾岩油藏泥质含量计算方法 | |
| Van den Hof et al. | Recent developments in model-based optimization and control of subsurface flow in oil reservoirs | |
| CN104272140A (zh) | 用于校准在储藏层建模中使用的渗透性的系统和方法 | |
| CN111027882A (zh) | 一种基于高阶神经网络利用常规测井资料评价脆性指数的方法 | |
| CN105317429A (zh) | 一种确定页岩地层干酪根含量的方法 | |
| CN107808032A (zh) | 页岩气储层含水饱和度计算方法及系统 | |
| Zhang et al. | Characterization of in situ stress state and joint properties from extended leak-off tests in fractured reservoirs | |
| Xiong et al. | Data‐driven design of wave‐propagation models for shale‐oil reservoirs based on machine learning | |
| CN113123783B (zh) | 一种精细评价地层体积含量的方法 | |
| CN115267935A (zh) | 数据与物理模型联合驱动的测井储层评价方法及装置 | |
| Xue et al. | An efficient automatic history matching method through the probabilistic collocation based particle filter for shale gas reservoir | |
| Zhao et al. | Determination of the geomechanical parameters and associated uncertainties in hydraulic fracturing by hybrid probabilistic inverse analysis | |
| Shchipanov et al. | A new approach to deformable fractured reservoir characterization: case study of the Ekofisk field | |
| CN113123782B (zh) | 一种利用线性规划评价地层组分的方法 | |
| Fang et al. | Brittleness index prediction method of tight reservoir based on grey correlation and analytic hierarchical process | |
| Delay et al. | Analytical 2D model to invert hydraulic pumping tests in fractured rocks with fractal behavior | |
| Hubinger et al. | A new equation solver for modeling turbulent flow in coupled matrix‐conduit flow models | |
| Wang et al. | Geophysical and production data history matching based on ensemble smoother with multiple data assimilation | |
| Hui et al. | An Integrated EDFM+ DPDK Hybrid Workflow for Hierarchical Treatment of Fractures in Practical Field Studies | |
| Li et al. | Correlations of permeability and geological characteristics based on mercury intrusion data and hierarchical statistical models: A case study |
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 |