張魯川，盧雙舫，肖佃師，黃文彪，李進步，郭思祺

（中國石油大學（華東）非常規(guī)油氣與新能源研究院，山東 青島 266580）

0 引言

美國頁巖油藏開發(fā)成功促進世界范圍內(nèi)頁巖油勘探開發(fā)進程，以頁巖油為代表的非常規(guī)油氣資源成為全球油氣勘探開發(fā)的新亮點[1－4].有機質(zhì)在頁巖油的生成和聚集中發(fā)揮重要作用，富有機質(zhì)泥頁巖層中含有海量油氣資源，它不僅是頁巖油生成的基礎(chǔ)，而且生烴過程中有機質(zhì)產(chǎn)生的有機微（納）米孔隙—微裂縫系統(tǒng)也是頁巖油重要的儲集空間和吸附載體[5].此外，頁巖油的分布不受構(gòu)造控制，有機碳質(zhì)量分數(shù)（TOC）成為頁巖油富集程度的主控因素，如北美Bakken組、Eagle Ford組頁巖油層的有機碳質(zhì)量分數(shù)普遍大于2.0%，最高達12.1%[6－7].因此，準確刻畫泥頁巖中有機碳質(zhì)量分數(shù)的空間分布成為頁巖油勘探開發(fā)中優(yōu)選有利區(qū)的關(guān)鍵.

目前，泥頁巖有機碳質(zhì)量分數(shù)評價主要有實驗分析和測井評價兩種手段，各有局限性.實驗分析價格昂貴但數(shù)據(jù)量小，難以評價有機碳質(zhì)量分數(shù)在井剖面上的變化；利用測井—地化響應(yīng)模型可方便、系統(tǒng)地確定目的層有機碳質(zhì)量分數(shù)的縱向變化[8－10]，但橫向預(yù)測能力不足，不能滿足泥頁巖強有機非均質(zhì)背景下頁巖油資源評價的需求.以松遼盆地南部小城子地區(qū)青山口組為例，筆者利用井震聯(lián)合反演技術(shù)刻畫泥頁巖有機碳質(zhì)量分數(shù)的空間分布，分析頁巖油資源空間分布特征，對松遼盆地南部頁巖油的勘探開發(fā)具有指導(dǎo)意義.

1 區(qū)域地質(zhì)概況

小城子地區(qū)位于松遼盆地東南隆起區(qū)王府凹陷南部，緊鄰釣魚臺隆起和青山口背斜，三者呈現(xiàn)“兩隆夾一凹”的構(gòu)造格局.該地區(qū)構(gòu)造演化主要經(jīng)歷隆起—裂谷—坳陷—抬升階段，呈現(xiàn)下斷上坳的雙層沉積結(jié)構(gòu).鉆井資料揭示，該地區(qū)自下而上沉積下白堊統(tǒng)沙河子組、營城組、登樓庫組和泉頭組，以及上白堊統(tǒng)青山口組、姚家組、嫩江組、四方臺組和明水組等地層[11].

青山口組是王府凹陷重要的烴源巖層，整體發(fā)育一套深湖—半深湖相泥巖、粉砂質(zhì)泥巖夾油頁巖地層.有機質(zhì)類型主要為Ⅱ—Ⅰ型，鏡質(zhì)體反射率Ro為0.4%～0.8%，屬低熟—成熟演化階段，生氣量不大，以生油為主[12].城深1和城4井地化分析資料顯示，該地區(qū)青山口組熱解烴含量（S1）在0.01～0.89 mg／g之間，平均為0.38mg／g；有機碳質(zhì)量分數(shù)在0.16%～5.91%之間，平均為2.02%.

2 單井評價方法

2.1 測井—地化響應(yīng)模型

