?

川中古隆起北斜坡蓬萊氣區上震旦統—下寒武統海相碳酸鹽巖天然氣多層系立體成藏

2024-01-11 11:04張本健宋澤章陶佳麗田興旺丁孝恒楊岱林張自力孫奕婷黃茂軒
海相油氣地質 2023年4期
關鍵詞:東壩龍王廟滄浪

張本健,周 剛,宋澤章,嚴 威,汪 華,陶佳麗,田興旺,丁孝恒,鐘 原,馬 奎,楊岱林,李 勇,張自力,陳 曦,孫奕婷,向 柱,黃茂軒

1 中國石油西南油氣田公司勘探開發研究院;2 油氣資源與探測國家重點實驗室;3 中國石油大學(北京)地球科學學院

0 前 言

近三年來,中國石油西南油氣田公司(簡稱西南油氣田)持續深化川中古隆起北斜坡(安岳氣田北側)深層(4 500~6 000 m)—超深層(>6 000 m)天然氣勘探。截至2021年底,公司已相繼在震旦系燈影組二段[1]和四段[2]、下寒武統滄浪鋪組一段[3]、二疊系茅口組二段[4]取得勘探突破,落實了萬億方新增地質儲量,形成了繼安岳氣區之外的另一個萬億方級大氣區——蓬萊氣區[2,5]。2022 年9 月,東壩1井在下寒武統龍王廟組測試獲氣20.28×104m3/d,首次在川中古隆起北斜坡龍王廟組獲工業性氣流,深層天然氣勘探再次獲重大戰略突破。

回顧西南油氣田在川中古隆起北斜坡上震旦統—下寒武統的油氣勘探與突破,大致經歷了以下幾個階段:

(1)燈影組二段(燈二段)勘探突破

2020 年5 月,風險探井蓬探1 井在燈二段獲121.98×104m3/d 的高產氣流[3],反映出燈二段良好的含氣性和勘探潛力[6]。此外,蓬探1 井西北側的中江2 井在燈二段同樣揭示了良好的含氣儲層(儲層總厚度100.6 m,測井解釋氣層33.6 m)及油氣顯示(測試獲氣3.36×104m3/d)[7]。

(2)滄浪鋪組勘探突破

2020 年,風險探井角探1 井在滄浪鋪組一段揭示了厚度約30 m 的鮞粒白云巖儲層(測井解釋孔隙度約為3.8%),初試獲51.62×104m3/d 的工業氣流,首次取得滄浪鋪組碳酸鹽巖儲層的勘探突破,表明北斜坡滄浪鋪組具備形成規模氣藏的條件[8]。2020年底開始鉆進的風險探井東壩1井未揭示滄浪鋪組儲層,這為界定滄浪鋪組儲層的展布范圍提供了依據。

(3)燈影組四段(燈四段)勘探突破

2021年8月,部署在射洪—鹽亭區塊東壩1井區的蓬深1 井在燈四段測試獲得天然氣3.4×104m3/d,證實了蓬萊氣區燈四段臺緣帶發育大型單斜構造背景下的地層-巖性氣藏。東壩1 井區燈四段不僅整體儲集條件優越(地層厚度達350 m,平均孔隙度為3.5%),而且含氣性好,是規模增儲的現實區塊。2022 年9 月,東壩1 井區蓬深13 井燈四上亞段(6 655.0~6 680.0 m、6 740.0~6 775.0 m)酸化后初測產氣10.56×104m3/d;2022年2月和4月,分別對東壩1 井燈四下亞段、上亞段進行酸化改造和測試求產,獲氣28.54×104m3/d(5月)。

(4)龍王廟組勘探突破

西南油氣田對斜坡區龍王廟組天然氣的勘探潛力進行了多年的論證。早在2015 年12 月,部署于磨溪地區北部(緊鄰北斜坡)的磨溪52 井在龍王廟組測試平均產氣18.42×104m3/d,推測斜坡區龍王廟組同樣具有良好的勘探潛力。之后,位于北斜坡南部的南充1 井和立探1 井分別揭示龍王廟組有效儲層厚度約為24 m和27 m,而且儲層固體瀝青充填嚴重,證實了北斜坡龍王廟組儲層中曾發生大規模油氣充注,形成了古油藏且發生了古油藏的裂解[9]。2022 年9 月,東壩1 井在龍王廟組測試獲工業性氣流(20.28×104m3/d),標志著斜坡區龍王廟組勘探取得重大突破。至此,西南油氣田在川中古隆起北斜坡上震旦統—下寒武統實現了多層系勘探突破。

