基于孔隙結構表征的低滲透砂巖流體賦存特征及滲透率評價

摘要：以渤海灣盆地G油田古近紀低滲透砂巖儲層為研究對象，綜合鑄體薄片、掃描電鏡（SEM）、高壓壓汞（MICP）及核磁共振（NMR）分析探討了孔隙結構及其對儲層品質和流體賦存的影響。低滲透砂巖儲層的孔喉大小分布范圍廣，微、納米級孔喉共存，且普遍呈現雙峰分布的特征。大孔喉系統主要由殘余粒間孔和溶蝕擴大孔組成；而微小孔喉系統則主要由粒內溶孔和晶間孔構成。研究結果表明，孔喉連通能力及連通孔喉的體積分數是影響低滲透砂巖儲層儲集物性和流體賦存特征的關鍵因素。物性越好的砂巖儲層中，連通性好的大孔喉（右峰）所占體積分數越高，儲層滲流能力越強；反之，物性越差的儲層孔隙系統則主要由連通性差的微小孔喉（左峰）主導，儲層滲流能力較弱。束縛流體主要賦存在孔喉半徑較小、微觀連通性較差、分形維數較大的小孔喉系統中；而中、大孔隙系統的微觀連通性好、分形維數較小，主要富集可動流體。研究還明確了半徑最大的孔喉是影響低滲透砂巖儲層滲透率和流體賦存特征的關鍵因素，并根據毛管壓力曲線和核磁共振T2（橫向弛豫時間）譜參數建立了兩套適用于低滲透砂巖儲層的滲透率評價模型。與經典模型（如Coates和Schlumberger Doll research（SDR））相比，新模型具有更好的有效性和適用性。

關鍵詞：渤海灣盆地；低滲透砂巖；孔隙結構；儲層質量；可動流體

doi：10.13278/j.cnki.jjuese.20230056

中圖分類號：TE122.2

文獻標志碼：A

Supported by the National Science and Technology Major Project Foundation of China （2016ZX05058001）

Fluid Occurrence State and Permeability Evaluation of Low-Permeability Sandstone Based on Pore Structure Characterization

Liu Zongbin， Li Chao， Lu Yan， Wang Ya， Huang Jianting

Tianjin Branch of CNOOC China Limited，Tianjin 300452，China

Abstract： This paper selected the Eocene low-permeability sandstones in the G oilfield， Bohai Bay basin as an example， and the influences of pore structure on reservoir quality and fluid occurrence were investigated by cast thin sections， scanning electron microscope （SEM）， mercury injection capillary pressure （MICP）， and nuclear magnetic resonance （NMR）. Due to a highly variable pore size distribution， micro-scale and nano-scale pore throats coexist. The pore-throat size distributions generally present bimodal behaviors. The large pore systems are mainly composed of residual intergranular pores and enhanced intergranular dissolution pores，while the micropore systems are dominated by intragranular dissolution pores and intergranular pores. The results indicate that pore connectivity is the key factor affecting the physical properties and fluid occurrence state of the low permeability sandstone reservoirs. The sandstones with better petrophysical properties have a high proportion of macroscopic pore throats （right peak） with good connectivity， and the permeability capacity of the sandstones is strong. In contrast， the pore systems of the sandstones with poor petrophysical properties are dominated by the microscopic pore throats （left peak） with poor pore connectivity， and the permeability capacity of the sandstones is weak. Bound fluids mainly exist in small pore systems with poor pore connectivity and large fractal dimensions， while large pore systems have good pore connectivity and small fractal dimensions and mainly contain movable fluids. The study also clarified that the permeability and fluid occurrence state of low-permeability sandstone reservoirs is mainly controlled by the relatively large pore throats in the pore network. It is revealed that the permeability prediction model based on MICP data and porosity derived from the “effective movable fluid method” provides a better estimation eect compared with the classic models （e.g.， Coates and SDR）.

Key words： Bohai Bay basin; low-permeability sandstones; pore structure; reservoir quality; movable fluids

0 引言

古近紀低滲透砂巖儲層作為渤海灣盆地G油田最重要的儲層類型之一，具有巨大的油氣勘探開發潛力[12]。低滲透砂巖儲層孔隙結構復雜、孔喉大小分布多變，微觀非均質性強，導致儲層評價及宏觀預測難度大[3]。另外，復雜的孔隙結構極大程度地控制著流體在砂巖中的賦存狀態及滲流規律[46]。因此，微觀孔隙結構對油氣藏評價及開發影響重大，決定著油氣藏產能的差異分布及油田的開發效果。準確評價低滲透砂巖儲層復雜孔隙結構，明確復雜孔隙結構對儲層品質、流體賦存的控制作用，探知不同孔隙系統控制下的滲流特征及滲流規律，對儲層精細評價、預測及提高采收率具有重大的理論和實際意義。