氣候變化和構(gòu)造運動等因素控制沉積環(huán)境變化，并影響泥頁巖中有機質(zhì)的分布和富集，使有機碳質(zhì)量分數(shù)分布在空間上表現(xiàn)為強非均質(zhì)性[13]，在測井響應(yīng)上，富有機質(zhì)泥頁巖常表現(xiàn)出高聲波時差、高電阻率、高自然伽馬、低密度等特征.測井資料具有縱向連續(xù)和高分辨率的優(yōu)勢，可以利用測井資料評價泥頁巖中有機碳質(zhì)量分數(shù)在井剖面上的變化.由EXXON／ESSO石油公司推導(dǎo)和實驗得出的ΔlgR模型（測井—地化響應(yīng)模型）[9]在烴源巖有機碳質(zhì)量分數(shù)評價中得到廣泛應(yīng)用.該技術(shù)利用聲波時差和電阻率曲線之間的幅度差（ΔlgR）刻畫泥頁巖中有機碳質(zhì)量分數(shù)的變化.隨著有機碳質(zhì)量分數(shù)的增大，聲波時差和電阻率呈現(xiàn)增大趨勢，且兩者的幅度差也增大，有機碳質(zhì)量分數(shù)TOC的計算公式為

式中：R、Rj為實測電阻率和基線對應(yīng)的電阻率；Δt、Δtj為實測聲波時差和基線對應(yīng)的聲波時差；K0為疊合因數(shù)，K0＝0.02；LOM為有機質(zhì)成熟度，由樣品分析得到.

2.2 改進模型

ΔlgR模型存在缺陷：首先，ΔlgR的計算依賴于基線值，而人為確定基線的過程繁瑣、主觀性強；其次，Passey Q R等[9]在推導(dǎo)疊合因數(shù)K的過程中應(yīng)用大量的經(jīng)驗公式和系數(shù)，僅具區(qū)域代表性，普適性差，利用疊合因數(shù)計算不同地區(qū)的有機碳質(zhì)量分數(shù)誤差較大；最后，對于地化分析資料不足的地區(qū)，LOM參數(shù)不易選取，導(dǎo)致有機碳質(zhì)量分數(shù)計算出現(xiàn)整體誤差.

對ΔlgR模型進行改進，建立一個操作簡便且適用于研究區(qū)的計算模型.由式（1）可知，疊合因數(shù)K具有兩方面意義：物理意義表示單位聲波時差所對應(yīng)的對數(shù)坐標下電阻率的單位數(shù)；地質(zhì)意義表示將有量綱的（Δt－Δtj）轉(zhuǎn)化為量綱為一的數(shù)，使之與lg（R／Rj）量級相當，共同構(gòu)成兩條曲線之間的幅度差ΔlgR.因此，動態(tài)疊合因數(shù)K表示為

假定基線后，可得

式中：Rmax、Rmin為疊合時電阻率刻度的最大和最小值；Δtmax、Δtmin為疊合時聲波時差刻度的最大和最小值.

改進后的模型表示為

式中：A為調(diào)節(jié)因子；B為有機碳背景值，即ΔlgR為0時的有機碳質(zhì)量分數(shù).

優(yōu)選動態(tài)疊合因數(shù)K是ΔlgR模型改進的關(guān)鍵，它制約烴類流體與干酪根的相對識別能力，并且在一定程度上抵消孔隙度對有機質(zhì)測井響應(yīng)的影響[14].

研究區(qū)共有鉆探井17口，對城深1井和城4井進行地化分析測試，其中城深1井測井響應(yīng)質(zhì)量好、實測地化數(shù)據(jù)較多.利用城深1井改進ΔlgR模型計算有機碳質(zhì)量分數(shù)結(jié)果，與實測數(shù)據(jù)的相關(guān)因數(shù)隨疊合因數(shù)變化曲線見圖1，隨著K增大，有機碳質(zhì)量分數(shù)計算和實測值的相關(guān)性呈現(xiàn)先增大后減小的趨勢，其中相關(guān)因數(shù)最大值對應(yīng)的K值即為最優(yōu)疊合因數(shù).

