白云巖化作用與白云巖孔隙的形成

摘 要 【意義】白云巖是一種常見的碳酸鹽巖和重要的油氣儲集層，四川盆地海相碳酸鹽巖層系中96%的探明地質儲量來自于白云巖儲層。因此，白云巖化作用和白云巖孔隙成因的研究引起了極大的關注，并取得了顯著進展，然而有關白云巖化作用對孔隙的影響仍存在較大爭議。【方法】通過歸納總結白云巖孔隙成因的主要理論觀點，結合多種類型交代反應的實驗模擬結果，闡明白云巖化過程中孔隙的形成與演化，并對白云巖孔隙成因的研究提出了展望。【進展】已有理論表明，白云巖的孔隙可以繼承于原巖，也可形成于等摩爾交代、白云巖化過程中的方解石溶解或者白云巖形成之后的溶蝕作用；白云巖化作用可以增加、保持或者降低儲層孔隙，取決于交代反應的環境、時間和流體性質。【結論與展望】綜合實驗模擬結果的共性特征發現，交代反應是一個溶解—運移—沉淀的微過程；新孔隙的形成貫穿交代反應的整個過程，并提供了流體、離子交換的通道，以保持交代反應的持續進行；新孔隙的形成并非受礦物摩爾體積控制，而是受控于相對溶解度，即反應物的溶解量大于生成物的沉淀量。交代過程中產生的新孔隙以微孔隙為主，對儲層有效孔隙度的貢獻仍待進一步評價。

關鍵詞 白云巖化；孔隙演化；實驗模擬；交代作用

第一作者簡介 王廣偉，男，1986年出生，博士，副教授，儲層沉積學，E-mail： wanggw@upc.edu.cn

中圖分類號 P618.13 文獻標志碼 A

0 引言

白云巖化作用是白云石礦物逐漸取代方解石或文石礦物的一種交代反應，微觀上表現為方解石、文石等反應物溶解、同時伴隨著白云石沉淀的耦合過程[1]。白云巖在全球范圍內廣泛分布[2?6]，常發育不同類型的儲集空間而成為優質的油氣儲集層[7?11]。據統計，四川盆地海相碳酸鹽巖層系已獲探明地質儲量2.34×1012 m3天然氣[12]，其中白云巖儲層中的探明儲量占總探明儲量的96.2%。因此，有關白云巖孔隙成因的研究引起了極大的關注，并取得了顯著的進展[13?19]。白云巖的孔隙可以形成于沉積期、大氣淡水暴露溶蝕[20?22]、白云巖化作用[23?24]或者埋藏溶蝕作用[25?28]等。其中，白云巖化作用對白云巖孔隙的影響仍存在較大爭議。基于孔隙與白云巖形成的相對時間，白云巖的孔隙可劃分為同生孔隙和后生孔隙[23]；同生孔隙形成時間與白云巖形成時間一致，是繼承/轉化的原巖孔隙，抑或是白云巖化過程中新形成的孔隙；后生孔隙則是白云巖形成之后形成的所有孔隙。

國內外學者對白云巖化作用本身對孔隙的影響（即同生孔隙的成因）已作出較多理論和實例分析，可總結為三個觀點[23，29]：（1）白云巖化作用可以增加儲層的孔隙度，方解石、文石轉變為白云石的過程中，因摩爾體積的差異，生成物沉淀的總體積小于反應物消耗的總體積，導致交代體系中固相體積收縮和孔隙體積增大；（2）白云巖化作用對儲層的孔隙度沒有影響，白云巖的孔隙主要繼承于灰巖原巖的孔隙；（3）白云巖化作用降低了儲層的物性，如過白云巖化作用。因此，白云巖化過程中孔隙的增加、保持或破壞似乎與交代反應的環境、流體性質和反應時間密切相關[29]。地質環境下的交代作用類型豐富，是一類流體參與的溶解—沉淀相耦合的微過程[30]。本文歸納總結了白云巖化過程中孔隙形成與演化的主要理論觀點，并列舉了四類典型的交代反應實驗模擬，通過總結和對比實驗的觀測結果，識別出交代反應的一些共性特征，為理解白云巖化作用對孔隙的影響提供不同的視角和理論借鑒。

1 白云巖孔隙成因研究現狀

碳酸鹽巖層系的油氣勘探與實踐表明，白云巖物性常優于同時期灰巖，尤其在古老地層中更為明顯[2，6]。針對這種現象的成因，前人做了大量地質實例分析和實驗研究，以下對白云巖孔隙成因的解釋和理論模式進行系統的梳理和總結。

1.1 白云巖化過程中的摩爾體積變化

受化學動力學障礙的限制，白云石極難從地表常溫、常壓環境的水體中直接沉淀[30?32]。規模性發育的白云巖主要為次生交代成因，反應過程可用公式（1）表示[33]：

2CaCO3+Mg2+→CaMg（CO3） 2+Ca2+ （1）