國內外學者[79]圍繞微觀孔隙結構特征，如孔喉幾何形態、類型、大小分布、連通性及潤濕性等開展了諸多研究工作。滲透率作為微觀孔隙結構的宏觀體現，也是評價儲層質量和控制流體在儲層中滲流行為的一個關鍵參數[1011]。目前，業內學者在滲透率評價方面開展了一些研究工作，并建立了用于滲透率評價的經驗方程、模型等，主要包括Kozeny-Carman （KC）方程、PaRiS方程、Winland模型、Schlumberger Doll research（SDR）模型以及Coates模型等[1216]。這些模型在孔隙結構簡單的中、高滲儲層的滲透率評價方面取得了較好的應用效果。但對于低滲透砂巖儲層，其復雜的孔隙結構系統和孔喉連通性是導致低滲透特性的重要因素，也是對滲透率評價起決定性作用的因素。以往的滲透率解釋模型難以有效地刻畫低滲透砂巖不同尺度孔喉及其連通性控制下的滲流特征差異。

本文綜合鑄體薄片、掃描電鏡（scanning electron microscope， SEM）、高壓壓汞（mercury injection capillary pressure， MICP）及核磁共振（nuclear magnetic resonance， NMR）精細刻畫了低滲透砂巖的微觀孔隙結構特征。在此基礎上，探究微觀孔隙結構特征對儲層品質和流體賦存的影響，建立適用于低滲透砂巖的滲透率預測模型，對低滲透砂巖儲層流體流動機制研究、儲層質量及非均質性評價具有重要意義。

1 區域地質概況

G油田位于中國渤海灣盆地的北東部（圖1）。古近系沙河街組是G油田的重要含油氣地層單元之一，自下而上可細分為4個沉積單元，即沙四段（Es4）、沙三段（Es3）、沙二段（Es2）、沙一段（Es1）（圖2）。沙四上亞段（Es4s）是本研究的目標層段，其地層埋深為2 160 ～3 563 m，發育廣泛濱淺湖灘壩

砂沉積體[1315]。巖性類型主要為細砂巖、粉砂巖和泥質粉砂巖（圖2）。由于受古地貌、古水流、古沉積物供應和古基準面變化的影響，砂體表現出復雜的橫向遷移和縱向變化特征，儲層宏觀和微觀非均質性強，而儲層微觀非均質性在很大程度上受控于微觀孔隙結構。

2 儲集空間與孔隙結構特征

2.1 儲集空間特征

根據鑄體薄片和掃描電鏡的巖相學分析可知，低滲透砂巖的孔隙系統主要由次生溶孔、殘余粒間孔隙、微裂縫和微孔隙組成，孔隙空間大小呈納米級—微米級不等，相應的與孔隙連接的喉道類型也復雜多樣。從鑄體薄片巖相學觀察結果來看（鑄體薄片無法獲取微孔隙體積分數），細砂巖平均總面孔率最高，孔隙系統主要由原生孔隙和次生溶孔組成；粉砂巖平均總面孔率次之，孔隙系統主要由次生溶孔組成；泥質粉砂巖的平均總面孔率最低，其孔隙系統主要由原生孔隙和次生溶孔組成（圖3）。

由不穩定組分溶蝕產生的次生溶孔是研究區主導的孔隙類型。部分次生溶孔與殘余粒間孔隙組成了擴大孔，顯著地改善了孔隙系統的儲集能力（圖4a、b）。與殘余粒間孔隙組成的擴大孔相連的喉道主要為孔隙縮小型和縮頸狀喉道，而局部不完全溶蝕作用形成的粒內次生溶孔則多被微孔喉構成的管束狀喉道連接（圖4c），在鑄體薄片的分辨率下呈孤立的死孔。

殘余粒間孔隙主要分布于粒徑較粗的細砂巖中，而在粉砂巖、泥質粉砂巖中少見。連接殘余粒間孔隙的喉道以孔隙縮小型、縮頸狀和片狀喉道為主?？紫犊s小型喉道主要發育在顆粒支撐的細砂巖中局部顆粒之間不接觸的位置（圖4d）；縮頸狀喉道多見于點接觸、點線接觸的砂巖中（圖4e）；片狀喉道常見于點接觸、點線接觸的砂巖中，喉道極細，是低滲透砂巖中最常見的喉道類型（圖4f）。