利用改進的ΔlgR模型評價有機碳質(zhì)量分數(shù)在縱向上的分布（見圖2）.由圖1選取最優(yōu)疊合因數(shù)K＝0.032 8，構(gòu)建有機碳質(zhì)量分數(shù)與ΔlgR的關(guān)系式為

有機碳質(zhì)量分數(shù)計算結(jié)果與實測數(shù)據(jù)之間的相關(guān)因數(shù)為0.86，遠大于K＝0.02時兩者的相關(guān)因數(shù)0.59.

圖1 城深1井改進模型有機碳質(zhì)量分數(shù)計算結(jié)果與實測數(shù)據(jù)的相關(guān)因數(shù)隨疊合因數(shù)變化曲線Fig.1 Relevance coefficient R2 between the calculated and measured data and its vary with superposition coefficient Kof well Chengshen

圖2 城深1井青山口組測井地化結(jié)果Fig.2 Logging geochemistry results of the Qingshankou formation in well Chengshen

3 有機碳質(zhì)量分數(shù)預(yù)測

采用改進的ΔlgR模型評價單井有機碳質(zhì)量分數(shù)，并通過內(nèi)插法預(yù)測有機碳質(zhì)量分數(shù)平面分布，易忽略泥頁巖有機非均質(zhì)性的影響，導(dǎo)致井間距較大或無井控制區(qū)域的有機碳質(zhì)量分數(shù)平面預(yù)測誤差較大，難以滿足復(fù)雜地質(zhì)背景下頁巖油勘探的需要[15].地震資料含有豐富的地層信息，橫向連續(xù)性好，但縱向識別能力有限；測井資料具有較高的縱向分辨率.結(jié)合兩者優(yōu)勢，利用井震聯(lián)合反演技術(shù)進行有機碳質(zhì)量分數(shù)空間預(yù)測，能夠提高有機碳質(zhì)量分數(shù)縱向預(yù)測精度和平面預(yù)測可信度.

3.1 有機碳質(zhì)量分數(shù)—波阻抗相關(guān)性分析

文中地質(zhì)統(tǒng)計學反演中將有機碳質(zhì)量分數(shù)參數(shù)視為區(qū)域化變量，需要進行有機碳質(zhì)量分數(shù)與地層縱波阻抗相關(guān)性分析.有機質(zhì)在測井響應(yīng)上表現(xiàn)為高聲波時差、低密度特征，泥頁巖有機碳質(zhì)量分數(shù)越高，其縱波阻抗越低.對松南小城子地區(qū)青山口組測井評價有機碳質(zhì)量分數(shù)與縱波阻抗進行交會分析（見圖3），顯示青山口組有機碳質(zhì)量分數(shù)隨地層波阻抗增大而減小，兩者的相關(guān)因數(shù)為0.81，滿足利用波阻抗對有機碳質(zhì)量分數(shù)進行地質(zhì)統(tǒng)計學反演的條件.

圖3 松南小城子地區(qū)青山口組有機碳質(zhì)量分數(shù)—縱波阻抗交會圖Fig.3 Crossplot between P－impedance and TOC of K1qnin Xiaochengzi area，southern Songliao basin

3.2 地質(zhì)統(tǒng)計學反演

地質(zhì)統(tǒng)計學反演是一種基于模型的反演方法，結(jié)合測井高縱向分辨和地震橫向連續(xù)的優(yōu)勢，具有在井點和井間分別符合測井解釋數(shù)據(jù)及原始地震資料的特點，反演結(jié)果具有較高的垂向分辨率，且更符合地質(zhì)規(guī)律[16－17].地質(zhì)統(tǒng)計學隨機反演易產(chǎn)生多個等概率數(shù)據(jù)體，結(jié)果多解性強，因此，需要通過優(yōu)選合理的變差函數(shù)，以及在反演過程中融合稀疏脈沖確定性反演，提高隨機反演的精度.

3.2.1 變差函數(shù)分析

