專利名稱：利用井眼圖像、數(shù)字巖石樣品以及多點(diǎn)統(tǒng)計(jì)算法生成數(shù)值假巖心的方法
技術(shù)領(lǐng)域：
本發(fā)明總體涉及利用基本多點(diǎn)模擬統(tǒng)計(jì)(MPS)算法從數(shù)字巖石或巖心樣品和井 眼成像測井圖生成數(shù)值假巖心的方法。更具體地，本專利說明書涉及由計(jì)算的X射線層析 成像(CT掃描)和地層微成像(FMI)測井圖生成3D數(shù)值巖心，并且在這些數(shù)值巖心中執(zhí)行 流動(dòng)模擬以理解所選擇的儲(chǔ)層內(nèi)的流體流動(dòng)路徑和采收率。
背景技術(shù)：
電井眼成像測井儀和聲波井眼成像測井儀廣泛用于對(duì)地下井眼進(jìn)行測井，以定位 和繪制例如地層界面的巖石層之間的邊界，并且目測和定向裂縫和斷裂。由于電測井儀是 帶有固定電極陣列的極板式裝置，所以通常在極板之間具有含漏失信息的間隙。電測井和 聲波測井由于非功能性電極、不足的極板壓力、井眼的不規(guī)律性、巖屑、測井儀偏心、或是差 聲波反射而通常具有帶不良數(shù)據(jù)質(zhì)量的間隔。數(shù)字巖石模擬由2D薄截面、掃描電子顯微鏡(SEM)圖像、計(jì)算機(jī)生成的球形單元、 激光掃描共聚焦顯微鏡圖像、以及各種類型的CT掃描(例如，傳統(tǒng)的微型CT，和同步計(jì)算的 顯微層析成像)構(gòu)造而成。CT掃描是應(yīng)用最廣泛的方法。CT掃描是由圍繞樣品旋轉(zhuǎn)的X 射線源生成的二維(2D)橫截面。由X射線衰減系數(shù)計(jì)算密度。對(duì)連續(xù)橫截面進(jìn)行掃描用 于構(gòu)造樣品的3D圖像。由于密度差在巖石和填充流體的孔隙之間較高，因此CT圖像可用 于顯像巖石孔隙系統(tǒng)。根據(jù)所用的裝置，分辨率在亞毫米到微米的尺度間。對(duì)于地下儲(chǔ)層模擬來說，多點(diǎn)統(tǒng)計(jì)(MPS)用于產(chǎn)生空間地質(zhì)和儲(chǔ)層特性區(qū)域的模 擬。這些方法是用諸如在井眼中測量的已知結(jié)果作為固定或“硬”數(shù)據(jù)的條件模擬，所述固 定或“硬”數(shù)據(jù)在模擬期間被絕對(duì)遵守。MPS利用一維、二維、或三維“訓(xùn)練圖像”作為對(duì)模 型地下特性區(qū)的定量模板。MPS模型從訓(xùn)練圖像獲取地質(zhì)結(jié)構(gòu)并將所述地質(zhì)結(jié)構(gòu)固定到數(shù) 據(jù)位置。這些結(jié)構(gòu)可以為先驗(yàn)地質(zhì)解釋或?yàn)楦拍钅Ｐ汀６帱c(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)(MPS)是一種新的前沿的地質(zhì)統(tǒng)計(jì)學(xué)方法。多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)允許 油藏模擬者通過訓(xùn)練圖像將他們的先驗(yàn)知識(shí)、解釋、或概念模型并入到油藏模擬過程中。這 些訓(xùn)練圖像是被認(rèn)為存在于研究中的儲(chǔ)層內(nèi)的結(jié)構(gòu)/特征的數(shù)值表征。一旦我們具有訓(xùn)練 圖像，MPS可以從訓(xùn)練圖像提取曲線結(jié)構(gòu)或復(fù)雜特征，并將所述曲線結(jié)構(gòu)或復(fù)雜特征賦給采 集樣品/觀察結(jié)果的儲(chǔ)層位置，從而產(chǎn)生更實(shí)際的油藏模型。將訓(xùn)練圖像引入到油藏模型 是一個(gè)里程碑。要注意的是在使用MPS時(shí)具有兩個(gè)要素訓(xùn)練圖像(概念模型)和實(shí)際數(shù) 據(jù)(觀察結(jié)果)。這兩個(gè)要素通常分開。然而，在實(shí)際應(yīng)用中，尤其是生成典型的3D訓(xùn)練圖像已經(jīng)被證明是MPS應(yīng)用中的瓶頸。生成連續(xù)變量訓(xùn)練圖像比生成分類訓(xùn)練圖像甚至更難。具有不同類型的電井眼成像測井儀和聲波井眼成像測井儀用于對(duì)地下井眼進(jìn)行測井以定位和繪制諸如地層界面的巖石層之間的邊界，并且顯像和定向裂縫和斷裂。例如，諸如Schlumberger公司的FMI (地層微成像儀)測井儀的電井眼成像儀可以在水基(導(dǎo)電)泥漿中運(yùn)行，所述測井儀基于50年代已經(jīng)在市場上可買得到的傾角測量儀技術(shù)。電井眼成像測井儀本質(zhì)上是成熟的地層傾角儀。成像測井儀具有在井筒周圍布置在極板上的微電阻率電極，所述極板壓靠在井壁上。從傾角測量儀到井眼成像儀的發(fā)展趨勢已經(jīng)從幾個(gè)電極發(fā)展到多個(gè)極板的復(fù)雜電極排列。參見Hurley，N. F. , 2004,Borehole Images, in Asquith,G.禾口Krygowski,D ；以及參見Basic Well LogAnalysis, 2nd Edition, AAPG Methods in Exploration Series No. 16，p. 151-164。首先將測井儀下入到具有閉合極板的井眼內(nèi)。在開始進(jìn)行測井時(shí)，將任意四個(gè)、六個(gè)、或八個(gè)極板壓靠在井壁上。極板的數(shù)量取決于測井裝置。電流通過電極被迫使進(jìn)入到巖石內(nèi)，并且在所述電流與地層相互作用之后傳感器測量所述電流。原始數(shù)據(jù)包括來自單個(gè)極板或極板對(duì)的多個(gè)電極讀數(shù)、井徑儀讀數(shù)和χ軸、y軸以及Z軸加速儀和磁力儀讀數(shù)。井斜和極板1(測井儀)方位由磁力儀 確定。用于電極和加速儀數(shù)據(jù)的采樣率非常高，通常大約120個(gè)樣品/ft(400個(gè)樣品/m)。井壁的覆蓋區(qū)是電極陣列的寬度、極板數(shù)量和井徑的函數(shù)。通常，在典型的井眼中，40% -80%的井壁被成像。井眼的非成像部分作為極板之間的空帶出現(xiàn)。井眼圖像通過將色彩圖分配給電阻率值的不同面元或范圍而生成。然后，將色彩像素布置在井眼周圍的所述色彩像素適當(dāng)?shù)膸缀挝恢弥?。按照慣例，諸如頁巖或被流體填充的裂縫的低電阻率特征顯示為深色。諸如砂巖和石灰?guī)r的高電阻率特征顯示為棕色、黃色和白色色調(diào)?？梢缘玫絻煞N主要類型的處理過的井眼圖像靜態(tài)和動(dòng)態(tài)。靜態(tài)圖象是已經(jīng)具有應(yīng)用到整個(gè)井的一個(gè)對(duì)比設(shè)定(contrast setting)的靜態(tài)圖象。所述靜態(tài)圖象提供整個(gè)井眼的巖石電阻率的相對(duì)變化的有用視圖。已經(jīng)具有施加在滑動(dòng)窗口內(nèi)的可變對(duì)比度的動(dòng)態(tài)圖象提供諸如孔洞、裂縫和地層界面的提高的特征視圖。動(dòng)態(tài)圖象提供具有極低電阻率(例如，頁巖)和極高電阻率(例如，碳酸鹽和結(jié)晶巖)的巖石內(nèi)的細(xì)微特征。電井眼成像儀的另一個(gè)示例可以在油基(不導(dǎo)電)泥漿中運(yùn)行，具體地，典型的油基泥漿的高泥漿電阻率(大于50ohm-m)不適于大多數(shù)電井眼成像儀。自從2001年開始，斯倫貝謝的OBMI (油基微成像儀)已經(jīng)可用于油基泥漿。這種測井儀通過使電流從每一個(gè)極板上的兩個(gè)較大電極通過進(jìn)入到地層而生成井眼圖像，所述每一個(gè)極板處于高電壓(大約300V)下。具有一系列緊密間隔的按鈕，所述按鈕以兩排每一排5個(gè)的方式位于4個(gè)極板中的每一個(gè)上。井眼圖像由緊密間隔的電極之間的電勢差(電壓降)生成。與井眼的非成像部分相對(duì)應(yīng)的寬間隙通常在在極板之間。另一方面，井眼圖像可以在鉆進(jìn)(例如，隨鉆測井或此后被稱作為“LWD”)期間獲得。斯倫貝謝測井儀的示例是可視地層電阻率(GVR) 測井儀和方位密度中子(ADN)測井儀。GVR使用旋轉(zhuǎn)電極，并且在水基泥漿中工作。ADN由方位密度讀數(shù)生成圖像，并且在任何泥漿中工作。當(dāng)在鉆井期間測井儀旋轉(zhuǎn)時(shí)，井眼被完全覆蓋，且沒有間隙。另一方面，井眼成像儀可以是被公知為井眼井下電視的聲波井眼成像儀，所述聲波井眼成像儀基于60年代首次研發(fā)的技術(shù)(1970年Geophysics，v. 35，p. 254-269中 Zemanek, J. ,Glenn,Ε. Ε.，Norton，L J.，禾口 Caldwell，R. L，1970 的 Formation evaluation by inspection with the borehoIeteIeviewer)0UBI (超聲波井眼成像儀)是斯倫貝謝的用于裸眼井應(yīng)用的主要聲波測井儀。位于井中心的UBI測井儀具有發(fā)射并記錄從井壁彈回的聲波的旋轉(zhuǎn)轉(zhuǎn)換器。兩個(gè)聲波振幅和傳 播時(shí)間被記錄并被處理成圖像。通常，井眼覆蓋率是100%，且在圖像中沒有間隙。然而，當(dāng) 測井儀偏心或井壁不規(guī)則時(shí)，可能產(chǎn)生質(zhì)量差的圖像。其中，巖石物理相可以被認(rèn)為是井眼圖像測井圖上的諸如孔洞、電阻率和導(dǎo)電塊 的特征標(biāo)記。1999年Dehghani等人的具體觀點(diǎn)認(rèn)為在孔洞附近存在增強(qiáng)的孔隙度和滲透 率的區(qū)域。Dehghani，K. ，Harris, P. Μ. ，Edwards, K. Α.禾口 Dees，W. Τ. ，1999，Modeling a vuggy carbonate reservoir :AAPG Bulletin, v. 83,p.19—42。Dehghani等人(1999)利用薄截面、SEM圖像、以及微滲透率測量值證明他們的 觀點(diǎn)。Schindler (2005)和Tanprasat (2005)利用標(biāo)記熒光的巖心圖片的分析顯示密集 的小孔洞優(yōu)先存在于大孔洞附近。參見Schindler，J.，2005，Quantification of vuggy porosity, Indian Basin field, NewMexico Unpublished M. S. thesis, Colorado School of Mines, Golden, CO. ； IMM Tanprasat, S. , 2005, Petrophysical analysis of vuggy porosity in theShu' aiba Formation of the United Arab Emirates :Unpublished M. S. thesis, Colorado School of Mines，Golden，CO.。這種小孔洞在井眼成像測井儀 的分辨率以下，所以所述小孔洞顯現(xiàn)為黑色區(qū)域，而不是在圖像測井圖中顯現(xiàn)為離散孔 洞。如果對(duì)于孔洞型碳酸鹽巖來說這是普遍情況，則電井眼圖像和聲波井眼圖像應(yīng)該在 孔洞附近具有高電導(dǎo)率或低振幅(黑色)區(qū)或暈圈。實(shí)際上，這種特征通常都能觀測到， 例如如圖3中所示。包圍孔洞的高電導(dǎo)率區(qū)和被公知為導(dǎo)電塊的增加的小尺度孔隙度形 成Schlumberger公司BorTex軟件的一部分的基礎(chǔ)，被Russell等人(2002)以及Hassall 等人(2004)引用。參見 Russell, S. D.，Akbar, M.，Vissapragada, B.，和 Walkden，G. Μ.， 2002, Rock types and permeability prediction fromdipmeter and image logs Shuaiba reservoir (Aptian), Abu Dhabi =AAPGBulletin, v.86, p.1709—1732 ；以及參 見 Hassal 1, J. K. ，F(xiàn)erraris, P. ，Al-Raisi Μ. ，Hurley, N. F. ，Boyd, Α.禾口 Allen，D. F.， 2004, Comparison of permeabilitypredictors from NMR, formation image and other logs in a carbonatereservoir :SPE preprint 88683, presented at the Ilth Abu Dhabi InternationalPetroleum Exhibition and Conference, Abu Dhabi, U. Α. Ε. ,10-13 October。Delhomme (1992)論證了在井眼圖像中繪制電阻率和非電阻塊的重要性。 Delhomme,J. P. ,1992,A quantitative characterization of formationheterogeneities based on borehole image analysis :Trans. 33rd SymposiumSPWLA, Paper T0 然而，由于 極板之間的間隙使得他的方法應(yīng)用受限。由于形狀的不確定性，Delhomme在高電阻區(qū)或低 電阻區(qū)周圍不能繪制閉合等高線。全井眼圖像(圖4)卻允許在井眼圖像中的電阻和/或 非電阻區(qū)周圍繪制閉合等高線。這種區(qū)域尤其在碳酸鹽巖儲(chǔ)層中可以提供多相性的重要測 量。這些區(qū)域通常比例如從巖石的CT掃描獲得的數(shù)字巖石或巖心樣品大得多。因此，如果 想在流動(dòng)模型中捕獲分米到米量級(jí)的多相性，則需要井眼圖像來進(jìn)行識(shí)別。
在井眼圖像測井圖中帶有特征標(biāo)記的區(qū)域(例如孔洞，電阻和導(dǎo)電塊)在這被 作為巖石物理相。其它作者，例如Leduc等(2002)和Mathis等(2003)將此結(jié)構(gòu)區(qū)域稱 作測井相(或電相)ο 參見 Leduc, J. P.，Delhaye-Prat, V.，Zaugg, P.并見 Mathis，B.， 2002，F(xiàn)MI*based sedimentaryfacies modelling，Surmont Lease (Athabasca，Canada) (abs.) :CSPG AnnualConvention, Calgary, Alberta，Canada，IOp.；以及參見 Mathis, B.，Leduc, J. P.禾口 Vandenabeele，Τ.，2003，F(xiàn)rom the geologists' eyes to synthetic coredescripti ons :Geological log modeling using well-log data (abs.) AAPGAnnual Meeting，Salt Lake City，UT，7p。由諸如黑色、棕色、和白色(圖4)的不同顏色表征的紋理可用作定義巖石物理相。 這種相具有復(fù)雜的三維形狀。如果是增強(qiáng)的孔隙度和滲透率區(qū)，則導(dǎo)電塊提供孔洞之間的 流動(dòng)連續(xù)區(qū)。公開的文獻(xiàn)具有使用以下技術(shù)構(gòu)建的多個(gè)數(shù)字巖石(或巖石和孔隙的數(shù)子巖石 模型)的示例，所述技術(shù)包括由二維薄截面或掃描電子顯微鏡(SEM)圖像、從測井圖解釋的 測井相、計(jì)算機(jī)生成的球形單元、激光掃描共聚焦顯微鏡、以及各種類型的CT掃描(傳統(tǒng) 的，微型CT，和同步計(jì)算的微層析成像)進(jìn)行的重構(gòu)。Bakke 和 Oren (1997)，Oren 等人(1998)，以及 Oren 和 Bakke (2002)發(fā)展了一 種由二維薄截面構(gòu)建三維孔隙網(wǎng)模型的技術(shù)。數(shù)值巖石(httpV/www.numericalrocks. com/)由二維薄截面計(jì)算三維孔隙模型。參見Bakke，S.和Oren，P. -Ε.，1997，3-D pore-scale modeling of sandstones andflow simulations in the pore networks SPE preprint 35479, European 3-DReservoir Modeling Conference, Stavanger, Norway, April 16—17，p. 136—149 ；Oren, P.-Ε.，Bakke, S.禾口 Arntzen，0.J.，1998， Extending predictivecapabilities to network models :SPE Journal, v. 3, p.324 ； 以及 Oren, P. -Ε.，禾口 Bakke，S. , 2002, Process based reconstruction of sandstones and predictionof transport properties Transport in Porous Media, v. 46, p. 311-343.此公司也使用由微CT掃描構(gòu)建的孔隙模型。Bakke等人(2002)將此技 術(shù)成功地應(yīng)用于砂粒狀白云巖。Duey (2008)和Suicmez以及Touai (2008)的文章總 結(jié)了由數(shù)值巖石處理的各種沙巖孔隙網(wǎng)的結(jié)果。參見Duey，R. ,2008, Quick analysis answers Heidrun question :Hart Energy Publishing， LP， accessed online athttp:// www.eandp.info/index2.php ？ area = article&articleld = 767， March 27，2008; 以及 Suicmez, V. S.，and Touai, Μ.，2008，Pore network modeling :Anew technology for SCAL predictions and interpretations Saudi Arabia Oiland Gas，Issue 5， p. 64-70。Wu等人(2006)介紹了一種使用三階馬爾可夫網(wǎng)格由二維薄截面生成三維數(shù) 值巖石模型的方法。參見 Wu，K.，Van Dijke，Μ. I. J.，Couples, G. D.，Jiang, Ζ.，Ma, J.， Sorbie, K.S., Crawford, J. , Young, L 禾口 Zhang, Χ. ,2006,3D stochastic modeling of heterogeneous porousmedia__Applications to reservoir rocks -Transport in Porous Media, v. 65，p. 443-467。Awwiller (2007)研發(fā)了一 種模擬比 Oren 和 Bakke (2002)描 述的砂巖的更復(fù)雜的砂巖的技術(shù)。Awwiller(2007)的專利申請US2007/0203677A1 (下 文)涉及此方面。Okabe和Blimt (2004，2005)利用多點(diǎn)模擬統(tǒng)計(jì)由二維薄截面得到三維 圖像。參見 0kabe，H.禾口 Blunt，Μ· J. ，2004，Prediction of permeability for porousmedia reconstructed usingmultiple—point statistics :Physical Review E, v. 70, p. 066135-1-10 ；以及參見Okabe, H.和 Blurt, Μ. J.，2005，Pore space reconstruction using multiple-pointstatisticsJournal of Petroleum Science and Engineering, v. 46，p. 121-137。Tomutsa和Radmilovic (2003)利用稀薄的離子束產(chǎn)生用于構(gòu)建亞微級(jí)孔 隙的三維模型的多個(gè)二維連續(xù)截面。參見Tomutsa，L.和Radmilovic，V.，2003，F(xiàn)ocused ion beam assisted three-dimensional rock imaging atsubmicron scale International Symposium of the Soc. Of Core Analysts, Pau, France, September 21—24, Paper SCA2003-47。Dvorkin等人(2003)介紹了物理數(shù)字巖石技術(shù)，其包括來由以下得出的孔隙級(jí) 數(shù)字模擬(a) 二維薄截面和統(tǒng)計(jì)指標(biāo)模擬，或(b)CT掃描。參見Dvorkin，J.，Kameda, Α. , Nur, Α. , Mese, Α.禾口 Tutuncu, Α. N. ,2003, Realtime monitoring of permeability, elastic moduli and strenth in sands andshales using Digital Rock Physics :SPE preprint 82246, presented at the SPEEuropean Formation Damage Conference, The Hague, Netherlands, Mayl3_14，7p。Dvorkin等人構(gòu)建了虛擬巖石的三維模型，并且使用 lattice-Boltzmarm方法進(jìn)行了流動(dòng)模擬。美國專利6516080 (下文)涉及此方面內(nèi)容。Leduc等人(2003)和Mathis等人(2003)(上文均提到)提出了由有限數(shù)量的 所述巖心、傳統(tǒng)的裸眼井測井圖以及井眼圖像測井圖生成“人造巖心”。在傳統(tǒng)的裸眼井測 井圖上使用點(diǎn)群分析。使用傳統(tǒng)的裸眼井測井圖和井眼圖像的結(jié)構(gòu)分析定義基于類似巖 性的深度間隔測井的“測井相”。使用“列聯(lián)表”在沒有巖心的井中計(jì)算假巖心。美國專利 6011557 (下文)與此內(nèi)容相關(guān)。Vahrenkamp等人(2008)描述了微型模型，即，尺寸小于1. Om3并在儲(chǔ)層尺度模 型中提供偽特性給體積單元的儲(chǔ)層模型。參見Vahrenkamp，V. C.，Creusen, A.，Tull, S.， Farmer,Α. ,Mookerjee,Α.禾口 Al Bahry,A. ,2008,Multi-scale heterogeneity modelling in a giant carbonate field, northern Oman(abs.) :GeoArabia, v. 13, No. 1, p. 248。微 型模型由于“基本巖石類型”(PRT)而增加，所述基本巖石類型“涵蓋以及分類整個(gè)范圍的 孔隙類型、尺寸、孔喉尺寸分布、毛細(xì)管入口壓力，相對(duì)滲透率等”?！癙RT”被組織成為“巖石 類型結(jié)合”(PTA)，其基于由井眼圖像測井圖確定的“沉積結(jié)構(gòu)”。利用井眼圖像測井圖在儲(chǔ) 層中分配RTA，并且觀測層模型、相模型、以及地震數(shù)據(jù)。Bryant等人(1993)和Behseresht等人(2007)說明了計(jì)算機(jī)生成的球形密度隨 機(jī)周期單元的數(shù)字巖石模型。參見Bryant，S.，Mellor，D.和Cade, C.，1993，Physically representative network models of transport in porousmedia American Institute of Chemical Engineers Journal,v. 39,No. 3,p. 387-396 ^^jALBehseresht, J. ,Bryant, S. L.禾口 Sepehrnoori, K. ,2007,Infinite-acting physically representative networks for capillarity-controlleddisplacements :SPE preprint 110581, presented at the SPE Annual TechnicalConference and Exhibition, Anaheim, CA, November 11—14,15p。 諸如Bosl等人(1998)和Holt(2001)的其它研究者對(duì)于流動(dòng)實(shí)驗(yàn)已經(jīng)生成了類似的數(shù) 字巖石模型。參見 Bosl，W. J, Dvorkin, J.和 Nur，Α.，1998，A study ofporosity and permeability using a lattice Boltzmann simulation :GeophysicalResearch Letters, v. 25，p. 1475-1478 ；并參見 Holt，R.M.，2001，Particle vs. laboratory modelling inin situ compaction :Physics and Chemistry of theEarth，Part A Solid Earth and Geodesy, v. 26，Issue 1—2，p.89—93。Fredrich等人(1995)和Fredrich (1999)提出應(yīng)用激光掃描共聚焦顯微鏡生成三維巖石圖像。參見 Fredrich, J. T.，Menendez, B.禾口 Wong，Τ. F.，1995，Imaging the pore structure of geomaterials Science, v. 268, p. 276-279 ；及參見Fredrich，J· Τ· , 1999, 3D imaging of porous media usinglaser scanning confocal microscop with application to microscale transportprocesses :Physics and Chemistry of the Earth， Part A Solid Earth andGeodesy, v. 24, Issue 7，p. 551-561。Connor 禾口 Fredrich (1999)運(yùn)用 Iattice-Boltzmann方法對(duì)這些數(shù)值巖石進(jìn)行流動(dòng)實(shí)驗(yàn)。參見0’Connor，R. Μ.和Fredrich， J. Τ. ,1999,Microscale flow modeling in geologic materials :Physics and Chemistry of the Earth，Part A :Solid Earth and Geodesy, v. 24,Issue 7，p.611—616。生成孔隙網(wǎng)的最普遍的方法是利用多種類型的CT掃描。Vinegar (1986)， Wellington和Vinegar (1987)，以及Withjack等人(2003)總結(jié)了這些技術(shù)并對(duì)計(jì)算的X射 線層析成像的多種應(yīng)用進(jìn)行了討論。參見Vinegar，H. J.，1986，X-ray CT md NMR imaging of rocks :JPT, p.257-259 ； # B Wellington S.L. ， and Vinegar, H. J. ，1987， X-ray computer izedtomography :JPT，p. 885-898 ；及參見 With jack，Ε. M.，Devier，C.禾口 Michael， G·，2003， The role of X-ray computed tomography in core analysis :SPEpreprint 83467， presented at the Western Region/AAPG Pacific Section JointMeeting, Long Beach, CA, May 19-24，2003，12p。Siddiqui 和 Khamees (2005)，和 Siddiqui 等人(2005)強(qiáng) 調(diào)了巖心三維圖像以及來自傳統(tǒng)和微CT掃描的片段(cutting)的應(yīng)用。參見Siddiqui， S.禾口Khamees，A. A.，2005，Data visualization challenges for displaying laboratory core and flow data inthree—dimensions :SPE preprint 106334， presented at the SPE TechnicalSymposium of Saudi Arabia，May 14-16,9p.；及參見 Siddiqui，S.禾口 Khamees, Α. Α.，2005，Data visualization challenges for displaying laboratory core andflow data in three-dimensions :SPE preprint 106334， presented at the SPETechnical Symposium of Saudi ArabiaiMay 14_16，9p0 Coles 等人(1996)、Fredrich 等人(2006)、以及Fredrich等人(2007)應(yīng)用同步計(jì)算的微層析成像在天然和合成沙巖中 構(gòu)建孔隙網(wǎng)的數(shù)值三維模型。參見 Coles，M.E.，Hazlett，R.D. ,MueggeiR. L.，JonesX W.， Andrews, B. Dowd, B. Siddons, P.，Peskin，A.，Spanne, P.禾口 Soll，W. Ε.，1996，Developments in synchrotronX-ray microtomography with applications to flow in porous media SPEpreprint 36531，presented at the SPE Annual Technical Conference andExhibition，Denver, C0，p. 413—424 ；參見 Fredrich，J. T. ，DiGiovanni，A. A.禾口 Noble， D.R，，2006， Predicting macroscopic transport properties usingmicroscopic image data Journal of Geophysical Research B Solid Earth，v.111，Issue 3 ；及參見 Fredrich,J. Τ. ,Haney,M. M.禾口 White，J，A.，2007，Predicting petrophysical properties using 3D image data (abs.) :MPG AnnualConvention, downloaded at http://www.aapg. org。他們使用Iattice-Boltzmann方法以模型滲透率。多點(diǎn)(或多重點(diǎn))統(tǒng)計(jì)方法(MPS)是一門在二十世紀(jì)90年代提出的用于生 成諸如地質(zhì)相的離散變化區(qū)域的條件模擬的空間統(tǒng)計(jì)插入算法方法。參見Guardimo，F(xiàn).禾口 Srivastava，R. Μ.，1993，Multivariate geostatistics :beyond bivariate moments :Geostatistics_Troia，A.Soares. Dordrecht, Netherlands，Kluwer Academic Publications, v. 1，p. 133-144。訓(xùn)練圖像是數(shù)值先驗(yàn)地質(zhì)模型，所述地質(zhì)模型包括相結(jié)構(gòu) 和被認(rèn)為存在于現(xiàn)實(shí)儲(chǔ)層中的關(guān)系。訓(xùn)練圖像本身是一種概念并且可以和手繪圖一樣簡 單，或者所述訓(xùn)練圖像可以由計(jì)算機(jī)工具生成。由Guardiano和Srivastava(1993)提出 的原始MPS算法通過對(duì)每個(gè)模擬點(diǎn)進(jìn)行再次掃描訓(xùn)練圖像來構(gòu)建多點(diǎn)條件概率分布函數(shù) (CPDF)。由于計(jì)算機(jī)中央處理單元(CPU)的限制，這種消耗時(shí)間的算法在當(dāng)時(shí)是不實(shí)用的。 Strebelle (2002)介紹了搜索樹的概念，所述搜索樹存儲(chǔ)訓(xùn)練圖像上的模板中得 到的圖案的所有復(fù)制。參見 Strebelle，S. ,2002, Conditionalsimulation of complex geological structures using multiple point statistics. Mathematical Geology, v. 34, p. 1-22。Strebelle (2002)的還被稱為SNESIM基本MPS算法已經(jīng)在用于儲(chǔ)層模 擬的多種應(yīng)用中使用，并且當(dāng)結(jié)合旋轉(zhuǎn)以及相似變換時(shí)已經(jīng)成為用于模擬河道沉積的參 考工具(Zhang，2002 ；Caers 和 Zhang, 2004 ；Strebelle 和 Zhang, 2004)。參見 Zhang， Τ.,2002, Multiple-point simulation of multiple reservoir facies :Unpublished M. S. thesis, Stanford University, California, 163p. ;#JAL Caers, J.禾口 Zhang，T.，2004， Multiple-point geostatistics :A quantitative vehicle for integration ofgeologic analogs into multiple reservoir models, in M. Grammer, P.M.Harrisand G.P. Eberli, eds. Integration of Outcrop and Modern Analogs inReservoir Modeling, AAPG. Memoir 80，p. 383-394.；及參見 Strebelle, S.禾口 Zhang，Τ. ，2004，Non-stationary multiple-point geostatistical models, inLeuangthong, 0. and Deutsch, C. V.，eds. Geostatistics, v. 1,p.235—244?；綧PS算法的速度遠(yuǎn)遠(yuǎn)快于Guardiano和Srivastava(1993)建立的原始算 法，但是尤其是對(duì)于較大的訓(xùn)練圖像進(jìn)行三維模擬時(shí)對(duì)隨機(jī)存取存儲(chǔ)器(RAM)有要求。 參見 Guardiano, F.禾口 Srivastava, R. Μ. ,1993, Multivariate geostatistics :beyond bivariate moments :Geostatistics_Troia, A.Soares.Dordrecht, Netherlands, Kluwer Academic Publications, v. 1，p. 133-144。這種對(duì)三維的RAM局限性需要做出可能導(dǎo)致三 維物體不恰當(dāng)?shù)男螤钤佻F(xiàn)的妥協(xié)。RAM的局限性還防止共同考慮太多的分類或類別，因此限 制了基本MPS算法對(duì)分類變量的模擬。為了處理分類及連續(xù)變化的訓(xùn)練圖像并減少RAM成本和改進(jìn)三維應(yīng)用時(shí)形狀再 現(xiàn)，提出 了諸如 FILTERSIM 的MPS算法(Zhang 2006a)。參見Zhang，Τ.，2006a，F(xiàn)ilter-based training image pattern classification for spatialpattern simulation Unpublished Ph. D. dissertation, Stanford University, California,153p。FILTERSIM 算法將一組局部過濾器施加到可以是分類或連續(xù)的訓(xùn)練圖像，以將局部圖案分組成圖案類 另O。然后，F(xiàn)ILTERSIM算法對(duì)基于所述分類的圖案進(jìn)行模擬。過濾器是具有與模板的每一 個(gè)像素位置相關(guān)聯(lián)的權(quán)值的集合的局部模板(窗口)。將過濾器施加到局部圖案產(chǎn)生過濾 器得分，所述得分被視作為所述局部圖案的數(shù)值總和。一組默認(rèn)或用戶定義的過濾器被設(shè) 計(jì)成使得每一個(gè)過濾器可以記錄在模板內(nèi)所看到的訓(xùn)練圖案的不同方面。這些過濾器用于 將訓(xùn)練圖案轉(zhuǎn)換成過濾器得分空間(score space)。這種圖案得分使圖案的尺寸顯著減小。 通過分隔所述有限尺寸的得分空間，根據(jù)訓(xùn)練圖案的過濾器得分對(duì)類似的訓(xùn)練圖案進(jìn)行分類。基本MPS算法的速度遠(yuǎn)遠(yuǎn)快于Guardiano和Srivastava (1993)建立的原始算法，但是尤其是對(duì)于較大的訓(xùn)練圖像進(jìn)行三維模擬時(shí)對(duì)計(jì)算機(jī)隨機(jī)存取存儲(chǔ)器(RAM)有要求。 這種對(duì)三維的RAM局限性需要做出可能導(dǎo)致三維物體不恰當(dāng)?shù)男螤钤佻F(xiàn)的妥協(xié)。RAM的局 限性還防止共同考慮太多的分類或類別，因此限制了基本MPS算法對(duì)分類變量的模擬?；?本MPS算法搜索條件數(shù)據(jù)事件的精確復(fù)制，每次為儲(chǔ)層模型構(gòu)建一個(gè)像素，以多點(diǎn)數(shù)據(jù)事 件為條件，并且不允許對(duì)在訓(xùn)練圖像中得到的圖案進(jìn)行任何過濾或平均。為了處理分類及連續(xù)變化的訓(xùn)練圖像并減少RAM成本和改進(jìn)三維應(yīng)用時(shí)形狀再現(xiàn)，Zhang(2006a)提出了一種新的稱作為FILTERSIM的MPS算法。FILTERSIM算法將一組 局部過濾器施加到可以是分類或連續(xù)的訓(xùn)練圖像，以將局部圖案分組成圖案類別。然后， FILTERSIM算法對(duì)基于所述分類的圖案進(jìn)行模擬。過濾器是具有與模板的每一個(gè)像素位置 相關(guān)聯(lián)的權(quán)值的集合的局部模板(窗口)。將過濾器施加到局部圖案產(chǎn)生過濾器得分，所述 得分被視作為所述局部圖案的數(shù)值總和。一組默認(rèn)或用戶定義的過濾器被設(shè)計(jì)成使得每一 個(gè)過濾器可以記錄在模板內(nèi)所看到的訓(xùn)練圖案的不同方面。這些過濾器用于將訓(xùn)練圖案轉(zhuǎn) 換成過濾器得分空間(score space)。這種圖案得分使圖案的尺寸顯著減小。通過分隔所 述有限尺寸的得分空間，根據(jù)訓(xùn)練圖案的過濾器得分對(duì)類似的訓(xùn)練圖案進(jìn)行分類。FILTERSIM算法開始于減小尺寸的過濾器得分空間內(nèi)的局部訓(xùn)練圖案分類。通過 確定哪一個(gè)圖案類別是最類似于局部條件數(shù)據(jù)事件、從圖案類別對(duì)具體圖案進(jìn)行采樣、然 后將采樣的圖案修補(bǔ)到模擬位置處的圖像上，模擬沿著通過模擬空間的序貫路徑進(jìn)行。模 擬隨機(jī)路徑和來自圖案類別的圖案取樣允許不同的模擬再現(xiàn)，然而所有都以相同的原始數(shù) 據(jù)為條件。由于任何圖案中的過濾器總和使尺寸減少，并且因?yàn)閳D案被分組到類別中，因此 算法根據(jù)RAM的請求而較快并且合理?；綧PS算法和FILTERSIM算法能夠遵守(honour)來自井或露頭內(nèi)獲得的數(shù)據(jù) 的絕對(duì)或所謂的“硬”約束，和來自地震數(shù)據(jù)、相概率場、以及旋轉(zhuǎn)和親力(affinity)(或 尺度)限制網(wǎng)格的條件或“軟”約束。所有這些數(shù)據(jù)在隨機(jī)模擬處理中使用，以生成地質(zhì) 相或巖石特性的1D、2D或3D圖像。因?yàn)樵谒婕暗腗PS模擬中具有隨機(jī)分量，所以由MPS 算法生成的特性場的單獨(dú)再現(xiàn)不同，但是整體再現(xiàn)為地學(xué)科學(xué)家和油藏工程師提供所模 擬的儲(chǔ)層體積內(nèi)的地質(zhì)相的空間分布和不確定性的定量估計(jì)。此外，這些算法遵守硬輸 入數(shù)據(jù)和軟輸入數(shù)據(jù)約束(Zhang, 2006a)參見 Zhang，Τ.，Switzer P.，and Journel Α.， 2006b,Filter-based classification oftraining image patterns for spatial pattern simulation !Mathematical Geology, v. 38,p.63—80。通常在FILTERSIM算法中使用6個(gè)定向2D默認(rèn)過濾器(Zhang，2006a ；Zhang等， 2006b)(這兩篇文章在以上都進(jìn)行了說明)。具有三種類型的過濾器平均過濾器、梯度過 濾器和曲率過濾器，并且每一種類型的過濾器都用于水平和垂直方向。平均過濾器的目的 在于局部化特征；梯度過濾器用于通過使不同特征(一階差)的對(duì)比顯著而檢測特征邊界； 曲率過濾器提供特征的二階差。為了反映大尺度結(jié)構(gòu)，使用多網(wǎng)格模擬，所述多網(wǎng)格模擬利用細(xì)網(wǎng)格模擬從粗網(wǎng) 格到細(xì)網(wǎng)格逐漸模擬多個(gè)網(wǎng)格的每一個(gè)水平，所述細(xì)網(wǎng)格由在粗網(wǎng)格處先前模擬的值約 束。在模擬的每一個(gè)水平處，重定尺度的過濾器施加在相應(yīng)的網(wǎng)格上(Zhang，2006a)。
具有兩種類型的訓(xùn)練圖像一種具有非常有限數(shù)量的分類，而另一種用于諸如 油藏巖石物理性質(zhì)的連續(xù)變量。多點(diǎn)地質(zhì)統(tǒng)計(jì)方法要求訓(xùn)練圖像的1D、2D或3D網(wǎng)格作 為在先驗(yàn)概念地質(zhì)模型，所述先驗(yàn)概念地質(zhì)模型包括研究中的空間屬性的圖案。在圖像 上顯現(xiàn)的不同特征的形狀被認(rèn)為是表征實(shí)際地質(zhì)特征的模型，且每一種分類通常表征不 同的地質(zhì)相或不同種類的地質(zhì)體。訓(xùn)練圖像通常要求含有“平穩(wěn)”圖案，即，所述圖案必 須不依賴于所述圖案在空間內(nèi)的位置(不根據(jù)任何平移而變化)并且必須在訓(xùn)練圖像 區(qū)域上重復(fù)。在訓(xùn)練圖像用于地質(zhì)模型的情況下，這種平穩(wěn)性可以包括但不限于地質(zhì)對(duì) 象定向平穩(wěn)性(其中，方向?qū)ο?特征在圖像上不旋轉(zhuǎn))和地質(zhì)尺度平穩(wěn)性(其中，圖 像上的對(duì)象/特征的尺寸在圖像上不改變)。(Caers, J.和Zhang，2004)。參見Caers， J. and Zhang, Τ. ,2004, Multiple-point geostatistics :A quantitative vehicle for integration of geologicanalogs into multiple reservoir models, in M. Grammer, P. M Harris and G. P. Eberli, eds. integration of Outcrop and Modern Analogs in ReservoirModeling, AAPG. Memoir 80, p.383—394。當(dāng)前的MPS算法隱含的問題是如何生成訓(xùn)練圖像。訓(xùn)練圖像被認(rèn)為是模擬或再現(xiàn)實(shí)際地質(zhì)特征，并且應(yīng)該盡可能地由現(xiàn)有的在地質(zhì)方面有意義的圖像得出。這需要對(duì)將 允許使用來自任何源(例如，手繪草圖、航空照片、衛(wèi)星圖象、地震體積、地質(zhì)對(duì)象模型、物 理尺度模型、或正演地質(zhì)處理模型)的圖像的統(tǒng)計(jì)和圖像處理方法進(jìn)行研究。與連續(xù)變化 的訓(xùn)練圖像的生成相比較，生成分類變化的訓(xùn)練圖像比較容易。基于對(duì)象的方法通常用 于生成具有分類變量的訓(xùn)練圖像。與增加的期望約束結(jié)合的基于區(qū)域的方法可以用于生 成連續(xù)變化的訓(xùn)練圖像(Zhang 等人，2005)。參見 Zhang, L.，Nair, N.，Jennings, J. W.， and Bryant, S. L. , 2005, Models andmethods for determining transport properties of touching-vug carbonates SPEpreprint 96027, presented at the SPE Annual Technical Conference andExhibition, Dallas, TX, October 9—12，9p。為了執(zhí)行MPS的模擬，訓(xùn)練圖像必須平穩(wěn)。然而，在大多數(shù)儲(chǔ)層模擬應(yīng)用中，地質(zhì) 沉積顯現(xiàn)出非平穩(wěn)的形狀/特征，所述非平穩(wěn)的形狀/特征反映儲(chǔ)層多相性和沉積的各相 異性。