微裂縫對砂巖儲集空間的貢獻較?。▓D4g），但微裂縫能夠顯著改善孔隙系統的連通性和滲流能力，也能夠通過提供有效的孔隙流體滲流通道，促進溶蝕反應的進行，從而對儲層品質的改造具有重要意義。

掃描電鏡顯示，微孔隙主要包括自生黏土礦物晶間孔（圖4h）及長石、巖屑等碎屑顆粒不完全溶蝕形成的微孔隙（圖4i）。自生黏土礦物晶間孔和碎屑顆粒溶蝕形成的微孔隙對提高儲層孔隙度（φ）具有一定的促進作用，但由于微孔隙的孔徑通常小于1 μm，對儲層滲流性能基本沒有貢獻。

2.2 孔隙結構特征

2.2.1 孔喉大小及分布

研究區低滲透砂巖儲層的孔隙結構復雜、類型多樣，孔喉大小分布范圍廣，自微米級到納米級均有發育（圖5）。按照其在低滲透砂巖儲層中的占比情況，可大致劃分三類。

占比最多的孔隙結構毛管壓力曲線在進汞階段形態較陡峭，平緩段不發育，且排驅壓力較高（圖5a）。該類孔隙結構主要發育在泥質粉砂巖和粉砂巖中，所對應的儲層儲集物性較差（滲透率（k）多分布在1×10-3 μm2以下）。喉道半徑大多小于1.0 μm（圖5b），孔隙系統以小孔喉為主，大孔喉發

育較少，小孔喉、中孔喉和大孔喉所占的平均體積分數分別為77.94%，20.51%和1.55%（圖5c）。

占比次多的孔隙結構的毛管壓力曲線在進汞階段曲線形態較平緩，同時表現出較高的排驅壓力（圖5d），表明孔隙系統以小孔喉為主?？讖椒植挤秶^窄，孔喉半徑基本在5.0 μm以下（圖5e），且小孔喉、中孔喉和大孔喉所占的平均體積分數分別為29.03%，37.79%和33.18%（圖5f）。該類孔隙結構主導的儲層儲集物性較好（滲透率多分布在1×10-3～10×10-3 μm2之間），主要發育在粉砂巖和細砂巖中。

占比最少的孔隙結構的毛管壓力曲線在進汞階段表現出一個較長的平緩段，具有較低的排驅壓力（圖5g），表明孔喉連通性較好。該類孔隙結構主要發育于細砂巖和粉砂巖中，所對應的儲層儲集物性最好（滲透率多分布在10×10-3 μm2以上），孔隙系統主要由大孔喉主導（圖5h）。據統計分析，小孔喉、中孔喉和大孔喉所占的平均體積分數分別為24.22%，24.43%和51.35%（圖5i）。

聯合壓汞核磁標定的方式獲取了低滲透砂巖全孔徑范圍內的孔隙大小分布[1617]。如圖6所示，基于核磁共振數據獲取的孔隙大小分布特征與壓汞獲取的喉道大小分布特征整體呈現相同的趨勢。當砂巖物性致密時（滲透率lt;0.1×10-3 μm2），孔隙和喉道均呈現單峰分布特征，以半徑小于0.1 μm的小孔喉占絕對主導（圖6a），儲集空間主要由晶間孔和微孔隙組成（圖4h、i）；隨著滲透率增加，砂巖孔喉開始呈現雙峰多峰分布特征（圖6b、c），主要由殘余粒間孔隙構成的半徑大于1.0 μm的大孔喉也開始發育（圖4a），但砂巖滲透率在小于5.0×10-3 μm2時，其孔喉系統仍是由小孔喉和中孔喉主導，其儲集空間仍主要由晶間孔和粒內溶孔組成；隨著滲透率增加，大孔喉開始占據主導地位（圖6d、e），殘余粒間孔（圖4d）和粒間溶蝕擴大孔發育，孔喉分布開始向大孔喉絕對主導的方向演化（圖6f）。

2.2.2 孔隙結構分形特征及表征

分形維數的概念是Mandelbrot于1977年首次提出的，其核心內涵是分形物體的自相似性不隨著放大倍數或者尺度的變化而改變[18]。通過計算核磁共振和高壓壓汞數據的分形維數，可以定量地表征孔隙結構在幾何學方面的不規則性和復雜程度，進而定量地評價孔隙結構的幾何非均質性[1920]。