變差函數(shù)是指區(qū)域化變量Z（x）在相距為Δx的兩個空間點x和x＋Δx的參數(shù)值Z（x）與Z（x＋Δx）的方差均值.作為刻畫區(qū)域化變量Z（x）空間差異性的基本工具，變差函數(shù)反映區(qū)域化變量空間變異程度隨距離變化的特征.通過計算不同方向的變差函數(shù)，表征有機碳空間分布的非均質(zhì)性，發(fā)揮測井數(shù)據(jù)和稀疏脈沖約束反演波阻抗數(shù)據(jù)體在垂向和橫向分辨率的優(yōu)勢.

求取變差函數(shù)是地質(zhì)統(tǒng)計學反演的關(guān)鍵步驟之一.在垂直方向上，可以利用測井資料采樣密集的特點獲取高分辨率的變差函數(shù).在水平方向上，由于井間距往往較大，可供用于擬合的數(shù)據(jù)點稀少、點距大，擬合的變差函數(shù)精度低，根據(jù)稀疏脈沖約束反演波阻抗數(shù)據(jù)體計算水平變差函數(shù)[18].

研究區(qū)青山口組有機碳質(zhì)量分數(shù)的垂向和水平方向變差函數(shù)曲線見圖4.由圖4（a）可見，垂向（z方向）變程僅為4.9m，反映青山口組有機碳質(zhì)量分數(shù)具有縱向變化快、非均質(zhì)性強的特點.由圖4（b）可見，不同水平方向變差函數(shù)曲線變程不同，反映青山口組有機碳質(zhì)量分數(shù)平面分布的差異性.結(jié)合目的層稀疏脈沖反演波阻抗沿層切片，確定x、y方向的變程分別為400.0和1 000.0m.

3.2.2 反演流程

泥頁巖有機碳質(zhì)量分數(shù)地質(zhì)統(tǒng)計學反演流程包括：（1）在測井曲線標準化[19]、層位標定及目的層低頻模型建立的基礎(chǔ)上，通過目的層稀疏脈沖約束波阻抗反演，得到反映泥頁巖變化趨勢的波阻抗反演數(shù)據(jù)體.（2）分析有機碳質(zhì)量分數(shù)的垂向、水平方向變差函數(shù)，優(yōu)選隨機反演參數(shù)，包括變程、塊金效應(yīng)及基臺值[20]等.（3）進行地質(zhì)統(tǒng)計學反演，包括隨機模擬和隨機反演過程，其中隨機模擬采用序貫高斯配置協(xié)模擬[21]，隨機反演采用模擬退火算法[16]，得到地質(zhì)意義上多個有機碳質(zhì)量分數(shù)數(shù)據(jù)體.（4）優(yōu)選有機碳質(zhì)量分數(shù)數(shù)據(jù)體，遵循原則：剖面上有機碳變化趨勢與地震資料的相符；平面上有機碳變化趨勢符合地質(zhì)沉積規(guī)律，與地震屬性反映的波阻抗趨勢基本一致；后驗井檢驗符合率較高.

圖4 研究區(qū)青山口組有機碳質(zhì)量分數(shù)垂向與水平方向變差函數(shù)Fig.4 Vertical and horizontal variogram of TOC of K1qnin the study area

由反演剖面（見圖5）可以看出，研究區(qū)青山口組地質(zhì)統(tǒng)計學反演（見圖5（b））較稀疏脈沖約束反演（見圖5（a））垂向分辨率明顯提高，薄層狀有機碳反映較好；有機碳主要呈層狀富集，其中青一段為有機碳主要富集分布層，有機碳質(zhì)量分數(shù)最高為9.20%，該段也是頁巖油資源主要發(fā)育段；有機碳質(zhì)量分數(shù)縱向變化快，橫向變化趨勢與地震剖面上的基本一致，空間非均質(zhì)性強.由平面（見圖6）可以看出，研究區(qū)西部和城10井附近（見圖6（a））有機碳質(zhì)量分數(shù)相對較高，與研究區(qū)構(gòu)造沉積規(guī)律一致.通過3口后驗井進行反演質(zhì)量檢驗，測井評價有機碳質(zhì)量分數(shù)數(shù)據(jù)與反演結(jié)果符合率為82.7%，反演效果較好.