該過程會消耗白云巖化流體中1 mol的Mg2+同時釋放1 mol的Ca2+至流體中，因此稱為等摩爾交代作用；同時，生成1 mol白云石會消耗2 mol方解石或文石。已知白云石的摩爾體積為64.365 cm3/mol，方解石的摩爾體積為36.934 cm3/mol，文石的摩爾體積為34.150 cm3/mol[34]。由于反應物與生成物摩爾體積的差異，反應前后體系中的固體體積會發生變化，孔隙體積也相應的發生改變。理論計算表明，若白云石交代方解石，則會造成礦物總體積縮小（36.934 cm3/mol×2 mol）-（64.365 cm3/mol×1 mol）=9.503 cm3，相應的 孔 隙 度 增 加 值 為 9.503 cm3（/ 36.934 cm3/mol×2 mol）×100%=12.86%；若白云石交代文石，固相體積會縮小（34.150 cm3/mol×2 mol）-（64.365 cm3/mol×1 mol）=3. 935 cm3，相應的孔隙度增加5.76%。由此，前人認為白云巖化過程中由于摩爾體積的變化會形成新的孔隙[33]，這便是經典的白云巖孔隙成因理論，并用于解釋碳酸鹽巖地層中白云巖孔隙度一般好于同時期灰巖的地質現象。等摩爾交代理論自提出后一直備受質疑，關鍵在于白云巖化過程中是否有外源CO23 -的混入。Weyl[35]研究認為，相對于Ca2+、Mg2+濃度，大多數地層水中CO23 -含量較低，并利用質量守恒原理證明白云巖化過程遵循等摩爾交代原理。基于對意大利Venetian Alps 地區中侏羅統的研究，Zempolich et al.[36]發現白云巖的面孔率多介于10%～15%，與理論增加值相似，因此認為白云巖化作用是一個等摩爾交代的過程。

然而，Lucia et al.[37]研究指出，全新世白云巖孔隙度與灰巖相似，甚至在波內爾島上新世—更新世碳酸鹽巖地層的白云巖孔隙度顯著低于同期灰巖；此外，鮞粒、生屑等顆粒的原始微結構在白云巖化之后被完美地保存下來[38]，這些均表明白云巖化作用能夠保持固相的原始格架，外源CO23 -的流入應該是一種普遍現象，白云巖化作用并非是等摩爾交代而是等體積交代的過程，可以用公式（2）和公式（3）表達：

1.743CaCO3（方解石）+Mg2++0.257CO23 -→CaMg（CO3）2 +0.743Ca2+ （2）

或

1.885CaCO3（ 文 石 ）+Mg2++0.115CO23 -→CaMg（CO3）2 +0.115Ca2+ （3）

交代反應體系中，生成1 mol 白云石僅消耗1.743 mol的方解石或1.855 mol的文石，固相體積的損失量與白云石的沉淀量相同，反應前后孔隙體積保持不變。等體積交代理論認為，白云巖化過程中沒有新孔隙的生成，白云巖的孔隙只是白云巖化過程中灰巖原巖孔隙的轉變，為繼承性孔隙；因此，白云巖的孔隙度取決于原巖。古代白云巖孔隙度大于相鄰灰巖的主要原因是白云巖具有更強的抗壓實—壓溶能力，在埋藏環境中能夠保存更多的孔隙[39]。四川盆地三疊系潮坪相白云巖形成于準同生—早成巖期，在近于等體積交代過程中，灰泥微孔隙轉變為白云巖晶間孔并在埋藏環境中保存下來，而未發生白云巖化的灰泥則被壓實成致密的泥晶灰巖[38]。Lucia[39]通過對比不同時代（全新世、更新世、新近紀和侏羅紀）灰巖與白云巖的孔隙特征指出，白云巖孔隙的形成與白云巖化作用本身無關，白云巖孔隙度即可以大于、等于或小于同時代的灰巖，因此認為白云巖孔隙為繼承性孔隙。

綜上所述，白云巖化流體在提供Mg2+離子的同時也提供了CO23 -離子，導致白云巖化反應路徑的差異，從而引起孔隙體積的變化，主要取決于白云巖化流體的CO23 -濃度。因此，根據流體提供CO23 -的含量，白云巖化過程可以用綜合反應模式（公式4）表示為[40]：

（2-x）CaCO3+Mg2++xCO23 -→CaMg（CO3）2+（1-x）Ca2+ （4）

x為白云巖化過程實際消耗流體中CO23 -的濃度，對于方解石和文石，當x 分別取值為0.257 和0.115時，反應前后礦物體積保持不變。以白云石交代方解石為例，當x等于0，即白云巖化發生在封閉的成巖體系，無外源CO23 -的混入，則白云巖化遵循等摩爾交代作用（公式1），孔隙度理論增加值為12.86%；當00.257時，則白云巖化后固相體積增加，存在過白云巖作用，孔隙度降低。白云巖化流體的性質最終決定了白云巖化作用是增加、降低還是保持孔隙。因此，在近地表環境下，白云巖化流體常富含外源CO23 -，交代反應傾向于等體積交代為主，在靠近鹵水源的位置當方解石被完全交代后可能發生過白云巖化；相比，在相對封閉的埋藏環境，云化流體的外源CO23 -濃度較低，白云巖化可能更傾向于等摩爾交代。

1.2 白云巖化過程中的溶解作用

白云巖化作用是碳酸鈣溶解—白云石再沉淀的微觀耦合過程。白云巖化流體主要為不同程度的蒸發海水[41?42]，并且隨著蒸發作用的增強，海水鹽度逐漸增加、pH值逐漸降低[43]。因此，蒸發性的海源流體作為白云巖化流體，對灰巖具有一定的溶解能力，且與灰巖地層的原始孔隙水混合，會進一步促進原巖的溶解[44?45]。

另一種孔隙成因理論為白云巖化晚期殘余方解石的溶蝕。在白云巖化過程中的晚期，地層流體受白云巖化反應的改造，可能對殘留的灰質成分有較強的溶解能力，導致殘余方解石的溶解速率大于白云石的沉淀速率[46]，從而形成孔隙。時間上，孔隙形成于白云巖化過程的晚期，屬于白云巖化的一部分。這一理論首先由Landes[47]提出，并用于解釋密歇根盆地局部高孔隙性白云巖儲層成因。Wendte[48]系統研究了不同白云巖化程度的儲層特征，發現白云巖的晶間孔實際為原巖的粒內孔和鑄模孔，且在厚層白云巖中部最發育，認為白云石首先交代原始灰巖顆粒的外部，內部殘存的方解石或文石在白云巖化晚期被溶蝕形成粒內孔。