通過繪制進汞飽和度（SHg）和毛管壓力（pc）的雙對數交會圖即可得到基于高壓壓汞數據計算的分形維數（D），lg（SHg）和lg（pc）交會圖中的轉折點實際上對應于Pittman交會圖的apex點[13]（圖7）。因此，rapex為區分較大孔喉系統與較小孔喉系統的孔徑界限，將孔喉系統分割為分形特征截然不同的兩部分[21]。根據毛管壓力曲線計算的分形結果可知[1920]，較小孔喉系統（孔喉半徑（r）＜rapex）的分形維數D2為2.1～2.8，平均值為2.3，而較大孔喉系統（r≥rapex）的分形維數D1變化較大，為3.4～6.5，平均值為4.7。因此，較大的孔喉系統要比較小的孔隙系統更為復雜，非均質性也更強，原因可能為大孔喉系統成因和幾何形狀更為復雜，而小孔喉系統通常為微孔喉，成因相對簡單，幾何形狀相對規則。

從基于核磁共振實驗獲取的累積孔隙體積分數（PV）和橫向弛豫時間（T2）的雙對數交會曲線（圖8）可以看出，隨著橫軸lg（T2）數值的增加，曲線的斜率逐漸減小，即孔隙半徑增大，分形維數增大，這也與高壓壓汞實驗得到的結論一致，即大孔喉相對于小孔喉具有更加復雜和非均質的幾何學特征，從而呈現出更高的分形維數。

3 孔隙結構對儲層品質的控制作用

孔隙度和滲透率是表征儲層巖石物理性質的關鍵參數，分別反映了儲層的儲集能力和滲流能力?？紫抖群蜐B透率的大小及分布與孔隙結構密切相關，其中微觀孔隙結構特征顯著控制了滲透率的大小及其與孔隙度的對應關系[21]。由于壓汞實驗能夠將儲層微觀孔隙結構參數與儲層宏觀物性參數聯系起來，因此主要采用壓汞參數來探討孔喉連通性和孔喉大小分布對儲集物性的影響。

3.1 孔喉連通性對儲層物性的影響

定義Papex為較大有效連通孔喉（r≥rapex）所占的體積分數[21]。這些孔隙體積均由較大的連通孔喉組成，也對滲透率起主要貢獻作用?；貧w分析表明，滲透率相同的砂巖樣品的Papex變化范圍較大，即滲透率不僅取決于連通孔隙的體積，還受孔喉連通能力的控制（圖9a）。一般而言，孔喉半徑越大，孔喉的連通能力越強。大孔喉和這些大孔喉連通的孔隙體積對控制儲層質量起著重要作用。此外，滲透率與排驅壓力和之間呈現出良好的負相關關系，即排驅壓力越大對應最大孔喉半徑越小，孔喉連通性越差，砂巖的滲流能力越差（圖9b）。

為了更加定量、深入地表征具有差異連通性的不同半徑的孔喉對儲層物性的影響，對半徑小于0.1 μm的孔喉和半徑大于等于1.0 μm的孔喉所連接的孔隙體積占總孔隙空間的體積分數進行統計分析。在孔喉半徑小于0.1 μm的孔喉主導的孔喉系統中，由于孔隙系統連通性差，儲層滲透率較低，一般小于1×10-3 μm2，且在同等孔隙度下，呈現出的滲透率顯著低于其他兩種孔隙系統，常出現“相對高孔低滲”的現象；而在以半徑大于等于1.0 μm的孔喉主導的孔隙系統中，孔隙空間連通性好，儲層滲透率較高，多大于10×10-3 μm2，且在同等孔隙度條件下顯示出最高的滲透率，常出現“相對低孔高滲”的現象（圖10）。

3.2 孔喉大小分布對儲層物性的影響

砂巖的最大孔喉半徑（rmax）和中值孔喉半徑（r50）與滲透率之間均存在很強的指數相關性（R2≥0.77，圖11），表明孔喉大小分布在控制砂巖儲層滲流能力方面起關鍵作用。泥質粉砂巖的最大孔喉半徑基本小于1.0 μm（圖11a），表明其孔隙系統由微孔喉主導。微孔喉在孔隙網絡中普遍不連通或連通性很差，對滲透率貢獻微乎其微，這是導致泥質粉砂巖樣品出現“相對高孔低滲”的主要原因。此外，微孔喉通常賦存在黏土礦物中，而黏土礦物的結晶習性通常導致孔隙空間具有很高的迂曲度，不利于巖