本文以東壩1 井龍王廟組重大勘探突破為契機,對蓬萊氣區特別是東壩1 井區上震旦統—下寒武統的沉積相、巖相、儲層進行分析,在此基礎上,從烴源巖及源-儲配置、沉積相及成巖作用、多層系含氣特征及圈閉等3個方面著手,對斜坡區多層系、立體含氣的成藏條件進行解剖,旨在深化川中古隆起北斜坡上震旦統—下寒武統多層系天然氣立體成藏的認識。

1 東壩1井震旦系—寒武系地層特征

蓬萊氣區面積約為2×104km2,構造上位于四川盆地中部,其南部緊鄰安岳氣田,西側緊鄰德陽—安岳裂陷槽,北側以九龍山為界;地理上位于遂寧以北,東靠南充,西臨德陽,北至綿陽。東壩1 井區位于蓬萊氣區西南部大型斜坡區(圖1a)[10]。

東壩1 井于2020 年12 月17 日開鉆,鉆井深度為6 720 m,自上而下揭示侏羅系蓬萊鎮組、遂寧組、沙溪廟組、自流井組,三疊系須家河組、雷口坡組、嘉陵江組、飛仙關組,二疊系長興組、龍潭組、茅口組、棲霞組、梁山組,下寒武統龍王廟組、滄浪鋪組、筇竹寺組,上震旦統燈影組四段和三段。對東壩1井上震旦統—下寒武統(圖1b)分述如下:

燈影組 厚度為434 m,其中,燈三段厚85 m,燈四段厚349 m,測井響應以低自然伽馬、高電阻率為特征。燈三段巖性以灰巖為主;燈四段可細分為2 段:燈四上亞段(6 286 m~6 485 m)以藻白云巖為主,僅在頂部發育少量灰巖;燈四下亞段(6 485 m~6 635 m)以泥粉晶白云巖為主。

筇竹寺組 上部主要發育灰色、深灰色泥巖夾少量粉砂質泥巖;下部發育灰色、深灰色頁巖。測井響應以典型的高自然伽馬、低電阻率為特征。烴源巖(含麥地坪組)厚度大(329 m)。

滄浪鋪組 厚度為146 m,可進一步細分為2段:滄浪鋪組一段(滄一段)厚59 m,底部巖性以碎屑巖為主,中上部以鮞粒白云巖為主,夾雜砂質白云巖、碎屑白云巖;滄浪鋪組二段(滄二段)厚87 m,中下部以白云巖為主,夾碎屑巖,上部為碎屑巖。

龍王廟組 頂部與下二疊統梁山組泥巖呈整合接觸,底部與滄浪鋪組砂巖呈整合接觸,地層總厚度為86 m,發育有效儲層28.3 m。龍王廟組巖性以泥粉晶白云巖為主,夾砂屑白云巖。

2 沉積相及儲層

2.1 沉積相

燈影組沉積時期,蓬萊氣區以典型的臺地相沉積為主(圖2a)。燈四段沉積早期,水體整體貧藻,蓬萊氣區的沉積環境主要為灘間海,多沉積泥晶白云巖;燈四段沉積中、晚期,水體變淺,各類菌藻大量生長,為微生物白云巖的沉積創造了有利條件,因此燈四段藻白云巖十分發育[11]。東壩1井位于燈四段臺緣帶丘灘相發育有利區。

圖2 蓬萊氣區上震旦統—下寒武統沉積微相展布圖Fig.2 Sedimentary microfacies plans of Upper Sinian-Lower Cambrian in Penglai gas area

蓬萊氣區滄一段主體沉積環境為清水、淺水陸棚(圖2b)。滄一段沉積早期,氣區西部為濱岸-混積陸棚相,以碎屑巖沉積為主;東部則主要為清水陸棚相,沉積碳酸鹽,陸棚邊緣發育顆粒灘相白云巖。滄一段沉積晚期,發生大規模海退,形成混積淺水陸棚沉積[12]。東壩1 井位于混積淺水陸棚相區。