2 交代反應的實驗模擬和結果

2.1 白云石[CaMg（CO3）2]和菱鎂礦（MgCO3）交代單晶方解石（CaCO3）

將邊長約為2 mm的立方體方解石單晶和1 mL濃度為1 mol/L的MgCl2溶液放置于高壓反應釜中，在200 ℃、飽和蒸氣壓（約16 bar）條件下反應。共計5組平行試驗，反應時間分別為1 d、3 d、7 d、14 d 和28 d，分別對每組試驗的產物進行觀察[49]，觀測結果如圖1所示。

交代反應顯示出方解石單晶被白云石和菱鎂礦自外向內逐漸交代（圖1a）。反應初期（1 d 和3d），方解石晶體外層首先被菱鎂礦交代，形成一個反應邊（reaction rim）；隨著反應的進行（14 d），白云石出現在方解石和菱鎂礦層之間（圖1b，c）；至反應的最后階段（28 d），生成的白云石又被菱鎂礦交代，反應邊幾乎全由菱鎂礦組成。整體而言，交代反應過程中，原始晶體的外形和大小未發生明顯變化（圖1a）。原始方解石單晶幾乎不發育孔隙，被白云石和菱鎂礦交代后，形成一個孔隙性的反應邊，發育大量的晶間微孔（圖1d）；反應邊與未反應的方解石之間存在明顯的間隙（圖1c），并隨交代時間的增加，反應邊和間隙的寬度均增大。反應后溶液中Ca離子濃度增加，Mg 離子濃度降低，Ca/Mg 值增大（表1）。

2.2 白云石（CaMg（CO3）2）交代文石鮞粒（CaCO3）

將50 mg 文石質鮞粒放入10 mL 含0.5 M 的MgCl2和0.4 M的CaCl2的溶液中，反應溫度為200 ℃、平衡蒸氣壓[50]。共計9組平行實驗，反應時長分別為6 h、24 h、36 h、48 h、60 h、72 h、96 h、168 h和180 h，產物過濾并用去離子水沖洗，干燥后用于觀測。

模擬結果顯示（圖2），白云巖化作用自鮞粒外層向內部逐漸延伸。反應6 h后，鮞粒表層可見零星分布的半自形白云石晶體（圖2a）。反應至24 h，白云石晶體趨于自形，且逐漸連片形成白云石包殼（圖2b）。隨著反應時間的增加（24～96 h），鮞粒內部圈層開始發育半自形—自形白云石（圖2c，d），白云石首先交代藻鉆孔和孔隙性紋層，然后再交代致密的部分；生成的微晶白云石顯示出與原始文石質鮞粒相似的紋層結構（圖2e，f），紋層厚度一般為數十微米。總體而言，交代反應形成的白云石主要沉淀在原始鮞粒輪廓范圍內，反應前后總體積不變，孔隙主要包括鮞粒紋層中的晶間微孔隙和少量的鑄模孔。白云石包殼和紋層的產生表明，大部分白云石沉淀發生在文石溶解位置附近的幾十微米范圍內，而少量鑄模孔的發育則表明，一些沉淀則發生在溶解位置數百微米范圍內[50]，反應物溶解的位置與生成物沉淀的位置存在一定距離的偏移。

2.3 方沸石（NaAlSi2O6·H2O） 交代白榴石（KAl?Si2O6）

在150 ℃溫度條件下，將白榴石晶體置于含3.5%的NaCl溶液的高壓釜中進行方沸石交代白榴石的實驗[51?52]。交代反應式如公式（5）所示：

KAlSi2O6+Na++H2O→NaAlSi2O6·H2O+K+ （5）

該過程消耗流體中1 mol的Na+同時釋放1 mol的K+至流體中，類似于白云巖化中的等摩爾交代作用。同時，生成1 mol 方沸石會消耗1 mol 白榴石。反應物白榴石的摩爾體積為89.302 cm3/mol，生成物方沸石的摩爾體積為97.574 cm3/mol[34]，生成物的摩爾體積大于反應物。與白云巖化作用的摩爾體積減少的情況相反，理論計算表明，1 mol 方沸石交代1 mol白榴石后，體系的固相體積應該增加9.26%，孔隙度相應的降低9.26%。

然而，實驗結果顯示，反應物白榴石為雙晶結構，孔隙不發育（圖3a～c），交代產物方沸石沿著白榴石晶體邊界向內部交代（圖3d）。掃描電鏡觀察顯示，方沸石由有序的微晶組成（圖3e），層間發育大量的微孔隙（圖3f）。這與傳統理論認為摩爾體積的增加降低孔隙的觀點相矛盾。由此表明，生成物的摩爾體積增加不會影響交代過程中新孔隙的形成，反應物中存在微孔隙似乎是一種普遍現象。新孔隙的形成表明交代過程中方沸石的沉淀量小于白榴石的溶解量，并超過摩爾體積增加的影響。

2.4 氯化鉀（KCl）交代溴化鉀（KBr）

室溫條件下，將大小為2 mm×2 mm×1 mm 的溴化鉀晶體放入0.25 mL的KCl飽和溶液中，KBr在溶液中發生溶解，增加了溶液中K+離子的含量，從而導致KCl過飽和而發生沉淀，即宏觀上表現的KCl交代KBr晶體，反應時長從1 min至12 d[53]。交代反應過程如公式（6）所示：

KBr+Cl-→KCl+Br- （6）

該過程消耗流體中1 mol的Cl-同時釋放1 mol的Br-至流體中，同樣類似于白云巖化的等摩爾交代作用。同時，生成1 mol 氯化鉀會消耗1 mol 溴化鉀。反應物溴化鉀的摩爾體積為43.273 cm3/mol，生成物氯化鉀的摩爾體積為37.519 cm3/mol[34]，生成物的摩爾體積小于反應物。理論計算表明，1 mol氯化鉀交代1 mol 溴化鉀后，體系的固相體積應該減少13.30%。