石中的流體流動。與之相反，細砂巖往往具有較大的孔喉系統，細砂巖的中值孔喉半徑一般大于0.5 μm（圖11b），最大孔喉半徑基本在1.0 μm以上（圖11a）。細砂巖中的大孔喉系統通常與埋藏過程中保存較好的原生孔隙有關，這些原生孔隙孔徑較大，更有可能被相鄰的喉道連通，并在砂巖中形成連通的路徑。即使這些較大的孔喉只占孔隙空間的一小部分，它們也對滲透率具有顯著的貢獻。

通過統計不同滲透率級別砂巖樣品孔喉半徑與絕對滲透率的關系可知：不同半徑的孔喉對滲透率的貢獻差異顯著，隨著孔喉半徑的增大，其對滲透率的絕對貢獻升高（圖12）。對于滲透率大于1×10-3 μm2的砂巖，孔隙系統中需存在半徑大于1.0 μm的孔喉，而滲透率大于10×10-3 μm2時，砂巖的孔隙系統中需存在半徑大于4.0 μm的孔喉（圖12），這也正是低滲透砂巖孔隙系統最大孔喉半徑與滲透率存在強正相關性（圖11a）的主要原因。因此，總的來說，連通的孔隙體積以及大孔喉的半徑及體積分數控制了儲層的滲流能力。

4 孔隙結構對流體賦存狀態的影響

4.1 可動流體與束縛流體分析

核磁共振實驗過程采用鹽度為12 000 mg/L的氯化鈉溶液對巖樣進行飽含水處理，并采集100%飽含水狀態下的核磁共振T2譜。以6 000 r/min的轉速對飽含水的巖樣進行離心處理，去除巖心柱塞中的自由水，然后再次測量以獲得離心狀態下的核磁共振T2譜。其中，飽含水的巖樣中存在兩種類型的流體[22]，即可動流體和束縛流體。離心后仍保留在黏土礦物微孔隙中或受毛細管力束縛而難以流動的流體稱為束縛流體，其賦存的孔隙稱為束縛流體孔隙[2324]。中、大孔隙內的流體由于受孔隙壁面巖石骨架作用力的影響較小，毛細管力對其束縛作用不大，主要富集可動流體。低滲透砂巖儲層中的流體賦存狀態與孔喉大小分布及連通性緊密相關：可動流體總是傾向于賦存在核磁獲取的孔隙大小分布的右峰，即連通性好、較大的孔隙系統內，并且隨著滲透率增加，右峰的可動流體孔隙占比顯著增加，左峰的束縛流體孔隙占比逐漸增大（圖13）。

4.2 可動流體賦存影響因素分析

低滲透砂巖儲層強烈的微觀非均質性，復雜多變的孔喉網絡極大程度地控制著砂巖儲層中流體的賦存狀態，是影響砂巖儲層中滲流行為的關鍵因素。

由于大孔隙通常具有更為復雜的幾何形態，因而在核磁共振實驗中，可動流體孔隙的分形維數（2.6～3.0）也明顯大于束縛流體孔隙的分形維數（1.2～2.4）（圖14a、b）。從核磁共振實驗獲取的分形維數來看：可動流體對應的大孔隙部分具有分形特征，且與滲透率具有較好的負相關性（圖14a），表明可動流體孔隙的幾何形態越復雜，越不利于流體在孔隙網絡中的滲流；而束縛流體孔隙則不具有分形特征，其分形維數值與滲透率大小基本沒有相關性（圖14b）。此外，作為微觀孔隙結構連通性的宏觀表現，低滲透砂巖儲層的滲透率和儲層品質指數都與可動流體飽和度呈現出較好的正相關性，表明隨著儲層品質變好和儲層滲透率的增大，儲層孔隙系統連通性變好，流體的可動性也隨之增大（圖14c、d）。可動流體參數能夠表征儲層的儲集能力和流體在巖石孔隙中的賦存特征，是評價儲層儲集、滲流性能的重要參數。

5 基于孔隙結構的滲透率評價

5.1 基于改進的Winland方程的滲透率評價模型

由于滲透率與孔喉大小分布之間存在很強的相關性，因此基于孔喉半徑的滲透率預測模型通常表現出較好的預測效果[8]。Winland模型描述了滲透率、孔喉半徑和孔隙度之間的關系，并指出常規儲層的滲透率主要受r35和孔隙度的控制[9]。有關學者[10，12]對Winland方程進行了擴展和重新擬合，發現與滲透率擬合度最高的孔喉半徑隨著儲層類型改變而變化。隨著儲層變得更加致密，最佳擬合半徑向大孔喉方向轉變。本研究通過對滲透率、孔隙度和間隔孔喉半徑（r5～r50）的多元回歸分析，建立了滲透率與間隔孔喉半徑的多元回歸方程。結果表明（表1），研究區低滲透砂巖儲層的滲透率與r10的相關性最好。此外，研究還發現沙四上亞段儲層的r10主要集中在0.1～5.0 μm之間，表明低滲透砂巖儲層的滲透率主要受大孔喉系統控制。