龍王廟組沉積時期,蓬萊氣區為局限臺地沉積環境,發育臺內灘、潮坪、潟湖等沉積亞相(圖2c),厚層顆粒灘是優質儲層發育的有利相帶。龍王廟組縱向上發育完整的海侵海退旋回(圖1b),海平面升降控制了高能灘體的發育。龍王廟組沉積早期,伴隨海平面上升,水體能量逐漸減弱,陸源供給減少,發育了第一套灘體,巖性主要為砂屑白云巖和鮞粒白云巖;龍王廟組沉積中、晚期,伴隨海平面下降,水體能量逐漸增強,發育了第二套灘體。受控于西高東低的古地貌格局,龍王廟組總體表現出西薄東厚的沉積特征,自西向東由混積潮坪相過渡至粒屑灘-砂屑灘-灘間海相,顆粒灘廣泛發育,為優質儲層發育奠定了良好基礎[13]。東壩1井位于粒屑灘相,靠近龍王廟組剝蝕邊界,巖溶作用強烈。

2.2 儲層

2.2.1 儲層巖性

蓬萊氣區燈四段儲層巖性主要為白灰色夾黑色藻白云巖、泥晶白云巖和晶粒白云巖(圖1b)。燈四上亞段丘灘相發育,以藻凝塊、藻疊層白云巖為主;燈四下亞段丘灘體欠發育,以泥粉晶白云巖為主,頂部為硅質白云巖和泥晶白云巖。受丘灘體分布控制,燈四上亞段儲層較發育,是主力儲層段[14]。蓬萊氣區燈四段臺緣帶被灘間海分隔成多個獨立的丘灘體,且與磨溪臺緣帶分隔開來。東壩1 井所在丘灘體沿臺緣呈近南北向分布[15-16](圖2a)。由連井儲層對比剖面(圖3)可見:橫向上燈四上亞段丘灘體儲層(孔隙度>3%)厚度從東壩1 井區(164 m)向角探1 井區(80 m)逐漸變小,下亞段丘灘體儲層厚度從東壩1 井區(15 m)向角探1 井區(45 m)逐漸變大。蓬萊氣區燈四段孔洞型優質儲層發育,儲集條件優越:儲集空間主要為溶蝕孔、洞,并以中、小型孔洞為主,常見白云石、石英或瀝青充填(圖4a,多為半充填),局部發育較大溶洞(直徑最大可達3 cm),偶見溶蝕縫。鑄體薄片鏡下顯示,燈四段溶蝕孔以粒間溶孔和晶間溶孔為主(圖4b)。

圖3 蓬萊氣區燈四段儲層連井對比剖面(剖面位置見圖1a)Fig.3 Inter-well reservoir correlation profile of the Dengying Member 4 in Penglai gas area(location is shown in Fig.1a)

圖4 蓬萊氣區上震旦統—下寒武統儲層巖性及孔隙結構Fig.4 Lithology and pore structure of Upper Sinian-Lower Cambrian reservoirs in Penglai gas area

蓬萊氣區滄浪鋪組儲層主要發育于滄一段[17],主要為灰黑色鮞粒白云巖夾砂質白云巖、砂屑白云巖(圖4c),多見水平層理,裂縫不發育;儲集空間以溶蝕孔、洞為主,裂縫不發育,瀝青沿孔隙周邊填充(圖4d),說明固體瀝青的填充晚于溶蝕孔洞的形成[18]。平面上,滄一段碳酸鹽巖厚值區主要分布于射洪西南部及北部。東壩1井揭示滄浪鋪組厚度為146 m,測井解釋為干層,巖心幾乎未見孔隙。

受古地貌控制,蓬萊氣區龍王廟組廣泛發育高能顆粒灘相,為優質儲層發育奠定了良好的基礎?;趦舆B井對比剖面(圖5)分析:縱向上,高能顆粒灘體主要發育在龍王廟組中上部;橫向上,高能灘相白云巖儲層(孔隙度>3%)穩定發育,厚度分布在12~40 m 之間(東壩1 井處厚28.1 m)。儲層巖性以砂屑白云巖、細—中晶白云巖為主。東壩1 井龍王廟組主要為粉晶白云巖與砂屑白云巖互層,發育4 套砂屑白云巖。龍王廟組儲集條件優越,儲集空間以粒間(溶)孔、溶洞為主(圖4e),固體瀝青在孔隙網絡中廣泛分布,呈半充填—全充填(圖4f),指示古油藏的存在。