實驗結果顯示，反應物溴化鉀晶體干凈透明，孔隙不發育（圖4a），隨著交代反應由外向內進行，乳白色的氯化鉀生成物逐漸增大（圖4b，c），反應時間至120 min后完全形成氯化鉀晶體（圖4d）。掃描電鏡觀察顯示，原始溴化鉀的孔隙度幾乎為零（圖4e），而由氯化鉀組成的反應邊發育大量的微孔隙（圖4e～g）。然而，由于微晶和微孔隙具有較大的表面能，形成的氯化鉀晶體的微結構仍會進一步演化和再平衡，導致晶體外邊緣微孔隙消失，形成一個顯著的亮邊結構（圖4h）。隨著平衡時間的增加，亮邊結構更趨明顯（圖4i），內部微孔隙完全被亮邊所包裹。這種結構與實際的白云石[54]地質樣品中的霧心—亮邊特征具有很好的相似性（圖4j）。

3 交代作用與孔隙的演化

交代反應存在反應物的溶解與生成物的沉淀，溶解量與沉淀量的相對變化必然引起固體格架的變化以及孔隙形態和大小的改變。基于不同類型實驗模擬的觀測結果，識別出交代反應的一些共性特征，為理解白云巖化過程中孔隙的演化提供借鑒。

綜合四組實驗模擬結果表明，交代反應基本保持了被交代礦物或顆粒的外部輪廓，生成物普遍發育微孔隙；微孔隙提供了流體和離子交換的通道，從而保持交代反應自外向內地持續進行。在三組交代實驗中，反應物方解石、溴化鉀和白榴石三種單晶礦物均不發育孔隙，即晶體的原始孔隙度幾乎為0，但是在生成物中均可觀察到大量的晶間微孔隙，表明交代反應確實可以形成新的孔隙。氯化鉀交代溴化鉀、白云石和菱鎂礦交代方解石、白云石交代文石鮞粒這三組反應，生成物的摩爾體積均小于反應物，似乎驗證了等摩爾交代控制了孔隙的形成。然而，方沸石交代白榴石后，生成物的摩爾體積大于反應物白榴石，但是在方沸石反應邊中仍可見到大量的微孔隙（圖3e，f），與理論計算的固相體積增加相矛盾。由此表明，交代過程中微孔隙的形成并非受礦物摩爾體積的控制，而是由于反應物的溶解量大于生成物的沉淀量造成固相體積的虧損。總體而言，交代過程中反應物首先發生溶解（如KBr溶解），向流體中釋放各種離子（如Br-和K+），造成流體對另一種礦物（如KCl）過飽和；流體通過孔隙性反應邊向反應前緣提供交代所需的離子（如Cl-），同時帶出反應后釋放的離子（如Br+）；交代反應過程中反應物的溶解量總是大于生成物的沉淀量，造成孔隙體積的增加。

此外，在白云石和菱鎂礦交代單晶方解石的實驗中，除了生成物中大量發育的微孔隙，反應物和生成物之間也發育明顯的間隙（gap），表明反應物溶解位置和生成物沉淀位置存在一定距離的偏離，溶解的物質進入流體并發生一定距離的運移后才發生沉淀。因此，白云巖化作用在微觀上可以看作一種反應物的溶解、物質運移至成核點、生成物沉淀三個耦合的微過程，即溶解—運移—沉淀的微過程（圖4c）。白云巖化作用中反應物在流體—礦物界面處發生微觀溶解（如鹵水與原始地層水混合造成對方解石、文石等被交代物的不飽和[44]），向邊界流體中釋放Ca2+和CO23 -離子，造成邊界流體對白云石礦物的過飽和；運移則涉及流體通過孔隙性反應邊向反應前緣提供交代所需的Mg2+離子，同時帶出反應后釋放的Ca2+離子。交代反應過程總體顯示出溶解量大于沉淀量、溶解位置與沉淀位置存在一定偏離的特征。從溶解—運移—沉淀三個微過程的角度分析交代過程及其產物特征，可以對白云巖化的一些現象作出較為合理的解釋。運移距離的大小影響了原始微結構的保存程度。只有運移距離小，溶解—沉淀相耦合的交代反應才能完美保存原始顆粒結構特征，如鮞粒內部的圈層結構。當運移距離較大，即物質溶解的部位和沉淀的部位距離較遠時，則溶解的部位可能演變為較大的孔隙，如實驗中觀察到的間隙和鑄模孔，同時原始結構保存得較差。運移距離可能受成巖流體成分、流速和pH值控制[51]。地質歷史中海水的Mg/Ca值、鹽度和白云巖化流體性質可能是影響白云巖化過程中原始顆粒保存程度和孔隙形成的重要因素[55?56]。

然而，新形成的交代產物由大量的微晶和微孔隙組成，具有較大的表面能。在表面能的驅動下，生成物會發生結構平衡和重結晶，導致晶體周緣微孔隙的消失（圖4h），形成亮邊結構（圖4i）。同時晶體內部的微孔隙被亮邊所封閉而形成孤立的孔隙和流體包裹體（圖4h）。由此形成的霧心—亮邊結構與地質樣品中的白云石結構非常相似（圖4j）。交代作用過程中形成的一部分微孔隙可能為無效孔隙，對儲集層有效孔隙度的貢獻有限。

綜上所述，不管交代產物摩爾體積是減小還是增大，交代過程中均存在新孔隙的形成，以保障交代反應的持續進行，礦物摩爾體積的減少不是交代過程中孔隙形成的主要控制因素。因此，可以推斷，白云巖化作用中孔隙的形成貫穿整個交代過程，溶解量大于沉淀量是孔隙形成的主要原因，新形成的孔隙可能由于結構平衡和重結晶而成為無效孔隙；運移距離則決定了白云巖化后原始結構的保存程度。