如圖15所示，預測滲透率與實測滲透率的擬合度較高，驗證了該滲透率預測模型的有效性。此外，根據建立的經驗方程，還可以通過滲透率和孔隙度反向預測所對應的孔喉半徑。

5.2 基于可動流體孔隙度的滲透率評價模型

前文研究表明，砂巖中不連通或弱連通的孔隙對滲透率的貢獻很小。因此，與總孔隙度相比，可動流體孔隙度（φm）更加適用于滲透率預測[11]，φm可以通過式（1）計算得到。φm與滲透率的回歸分析表明，φm和滲透率之間存在很強的指數相關性（R2=0.92，表2）。因此，利用φm來預測滲透率是可行的，基于φm的滲透率預測模型如式（2）所示。

φm=φNSm；（1）

k=0.0423e0.7899φm。（2）

式中：φN為核磁共振孔隙度，%；Sm為可動流體飽和度，%。

如圖16所示，基于可動流體孔隙度模型的滲透率預測值與實測值具有較高的吻合度。此外，與經典的Coates和Schlumberger Doll research（SDR）模型相比，可動流體孔隙度模型提供了更好的預測效果（表2，圖16），表明新模型對于滲透率預測具有更好的有效性和適用性。

6 結論

1）低滲透砂巖儲層的孔隙結構復雜、微觀非均質性強，孔喉大小分布范圍廣，普遍呈現微、納米級孔喉共存的多尺度分布特征。孔隙系統主要由次生溶孔、殘余粒間孔隙、

微裂縫和微孔隙組成。其中，大的殘余粒間孔和溶蝕擴大孔主要通過孔隙縮小型和縮頸狀喉道相連，而粒內溶孔和微孔隙則主要通過管束狀喉道連接。

2）儲層的宏觀物性極大程度上受控于微觀孔隙結構。整體上，孔喉半徑越大、連通性越好的孔喉對滲透率的貢獻越顯著。孔喉半徑大于等于1.0 μm的孔喉所主導的孔隙系統對滲透率的貢獻作用顯著；而由不連通或弱連通的孔喉（rlt;0.1 μm）所主導的孔隙系統的滲透率值多小于1×10-3 μm2。束縛流體主要賦存在孔喉半徑較小、微觀連通性較差的小孔喉系統中；而中、大孔隙系統的微觀連通性好，且孔隙內的流體受孔隙壁面巖石骨架作用力的影響較小，毛細管力對其束縛作用較小，主要富集可動流體。

3）低滲透砂巖儲層的孔喉大小及連通性的非均質性強，滲透率主要受控于孔隙網絡中半徑最大的孔喉，r10是評價研究區低滲透砂巖儲層滲透率的最佳孔喉參數。此外，考慮到砂巖儲層的滲透率主要由大而有效連通的孔隙貢獻，不連通或弱連通的孔隙對滲透率的貢獻很小，構建了基于可動流體孔隙度的滲透率預測模型，相較于經典的Coates和SDR模型，展現出更好的預測性能。

參考文獻（References）：

[1] 賈承造， 鄭民， 張永峰. 中國非常規油氣資源與勘探開發前景[J]. 石油勘探與開發， 2012， 39（2）： 129136.

Jia Chengzao， Zheng Min， Zhang Yongfeng. Unconventional Hydrocarbon Resources in China and the Prospect of Exploration and Development[J]. Petroleum Exploration and Development， 2012， 39（2）： 129136.

[2] 吳松濤， 朱如凱， 李勛， 等. 致密儲層孔隙結構表征技術有效性評價與應用[J]. 地學前緣， 2018， 25（2）： 191203.

Wu Songtao， Zhu Rukai， Li Xun， et al. Evaluation and Application of Porous Structure Characterization Technologies in Unconventional Tight Reservoirs[J]. Earth Science Frontiers， 2018， 25（2）： 191203.

[3] 李彥澤， 商琳， 王群會， 等. 斷陷盆地強非均質低滲致密砂巖儲層綜合分類： 以渤海灣盆地南堡凹陷為例[J]. 吉林大學學報（地球科學版）， 2022， 52（6）： 18301843.

