[發(fā)明專利]一種基于測井信息對(duì)地層層位進(jìn)行智能化劃分的方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201711310475.0
• 申請(qǐng)日: 2017-12-11
• 公開（公告）號(hào): CN108562950B
• 公開（公告）日: 2020-01-03
• 發(fā)明（設(shè)計(jì)）人: 王勇軍;姜艷玲;羅利;趙中明;胡振平;陳邦定;陳海祥;黃宏;張樹東;熊宇
• 申請(qǐng)（專利權(quán)）人: 中國石油天然氣集團(tuán)公司;中國石油集團(tuán)川慶鉆探工程有限公司
• 主分類號(hào): G01V9/00
• 分類號(hào): G01V9/00
• 代理公司: 51211 成都天嘉專利事務(wù)所(普通合伙)
• 代理人: 冉鵬程

1.一種基于測井信息對(duì)地層層位進(jìn)行智能化劃分的方法，其特征在于，包括如下步驟：

a、找曲線波峰、波谷；

b、計(jì)算兩口井之間的相對(duì)深度差和測量坐標(biāo)差；

c、對(duì)輸入曲線進(jìn)行平滑；

d、單曲線相似度檢測；

e、多曲線相似度檢測；

f、建立基于單曲線地層分層模型；

g、建立基于多曲線地層分層模型；

所述步驟d中，單曲線相似度檢測方法如下：

假設(shè)參考模板R，是一個(gè)M維的向量，即R＝{R(1)，R(2)，……，R(m)，……，R(M)}，每個(gè)分量是一個(gè)數(shù)或者是一個(gè)更小的向量；有一被測試模板T，是一個(gè)N維向量，即T＝{T(1)，T(2)，……，T(n)，……，T(N)}同樣每個(gè)分量是一個(gè)數(shù)或者是一個(gè)更小的向量，M不一定等于N，每個(gè)分量的維數(shù)相同；

由于M不一定等于N，要計(jì)算R和T的相似度，首先計(jì)算R的每一個(gè)分量和T中的每一個(gè)分量之間的距離，形成一個(gè)M*N的矩陣；

所述步驟e中，多曲線相似度檢測方法如下：

(1)多維點(diǎn)的定義：

Ri這個(gè)多維點(diǎn)的分量為(a1,a2,…,an)，有n種分量；

Tj這個(gè)多維點(diǎn)的分量為(b1,b2,…,bn)，有n種分量；

其中，ai,bi(i＝1,..,n)屬于同一種參量；

(2)多維點(diǎn)的歸一化：

把各個(gè)參量歸一化，多維點(diǎn)的所有參數(shù)都是在0～1的范圍之間；

對(duì)每一維的參量序列，如c＝{ci}進(jìn)行如下歸一化：

(3)多維點(diǎn)之間距離的計(jì)算

d(i,j)表示以左下點(diǎn)為起始坐標(biāo)原點(diǎn)的i行j列值，也就是原始序列R的第i個(gè)點(diǎn)和測試模版T的第j個(gè)點(diǎn)之間的歐式距離，有

所述步驟f中，建立基于單曲線地層分層模型：

已知(1):A測井曲線：A＝{ai}＝{a1,a2,…a8},其中ai都是列向量；

已知(2):B測井曲線：B＝{bi},其中bi都是列向量；

已知(3):A的地質(zhì)分層與人工分層L_A＝{L_A1,L_A2,L_A3,L_A4,L_A5}，其中L_Ai都是列向量，分別代表：地質(zhì)分層名稱，地質(zhì)分層深度，地質(zhì)分層厚度，人工分層深度，人工分層厚度；

已知(4):B的地質(zhì)分層L_B＝{L_B1,L_B2,B₃}，其中L_Bi都是列向量，分別代表：地質(zhì)分層名稱，地質(zhì)分層深度，地質(zhì)分層厚度；

對(duì)A的任意某j層，通過查找L_A的第j行，得到該層的分層數(shù)據(jù)：

L_A(j，1)地質(zhì)分層名稱,L_A(j，2)地質(zhì)分層深度，L_A(j，3)地質(zhì)分層厚度,L_A(j，4)人工分層深度，L_A(j，5)人工分層厚度；

根據(jù)L_A(j，4)人工分層深度，以及a₁即A井的深度序列，N＊1矩陣，行數(shù)代表序號(hào)，數(shù)據(jù)代表具體以米為單位的深度，N是總的測井?dāng)?shù)據(jù)點(diǎn)個(gè)數(shù)，找到第j層在A中的測量行序號(hào)k；

在A上取k-w:step:k+w曲線離散點(diǎn)作為一個(gè)匹配模版pattenA；

找到B上具有與L_A(j，1)相同地質(zhì)分層名稱的地質(zhì)分層序號(hào)m，也就是B上的第m層，或B矩陣的第m行,進(jìn)一步找到B的地質(zhì)分層結(jié)果：深度

L_B(m，2)與厚度L_B(m，3)

根據(jù)L_B(m，2)地質(zhì)分層深度，在B上取深度為L_B(m，2)-w_S1到

L_B(m，2)+w_B2的一段測井曲線x作為尋找最佳匹配的范圍；

x這段曲線的下限不應(yīng)該低于第m-1層的深度，即

L_B(m，2)-w_S1＞L_B(m-1，2)；同樣，x這段曲線的上限不應(yīng)該大于第m+1層的深度，即L_B(m，2)+wS2＜L_B(m+1，2)；w_B1和w_B2一般取50米；

得到B井搜索區(qū)域x的深度范圍后，進(jìn)一步根據(jù)b₁得到X的序號(hào)范圍；

從x這段曲線，從開始到結(jié)尾的length(x)個(gè)點(diǎn)進(jìn)行遍歷，第i個(gè)點(diǎn)p＝xi處，以p為中心，左右各偏移wa的深度并以step跳躍取點(diǎn)，得到一個(gè)B上的匹配模版pattenB，然后計(jì)算pattenA和pattenB的相似距離，結(jié)果為d_i(i＝1，2，3，...，length(x))；

定位所有d_i中，d_k等于min({d_i})的那個(gè)序號(hào)k，那么，子曲線x第k個(gè)序號(hào)對(duì)應(yīng)的深度值，就是得到的B上的分層點(diǎn)，它和A上的第j個(gè)分層點(diǎn)對(duì)應(yīng)的深度值一一對(duì)應(yīng)；

設(shè)最小距離曲線的點(diǎn)有y個(gè)，取這些點(diǎn)前per％的最小距離點(diǎn)，形成新的序列z；再對(duì)z根據(jù)序號(hào)排序后，看z的第一列在數(shù)學(xué)整數(shù)上的連續(xù)情況，連續(xù)多的序號(hào)序列，包含真實(shí)分層點(diǎn)的可能性更大，如果取的per％越大，則點(diǎn)數(shù)更多，點(diǎn)的組數(shù)量也就可能更多，如果要給出K個(gè)參考分層點(diǎn)，就取連續(xù)點(diǎn)數(shù)量排在前K個(gè)的組平均點(diǎn)序號(hào)，得到序號(hào)后，比如序號(hào)為n，則B井的計(jì)算機(jī)自動(dòng)分層深度就是b1(n),也就是B井深度矩陣的第n行對(duì)應(yīng)的數(shù)字，單位為米；至此，完成了自動(dòng)分層；

所述步驟g中，建立基于多曲線地層分層模型：與單曲線地層分層模型建立一致，其差別僅在于匹配曲線與被匹配曲線不再是單一曲線，而是曲線組，因此，對(duì)pattenA和pattenB的相似距離計(jì)算采用多維計(jì)算。