圖5 松南小城子地區(qū)青山口組反演剖面Fig.5 Inversion profile of K1qnin Xiaochengzi area，southern Songliao basin

圖6 研究區(qū)青一段反演與測井評價頁巖油富集資源段厚度Fig.6 Enriched resources section thickness of shale oil according to inversion data and well logging evaluation data of K1qn1 in the study area

4 應(yīng)用實例

利用有機碳質(zhì)量分數(shù)地質(zhì)統(tǒng)計學反演數(shù)據(jù)體，參照頁巖油資源分級評價標準[22]，統(tǒng)計青山口組有機碳質(zhì)量分數(shù)分布（見表1），刻畫目的層頁巖油資源空間分布.由表1可知，青一段有機碳質(zhì)量分數(shù)以大于2.0%為主，屬頁巖油富集資源段，頁巖油資源潛力大，可以作為當前頁巖油勘探重點層段；青二、三段有機碳質(zhì)量分數(shù)分布較低，以小于2.0%為主，主要屬頁巖油分散或低效資源段，頁巖油勘探潛力較低.

表1 松南小城子地區(qū)青山口組有機碳質(zhì)量分數(shù)分布Table 1 Frequency distribution of TOC of K1qnin Xiaochengzi area，southern Songliao basin

基于地質(zhì)統(tǒng)計學反演有機碳質(zhì)量分數(shù)數(shù)據(jù)體，可以得到研究區(qū)青一段頁巖油富集資源泥頁巖厚度（見圖6（a）），與由測井評價數(shù)據(jù)得到的頁巖油富集資源段厚度（見圖6（b））相比，井震聯(lián)合反演方法更好地刻畫井間有機碳質(zhì)量分數(shù)的復(fù)雜變化，更接近真實地質(zhì)規(guī)律.頁巖油富集資源段泥頁巖厚度的復(fù)雜變化也證實該地區(qū)青山口組存在較強的有機非均質(zhì)性.由圖6（a）可看出，研究區(qū)西部和城10井區(qū)有機碳質(zhì)量分數(shù)相對較高，屬頁巖油富集資源分布區(qū).

5 結(jié)論

（1）針對原始ΔlgR（測井—地化響應(yīng)）模型的不足，通過自動選取基線、動態(tài)優(yōu)選疊合因數(shù)等優(yōu)化，改進ΔlgR模型，提高對泥頁巖有機碳質(zhì)量分數(shù)刻畫的精度和可操作性.

（2）井震聯(lián)合反演方法有效預(yù)測研究區(qū)青山口組有機碳質(zhì)量分數(shù)空間分布，通過后驗井檢驗，測井評價有機碳質(zhì)量分數(shù)數(shù)據(jù)與反演結(jié)果吻合率可達82.7%.該方法為少井或勘探初期地區(qū)泥頁巖有機碳質(zhì)量分數(shù)空間預(yù)測提供新的思路.

（3）根據(jù)井震聯(lián)合反演結(jié)果，參照頁巖油資源分級評價標準，研究區(qū)青山口組頁巖油富集資源主要分布于青一段；平面上，西部和城10井區(qū)為頁巖油富集分布區(qū).

（References）：

[1]鄒才能，楊智，崔景偉，等.頁巖油形成機制、地質(zhì)特征及發(fā)展對策[J].石油勘探與開發(fā)，2013，40（1）：14－26.Zou Caineng，Yang Zhi，Cui Jingwei，et al.Formation mechanism，geological characteristics and development strategy of nonmarine shale oil in China[J].Petroleum Exploration and Development，2013，40（1）：14－26.