2.2.2 儲層物性

東壩1井燈四段、滄浪鋪組、龍王廟組均為典型的致密碳酸鹽巖,屬低孔、超低滲儲層(圖6)。本次巖心分析測試結果表明:①燈四段儲層孔隙度變化范圍在2.1%~4.9%之間,平均值為2.7%,并以2%~3%為主(圖6a,占樣品總數的81%);滲透率變化范圍在(0.190~1.534)×10-3μm2之間,平均值為0.339×10-3μm2。東壩1井燈四段儲層非均質性相對較弱,孔隙度—滲透率的相關性較好(圖6b)。②滄浪鋪組儲層孔隙度變化范圍在2.3%~7.9%之間(蓬深1 井、蓬深7 井等井),平均值為2.7%;滲透率變化范圍在(0.002~2.675)×10-3μm2之間,平均值為0.254×10-3μm2。③龍王廟組儲層孔隙度變化范圍在2.9%~5.1%之間(東壩1井、蓬深1 井等井),平均值為3.8%,且主要分布于3%~4%(占樣品總數的56%);滲透率變化范圍在(0.001~0.910) ×10-3μm2之 間,平 均 值 為0.227 ×10-3μm2。與燈影組樣品相比,滄浪鋪組和龍王廟組樣品的孔隙更發育,但非均質性更強、滲透率更低。

圖6 蓬萊氣區上震旦統—下寒武統儲層孔隙度分布直方圖和孔滲關系圖Fig.6 Porosity histogram and porosity-permeability correlations of Upper Sinian-Lower Cambrian reservoirs in Penglai gas area

3 震旦系—寒武系成藏條件

蓬萊氣區上震旦統—下寒武統多層系(燈影組四段、滄浪鋪組、龍王廟組)天然氣富集的主控因素可以歸納為以下3 個方面:①烴源巖質量和規模均較好,近源成藏。②受古地貌控制,丘灘相廣泛發育,橫向穩定分布;溶蝕等后期成巖作用改造強度大,有利于孔洞型優質儲層發育,儲集條件優越。③目的層系長期穩定處于北斜坡部位,受構造和巖性控制,發育大型單斜構造背景下相互獨立的巖性圈閉和構造-巖性復合圈閉。

3.1 優質烴源巖分布廣泛,近源成藏

蓬萊氣區上震旦統—下寒武統共發育2套烴源巖:燈三段泥頁巖和筇竹寺組黑色頁巖[5]。與筇竹寺組相比,燈三段烴源巖厚度小,雖對天然氣成藏有一定貢獻,但不是主力烴源巖[19-20]。平面上,筇竹寺組烴源巖在蓬萊氣區分布范圍廣,厚度變化范圍大(120~650 m):西北部厚度最大(>600 m),向東南方向逐漸減?。▓D7a)。筇竹寺組烴源巖有機質豐度 較 高[7,21],TOC介 于0.35%~2.53%,平 均 值 為0.92%,TOC大于1%的實測數據點占56%(圖7b),主要集中在筇竹寺組下段。

圖7 蓬萊氣區筇竹寺組烴源巖分布及TOC直方圖Fig.7 Thickness contour and TOC histogram of the source rock of Qiongzhusi Formation in Penglai gas area

上震旦統—下寒武統存在3套源-儲組合,由筇竹寺組(源)分別與燈影組、滄浪鋪組和龍王廟組(儲)構成。筇竹寺組作為主力烴源巖,供烴能力強,為燈四段、滄浪鋪組和龍王廟組油氣成藏奠定了充足的資源基礎。筇竹寺組烴源巖層側向接觸或直接覆蓋在燈四段儲層之上,兩者之間所形成的半包裹型源-儲配置關系[22]更有利于油氣充注,成藏效率更高,在空間上表現為上生下儲、旁生側儲。滄浪鋪組儲層直接覆蓋在筇竹寺組之上,滄一段顆粒灘是筇竹寺組優質烴源巖之上的第一套優質儲層,具有直接充注、優先捕獲油氣的有利條件[8],筇竹寺組烴源巖排出的油氣可原地直接向上運移、聚集成藏。龍王廟組儲層雖然在垂向上距離筇竹寺組較遠,但川中地區發育多條基底斷裂,形成良好的油氣輸導通道,油氣可沿斷裂向上運移充注[5];同時,龍王廟組廣泛發育的顆粒灘相優質儲層對烴類在目的層內部的運移輸導十分有利。