Li Yanze， Shang Lin， Wang Qunhui， et al. Comprehensive Classification of Highly Heterogeneous Low Permeability-Tight Sandstone Reservoir in Faulted Basin： Taking Nanpu Sag of Bohai Bay Basin as an Example[J]. Journal of Jilin University （Earth Science Edition）， 2022， 52（6）： 18301843.

[4] 余義常， 徐懷民， 高興軍， 等. 海相碎屑巖儲層不同尺度微觀剩余油分布及賦存狀態： 以哈得遜油田東河砂巖為例[J]. 石油學報， 2018， 39（12）： 13971409.

Yu Yichang， Xu Huaimin， Gao Xingjun， et al. Distribution and Occurrence Status of Microscopic Remaining Oil at Different Scales in Marine Clastic Reservoirs： A Case Study of Donghe Sandstone in Hadeson Oilfield[J]. Acta Petrolei Sinica， 2018， 39（12）： 13971409.

[5] 方向， 郝翠果， 劉俊田， 等. 新疆三塘湖盆地蘆草溝組混積巖型致密油形成條件與富集因素[J]. 地質通報， 2023， 42（2/3）： 397410.

Fang Xiang，Hao Cuiguo， Liu Juntian， et al. Geological Conditions and Enrichment Factors of the Lucaogou Formation Diamictite Tight Oil in Santanghu Basin， Xinjiang[J]. Geological Bulletin of China， 2023， 42（2/3）： 397410.

[6] 黃興， 高輝， 竇亮彬. 致密砂巖油藏微觀孔隙結構及水驅油特征[J]. 中國石油大學學報（自然科學版）， 2020， 44（1）： 8088.

Huang Xing，Gao Hui， Dou Liangbin. Micro Pore Structure and Water-Flooding Characteristics on Tight Sandstone Reservoir[J]. Journal of China University of Petroleum （Edition of Natural Science）， 2020， 44（1）： 8088.

[7] 汪成， 陳冬霞， 李洪博， 等. 珠江口盆地陸豐凹陷深部儲層原生孔隙特征及發育機制[J]. 吉林大學學報（地球科學版）， 2023， 53（1）： 4459.

Wang Cheng， Chen Dongxia， Li Hongbo， et al. Characteristics and Development Mechanism of Primary Pores in Deep Burial Reservoirs in Lufeng Sag， Pearl River Month Basin[J]. Journal of Jilin University （Earth Science Edition）， 2023， 53（1）： 4459.

[8] Swanson B F. A Simple Correlation Between Permeabilities and Mercury Capillary Pressures[J]. Journal of Petroleum Technology， 1981， 33： 24982504.

[9] Kolodzie S J.Analysis of Pore Throat Size and Use of the Waxman-Smits Equation to Determine OOIP in Spindle Field， Colorado[C]//Proceedings of the SPE Annual Technical Conference and Exhibition. Dallas： SPE， 1980： 2124.

[10] 曹英權， 曲希玉， 王清斌， 等. 渤中凹陷CFD64油田東營組儲層特征及物性控制因素[J]. 地質通報， 2024， 43（5）： 701711.

Cao Yingquan， Qu Xiyu， Wang Qingbin， et al. Reservoir Characteristics and Controlling Factors of Physical Properties of the Dongying Formation of the CFD64 Oilfield in the Bozhong Sag[J]. Geological Bulletin of China， 2024， 43（5）： 701711.

