一種地球物理勘探利用測井曲線自動求取巖性曲線的方法
【技術領域】
[0001] 本發明屬于石油地球物理勘探領域，具體涉及一種地球物理勘探利用測井曲線自 動求取巖性曲線的方法，針對測井約束地震反演問題，涉及應用測井資料進行巖性分析識 別技術，可應用于石油地球物理勘探中的巖性反演。
【背景技術】
[0002] 隨著油田勘探開發的深入以及儲層描述技術的發展，大的整裝的構造油氣藏很難 發現，目前重點研究的是復雜隱蔽性的巖性油氣藏，地震反演技術在巖性油氣藏的發現中 發揮了較大的作用。隨著隱蔽巖性油氣藏開發的進程，對地震反演的需求也越來越高，常規 的波阻抗反演技術在勘探時期能夠滿足簡單的巖性識別的要求，但針對精細儲層預測的需 求，波阻抗反演技術很存在一定的多解性，需要開展巖性反演，直接反演巖性數據體，為地 質人員提供直觀的、多解性小的巖性數據體，提高巖性識別的精度。巖性反演效果非常好， 但進行巖性反演時，要求井中有巖性曲線，實際的生產中，只有小部分的探井才進行巖性錄 井和取芯工作，才能有巖性曲線，其他的開發井等都沒有沒有巖性曲線，因此，需要對應用 測井曲線進行巖性曲線求取的技術開展研究工作。
[0003] 現有的應用測井曲線劃分巖性的方法主要是應用伽馬曲線來識別砂泥巖，伽馬低 的認為是砂巖，伽馬高的認為是泥巖，更準確的做法是應用伽馬曲線重新構建一條泥質含 量曲線（V sh)，采用公式如下：（公式1)
[0004]
[0005] Bq :純地層背景值
[0006] B0 = P sd · GRsd (或者 P ls · GRls)
[0007] P b、P sh、P sd、P ls分別是目的層、泥巖層、純砂巖、純石灰巖的體積密度；
[0008] GR、GRsh、GRsd、GRls分別是目的層、泥巖層、純砂巖、純石灰巖的自然伽馬測井。
[0009] 求得泥質含量后，在通過取值來判斷砂巖和泥巖等，該取值一般根據經驗，比如泥 質含量小于50為砂巖，大于50為泥巖。
[0010] 該方法主要用伽馬和密度曲線來求取泥質含量曲線，然后利用泥質含量曲線求取 巖性曲線，經過與實際的巖性錄井資料對比分析，發現存在一些誤差，包括這些純泥巖和純 砂巖等的自然伽馬值的確定也有很大的人為性，因此求取的泥質含量曲線不是很精確，所 以有必要開展泥質含量曲線方面的精確求取研究，需要同時結合其他多條曲線進行巖性曲 線的求取
【發明內容】
[0011] 本發明的目的在于解決上述現有技術中存在的難題，提供一種地球物理勘探利用 測井曲線自動求取巖性曲線的方法，提高巖性曲線的精度。
[0012] 本發明是通過以下技術方案實現的：
[0013] 一種地球物理勘探利用測井曲線自動求取巖性曲線的方法，包括：
[0014] (1)加載測井曲線，包括伽馬、自然電位、電阻率、中子和中子壽命曲線；
[0015] (2)分別對這5條曲線進行校正，得到校正好的曲線；
[0016] (3)從校正好的曲線上分別讀取5條曲線的最大值GMAX1和最小值GMIN 1 ;
[0017] (4)求取曲線泥質含量相對值SH1 :
[0018]
[0019] 式2中的SHLG1為測井曲線，即為步驟（1)中的5條曲線，i = 1，2,3,4,5分別代 表伽馬、自然電位、電阻率、中子和中子壽命曲線，即采用同一個公式計算各曲線的泥質含 量；
[0020] (5)計算泥質含量：
[0021]
[0022] W > > J
[0023] G⑶R為與地層年代有關的經驗參數；
[0024] (6)求得最終的泥質含量 Vsh : Vsh = min (Vshl，Vsh2，…，Vsh5);
[0025] (7)應用步驟（6)求得的最終泥質含量Vsh和實測密度曲線Density聯合求取巖 性曲線；
[0026] (8)輸出巖性曲線。
[0027] 所述步驟⑵中的校正包括：
[0028] 奇異值消除校正和或深度校正和或環境校正。
[0029] 所述步驟（5)中，如果是新地層，則GCUR取3. 7,如果是老地層，則GCUR取2。
[0030] 所述步驟（7)是這樣實現的：
[0031] 采用邏輯判斷：
[0032] 當Vsh < 50且Density〈地層密度值（KG/M3)時候，巖性曲線的取值為數值1，代 表的巖性為砂巖；當V sh > 50且Density >地層密度值（KG/M3)時候，巖性曲線的取值為數 值〇,代表的巖性為泥巖；這樣就構建了一條由數值1和〇表不的巖性曲線，1代表砂巖，〇代 表泥巖；
[0033] 所述砂泥巖密度門檻值是根據具體地層綜合分析后確定的值。
[0034] 與現有技術相比，本發明的有益效果是：常規泥質含量曲線的計算方法，應用的測 井曲線種類少，容易受到特殊巖性段的影響，同時，純的巖性地層密度值得讀取存在人為因 素，導致求得的泥質含量曲線與實際巖性錄井資料存在較大誤差。
[0035] 本發明能夠應用多條測井曲線，包括伽馬、中子、自然電位、電阻率、熱中子等曲線 求取泥質含量曲線，然后再應用泥質含量曲線和密度曲線聯合求取巖性曲線，能夠同時綜 合6條測井曲線所反映的巖性信息，同時考慮了地層年代信息，因此該方法求得的巖性曲 線在精度方面得到了較大提高，能夠在沒有巖性錄井資料的井中開展應用，尤其在大量沒 有做錄井的開發井中應用，可以快速準確的求得巖性曲線。
【附圖說明】
[0036] 圖1某井用公式1方法求得的泥質含量曲線與伽馬對比圖，可以看到由于特殊巖 性段的存在，導致整體求得的泥質含量曲線偏低，也就是說砂巖成分偏高了，該方法計算的 泥質含量曲線存在較大的誤差。
[0037] 圖2某井多條測井曲線及用本發明方法求得的泥質含量曲線，該方法求得的泥質 含量曲線比圖1求得的泥質含量在精度上得到了提高。
[0038] 圖3某井用常規方法和本發明所用方法求得泥質含量曲線對比圖，泥質含量1為 常規方法求得，泥質含量2為本發明方法求得，可以看到泥質含量2曲線更精確，受特殊巖 性段的影響很少。
[0039] 圖4某井用步驟⑦求得的巖性曲線（右1列，黑色為砂巖），最終求得的砂巖與伽 馬曲線和泥質含量曲線有很好的對應關系（對應低伽馬和低泥質含量），同時也受密度曲 線的影響，能