3.2 高能顆粒灘是儲層發育的基礎,巖溶作用是改善儲層物性條件的關鍵

燈四段儲層的形成和演化受沉積作用和成巖作用的雙重控制。蓬萊氣區燈四段沉積時期,在海平面下降過程中,水體中藻類逐漸繁育,上亞段發育多套丘灘;隨著海平面持續下降,到燈四晚期逐漸轉變為淺水高能環境,發育厚層疊置丘灘體(圖2)。燈影組經歷了多期次巖溶作用,其中桐灣期表生巖溶對改善儲層起到主要作用[23-24]。地表垂向節理和裂縫較為發育,為流體下滲提供了通道,大氣淡水沿節理和裂縫不斷向下滲透、淋濾,導致了長時間、大規模的巖溶,最終形成大量具有儲集意義的溶蝕孔、洞和晶間孔(圖8a—8d)[25]。

滄一段灘相白云巖儲層主要受白云石化顆粒灘相分布及加里東末期巖溶作用的雙重影響[12]。蓬萊氣區處于裂陷槽東側古臺緣高部位,受海平面下降的影響,滄浪鋪組發育白云石化顆粒灘相,為后期灘相儲層的形成提供了良好的物質基礎[26]。川中古隆起在加里東運動時期持續抬升,導致滄浪鋪組遭受剝蝕,形成剝蝕區。根據陰極發光特征(圖8e,8f)分析,后期地表淡水對北斜坡區顆粒灘體進行表生溶蝕改造,沿剝蝕區潛流帶和垂直滲流帶進入滄一段中,進一步改善了滄一段顆粒白云巖儲層。

龍王廟組白云巖儲層的形成、發育受沉積相[27]與巖溶作用共同控制。蓬萊氣區龍王廟組為局限臺地沉積環境,發育粒屑灘、砂屑灘、灘間海、混積潮坪(圖2)。在海平面下降的過程中,形成臺內高能顆粒灘并發生白云石化,成為龍王廟組優質儲層發育的基礎。晚期在大氣淡水溶蝕的影響下,古隆起上的顆粒白云巖儲層進一步改善。東壩1 井區靠近剝蝕邊界,巖溶作用強烈,儲集空間以溶蝕孔、晶間孔、裂縫為主(圖8g—8i),多數被瀝青充填。

3.3 斜坡背景上發育多個相互獨立的巖性圈閉,形成巖性氣藏群

蓬萊氣區位于川中古隆起北斜坡,構造特征為向北傾斜的單斜構造,局部發育小型構造圈閉,但構造圈閉對氣藏不起主要控制作用。燈四段氣藏受斜坡構造背景和巖性圈閉共同控制??v向上,燈四上、下亞段受丘間巖性致密帶分隔,形成相對獨立、互不連通的大型天然氣藏(圖3)。東壩1井燈四上亞段為主要的含氣層段,測試不產地層水,產氣20.28×104m3/d(穩定油壓為24.51 MPa),氣藏壓力系數為1.05,為常壓氣藏;燈四下亞段儲層局部發育,橫向不連片,測試日產氣8.26×104m3(穩定油壓為2.97 MPa),日產水66.8×104m3。橫向上,蓬萊氣區不同井區燈四上亞段的氣層底界各不相同,角探1 井區、蓬深7 井區和東壩1 井區分別為-7 110 m、-6 860 m 和-6 350 m(圖9)。構造位置較低的角探1井區燈四上亞段氣藏底部見水,而構造高部位的東壩1井區燈四上亞段整體含氣,未發現氣水界面,這表明東壩1井區燈四上亞段藻丘灘巖性圈閉頂底板及上傾方向發育的巖性致密帶的封堵是成藏關鍵因素。

滄浪鋪組氣藏分布受巖性與構造共同控制,在二者的共同作用下形成了具有北低南高單斜構造背景的構造-巖性復合圈閉。滄浪鋪組上部發育的致密砂泥巖層和龍王廟組致密厚層碳酸鹽巖是直接蓋層(圖1b),相變的致密灰巖等碳酸鹽巖和砂泥巖形成側向遮擋。滄浪鋪組顆粒白云巖儲層含氣性好,為天然氣富集帶,橫向上形成相對獨立的巖性氣藏(圖9)。滄浪鋪組氣層底界各不相同,角探1井區、蓬深1井區氣層底界海拔分別為-6 600 m、-6 300 m。角探1井區滄一段氣層厚度為14.5 m,測試日產氣51.62×104m3。東壩1井則不在有利區里。