[2]林森虎，鄒才能，袁選俊，等.美國致密油開發(fā)現(xiàn)狀及啟示[J].巖性油氣藏，2011，23（4）：25－30.Lin Senhu，Zou Caineng，Yuan Xuanjun，et.al.Status quo of tight oil exploitation in the United States and its implication[J].Lithologic Reservoirs，2011，23（4）：25－30.

[3]康玉柱.中國非常規(guī)泥頁巖油氣藏特征及勘探前景展望[J].天然氣工業(yè)，2012，32（4）：1－5.Kang Yuzhu.Characteristics and exploration prospect of unconventional shale gas reservoirs in China[J].Natural Gas Industry，2012，32（4）：1－5.

[4]賈承造，鄭民，張永峰.中國非常規(guī)油氣資源與勘探開發(fā)前景[J].石油勘探與開發(fā)，2012，39（2）：129－136.Jia Chengzao，Zheng Min，Zhang Yongfeng.Unconventional hydrocarbon resourcesin China and the prospect of exploration and development[J].Petroleum Exploration and Development，2012，39（2）：129－136.

[5]吳偉，劉惟慶，唐玄，等.川西坳陷富有機質(zhì)頁巖孔隙特征[J].中國石油大學學報：自然科學版，2014，38（4）：1－8.Wu Wei，Liu Weiqing，Tang Xuan，et al.Organic rich shale pore characteristics of western Sichuan depression[J].Journal of China University of Petroleum：Natural Science Edition，2014，38（4）：1－8.

[6]Mullen J，Lowry J C，Nwabuoku K C.Lessons learned developingthe Eagle Ford Shale[C].SanAntonio：Tight Gas Completions Conference，2010.

[7]Schmoker J W，Hester T C.Organic carbon in Bakken formation，United States portion of Williston basin[J].AAPG Bulletin，1983，67（12）：2165－2174.

[8]張寒，朱光有.利用地震和測井信息預(yù)測和評價烴源巖——以渤海灣盆地富油凹陷為例[J].石油勘探與開發(fā)，2007，34（1）：55－59.Zhang Han，Zhu Guangyou.Using seismic and log information to predict and evaluate hydrocarbonsource rocks：An example from rich oil depressions in Bohai bay[J].Petroleum Exploration and Development，2007，34（1）：55－59.

[9]Passey Q R，Creaney S.A practical model for organic richness from porosity and resistivity logs[J].AAPG Bulletin，1990，74（12）：1777－1794.

[10]胡慧婷，盧雙舫，劉超，等.測井資料計算源巖有機碳含量模型對比及分析[J].沉積學報，2011，29（6）：1199－1205.Hu Huiting，Lu Shuangfang，Liu Chao，et al.Models for calculating organic carbon content from logging information：Comparison and analysis[J].Acta Sedimentologica Sinica，2011，29（6）：1199－1205.

[11]郭英海，李壯福，李熙哲，等.松遼盆地王府凹陷的沉積充填及生儲蓋組合[J].中國礦業(yè)大學學報，2001，30（1）：32－36.Guo Yinghai，Li Zhuangfu，Li Xizhe，et al.Sedimentary filling and source－reservoir－seal combination at Wangfu depression[J].Journal of China University of Mining ＆ Technology，2011，30（1）：32－36.

[12]趙佳楠，陳永進，姜文斌.松遼盆地南部白堊系青山口組頁巖氣儲層評價及生儲有利區(qū)預(yù)測[J].東北石油大學學報，2013，37（2）：26－36.Zhao Jia＇nan，Chen Yongjin，Jiang Wenbin.Shale gas reservoir evaluation and prediction of the shale gas reservoir＇s favorable zone in the cretaceous Qingshankou formation of the southern Songliao basin[J].Journal of Northeast Petroleum University，2013，37（2）：26－36.

[13]朱光有，金強.烴源巖的非均質(zhì)性研究——以東營凹陷牛38井為例[J].石油學報，2002，23（5）：34－39.Zhu Guangyou，Jin Qiang.Study on source rock heterogeneity：A case of Niu－38well in Dongying depression[J].Acta Petrolei Sinica，2002，23（5）：34－39.