4 存在的問題與發展趨勢

現今普遍認為，白云巖的成因類型多樣[57]，然而規模性白云巖主要為近地表環境下的次生交代成因[42]。雖然地質環境中的交代作用類型眾多[58]，存在共性特征和各自的特殊性，但是均表現出從顆粒邊緣向內部逐漸反應的過程，形成一個孔隙性反應邊。如何有效定量表征反應邊的孔隙含量，從而定量比較交代反應前后孔隙的變化，有助于深入理解交代反應過程中的孔隙演化。本文列舉的四組交代反應的實驗模擬，除了氯化鉀交代溴化鉀的實驗為常溫、常壓條件，其他三組均是高溫、短時間的模擬，與地質歷史時期白云巖化的環境存在顯著差異。此外，海水作為白云巖化流體的主要來源，海水成分在地質歷史中持續演變[56]。如何有效改進白云巖化實驗模擬的溫壓條件，使之更加符合白云巖形成的自然地質條件，成為深入理解白云巖孔隙成因和演化的關鍵。隨著實驗技術的提高，在考慮白云巖化流體性質的條件下，可以開展相對低溫—長時間的白云巖化實驗模擬，如Usdowski[59]開展了反應溫度介于60 ℃～90 ℃、反應時間長達7年的白云巖化實驗模擬。近二十年來，數值模擬技術逐漸應用于白云巖化過程中礦物和孔隙演化的動態定量分析[60?63]，有效彌補實驗模擬在時間尺度上的不足。相對低溫—長時間的實驗模擬結果可以為白云巖化數值模擬提供約束條件，形成實驗模擬和數值模擬相結合的方法，用于綜合分析白云巖化作用過程及其孔隙的形成與演化。

5 結論

（1） 綜合四種類型的交代反應模擬實驗表明，無論生成物礦物摩爾體積的減少或增加，交代反應過程中均存在微孔隙的形成，微孔隙提供了流體和離子的交換通道，保持了交代反應的持續進行。

（2） 白云巖化作用是溶解—運移—沉淀三個耦合的微過程。白云巖孔隙的形成并非等摩爾交代相關的礦物體積收縮，而是取決于反應物與生產物的相對溶解度，即溶解量大于沉淀量。交代過程中形成的孔隙主要為晶間微孔，隨著晶體結構的平衡和重結晶，部分微孔隙轉變為無效孔隙。

（3） 實際地質條件下，由于白云巖化流體性質、被交代礦物類型和反應時間的不同，形成的白云巖孔隙和結構存在差異。隨著實驗模擬和分析測試技術的提高，有待開展更多不同條件下白云巖化的實驗模擬，結合數值模擬方法，深化白云巖化過程和孔隙形成演化的認識。

致謝 感謝評審專家對本文提出的建設性修改意見，使得稿件質量得以提升。中國石油大學（華東） 郝芳教授對本文的撰寫提供了諸多指導和幫助。

參考文獻（References）

［1］ Warren J. Dolomite： Occurrence， evolution and economically important

associations[J]. Earth-Science Reviews， 2000， 52（1/2/3）：

1-81.

［2］ 陳代釗，錢一雄. 深層—超深層白云巖儲集層：機遇與挑戰[J].

古地理學報，2017，19（2）：187-196. ［Chen Daizhao， Qian

Yixiong. Deep or super-deep dolostone reservoirs： Opportunities

and challenges[J]. Journal of Palaeogeography， 2017， 19（2）：

187-196.］

［3］ 文華國，霍飛，郭佩，等. 白云巖—蒸發巖共生體系研究進展及

展望[J]. 沉積學報，2021，39（6）：1319-1343.［Wen Huaguo， Huo

Fei， Guo Pei， et al. Advances and prospects of dolostoneevaporite

paragenesis system[J]. Acta Sedimentologica Sinica，

2021， 39（6）： 1319-1343.］

［4］ 何治亮，馬永生，張軍濤，等. 中國的白云巖與白云巖儲層：分

布、成因與控制因素[J]. 石油與天然氣地質，2020，41（1）：1-14.

［He Zhiliang， Ma Yongsheng， Zhang Juntao， et al. Distribution，

genetic mechanism and control factors of dolomite and dolomite

reservoirs in China[J]. Oil amp; Gas Geology， 2020， 41（1）： 1-14.］

［5］ 馬鋒，楊柳明，顧家裕，等. 世界白云巖油氣田勘探綜述[J]. 沉積學

報，2011，29（5）：1010-1021.［Ma Feng， Yang Liuming， Gu Jiayu， et

al. The summary on exploration of the dolomite oilfields in the

world[J]. Acta Sedimentologica Sinica， 2011， 29（5）： 1010-1021.］

［6］ Allan J R， Wiggins W D. Dolomite reservoirs： Geochemical techniques

for evaluating origin and distribution[J]. American

Association of Petroleum Geologists， 1993， 36： 1-129.

［7］ Ma Y S， Zhang S C， Guo T L， et al. Petroleum geology of the Puguang

sour gas field in the Sichuan Basin， SW China[J]. Marine

and Petroleum Geology， 2008， 25（4/5）： 357-370.

［8］ 馬永生，蔡勛育，趙培榮，等. 四川盆地大中型天然氣田分布特

征與勘探方向[J]. 石油學報，2010，31（3）：347-354.［Ma Yongsheng，

Cai Xunyu， Zhao Peirong， et al. Distribution and further

exploration of the large-medium sized gas fields in Sichuan Basin

[J]. Acta Petrolei Sinica， 2010， 31（3）： 347-354.］