龍王廟組在斜坡背景上發育多個相互獨立的巖性圈閉,形成巖性氣藏群,構造高部位整體含氣(圖9)。龍王廟組氣層底界也各不相同,東壩1 井區、蓬深1 井區氣層底界的海拔分別為-5 731 m、-6 157 m。東壩1井區龍王廟組儲層厚度為28.1 m,主要為顆粒白云巖,在5 727.8~5 772.7 m 井段測井解釋有3 個氣層,測試日產氣20.28×104m3(穩定油壓為25.06 MPa)。

截至2022年底,西南油氣田在蓬萊氣區燈四段刻畫巖性圈閉面積累計約2 280 km2,落實資源量約5 400×108m3(其中,東壩1 井區燈四上亞段新增天然氣含氣面積807.3 km2,新增天然氣控制地質儲量2 386.53×108m3);在滄一段刻畫巖性圈閉7個,含氣面積為1 400 km2,預計儲量規模為2 800×108m3;在龍王廟組刻畫圈閉5個,總面積為2 050 km2,落實資源量4 000×108m3(其中,東壩1 井所在巖性圈閉含氣面積為238 km2,儲量規模達480×108m3)。東壩1井的突破,標志著龍王廟組成為蓬萊氣區震旦系—寒武系立體勘探、規模增儲的重要新層系。

4 勘探突破經驗

堅持不懈探索,突破地質認識,是川中古隆起斜坡區超深層天然氣勘探取得突破的關鍵。蓬萊氣區的勘探部署工作持續4 年時間,在地質上敢于突破傳統理論認識,在技術上綜合應用超深層信號采集技術、高精度測井技術及地震儲層預測技術,形成優化鉆完井及儲層改造技術系列,為超深層鉆探目標的落實及實施提供了有效支撐;圍繞氣藏描述和儲量參數開展了大量攻關研究,為儲量估算提供了堅實的地質認識和技術支撐;積極探索斜坡區臺緣帶大型巖性圈閉目標,在前期“構造-巖性控圈、立體成藏、復式聚集”認識的基礎上,持續深化巖性圈閉形成條件及成藏富集規律研究,有效支撐蓬萊氣區整體部署、立體勘探,最終實現了燈四段、滄浪鋪組及龍王廟組勘探重大突破??碧綄嵺`表明,四川盆地超深層油氣資源豐富,持續深化超深層成儲、成藏地質理論認識,注重地質理論和工程技術齊頭并進、相互促進,繼續向萬米超深層進軍,必將開啟優質規模增儲上產新領域。

5 結 論

東壩1 井下寒武統龍王廟組取得重大勘探突破,證實了川中古隆起北斜坡上震旦統—下寒武統超深層古老海相碳酸鹽巖地層多層系立體成藏:分布廣泛的下寒武統筇竹寺組厚層、優質烴源巖,為多層系立體成藏奠定了良好的物質基礎;良好的源-儲配置關系(上生下儲、旁生側儲)保證了斜坡區多層系天然氣高效充注;沉積微相(高能顆粒灘相)是多套優質儲層形成的基礎,而后期廣泛發育的巖溶作用進一步改善了儲層的物性條件,直接控制了優質儲層的形成和分布;大型單斜構造背景下相互獨立的巖性圈閉和構造-巖性復合圈閉廣泛發育,形成縱向上互相疊置、橫向上連片的巖性氣藏群。深化超深層成儲、成藏地質認識,將會推動持續開拓四川盆地超深層油氣增儲上產新領域。

猜你喜歡
東壩龍王廟滄浪
北京住總集團兩級工會慰問東壩區域項目
山西文明守望工程掠影
論魯敏小說的“輕”與“重”
《蘭心詩鈔》“均州滄浪”意象的“元典”元素分析
滄浪
奇怪森林
川西南龍王廟組油氣成藏控制因素
均州:滄浪之上,滄浪之下
只見滄浪不聞君
西南油氣田勘探獲得重要新發現
91香蕉高清国产线观看免费-97夜夜澡人人爽人人喊a-99久久久无码国产精品9-国产亚洲日韩欧美综合