[14]劉超，盧雙舫，黃文彪，等.ΔlogR技術(shù)改進及其在烴源巖評價中的應(yīng)用[J].大慶石油地質(zhì)與開發(fā)，2011，30（3）：27－31.Liu Chao，Lu Shuangfang，Huang Wenbiao，et al.Improvement ofΔlogRand its application in source rocks evalution[J].Petroleum Geology and Oilfield Development in Daqing，2011，30（3）：27－31.

[15]張營革.泥頁巖油氣藏有機碳含量地球物理預(yù)測方法——以濟陽坳陷羅家地區(qū)沙三段下亞段為例[J].油氣地質(zhì)與采收率，2013，20（5）：64－67.Zhang Yingge.Geophysical prediction method for organic carbon content in shale oil and gas reservoir：A case of lower Es3submember of Luojia area，Jiyang depression[J].Petroleum Geology and Recovery Efficiency，2013，20（5）：64－67.

[16]孫思敏，彭仕宓.基于模擬退火算法的地質(zhì)統(tǒng)計學反演方法[J].石油地球物理勘探，2007，42（1）：38－43.Sun Simin，Peng Shimi.Inversion of geostatistics based on simulated annealing algorithm[J].Oil Geophysical Prospecting，2007，42（1）：38－43.

[17]李方明，計智鋒，趙國良，等.地質(zhì)統(tǒng)計反演之隨機地震反演方法——以蘇丹 M盆地P油田為例[J].石油勘探與開發(fā)，2007，34（4）：451－455.Li Fangming，Ji Zhifeng，Zhao Guoliang，et al.Methodology and application of stochastic seismic inversion：A case from P oilfield，M basin，Sudan[J].Petroleum Exploration and Development，2007，34（4）：451－455.

[18]孫思敏，彭仕宓.地質(zhì)統(tǒng)計學反演方法及其在薄層砂體預(yù)測中的應(yīng)用[J].西安石油大學學報：自然科學版，2007，22（1）：41－48.Sun Simin，Peng Shimi.Geostatistical inversion method and its application in the prediction of thin reservoirs[J].Journal of Xi′an Shiyou University：Natural Science Edition，2007，22（1）：41－48.

[19]肖佃師，黃文彪，張小剛，等.王府凹陷青山口組含油氣泥頁巖層的測井曲線標準化[J].東北石油大學學報，2014，38（1）：46－53.Xiao Dianshi，Huang Wenbiao，Zhang Xiaogang，et al.Normalization of logging data of oil bearing mud shale of Qingshankou formation in Wangfu depression[J].Journal of Northeast Petroleum University，2014，38（1）：46－53.

[20]王云專，付檢剛，曲忠勇，等.變差函數(shù)在儲層地質(zhì)建模中的應(yīng)用[J].科學技術(shù)與工程，2010，10（29）：7147－7150.Wang Yunzhuan，F(xiàn)u Jiangang，Qu Zhongyong，et al.The application of variogram in reservoir geological modeling[J].Science Technology and Engineering，2010，10（29）：7147－7150.

[21]姜曉宇，計智鋒，毛鳳軍，等.隨機模擬技術(shù)在儲層物性研究中的應(yīng)用[J].地球物理學進展，2014，29（6）：2665－2668.Jiang Xiaoyu，Ji Zhifeng，Mao Fengjun，et al.Application of random simulation method in the study of physical property of reservoir[J].Progress in Geophysics，2014，29（6）：2665－2668.

[22]盧雙舫，黃文彪，陳方文，等.頁巖油氣資源分級評價標準探討[J].石油勘探與開發(fā)，2012，39（2）：249－256.Lu Shuangfang，Huang Wenbiao，Chen Fangwen，et al.Classification and evaluation criteria of shale oil and gas resources：Discussion and application[J].Petroleum Exploration and Development，2012，39（2）：249－256.