發(fā)明內(nèi)容
本發(fā)明涉及一種用于生成被稱為數(shù)值假巖心的巖石和孔隙的三維(3D)計(jì)算機(jī)模 型的方法。該技術(shù)使用全井眼圖像、數(shù)字巖石樣品、和多點(diǎn)統(tǒng)計(jì)(MPS)以對(duì)其中實(shí)際沒有采 集巖心但是具有已測錄井眼圖像的測井間隔重現(xiàn)3D假巖心。數(shù)字巖心樣品通過多點(diǎn)地質(zhì) 統(tǒng)計(jì)用于引導(dǎo)假巖心的3D特征再現(xiàn)，并且最終的假巖心由全井眼圖像限制。當(dāng)結(jié)合參照附圖時(shí)，本發(fā)明的進(jìn)一步特征和優(yōu)點(diǎn)將從以下詳細(xì)說明變得顯而易 見。


本申請文件包括至少一副彩圖。具有彩圖的本專利或?qū)＠暾埑霭嫖锏母北驹谡?求和支付必要費(fèi)用時(shí)將由辦事處來提供。以下以本發(fā)明的示例性實(shí)施例的非限制性示例的方式參照所示的多個(gè)附圖在以下詳細(xì)說明中進(jìn)一步描述本發(fā)明，其中，相同的附圖標(biāo)記在附圖的幾幅視圖中表示相似的 部件，其中圖Ia示出了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的被定義為在孔洞碳酸鹽巖地層中具 體地為單程地層微成像儀(FMI)的井眼圖像測井圖的3英尺(或1米)的間隔的訓(xùn)練圖 像；圖Ib圖示了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的使用與圖Ia相同的訓(xùn)練圖像的 FILTERISM的多點(diǎn)模擬算法(MPS)的結(jié)果，其依據(jù)本發(fā)明至少一個(gè)實(shí)施例；圖2a示出了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的與圖Ia相同的訓(xùn)練圖像；圖2b圖示了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的圖Ia的訓(xùn)練圖像的第一次再現(xiàn) (realization)的結(jié)構(gòu)，并與圖1中的孔洞地層的原始測量數(shù)據(jù)相匹配，使得所測量的數(shù)據(jù) 在第一再現(xiàn)(圖2a)、第二再現(xiàn)(圖2b)、以及第三再現(xiàn)(圖2c)被遵守，所述第一再現(xiàn)、第 二再現(xiàn)以及第三再現(xiàn)提供了被模擬區(qū)域的細(xì)微變化；圖2c圖示了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的圖Ia的訓(xùn)練圖像的第二再現(xiàn)的結(jié) 果，并與圖Ia中的孔洞地層的原始測量數(shù)據(jù)相匹配；圖2d圖示了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的圖Ia的訓(xùn)練圖像的第三再現(xiàn)的結(jié) 果，并與圖Ia中的孔洞地層的原始測量數(shù)據(jù)相匹配；圖3a示出了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的孔洞碳酸鹽巖的隨鉆測井(LWD) 圖像中顯示為黑色特征，使得導(dǎo)電塊以棕色示出，而非導(dǎo)電塊以白色示出的孔洞的靜態(tài)圖 像；圖3b示出了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的孔洞碳酸鹽巖的隨鉆測井(LWD) 圖像中顯示為黑色特征，使得導(dǎo)電塊以棕色示出，而非導(dǎo)電塊以白色示出的孔洞的動(dòng)態(tài)圖 像；圖4示出了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的顯示出畫出電圖像的較低電阻率區(qū) 域的等高線(綠線)的圖Ia的孔洞孔隙度的全井眼圖像；圖5a和5b示出了根據(jù)本發(fā)明的至少一個(gè)實(shí)施例的使用井眼圖像、數(shù)字巖石樣品 和多點(diǎn)統(tǒng)計(jì)生成數(shù)值假巖心的方法的示意性方框圖；圖6示出了根據(jù)本發(fā)明的實(shí)施例的在透明、板狀孔洞碳酸鹽巖樣品(約6英寸高 和直徑為4英寸)的CT掃描的三維視圖中圖Ia的訓(xùn)練圖像，使得所示孔洞為紅色到灰色， 而巖石基質(zhì)被示出為藍(lán)色；圖7示出了根據(jù)本發(fā)明的實(shí)施例的在具體深度處的被彎曲成與井眼直徑(井眼直 徑約為8. 5英寸，并且長度間隔為3英尺或1米)。圖6的CT掃描訓(xùn)練圖像精確地標(biāo)刻并且 位于井眼圖像中心的正確深度處；其中孔洞被示出為紅色，而巖石基質(zhì)被示出為淺藍(lán)色；圖8示出了根據(jù)本發(fā)明的實(shí)施例的由訓(xùn)練圖像和圖7相應(yīng)的全井眼圖像生成的數(shù) 值兩相假巖心的結(jié)果，使得孔洞被示出為紅色，而巖石基質(zhì)被示出為深藍(lán)色；圖9示出了根據(jù)本發(fā)明的實(shí)施例的由訓(xùn)練圖像和圖7相應(yīng)的全井眼圖像生成的數(shù)值假巖心的結(jié)果，使得孔洞不可見而巖石基質(zhì)被示出為灰色。任意懸掛部分(見圖9上部) 表示這是三維模型，且在每個(gè)部分帶有大量孔隙(間隔的長度約1英尺或0. 3米，并且井眼 直徑約為8. 5英寸)；圖10示出了根據(jù)本發(fā)明的實(shí)施例的由訓(xùn)練圖像和圖6相應(yīng)的全井眼圖像生成的數(shù)值假巖心的結(jié)果。導(dǎo)電塊(被示出為紅色)在如圖4所示的導(dǎo)電塊周圍具有類似等高 線。要注意的是示出了三個(gè)巖石物理相孔隙或孔洞(被示出為綠色)、導(dǎo)電塊(被示出為 紅色)、以及巖石基質(zhì)(被示出為藍(lán)色)；圖11示出了圖10的具有不可見巖石基質(zhì)的數(shù)值的假巖心，所述不可見巖石基質(zhì) 強(qiáng)調(diào)孔隙或孔洞(被示出為綠色)與導(dǎo)電塊(被示出為紅色)的三維互聯(lián)性質(zhì)；圖12示出了根據(jù)本發(fā)明的實(shí)施例的如圖10所示的數(shù)值假巖心，以及孔隙或孔洞 (被示出為綠色)的復(fù)雜三維形狀，導(dǎo)電塊和巖石基質(zhì)都不可見；圖13示出了根據(jù)本發(fā)明的實(shí)施例的基于在井眼圖像或其它測井圖中觀測的層可 以再取樣或者再分配到其中徑向網(wǎng)格可以被分層的徑向網(wǎng)格的數(shù)值假巖心；
圖14a示出了根據(jù)本發(fā)明的實(shí)施例的在流動(dòng)模擬中用于導(dǎo)電塊的相對(duì)滲透率曲 線；其中具有三條不同的相對(duì)滲透率曲線對(duì)應(yīng)于假巖心模型中諸如孔洞、導(dǎo)電塊和緊密巖 石基質(zhì)的三種巖石類型的每一個(gè)；以及圖14b示出了根據(jù)本發(fā)明的實(shí)施例的通過數(shù)值假巖心的流動(dòng)模擬結(jié)果，使得一連 串水微型注射器圍繞在直徑以外，一連串微型發(fā)生器包圍假巖心的內(nèi)直徑顏色代表油飽 和(So)，其中在此流動(dòng)模型中顯示了多相進(jìn)步和突破。
具體實(shí)施例方式這里示出的細(xì)節(jié)僅以示例的方式并且出于說明本發(fā)明的實(shí)施例的示例性說明，并 在提供什么是最有用并且易于理解本發(fā)明的原理和概念方面的說明時(shí)給出。在這點(diǎn)上，并 未試圖以比基本理解本發(fā)明所需更詳盡的方式顯示本發(fā)明的結(jié)構(gòu)細(xì)節(jié)，且結(jié)合附圖的說明 使本領(lǐng)域技術(shù)人員理解本發(fā)明是如何在實(shí)際中實(shí)施的幾個(gè)形式。此外，在各個(gè)圖中涉及的 相同的附圖標(biāo)記及名稱表示相同元件。此外，以下優(yōu)選實(shí)施例的細(xì)節(jié)描述，參考形成所述實(shí) 施例的一部分的附圖，對(duì)特定的實(shí)施例做出了使本發(fā)明更容易理解的圖解。應(yīng)了解其它應(yīng) 用的實(shí)施方式和結(jié)構(gòu)的變化并不脫離本發(fā)明的保護(hù)范圍。本發(fā)明涉及生成被稱作為數(shù)值假巖心的巖石和孔隙的三維(3D)計(jì)算機(jī)模型的方 法。該技術(shù)應(yīng)用全井眼圖像、數(shù)字巖石樣品，和多點(diǎn)模擬統(tǒng)計(jì)(MPS)以對(duì)于沒有實(shí)際采集巖 心但是具有井眼圖像的被測錄間隔重構(gòu)三維假巖心。數(shù)字巖心樣品通過多點(diǎn)地質(zhì)統(tǒng)用于引 導(dǎo)假巖心的三維特性再現(xiàn)，并且最終的假巖心由全井眼圖像限制。根據(jù)本發(fā)明的一個(gè)方面，基本思想在于使用訓(xùn)練圖像我們直接使用數(shù)據(jù)本身 [在全井眼生成期間，所述訓(xùn)練圖像是具有大于整個(gè)區(qū)域的覆蓋率60%的原始2D不完全圖 像(連續(xù)變量訓(xùn)練圖像)，而在假象重構(gòu)中，計(jì)算機(jī)輔助測試掃描數(shù)字巖芯直接用作隨后與 全井眼圖像數(shù)據(jù)結(jié)合的3D訓(xùn)練圖像]。因此，應(yīng)用MPS的整個(gè)過程變成為由數(shù)據(jù)驅(qū)動(dòng)。應(yīng) 該在我們的專利備忘錄和條款中強(qiáng)調(diào)此優(yōu)點(diǎn)。井眼圖像給出了與井眼碰撞的巖石和流體的定向電和聲波圖像。全井眼圖像是井 壁360度全視圖。通過設(shè)計(jì)，大多數(shù)聲波測井儀和隨鉆測井儀生成全井眼圖像。因?yàn)樵跇O板 間存在間隙，而且可能損壞一些電極或極板，因此大多數(shù)電阻率成像測井儀需要模擬全井 眼圖像。建模過程應(yīng)用MPS中提出的連續(xù)變量算法。全井眼圖像給出了特征的物理位置， 例如在井壁上的巖石內(nèi)的地層界面、孔隙、和導(dǎo)電和電阻塊。這些概括了公知為巖石物理相 的復(fù)雜三維體積。
數(shù)字巖石或巖心樣品可由二維薄截面、掃描電子顯微鏡圖像、共聚焦顯微鏡圖像、 或者計(jì)算機(jī)生成的球形單元生成。然而，大多數(shù)數(shù)字巖石或巖心樣品由對(duì)巖石和流體的計(jì) 算的層析成像掃描(CT掃描)生成。當(dāng)掃描器以不同的角度將X射線通過巖石發(fā)射到接收 器時(shí)，層析成像圖產(chǎn)生。X射線衰減轉(zhuǎn)化為密度，并且連續(xù)的截面構(gòu)建成巖石和孔隙的三維 視圖。根據(jù)本發(fā)明一個(gè)方面，諸如CT掃描的每個(gè)數(shù)字巖石樣品直接作為訓(xùn)練圖像。具有離 散變量的訓(xùn)練圖像，且在圖像的每一個(gè)像素處屬性為巖石(白色)或孔隙(黑色)。訓(xùn)練圖 像可以具有任意形狀的邊界或者包括多個(gè)不規(guī)則孔。數(shù)值假巖心利用MPS中的離散變量算法產(chǎn)生。整數(shù)值被分配給每個(gè)巖石物理相，例如巖石(0)、孔隙(1)、以及導(dǎo)電塊(2)。數(shù)字巖石或巖心樣品用作訓(xùn)練圖像，即，定量模板 用作模型特性區(qū)。全井眼圖像利用圓柱形包絡(luò)包圍數(shù)值假巖心，以滿足模型的條件。每個(gè) 數(shù)值假巖心絕對(duì)遵守?cái)?shù)字巖石或巖心樣品和全井眼圖像。數(shù)值假巖心可被網(wǎng)格化為適于流 體-流動(dòng)模擬的模型。毛細(xì)管壓力和相對(duì)滲透率曲線由概念模型、專門的巖心分析、或精密 尺度孔隙網(wǎng)模型建立的技術(shù)提供。圖Ia是被定義為在孔洞碳酸鹽巖中的單程地層微成像儀(FMI)的井眼圖像測井 圖的3英尺(1米)間隔的訓(xùn)練圖像?？锥词翘畛溆兴@井泥漿的孔隙，并且所述孔洞被 示出為黑色、低電阻率點(diǎn)。注意圖Ia極板間的間隙。如圖所示，基于像素的、自定義的3X3 的模板(底部)移動(dòng)通過訓(xùn)練圖像，從而檢測圖案并且對(duì)每個(gè)測量像素周圍的相鄰區(qū)域給 出過濾得分。這為利用FILTERSIM的MPS模擬提供了基礎(chǔ)。圖Ib中的全井眼圖像給出了應(yīng) 用圖Ia中的整個(gè)圖像作為訓(xùn)練圖像的FILTERSIM的再現(xiàn)。再現(xiàn)被條件化，使得所述再現(xiàn)與 原始測量數(shù)據(jù)完美匹配。沿著圖像的頂部給出了方位。縮寫N=北；E=東；S=南；以及W =西。垂直比例沒有放大。鉆頭的大小為8. 5in(21. 5cm)。具體地，圖Ib給出了利用圖Ia 的相同訓(xùn)練圖像的FILTERSIM的多點(diǎn)模擬統(tǒng)計(jì)(MPS)結(jié)果。因此，通過應(yīng)用FILTERSIM，我 們利用合適的基于像素的模板為每個(gè)訓(xùn)練圖像確定過濾器得分。如在測量數(shù)據(jù)中所觀察到 的，這些分值量化了圖案和所述圖案的概率，然后利用多點(diǎn)模擬統(tǒng)計(jì)(圖la)的FILTERSIM 算法生成全井眼圖像。圖2b、2c以及2d給出了在孔洞碳酸鹽巖地質(zhì)中的全井眼圖像的三種再現(xiàn)。此處 圖2a給出了與圖Ia相同的訓(xùn)練圖像。要注意的是測量數(shù)據(jù)在每一個(gè)再現(xiàn)中被遵守，并且 在模擬區(qū)域中具有細(xì)微變化。圖2b給出了圖2a訓(xùn)練圖像的第一再現(xiàn)的結(jié)果，并且與圖Ia 的孔洞性地層的原始測量數(shù)據(jù)相匹配。圖2c給出了圖2a訓(xùn)練圖像的第二再現(xiàn)的結(jié)果，并 且與圖Ia的孔洞性地層的原始測量數(shù)據(jù)相匹配。圖2d給出了圖2a訓(xùn)練圖像的第三再現(xiàn) 的結(jié)果，并且與圖Ia的孔洞性地層的原始測量數(shù)據(jù)相匹配。圖3a給出了在孔洞碳酸鹽巖的隨鉆測井(LWD)圖像中被示出為黑色特征的孔洞 的靜態(tài)圖像(參見 Xiao, L.，Jun, C.，Duo, Y. S.，Han, S. Y.，Xia, W. H.和 Xi, W. Y.，2007 Fully integrated solution for LffD resistivity imageapplication a case study from B eibu Gulf,China :1st SPffLA India RegionalConference,Formation Evaluation in Horizontal Wells, Mumbai, Marchl9_20，IOp.)，使得導(dǎo)電塊被示出為棕色，而非導(dǎo)電塊 被示出為白色。圖3b給出了在洞型碳酸鹽巖的隨鉆測井(LWD)圖像中被示出為黑色特征 的孔洞的動(dòng)態(tài)圖像，使得導(dǎo)電塊被示出為棕色，而非導(dǎo)電塊被示出為白色。Dehghani等人 (1999)介紹了增強(qiáng)孔隙度和滲透率的區(qū)域存在于孔洞附近。參見Dehghani，K.，Harris,P. Μ. , Edwards, K. Α.禾口 Dees, W. Τ. ,1999, Modeling a vuggy carbonate reservoir :AAPG Bulletin, v. 83，p.19-42。在圖4中，孔洞孔隙度的全井眼圖像(與圖1和2相同)顯示了勾勒電圖像中的 較低電阻區(qū)的輪廓的等高線(綠線)。沿著圖像的頂部給出了方位?？s寫N=北；E=東； S=南；以及W=西。垂直比例沒有放大。鉆頭的大小為8. 5in (21. 5cm)。圖5示出了應(yīng)用井眼圖像、數(shù)字巖石和多點(diǎn)統(tǒng)計(jì)生成數(shù)值假巖心的方法的示意性 方框圖。步驟100包括從可以用作訓(xùn)練圖像的儲(chǔ)層采集井眼成像測井圖的深度限定間隔。 訓(xùn)練圖像被定向；訓(xùn)練圖像是2維(二維)連續(xù)變化的數(shù)值的標(biāo)量陣列，且在極板之間具有 間隙和/或沒有需要修補(bǔ)的區(qū)域，使得極板表征測量值，而間隙是井眼中的非成像部分。仍然參照圖5中的步驟100，原始全井眼圖像通過諸如聲波裝置和隨鉆測井儀的 一些測井儀生成。在其它情況下，必須生成全井眼圖像。這對(duì)通常在極板之間具有間隙的 電阻率測井尤其正確。還參照圖5中的步驟100，根據(jù)發(fā)明的一個(gè)方面，選擇井眼圖像測井圖的深度限定 間隔，例如，間隔可以為所測深度的lft、3ft、或IOft (0. 3，1，或3m)測量深度。用戶可能想 要選擇厚間隔或薄間隔，這取決于分層、斷裂以及其它非均勻圖案的觀測量。全井眼圖像可利用如圖5的步驟100所示的MPS中的FILTERSIM算法生成。 FILTERSIM利用過濾器得分進(jìn)行分組然后模擬極板間的間隙內(nèi)的圖案，在所述極板間的間 隙中不存在測量數(shù)據(jù)。用于單個(gè)再現(xiàn)的方法將隨機(jī)占據(jù)像素點(diǎn)位置，并且根據(jù)一組過濾器 得分進(jìn)行繪圖以為隨機(jī)點(diǎn)選擇合適的圖案。因?yàn)檫@些都是條件模擬，所以測量的數(shù)據(jù)被遵 守。與極板的邊緣相鄰的圖案與實(shí)際極板上觀察的圖案相匹配。模擬的像素顏色、連續(xù)變 量的頻率分布與測量顏色的頻率分布完全匹配。。圖5中的步驟200包括從儲(chǔ)層中采集至少一個(gè)巖心樣品，然后從采集的巖心樣品 中生成數(shù)字巖心樣品。采集的巖心樣品類似地描繪了井眼和/或儲(chǔ)層中的一個(gè)或多個(gè)深度 限定間隔的代表性特征和/或結(jié)構(gòu)。例如，諸如CT掃描的數(shù)字巖石或巖心樣品可以用作訓(xùn) 練圖像。訓(xùn)練圖像是離散數(shù)值的三維陣列。在二-相模型中，例如，巖石具有數(shù)值0，而孔隙 具有數(shù)值1。在三-相模型中，巖石具有數(shù)值0，孔隙具有數(shù)值1，而導(dǎo)電塊具有數(shù)值2。對(duì) 相的數(shù)量沒有限制。各個(gè)相體(體積)的輪廓可具有任意形狀或尺寸。圖5的步驟300包括對(duì)采集的例如生成全井眼圖像的井眼成像測井?dāng)?shù)據(jù)進(jìn)行預(yù)建 模。圖5的步驟300a公開了對(duì)井眼圖像測井圖的深度限定間隔進(jìn)行選擇的步驟。全 井眼圖像利用多點(diǎn)模擬統(tǒng)計(jì)(MPS)中的FILTERSIM算法生成，其中MPS模擬從訓(xùn)練圖像捕 獲地質(zhì)結(jié)構(gòu)并且將所述地質(zhì)結(jié)構(gòu)固定到數(shù)據(jù)位置。FILTERSIM利用過濾器得分進(jìn)行分組，然后模擬極板間的間隙中的圖案，在所述極板間的間隙中不存在測量數(shù)據(jù)。圖5的步驟300b公開了截成為識(shí)別諸如孔隙或孔洞；在井壁上的巖石的導(dǎo)電和電 阻塊；以及地層界面或巖石基質(zhì)的巖石物理相的步驟。用于單個(gè)再現(xiàn)的方法將隨機(jī)地占據(jù) 像素位置，并且由一組過濾器得分進(jìn)行繪圖以為隨機(jī)點(diǎn)選擇合適的圖案。因?yàn)檫@些都是條件模擬，所以測量的數(shù)據(jù)被遵守。與極板的邊緣相鄰的圖案與實(shí)際極板上觀察的圖案相匹 配。模擬的像素顏色、連續(xù)變量的頻率分布與測量顏色的頻率分布完全匹配。圖5的步驟 300c公開了將全井眼圖像卷曲成帶尺度的圓柱形形狀的步驟。具體地，對(duì)于常規(guī)解釋來說，檢查三維井眼圖像是困難的。因此，通常是沿著真北方向?qū)⒕睬虚_，然后展開圓柱體直到所述井眼圖像變成2D視圖。在高度偏斜的水平井中，通常沿著鉆孔的頂部切開井眼圖像。 與圓柱形井筒相交的平面特征在2D視圖中表現(xiàn)為正弦波。為生成三維數(shù)值假巖心，需要將 二維全井眼圖像卷曲成其原始三維形狀。為了進(jìn)行此，必須已知井眼直徑?？梢杂刹捎迷?始井眼圖像的井陘測井圖來確定所述井眼直徑。圖像的比例必須是1 1，即，垂直比例沒 有放大。圖5的步驟400公開了對(duì)采集的數(shù)字巖心樣品進(jìn)行預(yù)建模的步驟。步驟400a將 由二-相轉(zhuǎn)變成三_相，使得數(shù)字巖石中的孔隙或孔洞、導(dǎo)電塊和巖石基質(zhì)已經(jīng)被識(shí)別。例 如，如果導(dǎo)電塊在數(shù)字巖石(訓(xùn)練圖像)中沒有被識(shí)別，則可以由固定數(shù)量的體素通過孔隙 相的體變(dilation)對(duì)所述導(dǎo)電塊進(jìn)行模擬。這樣，可以在數(shù)值模擬巖心中捕獲大尺度多 相性。典型地，數(shù)字巖石或巖心訓(xùn)練圖像的體素分辨率比全井眼圖像分辨率更高。為了確 保產(chǎn)生的假巖心模型與全井眼圖像具有一致的尺度特征，根據(jù)數(shù)字巖心分辨率與全井眼圖 像分辨率的比對(duì)訓(xùn)練圖像進(jìn)行粗糙取樣。圖5的步驟500公開了通過合并從采集的數(shù)值巖心數(shù)據(jù)和采集的井眼圖像測井圖 數(shù)據(jù)取出的結(jié)構(gòu)而生成數(shù)值假巖心的再現(xiàn)的步驟。步驟500a公開了運(yùn)用MPS SNESIN算法 生成2個(gè)或更多相的數(shù)值假巖心的再現(xiàn)，以使所述再現(xiàn)條件化，從而與在數(shù)字巖石的訓(xùn)練 圖像和全井眼圖像中觀測的相尺寸和形狀相匹配。在全井眼圖像和數(shù)字巖石上已經(jīng)繪制導(dǎo) 電塊(圖3和4)的情況下，可以對(duì)多個(gè)兩個(gè)相進(jìn)行建模。如果導(dǎo)電塊在數(shù)字巖石(訓(xùn)練圖 像)中沒有被識(shí)別，可以由固定數(shù)量的體素通過孔隙相的體變對(duì)所述導(dǎo)電塊進(jìn)行模擬。這 樣，可在數(shù)值模擬巖心中捕獲大尺度多相性。數(shù)值假巖心的半徑大小和高度僅僅受限于可 用的計(jì)算機(jī)內(nèi)存的限制。圖5的步驟600公開了將數(shù)值假巖心再取樣到徑向網(wǎng)格的步驟。為了實(shí)施數(shù)值假 巖心的流動(dòng)研究，將使用先前的步驟生成的笛卡爾坐標(biāo)系數(shù)值假巖心重新劃分到圓柱坐標(biāo) 形式的徑向網(wǎng)格中。在每個(gè)笛卡爾體素處，根據(jù)巖石類型(基質(zhì)，孔洞，或?qū)щ妷K)分配恒 定的孔隙度或滲透率。對(duì)于笛卡爾坐標(biāo)系的來每一個(gè)徑向網(wǎng)格單元說，由于所述單位包括 許多笛卡爾體素，因此生成了平均孔隙度和滲透率。通過對(duì)單元內(nèi)的笛卡爾體素的所有孔 隙度進(jìn)行算術(shù)平均而獲得平均孔隙度；通過執(zhí)行幾何平均得到滲透率。圖5的步驟700公開了執(zhí)行近井眼的流動(dòng)模擬的步驟。諸如注水的流體流動(dòng)數(shù)值 模擬在構(gòu)建假巖心時(shí)實(shí)施以評(píng)價(jià)重要的參數(shù)，例如，含水量、原油采收率、和采收率。對(duì)于 數(shù)值假巖心中的不同相來說，毛細(xì)管壓力和相對(duì)滲透率查找表提供了供給到諸如Ecllipse 的流動(dòng)模擬器的值。這是根據(jù)假巖心模型量化碳酸鹽巖石多相性對(duì)流體流動(dòng)影響的關(guān)鍵 步驟。毛細(xì)管壓力和相對(duì)滲透率可以由帶有相同巖石類型的巖心樣品的SCAL或MICP數(shù)據(jù) (如果可以獲得)獲得。更通常的情況是，孔隙網(wǎng)建模技術(shù)可以用于評(píng)價(jià)上述不同巖石類型 的參數(shù)。圖5的步驟800公開了確定用于提高評(píng)價(jià)諸如含水量、原油采收率、和采收率的參 數(shù)中的潤濕性作用的過程的步驟。圖6示出了板狀孔洞碳酸鹽巖樣品(約6英寸高和直徑為4英寸)的CT掃描的三 維視圖中圖Ia的訓(xùn)練圖像，使得孔洞被示出為紅色到灰色，而巖石基質(zhì)被示出為藍(lán)色。圖6 還顯示出所選擇的訓(xùn)練圖像，孔洞碳酸鹽巖的CT掃描(GoWelly，S.，2003，3D analysis ofvug connectivity,Indian Basinfield,New Mexico Unpublished M. S. thesis,Colorado School of Mines,Golden, CO) ο需要注意的是，圖la示出了來自相同深度和在相同的地層 的井眼圖像測井圖的3-ftl-m)間隔，并且孔洞是肉眼可見的大的、沒有規(guī)則的孔隙。更進(jìn) 一步，在井眼圖像測井圖中，因?yàn)榭紫短畛溆兴@井泥漿，并且所述孔隙導(dǎo)電，因此導(dǎo)電 孔洞可以被顯現(xiàn)為黑斑。圖7示出了在具體深度處的被彎曲成與井眼直徑(井眼直徑約為8. 5英寸，并且 長度間隔為3英尺或1米)。圖6的CT掃描訓(xùn)練圖像精確地標(biāo)刻并且位于井眼圖像中心的 正確深度處，其中孔洞被示出為紅色，而巖石基質(zhì)被示出為淺藍(lán)色。如上所述，井眼圖像卷 曲成三維圓柱形形狀，使得全井眼圖像并未示出，因?yàn)樗鼍蹐D像允許以相對(duì)尺度的訓(xùn) 練圖像(CT掃描)(位于圓柱體的中心)和測井圖像的視圖。為了觀察數(shù)值假巖心的模擬 (以下參見圖8)，需注意的是使用訓(xùn)練圖像和全井眼圖像。圖8圖示了由對(duì)應(yīng)圖7的訓(xùn)練圖像和全井眼圖像生成數(shù)值假巖心的結(jié)果，使得孔洞顯現(xiàn)為紅色而巖石基質(zhì)顯現(xiàn)為深藍(lán)色。要注意的是間隔的長度可以為3英尺(1米)，井 眼直徑可以為8. 5英寸(22厘米)。圖9示出了由對(duì)應(yīng)圖7的訓(xùn)練圖像和全井眼圖像生成數(shù)值假巖心的結(jié)果，使得孔洞不可見而巖石基質(zhì)顯現(xiàn)為灰色。進(jìn)一步地，任意懸浮的部分(見圖9上部)表示這是三維 模型，且在每個(gè)部分帶有大量孔(間隙長度約1英尺或0. 3米以及井眼直徑約8. 5英寸)。圖10示出了由對(duì)應(yīng)圖6的訓(xùn)練圖像和全井眼圖像生成的數(shù)值假巖心的結(jié)果，數(shù)值假巖心(參見圖9)示出了使用具有在每個(gè)孔周圍的8像素半徑的圓形體變。導(dǎo)電塊(顯 示為紅色)類似于如圖4所示的導(dǎo)電塊周圍的等高線。因此，此模型具有三個(gè)巖石物理相 孔隙或孔洞(顯現(xiàn)為綠色)、導(dǎo)電塊(顯現(xiàn)為紅色)、以及巖石基質(zhì)(顯現(xiàn)為藍(lán)色)。導(dǎo)電塊 提供孔隙或孔洞間的三維連通性，并且允許捕獲在大多數(shù)碳酸鹽巖石中所固有的多相性。圖11示出了圖10的具有不可見巖石基質(zhì)的數(shù)值的假巖心，所述不可見巖石基質(zhì)強(qiáng)調(diào)孔隙或孔洞(被示出為綠色)與導(dǎo)電塊(被示出為紅色)的三維互聯(lián)性質(zhì)。要注意的 是間隔長度可以是3英尺(1米)，井眼圖像直徑可為8. 5英寸(22厘米)。圖12示出了如圖10所示的數(shù)值假巖心，以及孔隙或孔洞(被示出為綠色)的復(fù)雜三維形狀，導(dǎo)電塊和巖石基質(zhì)都不可見。要注意的是間隔長度可以是3英尺(1米)，井眼 圖像直徑可為8. 5英寸(22厘米)。圖13示出了基于在井眼圖像或其它測井圖中觀測的層可以再取樣或者再劃分到其中徑向網(wǎng)格可以被分層的徑向網(wǎng)格的數(shù)值假巖心。圖13還示出了被生成以研究數(shù)值假 巖心模型的流動(dòng)特性的徑向網(wǎng)格。要注意的是原始數(shù)值假巖心內(nèi)部部分(直徑約4英寸) 被擴(kuò)孔，以允許布置微型發(fā)生器。微型注射器可位于假巖心的外邊界周圍。圖14a示出了在流動(dòng)模擬中用于導(dǎo)電塊的相對(duì)滲透率曲線；其中具有三條不同的相對(duì)滲透率曲線對(duì)應(yīng)于假巖心模型中諸如孔洞、導(dǎo)電塊和緊密巖石基質(zhì)的三種巖石類型的 每一個(gè)。具體地，圖14a示出了在流動(dòng)模擬中用于導(dǎo)電塊的相對(duì)滲透率曲線。注意圖14a 和圖14b中的假巖心均可為約1英尺(0. 3米)高，外徑約為8. 5英寸(22厘米)(或巖心 寬度為8. 5英寸(22厘米)以及內(nèi)徑約為4英寸)。圖14b示出了根據(jù)本發(fā)明的實(shí)施例的通過數(shù)值假巖心的流動(dòng)模擬結(jié)果，使得一連 串水微型注射器圍繞在直徑以外，一連串微型發(fā)生器包圍假巖心的內(nèi)徑顏色代表油飽和(So)， 其中在此流動(dòng)模型中顯示了多相進(jìn)步和突破。要注意的是圖14b顯示了注水?dāng)?shù)值假 巖心的油飽和度剖面圖。 進(jìn)一步地，雖然已經(jīng)參照示例性實(shí)施例說明了本發(fā)明，但是應(yīng)當(dāng)理解的是這里已 經(jīng)使用的詞匯是說明和示例性詞匯，而不是限制性性詞匯。在本發(fā)明的各個(gè)方面中在不背 離本發(fā)明的保護(hù)范圍和精神的情況下可以在所附權(quán)利要求的范圍內(nèi)如所述和所修改的做 各種改變。雖然在此已經(jīng)參照具體裝置、材料和實(shí)施例說明了本發(fā)明，但是本發(fā)明旨在不限 于這里所公開的細(xì)節(jié)；相反，本發(fā)明延伸到諸如在所附權(quán)利要求的保護(hù)范圍內(nèi)的所用功能 等價(jià)的結(jié)構(gòu)、方法和用途。
權(quán)利要求
一種用于生成數(shù)值假巖心模型的方法，包括以下步驟a)從包括儲(chǔ)層的深度限定間隔的所述儲(chǔ)層獲得測井?dāng)?shù)據(jù)，并且將所述測井?dāng)?shù)據(jù)處理成至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分，所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)具有未識(shí)別的井眼圖像數(shù)據(jù)；b)檢查所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分、其它已處理測井?dāng)?shù)據(jù)中的一個(gè)，或檢查所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分和其它已處理測井?dāng)?shù)據(jù)，以生成所述未識(shí)別的井眼圖像數(shù)據(jù)，然后將生成的所述未識(shí)別的井眼圖像數(shù)據(jù)處理成所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分，以生成至少一個(gè)卷曲的全井眼圖像數(shù)據(jù)；c)采集來自所述儲(chǔ)層的至少一個(gè)巖心、所述測井?dāng)?shù)據(jù)的一部分中的一個(gè)，或采集所述至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)的一部分，以及由采集的至少一個(gè)巖心、所述測井?dāng)?shù)據(jù)的一部分中的一個(gè)，或所述至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)的一部分生成數(shù)字巖心數(shù)據(jù)，使得生成的所述數(shù)字巖心數(shù)據(jù)表征所述儲(chǔ)層的一個(gè)或多個(gè)深度限定間隔的特征或結(jié)構(gòu)中的一個(gè)；以及d)處理生成的所述數(shù)字巖心數(shù)據(jù)的一部分、所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分或所述測井?dāng)?shù)據(jù)的一部分，以生成所述數(shù)值假巖心模型的再現(xiàn)。
2.如權(quán)利要求1所述的方法，其中，所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分 是連續(xù)變化的數(shù)值的二維(2D)標(biāo)量陣列的訓(xùn)練圖像。
3.如權(quán)利要求1所述方法，其中，所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分包 括已處理的原始數(shù)據(jù)，所述原始數(shù)據(jù)包括測量值和非測量值。
4.如權(quán)利要求1所述的方法，其中，所述未識(shí)別的井眼圖像數(shù)據(jù)表征所述至少一個(gè)能 夠解釋的井眼圖像數(shù)據(jù)的一部分的非測量值和所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的 一部分中的數(shù)據(jù)間隙中的一個(gè)。
5.如權(quán)利要求4所述的方法，其中所述數(shù)據(jù)間隙來自包括下列所述中的一個(gè)的組所 述儲(chǔ)層內(nèi)的至少一個(gè)損壞極板、所述儲(chǔ)層內(nèi)的至少一個(gè)損壞區(qū)域、相對(duì)于儲(chǔ)層內(nèi)的井壁具 有不足極板壓力的至少一個(gè)極板、被妨礙與儲(chǔ)層內(nèi)的所述井壁接觸的至少一個(gè)極板或所述 儲(chǔ)層內(nèi)的至少一個(gè)不能工作的極板。
6.如權(quán)利要求1所述的方法，其中，所述其它已處理的所述測井?dāng)?shù)據(jù)來自包括隨鉆測 井?dāng)?shù)據(jù)或電纜線測井?dāng)?shù)據(jù)、以及所述隨鉆測井?dāng)?shù)據(jù)和所述電纜線測井?dāng)?shù)據(jù)的一些組合中的 一個(gè)的組。
7.如權(quán)利要求6所述方法，其中，電纜測井?dāng)?shù)據(jù)包括伽馬射線測井圖、密度測井圖、聲 波測井圖、中子測井圖、井徑測井圖、和電阻率測井圖中的至少一個(gè)。
8.如權(quán)利要求1所述的方法，其中，所述至少一個(gè)卷曲的全井眼圖像數(shù)據(jù)是二維(2D)、 三維(3D)或二維和三維中的一種。
9.如權(quán)利要求1所述的方法，其中，所述生成至少一個(gè)全井眼圖像的步驟包括以下步驟1)從所述儲(chǔ)層的深度限定間隔選擇所述井眼圖像測井圖的深度限定間隔，使用多點(diǎn)統(tǒng) 計(jì)(MPS)算法處理所述井眼圖像測井圖的深度限定間隔，其中所述多點(diǎn)統(tǒng)計(jì)算法提供從所 述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分捕獲未識(shí)別的地質(zhì)結(jié)構(gòu)的模型，以識(shí)別所述 至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分中的數(shù)據(jù)位置；2)使用過濾器得分處理已識(shí)別的所述數(shù)據(jù)位置，以分組并然后模擬未識(shí)別所述井眼圖 像數(shù)據(jù)中的圖案，藉此，未識(shí)別的所述井眼圖像數(shù)據(jù)表征所述至少一個(gè)能夠解釋的井眼圖 像的一部分的非測量值和所述至少一個(gè)能夠解釋的井眼圖像的一部分中的數(shù)據(jù)間隙中的 一個(gè)；3)識(shí)別所述儲(chǔ)層的巖石物理相，其中所述巖石物理相來自包括下列所述中的一個(gè)的 組裂縫、孔隙、孔洞、所述儲(chǔ)層的井壁內(nèi)的巖石的導(dǎo)電塊、所述儲(chǔ)層的所述井壁內(nèi)的巖石的 電阻塊、地層界面和巖石基質(zhì)；以及4)將未識(shí)別的所述井眼圖像數(shù)據(jù)處理成所述至少一個(gè)能夠解釋的井眼圖像的一部分， 以生成所述至少一個(gè)全井眼圖像，然后利用已知的井眼直徑，將所述至少一個(gè)全井眼圖像 卷曲成原始三維(3D)形狀，藉此，數(shù)值假巖心被條件化，以與所述全井眼圖像的卷曲部分 和所述至少一個(gè)能夠解釋的井眼圖像的一部分中的一個(gè)相匹配。
10.如權(quán)利要求1所述的方法，其中，所述測井?dāng)?shù)據(jù)來自包括下列所述中的一個(gè)的組 具有多個(gè)探測深度的測井?dāng)?shù)據(jù)、來自一個(gè)或多個(gè)儲(chǔ)層的巖心數(shù)據(jù)的模擬模型、和所述測井 數(shù)據(jù)和所述模擬模型的組合中的一個(gè)。
11.如權(quán)利要求1所述的方法，其中，所述生成數(shù)字巖心數(shù)據(jù)的步驟包括以下步驟從采集的所述至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)的一部分中的一個(gè)識(shí)別巖石物理相，其中所述巖石物理相來自包括下列所述中的一個(gè)的組裂縫、孔隙、孔洞、所述儲(chǔ)層的井壁內(nèi)的 巖石的導(dǎo)電塊、所述儲(chǔ)層的井壁內(nèi)的巖石的電阻塊、地層界面和巖石基質(zhì)。
12.如權(quán)利要求11所述的方法，其中，識(shí)別所述巖石物理相的步驟包括以下步驟通過固定數(shù)量的體素由采集的所述至少一個(gè)巖心的其它巖石物理相的體變進(jìn)行模擬。
13.如權(quán)利要求1所述的方法，其中，所述至少一個(gè)巖心從至少一個(gè)其它儲(chǔ)層采集。
14.如權(quán)利要求1所述的方法，其中，來自所述至少一個(gè)巖心的第二巖心從一個(gè)或多個(gè)儲(chǔ)層獲得。
15.如權(quán)利要求1所述的方法，其中，步驟d)包括以下步驟將生成的所述數(shù)值假巖心模型的再現(xiàn)的數(shù)字文件繪制到數(shù)字介質(zhì)或硬拷貝介質(zhì)上。
16.如權(quán)利要求1所述的方法，其中，所述步驟d)的生成數(shù)值假巖心模型的再現(xiàn)包括以下步驟對(duì)兩個(gè)或更多個(gè)巖石物理相使用多點(diǎn)統(tǒng)計(jì)(MPS)算法生成數(shù)值假巖心的再現(xiàn)，以使所述再現(xiàn)條件化，從而與在生成的所述數(shù)字巖心數(shù)據(jù)的一部分和生成的所述至少一個(gè)卷曲的 全井眼圖像數(shù)據(jù)的至少一個(gè)能夠解釋的井眼圖像的一部分中觀測的相尺寸和形狀相匹配。
17.如權(quán)利要求16所述的方法，其中，兩個(gè)或更多個(gè)巖石物理相來自包括下列所述中的一個(gè)的組采集的所述至少一個(gè)巖心、測井?dāng)?shù)據(jù)的一部分、以及采集的所述至少一個(gè)巖心 和所述測井?dāng)?shù)據(jù)的部分，其中所述巖石物理相來自包括下列所述中的一個(gè)的組裂縫、孔 隙、孔洞、所述儲(chǔ)層的井壁內(nèi)的巖石的導(dǎo)電塊、所述儲(chǔ)層的井壁內(nèi)的巖石的電阻塊、地層界 面和巖石基質(zhì)。
18.如權(quán)利要求15所述的方法，進(jìn)一步包括以下步驟e)將步驟d)的所述數(shù)值假巖心再取樣到徑向網(wǎng)格。
19.如權(quán)利要求18所述的方法，其中，所述徑向網(wǎng)格提供所述數(shù)值假巖心模型的流動(dòng)研究，并且包括以下步驟1)構(gòu)建具有一個(gè)或多個(gè)形狀或一個(gè)或多個(gè)層中的一個(gè)的圓柱形網(wǎng)格；2)將來自笛卡爾坐標(biāo)的所述數(shù)值假巖心模型再采樣到圓柱坐標(biāo)，根據(jù)所述體素的相關(guān) 聯(lián)的巖石類型，所述笛卡爾坐標(biāo)的體素具有恒定孔隙度或恒定滲透率，其中每一個(gè)圓柱形 單元包括所述數(shù)值假巖心模型的多個(gè)笛卡爾體素；3)根據(jù)包括的所述多個(gè)笛卡爾體素，為每個(gè)圓柱形單元分配平均孔隙度和平均滲透 率；以及4)根據(jù)包括的所述多個(gè)笛卡爾體素的主要巖石類型，為每個(gè)圓柱形單元分配相對(duì)滲透 率和毛細(xì)管壓力曲線，對(duì)供給到模擬流動(dòng)器中的數(shù)值假巖心中的不同相使用毛細(xì)管壓力、 相對(duì)滲透率以及所述毛細(xì)管壓力和所述相對(duì)滲透率中的一個(gè)的參考表，由此根據(jù)所述假巖 心模型量化巖石多相性對(duì)流體流動(dòng)的影響，所述假巖心模型涉及具有相同巖石類型的一個(gè) 或多個(gè)巖石樣品的專門的巖心分析(SCAL)數(shù)據(jù)、壓汞毛細(xì)管壓力(MICP)數(shù)據(jù)、孔隙網(wǎng)結(jié)構(gòu) 模型、以及所述專門的巖心分析(SCAL)數(shù)據(jù)、所述壓汞毛細(xì)管壓力(MICP)數(shù)據(jù)和所述孔隙 網(wǎng)結(jié)構(gòu)模型的一些組合中的一個(gè)。
20.如權(quán)利要求19所述的方法，其中，所述相關(guān)聯(lián)的巖石類型來自包括下列所述中的 一個(gè)的組裂縫、孔隙、孔洞、所述儲(chǔ)層的井壁內(nèi)的巖石的導(dǎo)電塊、所述儲(chǔ)層的井壁內(nèi)的巖石 的電阻塊、地層界面和巖石基質(zhì)。
21.如權(quán)利要求18所述的方法，進(jìn)一步包括以下步驟f)執(zhí)行諸如井井眼區(qū)域的一個(gè)或多個(gè)地理相關(guān)聯(lián)的儲(chǔ)層的流動(dòng)模擬。
22.如權(quán)利要求21所述的方法，其中，所述執(zhí)行流動(dòng)模擬包括單相或多相流動(dòng)中的一 個(gè)，使得對(duì)步驟e)計(jì)算的所述數(shù)值假巖心模型執(zhí)行所述流動(dòng)模擬以評(píng)價(jià)一個(gè)或多個(gè)參數(shù)。
23.如權(quán)利要求22所述的方法，其中，所述一個(gè)或多個(gè)參數(shù)來自包括下列所述中的一 個(gè)的組孔隙度、滲透率、毛細(xì)管壓力、相對(duì)滲透率、含水量、原油采收率和采收率。
24.如權(quán)利要求21所述的方法，進(jìn)一步包括g)實(shí)施對(duì)與流動(dòng)相關(guān)的參數(shù)的敏感性分析，以確定用于注水、注氣、注蒸汽和提高采收 率(EOR)中的一個(gè)的過程。
25.如權(quán)利要求24所述的方法，其中，所述與流動(dòng)相關(guān)的參數(shù)來自包括下列所述中的 一個(gè)的組孔隙度、滲透率、毛細(xì)管壓力、相對(duì)滲透率、含水量、原油采收率和采收率。
26.如權(quán)利要求21所述的方法，進(jìn)一步包括以下步驟g)確定用于評(píng)價(jià)所估計(jì)的參數(shù)的潤濕性作用的過程，其所述估計(jì)的參數(shù)包括毛細(xì)管壓 力、相對(duì)滲透率、含水量、原油采收率和采收率中的一個(gè)。
27.一種用于生成數(shù)值假巖心模型的方法，包括以下步驟a)從包括儲(chǔ)層的深度限定間隔的所述儲(chǔ)層獲得測井?dāng)?shù)據(jù)，并且將所述測井?dāng)?shù)據(jù)處理成 至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分；b)檢查所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分、其它已處理測井?dāng)?shù)據(jù)中的一 個(gè)，或檢查所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分和其它已處理測井?dāng)?shù)據(jù)，以生 成至少一個(gè)卷曲的全井眼圖像數(shù)據(jù)；c)采集來自所述儲(chǔ)層的至少一個(gè)巖心、所述測井?dāng)?shù)據(jù)的一部分中的一個(gè)，或采集所述 至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)的一部分，使得生成的所述數(shù)字巖心數(shù)據(jù)表征所述儲(chǔ)層的一 個(gè)或多個(gè)深度限定間隔的特征或結(jié)構(gòu)中的一個(gè)；d)處理生成的所述數(shù)字巖心數(shù)據(jù)的一部分、所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的 一部分、所述測井?dāng)?shù)據(jù)的一部分，或所述數(shù)字巖心數(shù)據(jù)的一部分、所述至少一個(gè)能夠解釋的 井眼圖像數(shù)據(jù)的一部分和所述測井?dāng)?shù)據(jù)的一部分的一些組合，以生成所述數(shù)值假巖心模型 的再現(xiàn)。
28.如權(quán)利要求27所述的方法，其中，所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部 分是連續(xù)變化的數(shù)值的二維(2D)標(biāo)量陣列的訓(xùn)練圖像。
29.如權(quán)利要求27所述方法，其中，其它已處理的所述測井?dāng)?shù)據(jù)來自包括隨鉆測井?dāng)?shù) 據(jù)或電纜線測井?dāng)?shù)據(jù)、以及所述隨鉆測井?dāng)?shù)據(jù)和所述電纜線測井?dāng)?shù)據(jù)的一些組合中的一個(gè) 的組。
30.如權(quán)利要求27所述的方法，其中，所述測井?dāng)?shù)據(jù)來自包括下列所述中的一個(gè)的組 具有多個(gè)探測深度的測井?dāng)?shù)據(jù)、來自一個(gè)或多個(gè)儲(chǔ)層的巖心數(shù)據(jù)的模擬模型、以及所述測 井?dāng)?shù)據(jù)和所述模擬模型的組合。
31.如權(quán)利要求27所述的方法，其中，所述生成數(shù)字巖心數(shù)據(jù)的步驟包括以下步驟從采集的所述至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)的一部分中的一個(gè)識(shí)別巖石物理相，其中所述巖石物理相來自包括下列所述中的一個(gè)的組裂縫、孔隙、孔洞、所述儲(chǔ)層的井壁內(nèi)的 巖石的導(dǎo)電塊、所述儲(chǔ)層的井壁內(nèi)的巖石的電阻塊、地層界面和巖石基質(zhì)。
32.如權(quán)利要求27所述的方法，其中，識(shí)別所述巖石物理相的步驟包括以下步驟由固定數(shù)量的體素通過采集的所述至少一個(gè)巖心的其它巖石物理相的體變進(jìn)行模擬。
33.一種用于生成數(shù)值假巖心模型的系統(tǒng)，包括a)從包括儲(chǔ)層的深度限定間隔的所述儲(chǔ)層獲得測井?dāng)?shù)據(jù)，并且將所述測井?dāng)?shù)據(jù)處理成 至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分，所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)具有 未識(shí)別的井眼圖像數(shù)據(jù)；b)檢查所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分、其它已處理測井?dāng)?shù)據(jù)中的一 個(gè)，或檢查所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分和其它已處理測井?dāng)?shù)據(jù)，以生 成所述未識(shí)別的井眼圖像數(shù)據(jù)，然后將生成的所述未識(shí)別的井眼圖像數(shù)據(jù)處理成所述至少 一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分，以生成至少一個(gè)卷曲的全井眼圖像數(shù)據(jù)；c)采集來自所述儲(chǔ)層的至少一個(gè)巖心、所述測井?dāng)?shù)據(jù)的一部分中的一個(gè)，或采集所述 至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)的一部分，并且由采集的至少一個(gè)巖心、所述測井?dāng)?shù)據(jù)的一 部分中的一個(gè)或所述至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)的一部分生成數(shù)字巖心數(shù)據(jù)，使得生成 的所述數(shù)字巖心數(shù)據(jù)表征所述儲(chǔ)層的一個(gè)或多個(gè)深度限定間隔的特征或結(jié)構(gòu)中的一個(gè)；以 及d)處理生成的所述數(shù)字巖心數(shù)據(jù)的一部分、所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的 一部分和所述測井?dāng)?shù)據(jù)的一部分，或生成的所述數(shù)字巖心數(shù)據(jù)的一部分、所述至少一個(gè)能 夠解釋的井眼圖像數(shù)據(jù)的一部分和所述測井?dāng)?shù)據(jù)的一部分的一些組合以生成所述數(shù)值假 巖心模型的再現(xiàn)，其中，所述生成數(shù)值假巖心模型的再現(xiàn)包括以下步驟對(duì)兩個(gè)或更多個(gè)巖石物理相使用多點(diǎn)統(tǒng)計(jì)(MPS)算法生成數(shù)值假巖心的再現(xiàn)，以使所 述再現(xiàn)條件化，從而與在生成的所述數(shù)字巖心數(shù)據(jù)的一部分和生成的所述至少一個(gè)卷曲的 全井眼圖像數(shù)據(jù)的至少一個(gè)能夠解釋的井眼圖像的一部分中觀測的相尺寸和形狀相匹配。
34. 一種用于生成數(shù)值假巖心模型的系統(tǒng)，包括a)從包括儲(chǔ)層的深度限定間隔的所述儲(chǔ)層獲得測井?dāng)?shù)據(jù)，并且將所述測井?dāng)?shù)據(jù)處理成至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分；b)檢查所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分、其它已處理測井?dāng)?shù)據(jù)中的一個(gè)，或檢查所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分和其它已處理測井?dāng)?shù)據(jù)，以生 成至少一個(gè)卷曲的全井眼圖像數(shù)據(jù)；c)采集來自所述儲(chǔ)層的至少一個(gè)巖心、所述測井?dāng)?shù)據(jù)的一部分中的一個(gè)，或采集所述至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)的一部分，使得生成的所述數(shù)字巖心數(shù)據(jù)表征所述儲(chǔ)層的一 個(gè)或多個(gè)深度限定間隔的特征或結(jié)構(gòu)中的一個(gè)；以及d)處理生成的所述數(shù)字巖心數(shù)據(jù)的一部分、所述至少一個(gè)能夠解釋的井眼圖像數(shù)據(jù)的一部分、所述測井?dāng)?shù)據(jù)的一部分，或所述數(shù)字巖心數(shù)據(jù)的一部分、所述至少一個(gè)能夠解釋的 井眼圖像數(shù)據(jù)的一部分和所述測井?dāng)?shù)據(jù)的一部分的一些組合，以生成所述數(shù)值假巖心模型 的再現(xiàn)，其中，所述生成數(shù)字巖心數(shù)據(jù)的步驟包括以下步驟從采集的所述至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)的一部分中的一個(gè)識(shí)別相，其中所述相來自包括下列所述中的一個(gè)的組裂縫、孔隙、孔洞、所述儲(chǔ)層的井壁內(nèi)的巖石的導(dǎo)電塊、所述 儲(chǔ)層的井壁內(nèi)的巖石的電阻塊、地層界面和巖石基質(zhì)。
全文摘要
本發(fā)明公開了一種用于生成數(shù)值假巖心模型的方法和系統(tǒng)，包括a)從具有儲(chǔ)層的深度限定間隔的儲(chǔ)層獲得測井?dāng)?shù)據(jù)，并且將測井?dāng)?shù)據(jù)處理成至少一個(gè)可解釋的井眼圖像數(shù)據(jù)，所述至少一個(gè)可解釋的井眼圖像數(shù)據(jù)具有未識(shí)別的井眼圖像數(shù)據(jù)；b)檢查至少一個(gè)可解釋的井眼圖像數(shù)據(jù)、其它已處理測井?dāng)?shù)據(jù)中的一個(gè)，或檢查至少一個(gè)可解釋的井眼圖像數(shù)據(jù)和其它已處理數(shù)據(jù)，以生成未識(shí)別的井眼圖像數(shù)據(jù)，將生成的未識(shí)別的井眼圖像數(shù)據(jù)處理成至少一個(gè)可解釋的井眼圖像數(shù)據(jù)，以生成卷曲的全井眼圖像數(shù)據(jù)；c)采集來自儲(chǔ)層的至少一個(gè)巖心、測井?dāng)?shù)據(jù)中的一個(gè)，或采集所述至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)，并且由采集的至少一個(gè)巖心、測井?dāng)?shù)據(jù)中的一個(gè)或所述至少一個(gè)巖心和所述測井?dāng)?shù)據(jù)生成數(shù)字巖心數(shù)據(jù)，使得生成的所述數(shù)字巖心數(shù)據(jù)表征儲(chǔ)層的一個(gè)或多個(gè)深度限定間隔的特征或結(jié)構(gòu)中的一個(gè)；以及d)處理生成的所述數(shù)字巖心數(shù)據(jù)、至少一個(gè)可解釋的井眼圖像數(shù)據(jù)或所述測井?dāng)?shù)據(jù)，以生成數(shù)值假巖心模型的再現(xiàn)。
文檔編號(hào)G03B35/16GK101802649SQ200980000184
公開日2010年8月11日 申請日期2009年4月10日 優(yōu)先權(quán)日2008年4月10日
發(fā)明者尼爾·弗蘭西斯·赫爾利, 張團(tuán)峰, 趙偉書 申請人:普拉德研究及開發(fā)股份有限公司