［9］ 蘆飛凡，譚秀成，王利超，等. 川中地區中二疊統棲霞組灘控巖溶

型白云巖儲層特征及主控因素[J]. 沉積學報，2021，39（2）：456-

469.［Lu Feifan， Tan Xiucheng， Wang Lichao， et al. Characteristics

and controlling factors of dolomite reservoirs within shoalcontrolled

karst in the Middle Permian Qixia Formation， central

Sichuan Basin[J]. Acta Sedimentologica Sinica， 2021， 39（2）：

456-469.］

［10］ 田興旺，楊岱林，鐘佳倚，等. 基于CT成像技術的白云巖儲層

微觀表征：以川中磨溪—龍女寺臺內地區震旦系燈影組四段

為例[J]. 沉積學報，2021，39（5）：1264-1274.［Tian Xingwang，

Yang Dailin， Zhong Jiayi， et al. Microscopic characterization of

dolomite reservoirs by CT imaging： A case study of the Dengsi

Formation in Moxi-Longnvsi area， central Sichuan[J]. Acta Sedimentologica

Sinica， 2021， 39（5）： 1264-1274.］

［11］ 王勇，施澤進，孟興平，等. 川東南龍王廟組埋藏及混合水白

云巖化作用[J]. 沉積學報，2021，39（6）：1517-1531.［Wang

Yong， Shi Zejin， Meng Xingping， et al. Burial dolomitization

and mixed water dolomitization in Longwangmiao Formation，

southeastern Sichuan Basin[J]. Acta Sedimentologica Sinica，

2021， 39（6）： 1517-1531.］

［12］ 馬新華，楊雨，文龍，等. 四川盆地海相碳酸鹽巖大中型氣田

分布規律及勘探方向[J]. 石油勘探與開發，2019，46（1）：1-13.

［Ma Xinhua， Yang Yu， Wen Long， et al. Distribution and exploration

direction of medium-and large-sized marine carbonate gas

fields in Sichuan Basin， SW China[J]. Petroleum Exploration

and Development， 2019， 46（1）： 1-13.］

［13］ 馬永生，蔡勛育，趙培榮，等. 深層超深層碳酸鹽巖優質儲層

發育機理和“三元控儲”模式：以四川普光氣田為例[J]. 地質學

報，2010，84（8）：1087-1094.［Ma Yongsheng， Cai Xunyu， Zhao

Peirong， et al. Formation mechanism of deep-buried carbonate

reservoir and its model of three-element controlling reservoir： A

case study from the Puguang oilfield in Sichuan[J]. Acta Geologica

Sinica， 2010， 84（8）： 1087-1094.］

［14］ 趙文智，沈安江，鄭劍鋒，等. 塔里木、四川及鄂爾多斯盆地白

云巖儲層孔隙成因探討及對儲層預測的指導意義[J]. 中國科

學（D 輯）：地球科學，2014，44（9）：1925-1939.［Zhao Wenzhi，

Shen Anjiang， Zheng Jianfeng， et al. The porosity origin of

dolostone reservoirs in the Tarim， Sichuan and Ordos Basins and its

implication to reservoir prediction[J]. Science China （Seri. D）：

Earth Sciences， 2014， 44（9）： 1925-1939.］

［15］ 沈安江，趙文智，胡安平，等. 海相碳酸鹽巖儲集層發育主控

因素[J]. 石油勘探與開發，2015，42（5）：545-554.［Shen Anjiang，

Zhao Wenzhi， Hu Anping， et al. Major factors controlling

the development of marine carbonate reservoirs[J]. Petroleum

Exploration and Development， 2015， 42（5）： 545-554.］

［16］ Purser B H， Brown A， Aissaoui D M. Nature， origins and evolution

of porosity in dolomites [M]//Purser B， Tucker M， Zenger

D. Dolomites： A volume in honour of dolomieu. Oxford： Blackwell，

1994： 283-308.

［17］ Hao F， Zhang X F， Wang C W， et al. The fate of CO2 derived

from thermochemical sulfate reduction （TSR） and effect of TSR

on carbonate porosity and permeability， Sichuan Basin， China

[J]. Earth-Science Reviews， 2015， 141： 154-177.

［18］ 馬永生，蔡勛育，趙培榮. 深層、超深層碳酸鹽巖油氣儲層形

成機理研究綜述[J]. 地學前緣，2011，18（4）：181-192.［Ma

Yongsheng， Cai Xunyu， Zhao Peirong. The research status and

advances in porosity evolution and diagenesis of deep carbonate

reservoir[J]. Earth Science Frontiers， 2011， 18（4）： 181-192.］

［19］ 喬占峰，張哨楠，沈安江，等. 塔里木和四川盆地白云巖規模

優質儲層形成與發育控制因素[J]. 石油與天然氣地質，2022，

43（1）：92-104.［Qiao Zhanfeng， Zhang Shaonan， Shen Anjiang，

et al. Controls on formation and development of large-sized

high-quality dolomite reservoirs in the Tarim and Sichuan Basins

[J]. Oil amp; Gas Geology， 2022， 43（1）： 92-104.］

［20］ 潘立銀，劉占國，李昌，等. 四川盆地東部下三疊統飛仙關組白云

巖化作用及其與儲集層發育的關系[J]. 古地理學報，2012，14

（2）：176-186. ［Pan Liyin， Liu Zhanguo， Li Chang， et al.

Dolomitization and its relationship with reservoir development of

the Lower Triassic Feixianguan Formation in eastern Sichuan

Basin[J]. Journal of Palaeogeography， 2012， 14（2）： 176-186.］

［21］ 潘立銀，郝毅，梁峰，等. 白云巖儲層成因的激光原位U-Pb 定

年和同位素地球化學新證據：以四川盆地西北部中二疊統棲