[11] Zheng S J， Yao Y B， Liu D M， et al. Characterizations of Full-Scale Pore Size Distribution， Porosity and Permeability of Coals： A Novel Methodology by Nuclear Magnetic Resonance and Fractal Analysis"""" Theory[J]. International Journal of Coal Geology， 2018， 196： 148158.

[12] Rezaee R， Saeedi A， Clennell B. Tight Gas Sands Permeability Estimation from Mercury Injection Capillary Pressure and Nuclear Magnetic Resonance Data[J]. Journal of Petroleum Science and Engineering， 2012， 88/89： 9299.

[13] Wang Y， Yang S C， Zhang Y F， et al. Investigation of Pore Structure and Reservoir Quality of Eocene Beach-Bar Sandstones in Dongying Depression， Bohai Bay Basin， East China[J]. Journal of Petroleum Science and Engineering， 2020， 189： 106854.

[14] 葉興樹， 王偉鋒， 戴俊生， 等. 東營凹陷沙三—東營期斷裂活動特征[J]. 中國石油大學學報（自然科學版）， 2006， 30（4）： 711.

Ye Xingshu， Wang Weifeng， Dai Junsheng， et al. Characteristics of Fault Activities of Sha3 Member and Dongying Periods in Dongying Depression[J]. Journal of China University of Petroleum （Edition of Natural Science）， 2006， 30（4）： 711.

[15] 操應長， 楊田， 宋明水， 等. 陸相斷陷湖盆低滲透碎屑巖儲層特征及相對優質儲層成因： 以濟陽坳陷東營凹陷古近系為例[J]. 石油學報， 2018， 39（7）： 727743.

Cao Yingchang， Yang Tian， Song Mingshui， et al. Characteristics of Low-Permeability Clastic Reservoirs and Genesis of Relatively High-Quality Reservoirs in the Continental Rift Lake Basin： A Case Study of Paleogene in the Dongying Sag， Jiyang Depression[J]. Acta Petrolei Sinica， 2018， 39（7）： 727743.

[16] 李愛芬， 任曉霞， 王桂娟， 等. 核磁共振研究致密砂巖孔隙結構的方法及應用[J]. 中國石油大學學報（自然科學版）， 2015， 39（6）： 9298.

Li Aifen， Ren Xiaoxia， Wang Guijuan， et al. Characterization of Pore Structure of Low Permeability Reservoirs Using a Nuclear Magnetic Resonance Method[J]. Journal of China University of Petroleum （Edition of Natural Science）， 2015， 39（6）： 9298.

[17] 劉可禹， 路研. 砂巖孔徑分布測定方法及其應用： CN202011156323.1[P]. 20220617.

Liu Keyu， Lu Yan. Measurement Method of Sandstone Pore-Throat Size Distribution and Their Application： CN202011156323.1[P]. 20220617.

[18] Mandelbrot B B， Passoja D E， Paullay A J. Fractal Character of Fracture Surfaces of Metals[J]. Nature， 1984， 308： 721722.

[19] 趙越， 李磊， 司運航， 等. 淺層頁巖氣儲層孔隙分形特征及控制因素： 以云南昭通地區龍馬溪組為例[J]. 吉林大學學報（地球科學版）， 2022， 52（6）： 18131829.

Zhao Yue， Li Lei， Si Yunhang， et al. Fractal Characteristics and Controlling Factors of Pores in Shallow Shale Gas Reservoirs： A Case Study of Longmaxi Formation in Zhaotong Area， Yunnan Province[J]. Journal of Jilin University （Earth Science Edition）， 2022， 52（6）： 18131829.

[20] 張全培， 王海紅， 劉美榮， 等. 超低滲透儲層全孔徑分布及其分形特征研究[J]. 中國礦業大學學報， 2020， 49（6）： 11371149.

Zhang Quanpei， Wang Haihong， Liu Meirong， et al. Study of the Full Pore Size Distribution and Fractal Characteristics of Ultra-Low Permeability Reservoir[J]. Journal of China University of Mining amp; Technology， 2020， 49（6）： 11371149.

[21] Lai J， Wang G W， Wang Z Y， et al. A Review on Pore Structure Characterization in Tight Sandstones[J]. Earth Science Reviews， 2018， 177： 436457.

[22] 吳蒙， 秦勇， 申建， 等. 致密砂巖儲層束縛水飽和度影響因素： 以鄂爾多斯盆地臨興地區為例[J]. 吉林大學學報（地球科學版）， 2022， 52（1）： 6879.

Wu Meng， Qin Yong， Shen Jian， et al. Influencing Factors of Irreducible Water Saturation in Tight Sandstone Reservoirs： A Case Study of Linxing Area in Ordos Basin[J]. Journal of Jilin University （Earth Science Edition）， 2022， 52（1）： 6879.

[23] 房濤， 張立寬， 劉乃貴， 等. 核磁共振技術定量表征致密砂巖氣儲層孔隙結構： 以臨清坳陷東部石炭系—二疊系致密砂巖為例[J]. 石油學報， 2017，38（8）：902915.

Fang Tao， Zhang Likuan， Liu Naigui， et al. Quantitative Characterization of Pore Structure of Tight Sandstone Reservoirs by NMR T2 Spectrum Technology： A Case Study of Carboniferous-Permian Tight Sandstone Reservoir in Linqing Depression[J]. Acta Petrolei Sinica， 2017， 38（8）： 902915.

[24] Testamanti M N， Rezaee R. Determination of NMR T2cut-off for Clay Bound Water in Shales： A Case Study of Carynginia Formation， Perth Basin， Western Australia[J]. Journal of Petroleum Science and Engineering， 2017， 149： 497503.