霞組白云巖儲層為例[J]. 石油學報，2022，43（2）：223-233.

［Pan Liyin， Hao Yi， Liang Feng， et al. New evidence of laser insitu

U-Pb dating and isotopic geochemistry for the genesis of dolomite

reservoir： A case study of dolomite reservoir from Middle

Permian Qixia Formation in northwestern Sichuan Basin[J]. Acta

Petrolei Sinica， 2022， 43（2）： 223-233.］

［22］ 李明隆，譚秀成，羅冰，等. 四川盆地西北部中二疊統棲霞組

相控早期高頻暴露巖溶特征及啟示[J]. 中國石油勘探，2020，

25（3）：66-82.［Li Minglong， Tan Xiucheng， Luo Bing， et al.

Characteristics of facies-controlled and early high-frequency

exposed karstificationin the Qixia Formation of Middle Permian

in the northwest of Sichuan Basin and its significance[J]. China

Petroleum Exploration， 2020， 25（3）： 66-82.］

［23］ Wang G W， Li P P， Hao F， et al. Dolomitization process and its

implications for porosity development in dolostones： A case

study from the Lower Triassic Feixianguan Formation， Jiannan

area， eastern Sichuan Basin， China[J]. Journal of Petroleum

Science and Engineering， 2015， 131： 184-199.

［24］ Maliva R G， Budd D A， Clayton E A， et al. Insights into the

dolomitization process and porosity modification in sucrosic

dolostones， Avon Park Formation （Middle Eocene）， east-central

Florida， U. S. A. [J]. Journal of Sedimentary Research， 2011， 81

（3）： 218-232.

［25］ 朱光有，張水昌，梁英波，等. TSR對深部碳酸鹽巖儲層的溶蝕

改造：四川盆地深部碳酸鹽巖優質儲層形成的重要方式[J]. 巖

石學報，2006，22（8）：2182-2194. ［Zhu Guangyou， Zhang

Shuichang， Liang Yingbo， et al. Dissolution and alteration of the

deep carbonate reservoirs by TSR： An important type of deepburied

high-quality carbonate reservoirs in Sichuan Basin[J]. Acta

Petrologica Sinica， 2006， 22（8）： 2182-2194.］

［26］ 朱光有，張水昌，梁英波. 四川盆地深部海相優質儲集層的形

成機理及其分布預測[J]. 石油勘探與開發，2006，33（2）：161-

166.［Zhu Guangyou， Zhang Shuichang， Liang Yingbo. Formation

mechanism and distribution prediction of high-quality marine

reservoir in deeper Sichuan Basin[J]. Petroleum Exploration

and Development， 2006， 33（2）： 161-166.］

［27］ Zhu D Y， Meng Q Q， Jin Z J， et al. Formation mechanism of

deep Cambrian dolomite reservoirs in the Tarim Basin， northwestern

China[J]. Marine and Petroleum Geology， 2015， 59：

232-244.

［28］ 王興志，張帆，蔣志斌，等. 四川盆地東北部飛仙關組儲層研

究[J]. 地學前緣，2008，15（1）：117-122.［Wang Xingzhi， Zhang

Fan， Jiang Zhibin， et al. A study of Feixianguan reservoir in

northeast Sichuan Basin[J]. Earth Science Frontiers， 2008， 15

（1）： 117-122.］

［29］ 張學豐，劉波，蔡忠賢，等. 白云巖化作用與碳酸鹽巖儲層物

性[J]. 地質科技情報，2010，29（3）：79-85.［Zhang Xuefeng，

Liu Bo， Cai Zhongxian， et al. Dolomitization and carbonate

reservoir formation[J]. Geological Science and Technology

Information， 2010， 29（3）： 79-85.］

［30］ Lippmann F. Sedimentary carbonate minerals[M]. New York：

Springer， 1973： 1-228.

［31］ Land L S. Failure to precipitate dolomite at 25℃ from dilute solution

despite 1000-fold oversaturation after 32 years[J]. Aquatic

Geochemistry， 1998， 4（3/4）： 361-368.

［32］ Moore C H. Carbonate reservoirs： Porosity evolution and diagenesis

in a sequence stratigraphic framework[M]. New York：

Elsevier， 2001： 1-340.

［33］ Murray R C. Origin of porosity in carbonate rocks[J]. Journal

of Sedimentary Research， 1960， 30（1）： 59-84.

［34］ Blanc P， Lassin A， Piantone P， et al. Thermoddem： A geochemical

database focused on low temperature water/rock interactions

and waste materials[J]. Applied Geochemistry， 2012， 27（10）：

2107-2116.

［35］ Weyl P K. Porosity through dolomitization： Conservation-ofmass

requirements[J]. Journal of Sedimentary Research， 1960，

30（1）： 85-90.

［36］ Zempolich W G， Hardie L A. Geometry of dolomite bodies within

deep-water resedimented oolite of the Middle Jurassic Vajont limestone，

Venetian Alps， Italy： Analogs for hydrocarbon reservoirs

created through fault-related burial dolomitization[M]//Kupecz J

A， Gluyas J G， Bloch S. Reservoir quality prediction in sandstones

and carbonates. Tulsa， Okla： AAPG Memoirs， 1997： 127-162.

［37］ Lucia F J， Major R P. Porosity evolution through hypersaline reflux

dolomitization//Purser B， Tucker M， Zenger D. Dolomites：

A volume in honour of dolomieu. Oxford： Blackwell， 1994：

325-341.

［38］ 陳彥華，劉鶯，孫妥. 白云化過程中巖石孔隙體積的變化[J].

石油實驗地質，1985，7（1）：29-37.［Chen Yanhua， Liu Ying，

Sun Tuo. Change of pore volume in dolomitization[J]. Petroleum

Geology amp; Expeximent， 1985， 7（1）： 29-37.］

［39］ Lucia F J. Origin and petrophysics of dolostone pore space[M]//

Braithwaite C J R， Rizzi G， Darke G. The geometry and petrogenesis

of dolomite hydrocarbon reservoirs. London： Geological

Society， London， Special Publication， 2004， 235（1）： 141-155.

［40］ Morrow D W. Diagenesis 1. dolomite-Part 1： The chemistry of

dolomitization and dolomite precipitation[J]. Geoscience Canada，

1982， 9（1）： 5-13.

［41］ Hardie L A. Dolomitization： A critical view of some current

views[J]. Journal of Sedimentary Research， 1987， 57（1）：

166-183.

［42］ Land L S. The origin of massive dolomite[J]. Journal of Geological

Education， 1985， 33（2）： 112-125.

［43］ Lazar B， Starinsky A， Katz A， et al. The carbonate system in

hypersaline solutions： Alkalinity and CaCO3 solubility of

evaporated seawater[J]. Limnology and Oceanography， 1983， 28

（5）： 978-986.

［44］ Sun S Q. Skeletal aragonite dissolution from hypersaline seawater：

A hypothesis[J]. Sedimentary Geology， 1992， 77（3/4）：

249-257.

［45］ Sun S Q. Dolomite reservoirs： Porosity evolution and reservoir

characteristics[J]. AAPG Bulletin， 1995， 79（2）： 186-204.

［46］ Machel H G. Concepts and models of dolomitization： A critical

reappraisal[M]//Braithwaite C J R， Rizzi G， Darke G. The

geometry and petrogenesis of dolomite hydrocarbon reservoirs.

London： Geological Society， London， Special Publication， 2004，

235（1）： 7-63.

［47］ Landes K K. Porosity through dolomitization[J]. AAPG Bulletin，

1946， 30（3）： 305-318.

［48］ Wendte J. Origin of molds in dolostones formed by the dissolution

of calcitic grains： Evidence from the Swan Hills Formation

in west-central Alberta and other Devonian formations in Alberta

and northeastern British Columbia[J]. Bulletin of Canadian

Petroleum Geology， 2006， 54（2）： 91-109.

［49］ Jonas L， Müller T， Dohmen R， et al. Transport-controlled hydrothermal

replacement of calcite by Mg-carbonates[J]. Geology，

2015， 43（9）： 779-782.

［50］ Zempolich W G， Baker P A. Experimental and natural mimetic

dolomitization of aragonite ooids[J]. Journal of Sedimentary Research，

1993， 63（4）： 596-606.

［51］ Xia F， Brugger J， Ngothai Y， et al. Three-dimensional ordered arrays

of zeolite nanocrystals with uniform size and orientation by a

pseudomorphic coupled dissolution-reprecipitation replacement

route[J]. Crystal Growth amp; Design， 2009， 9（11）： 4902-4906.

［52］ Xia F， Brugger J， Chen G R， et al. Mechanism and kinetics of

pseudomorphic mineral replacement reactions： A case study of

the replacement of pentlandite by violarite[J]. Geochimica et

Cosmochimica Acta， 2009， 73（7）： 1945-1969.

［53］ Putnis C V， Tsukamoto K， Nishimura Y. Direct observations of

pseudomorphism： Compositional and textural evolution at a

fluid-solid interface[J]. American Mineralogist， 2005， 90（11/

12）： 1909-1912.

［54］ Budd D A， Mathias W D. Formation of lateral patterns in rock

properties by dolomitization： Evidence from a Miocene reaction

front （Bonaire， Netherlands Antilles）[J]. Journal of Sedimentary

Research， 2015， 85（9）： 1082-1101.

［55］ Sibley D F. Secular changes in the amount and texture of dolomite[

J]. Geology， 1991， 19（2）： 151-154.

［56］ Gabellone T， Whitaker F. Secular variations in seawater chemistry

controlling dolomitization in shallow reflux systems： Insights

from reactive transport modelling[J]. Sedimentology， 2016， 63

（5）： 1233-1259.

［57］ 黃擎宇，劉偉，張艷秋，等. 白云石化作用及白云巖儲層研究

進展[J]. 地球科學進展，2015，30（5）：539-551. ［Huang

Qingyu， Liu Wei， Zhang Yanqiu， et al. Progress of research on

dolomitization and dolomite reservoir[J]. Advances in Earth

Science， 2015， 30（5）： 539-551.］

［58］ Putnis A. Mineral replacement reactions[J]. Reviews in Mineralogy

and Geochemistry， 2009， 70（1）： 87-124.

［59］ Usdowski E. Synthesis of dolomite and geochemical implications[

M]//Purser B， Tucker M， Zenger D. Dolomites： A volume

in honour of dolomieu. Oxford： Blackwell Scientific， 1994：

345-360.

［60］ Jones G D， Xiao Y T. Dolomitization， anhydrite cementation，

and porosity evolution in a reflux system： Insights from reactive

transport models[J]. AAPG Bulletin， 2005， 89（5）： 577-601.

［61］ Al-Helal A B， Whitaker F F， Xiao Y T. Reactive transport modeling

of brine reflux： Dolomitization， anhydrite precipitation， and

porosity evolution[J]. Journal of Sedimentary Research， 2012， 82

（2）： 196-215.

［62］ Garcia-Fresca B， Lucia F J， Sharp J M， et al. Outcropconstrained

hydrogeological simulations of brine reflux and

early dolomitization of the Permian San Andres Formation[J].

AAPG Bulletin， 2012， 96（9）： 1757-1781.

［63］ Whitaker F F， Xiao Y T. Reactive transport modeling of early

burial dolomitization of carbonate platforms by geothermal convection[

J]. AAPG Bulletin， 2010， 94（6）： 889-917.