李 蕭，勞海港，胡秋媛，趙 利

(1.中國石油大學 勝利學院 油氣工程學院，山東 東營 257000； 2.華北理工大學 礦業(yè)工程學院，河北 唐山 063009；3. 中國石化 石油勘探開發(fā)研究院 無錫石油地質研究所，江蘇 無錫 214126)

?

準噶爾盆地西北緣烏夏斷裂帶構造特征及物理模擬實驗

李 蕭1，勞海港2，胡秋媛1，趙 利3

(1.中國石油大學 勝利學院 油氣工程學院，山東 東營 257000； 2.華北理工大學 礦業(yè)工程學院，河北 唐山 063009；3. 中國石化 石油勘探開發(fā)研究院 無錫石油地質研究所，江蘇 無錫 214126)

依據大量地震勘探資料，對準噶爾盆地西北緣烏夏斷裂帶的構造特征進行了系統(tǒng)研究，并將其構造演化劃分為碰撞前陸盆地、前陸盆地逆沖推覆、陸內拗陷逆沖褶皺、陸內拗陷抬升剝蝕和覆蓋定型5個階段。在此基礎上，采用構造物理模擬實驗方法對研究區(qū)典型剖面晚石炭世以來的構造演化進行了正演模擬，取得了與實際地質情況較高的吻合度。模擬結果表明，烏夏斷裂帶經歷了由前陸盆地“前展式”逆沖向陸內拗陷“后展式”逆沖的轉變。同時，實驗結果進一步證實，NW-SE向擠壓應力場的控制是烏夏斷裂帶形成、演化的根本原因，其主要受控于板塊活動及造山帶、邊界控盆斷裂的活動等關鍵因素。

構造特征；物理模擬實驗；烏夏斷裂帶；準噶爾盆地西北緣

前陸沖斷帶指處于盆地與造山帶間，由造山帶向盆地大規(guī)模逆掩推覆的沖斷系統(tǒng)[1-2]。近年來，前陸沖斷帶已成為我國中西部油氣勘探的重點，勘探力度逐步增大。烏夏斷裂帶地處準噶爾盆地西北緣前陸沖斷帶東段，為盆地西北緣中構造活動最強烈的地帶，具有典型的逆沖推覆性質。前人針對烏夏斷裂帶的構造特征、沉積模式、油氣聚集及成藏等方面進行了研究[3-7]，但由于前陸沖斷帶受多期造山運動改造，構造變形復雜，對其構造演化及形成機制的認識仍存在分歧，尤其缺少一定的實際變形支持。構造物理模擬實驗是真實再現(xiàn)盆內各種構造形跡的重要方法，可以直觀、有效地反映前陸沖斷帶的變形過程及構造發(fā)育特征，進而為其動力學機制的研究提供重要依據。因此，筆者基于烏夏斷裂帶構造特征研究的基礎上，結合構造物理模擬實驗，對其典型剖面的發(fā)育演化進行了正演模擬，以期獲得對該區(qū)現(xiàn)今構造發(fā)育演化的合理解釋。

1 區(qū)域地質概況

烏夏斷裂帶地處準噶爾盆地西北緣前陸沖斷帶東段，北抵哈拉阿拉特山，南接瑪湖凹陷(圖1)，自晚石炭世起，經歷海西晚期、印支、燕山、喜馬拉雅等多期構造運動持續(xù)發(fā)育，形成了獨具特色的構造形態(tài)。依據構造位置、基底特征的不同，可將烏夏斷裂帶自北向南劃分為山前隆起區(qū)、烏夏沖斷帶和南部斜坡帶3個構造單元，整體呈現(xiàn)“南北分帶”的構造格局(圖1)。其中，山前隆起區(qū)以擠壓逆沖為主，變形強度大，地層嚴重剝蝕，保存有多期不整合面；中部沖斷帶斷展褶皺發(fā)育，除頂部地層遭受剝蝕外，整體保存較好；南部斜坡帶變形微弱，地層發(fā)育齊全，石炭系—第四系保存相對完整。

2 烏夏斷裂帶構造特征

2.1 斷裂特征

烏夏斷裂帶是準噶爾盆地西北緣中典型的大型疊瓦沖斷系統(tǒng)，由多條近NE向和NEE向弧形逆斷層及逆掩斷層組成。大量地震勘探資料表明，烏夏斷裂帶平面上多呈NE向、NEE向延伸，均出現(xiàn)不同程度的弧形，主斷裂弧形幅度小，次級斷裂弧形幅度大，弧頂凸向盆內，斷面傾向盆外；剖面上傾向NW、NNW，沖斷面多為上陡下緩，常呈逆沖疊瓦狀向深部收斂；同時，上盤沉積厚度明顯小于相應的下盤，且上盤一般為連續(xù)沉積，表明斷裂活動具有典型的“同生性”。

依據研究區(qū)內斷裂規(guī)模及其控盆作用，可將烏夏斷裂帶自西向東劃分為烏爾禾斷裂、夏子街斷裂、夏紅斷裂3條一級斷裂，均為主要控盆斷裂，其下又包含一系列控制次級構造形成的二級斷裂。表1詳細討論了研究區(qū)3條一級斷裂的構造特征。

2.2 構造演化分析

烏夏斷裂帶既繼承了準噶爾盆地晚古生代以來疊合盆地的特點[8-9]，又有其獨特的構造變形特征。依據區(qū)域典型剖面的地震勘探成果，結合構造演化特征及前人研究成果，筆者認為，自晚石炭世之后，烏夏斷裂帶經歷了5個演化階段。

2.2.1 碰撞前陸盆地階段(C3)

晚石炭世，伴隨西準噶爾洋殼NW向俯沖、消減于哈薩克斯坦板塊，洋盆完全閉合，陸陸發(fā)生碰撞(圖3a)，產生達爾布特斷裂，盆地西北緣隆起成陸，開始出現(xiàn)強烈的碰撞造山運動。

2.2.2 前陸盆地逆沖推覆階段(P)

晚石炭世洋盆閉合，但地殼深部的洋殼仍繼續(xù)俯沖、消減。早二疊世，哈山作為烏夏斷裂帶北緣推覆特征最明顯的推覆構造之一，已初見雛形，為烏夏斷裂帶NW向擠壓變形提供了動力來源，同時伴隨廣泛的火山活動。

中二疊世，深部洋殼的俯沖、消減已趨于停止，火山活動趨于停息，地殼趨于穩(wěn)定[10-11]，研究區(qū)發(fā)生強烈逆沖推覆，前陸盆地持續(xù)發(fā)育至高峰，達爾布特斷裂為研究區(qū)在北緣的控盆斷裂。造山帶持續(xù)SE向強烈推進擠壓，在山前依次前展式發(fā)育了一系列NE向斷裂(圖3b)，如烏爾禾南斷裂、夏5井斷裂。

2.2.3 陸內拗陷逆沖褶皺階段(T)

進入三疊紀，準噶爾盆地整體抬升，盆地進入整體陸內拗陷階段，受印支運動影響，由NW-SE向擠壓轉為近SN向擠壓，并伴有左旋剪切，在烏夏斷裂帶自南向北后展式發(fā)育了一系列近EW向斷裂(圖3c)，如夏21井斷裂、夏10井斷裂東段、烏蘭林格斷裂等。同時，二疊紀形成的NE向斷裂持續(xù)逆沖推覆，并伴隨大規(guī)模褶皺發(fā)育。三疊紀末，受SN向強烈擠壓和重新活動的達爾布特斷裂影響，斜坡地帶持續(xù)隆起，哈山形成。

圖1 準噶爾盆地西北緣烏夏斷裂帶構造綱要Fig.1 Tectonic pattern of the Wuxia fault belt on the northwestern margin of Junggar Basin表1 準噶爾盆地西北緣烏夏斷裂帶各一級斷裂構造特征Table 1 Characteristics of first-order faults in the Wuxia fault belt on the northwestern margin of Junggar Basin

斷裂名稱次級斷裂斷裂產狀走向傾向傾角滑脫面平面形態(tài)剖面形態(tài)斷距主要活動時期烏爾禾斷裂 百烏斷裂、烏爾禾斷裂、烏南斷裂NENW上部50°～60°；下部10°～20° 呈勺形(圖2a，b) 斷裂呈弧形，整體近平行排列 上陡下緩呈鏟式，為逆掩斷層(圖2a，b) 剖面斷至三疊系，自下而上逐漸減小，下部斷距大者可達750m，上部僅50m晚二疊世—晚三疊世夏子街斷裂 夏10井斷裂、夏21井斷裂、夏紅北斷裂近EW近N35°～55° 呈勺形(圖2c)；石炭系和二疊系推移距離逐漸增大，可達15km 斷裂呈弧形，夏紅北斷裂為前鋒斷層，與烏南斷裂相接，整體近平行排列 剖面呈坡坪式，具平緩的“躺椅”形態(tài)；夏21井斷裂與夏紅北斷裂強烈擠壓破碎呈“反Y”字形，均為逆掩斷層(圖2c) 斷裂西部斷至二疊系，東部斷至三疊系，斷距下大上小，下部斷距大者可達500m，上部僅為75m三疊紀夏紅斷裂 夏紅南斷裂、夏59井斷裂、夏5井斷裂NEENNW25°～50° 呈舒緩波狀(圖2d)；石炭系—二疊系長距離逆掩推覆，平均推移距離達14km 斷裂呈弧形，整體近平行排列 整體以低角度逆沖為主，斷面前部陡、中部平、尾部向深處逐漸收斂為一條(圖2d) 斷裂西部斷至二疊系頂部，中部斷至侏羅系，東部斷至三疊系；二疊系斷距下大上小，1400～300m晚二疊世—晚三疊世

圖2 準噶爾盆地西北緣烏夏斷裂帶典型剖面構造特征 剖面位置見圖1。Fig.2 Characteristics of typical profiles of the Wuxia fault belt on the northwestern margin of Junggar Basin

2.2.4 陸內拗陷抬升剝蝕階段(J-K)

燕山運動早中期，為準噶爾西北緣陸內拗陷的消亡期，而烏夏斷裂帶仍有繼承性活動，夏59井斷裂在此階段發(fā)育。早—中侏羅世，除哈山部分區(qū)域處于較高地勢，全區(qū)穩(wěn)定下沉，具有明顯的超覆特征(圖3d)。白堊紀時，受燕山運動晚期影響，構造活動較弱，除哈山主體仍有抬升外，烏夏斷裂帶斷裂和褶皺基本停止。

圖3 準噶爾盆地西北緣烏夏斷裂帶構造演化示意Fig.3 Tectonic evolution of the Wuxia fault belt on the northwestern margin of Junggar Basin

2.2.5 覆蓋定型階段(N-Q)

至新生代時，整體平穩(wěn)充填，自SE向NW低角度超覆，由南向北逐漸尖滅，研究區(qū)構造覆蓋定型(圖3e)。

3 烏夏斷裂帶構造物理模擬實驗

3.1 實驗原理與設備

模擬實驗中，實驗模型與構造原型、模型材料與實際地層間均要遵循相似性原理[12-13]。本次實驗采用中國石油大學(華東)自行研制的DGU-Ⅰ型地質構造物理模擬實驗裝置，驅動方式為液壓驅動，注液方式為平流泵自動注液。

3.2 實驗條件與模型設計

依據前述的構造演化分析，將物理模擬實驗的條件設置為：地質時期——晚石炭世(C3)以來；動力來源——板塊活動及造山帶、邊界控盆斷裂活動產生的擠壓應力場；演化過程——晚石炭世以來的持續(xù)擠壓、隆升?？紤]上述地質要素基礎上，依據相似性原理，選擇研究區(qū)典型NW-SE向剖面(C-C’，見圖2)，設計構造物理模擬實驗如下。

圖4 構造物理模擬實驗條件示意Fig.4 Physical modeling for tectonic evolution

實驗砂體尺寸長80 cm，寬20 cm，高15 cm(圖4)。為便于觀察，用黃色天然石英砂模擬地層，染色的石英砂模擬地質界線，模型相似系數為3.6×10-6。與隆起區(qū)/造山帶相比，斷陷區(qū)/拗陷區(qū)具有“下強上弱”的流變學特征[14]?；诖耍瑢嶒災Ｐ妥髠蠕佋O15 cm厚的石英砂，右側底部鋪設8 cm厚的聚苯乙烯泡沫基底(預設一定角度以實現(xiàn)深部的俯沖、消減)，上部鋪設7 cm厚的石英砂(圖4)，通過兩側平流泵提供力源。

3.3 實驗過程與結果分析

3.3.1 實驗過程

本次構造物理模擬實驗是在著重考慮邊界擠壓、深部俯沖2大因素對逆沖構造形成、演化影響的基礎上而實現(xiàn)的。實驗結果需從統(tǒng)計的角度綜合評價，因此在不斷修正、改進過程中共進行了5組實驗，此處僅選擇擬合程度最高的一組分析、討論。

依據區(qū)域構造演化分析，實驗采用雙向擠壓的方式模擬烏夏斷裂帶的擠壓逆沖活動，施力大小通過平流泵中的蒸餾水流速控制。實驗歷時72 min，時間相似比為2.86 Ma/min。其中，兩側平流泵推進的速度依次為0.21 cm/min(0～7 min)，0.25 cm/min(8～23 min)，0.13 cm/min(24～41 min)，0.08 cm/min(42～60 min)，0.05 cm/min(61～72 min)。

3.3.2 實驗結果分析

根據地質模型的設計要求，烏夏斷裂帶典型剖面晚石炭世以來的構造演化過程在構造模擬實驗中均獲得再現(xiàn)[15]。圖5a為初始狀態(tài)，伴隨持續(xù)的擠壓作用，實驗進行至10 min時砂體右側開始發(fā)育逆沖斷層①，傾向NW，斷面上陡下緩(圖5b)，對應中晚二疊世研究區(qū)山前發(fā)育的夏10井斷裂西段；23 min時，逆沖斷層①持續(xù)活動，同時在其北西側發(fā)育逆沖斷層②，產狀與斷層①基本一致(圖5c)，對應研究區(qū)NE向發(fā)育的夏紅北斷裂。伴隨斷層的持續(xù)逆沖活動，研究區(qū)地層開始出現(xiàn)大規(guī)模褶皺變形，30 min時，斷層②北西側地層大幅度褶皺隆升(圖5d)。實驗進行至41 min，斷層①與斷層②之間發(fā)育逆沖斷層③，對應三疊紀開始活動的夏21井斷裂；同時，褶皺地帶持續(xù)擠壓、隆升，哈山形成。實驗持續(xù)進行至72 min，除隆起部位有少量抬升外，剖面基本保持了前一階段的構造特征(圖5e)。

圖5 準噶爾盆地西北緣 烏夏斷裂帶構造物理模擬實驗過程Fig.5 Simulation process of tectonic evolution of the Wuxia fault belt on the northwestern margin of Junggar Basin

已有研究資料表明，中寒武世—早石炭世，準噶爾盆地西北緣位于準噶爾古地塊西北部的被動大陸邊緣。伴隨西準噶爾洋殼NW向俯沖、消減，陸陸碰撞，洋盆閉合，深部洋殼俯沖持續(xù)進行[16]；直至中晚二疊世，達爾布特斷裂成為研究區(qū)主要的控盆斷裂，造山帶持續(xù)向盆內強烈擠壓，這些均為研究區(qū)擠壓變形提供了強大的NW-SE向動力來源。構造物理模擬實驗證實了上述因素的主控作用，中晚二疊世山前發(fā)育的逆沖斷層①——夏10井斷裂西段，到其北西側發(fā)育的逆沖斷層②——夏紅北斷裂，研究區(qū)處于前陸盆地的前展式逆沖階段。進入三疊紀，受印支運動影響產生的近SN向擠壓應力場，一方面派生出NW-SE向的擠壓應力場，促使前期發(fā)育的逆沖斷層①、②持續(xù)活動；另一方面導致逆沖斷層①、②間發(fā)育逆沖斷層③——夏21井斷裂。至此，研究區(qū)由前陸盆地的“前展式”逆沖階段向陸內拗陷的“后展式”逆沖階段轉變。進入侏羅紀，研究區(qū)構造活動強度明顯減弱，構造變形基本定型。

4 結論

(1)烏夏斷裂帶是一組具有逆沖推覆性質的斷裂帶，由多條弧形逆斷層及逆掩斷層組成。平面上斷裂多呈NE向、NEE向延伸，均呈不同程度的弧形，近平行排列；剖面上傾向NW、NNW，沖斷面為勺形或舒緩波狀，常呈逆沖疊瓦狀向深部收斂。斷裂構造活動強烈，以晚二疊世—晚三疊世活動為主，多表現(xiàn)為逆沖/逆掩斷層性質，且斷層相關褶皺發(fā)育。

(2)烏夏斷裂帶自晚石炭世構造演化可分為5個階段，包括晚石炭世后碰撞前陸盆地階段、二疊紀前陸盆地逆沖推覆階段、三疊紀陸內拗陷逆沖褶皺階段、侏羅紀—白堊紀陸內拗陷抬升剝蝕階段和新生代覆蓋定型階段。

(3)烏夏斷裂帶典型剖面晚石炭世以來的構造演化過程，在構造物理模擬實驗結果中均獲得再現(xiàn)。斷層的發(fā)育演化序列表明，烏夏斷裂帶經歷了由前陸盆地“前展式”逆沖向陸內拗陷“后展式”逆沖的轉變。同時，模擬結果進一步證實，北西—南東向擠壓應力場為烏夏斷裂帶形成、演化的主控因素，板塊活動及造山帶、邊界控盆斷裂的活動為其提供了重要的動力來源。

[1] Morley C K.A classification of thrust fronts[J].AAPG Bulletin,1986,70(1):12-25.

[2] Vann I R,Graham R H,Hayward A B.The structure of mountain fronts[J].Journal of Structural Geology,1986,8(3/4):215-227.

[3] 蔚遠江,何登發(fā),雷振宇,等.準噶爾盆地西北緣前陸沖斷帶二疊紀逆沖斷裂活動的沉積響應[J].地質學報,2004,78(5):612-625.

Yu Yuanjiang,He Dengfa,Lei Zhenyu,et al.Sedimentary response to the activity of the Permian thrusting fault in the foreland thrust belt of the northwestern Junggar Basin[J].Acta Geologica Sinica,2004,78(5):612-625.

[4] 王穎,王英民,趙錫奎.構造模擬實驗在構造研究中的應用:以樁西潛山為例[J].石油實驗地質,2004,26(3):308-312.

Wang Ying,Wang Yingmin,Zhao Xikui.Application of simulation experiment to the study of structural evolution:An example of the Zhuangxi buried hill[J].Petroleum Geology & Experiment,2004,26(3):308-312.

[5] 白振華,孫培元,孫寶才.準噶爾盆地烏夏斷裂帶構造特征及其油氣成藏意義[J].石油地質與工程,2011,25(3):19-22.

Bai Zhenhua,Sun Peiyuan,Sun Baocai.The structural characters and its implications for hydrocarbon accumulation in Wu-Xia fault belt,Junggar Basin[J].Petroleum Geology and Engineering,2011,25(3):19-22.

[6] 馮建偉,戴俊生,葛盛權.準噶爾盆地烏夏斷裂帶構造演化及油氣聚集[J].中國石油大學學報(自然科學版),2008,32(3):23-29.

Feng Jianwei,Dai Junsheng,Ge Shengquan.Structural evolution and pool-forming in Wuxia fault belt of Junggar Basin[J].Journal of China University of Petroleum,2008,32(3):23-29.

[7] 宋永東,戴俊生,吳孔友.準噶爾盆地西北緣烏夏斷裂帶構造特征與油氣成藏模式[J].西安石油大學學報(自然科學版),2009,24(3):17-20.

Song Yongdong,Dai Junsheng,Wu Kongyou.Study on the structural features and hydrocarbon accumulation modes of Wuxia fault belt in the northwest margin of Junggar Basin[J].Journal of Xi’an Shiyou University (Natural Science Edition),2009,24(3):17-20.

[8] 劉玉萍,尹宏偉,張潔,等.褶皺—沖斷體系雙層滑脫構造變形物理模擬實驗[J].石油實驗地質,2008,30(4):424-428.

Liu Yuping,Yin Hongwei,Zhang Jie,et al.Physical modeling on double decollements structure in fold-and-thrust belt[J].Petro-leum Geology & Experiment,2008,30(4):424-428.

[9] 鄭孟林,金之鈞.準噶爾盆地北東向構造及其油氣地質意義[J].地球科學(中國地質大學學報),2004,29(4):467-472.

Zheng Menglin,Jin Zhijun.NE-trending structures and their significance on petroleum geology in Junggar Basin[J].Earth Science(Journal of China University of Geosciences),2004,29(4):467-472.

[10] 陳新,盧華復,舒良樹,等.準噶爾盆地構造演化分析新進展[J].高校地質學報,2002,8(3):257-267.

Chen Xin,Lu Huafu,Shu Liangshu,et al.Study on tectonic evolution of Junggar Basin[J].Geological Journal of China Universities,2002,8(3):257-267.

[11] 王圣柱,林會喜,張奎華,等.準噶爾盆地北緣哈山構造帶油氣輸導系統(tǒng)研究[J].特種油氣藏,2015,22(6):30-34.

Wang Shengzhu,Lin Huixi,Zhang Kuihua,et al.Study on hydrocarbon transport system in Kazan tectonic zone at north margin of Junggar Basin[J].Special Oil & Gas Reservoirs,2015,22(6):30-34.

[12] Ramberg H.Gravity,Deformation and the earth’s crust in theory,experiment,and geological application[M].London:Academic Press,1981.

[13] 戈紅星,Vendeville B C,Jackson M P A.前陸褶皺沖斷帶厚皮縮短鹽構造運動的物理模擬[J].高校地質學報,2004,10(1):39-49.

Ge Hongxing,Vendeville B C,Jackson M P A.Physical models of thick-skinned contractional salt tectonics in a foreland fold-and-thrust belt[J].Geological Journal of China Universities,2004,10(1):39-49.

[14] 胡秋媛,李理.魯西地區(qū)晚中生代—古近紀伸展構造的應力場數值模擬[J].石油實驗地質,2015,37(2):259-266.

Hu Qiuyuan,Li Li.Numerical simulations of tectonic stress fields for Late Mesozoic-Paleogene extensional tectonics in western Shandong[J].Petroleum Geology & Experiment,2015,37(2):259-266.

[15] 胡秋媛,董大偉,趙利,等.準噶爾盆地車排子凸起構造演化特征及其成因[J].石油與天然氣地質,2016,37(4):556-564.

Hu Qiuyuan,Dong Dawei,Zhao Li,et al.Tectonic evolutionary characteristics and their causes of Chepaizi Uplift in Junggar Basin[J].Oil & Gas Geology,2016,37(4):556-564.

[16] 何登發(fā),尹成,杜社寬,等.前陸沖斷帶構造分段特征:以準噶爾盆地西北緣斷裂構造帶為例[J].地學前緣,2004,11(3):91-101.

He Dengfa,Yin Cheng,Du Shekuan,et al.Characteristics of structural segmentation of foreland thrust belts:A case study of the fault belts in the northwestern margin of Junggar Basin[J].Earth Science Frontiers,2004,11(3):91-101.

(編輯 徐文明)

Tectonic evolution and its physical simulation of Wuxia fault belt in the northwestern margin of Junggar Basin

Li Xiao1, Lao Haigang2, Hu Qiuyuan1, Zhao Li3

(1.DepartmentofOil&GasEngineering,ShengliCollegeofChinaUniversityofPetroleum,Dongying,Shandong257000,China; 2.CollegeofMiningEngineering,NorthChinaUniversityofScienceandTechnology,Tangshan,Hebei063009,China; 3.WuxiResearchInstituteofPetroleumGeology,SINOPEC,Wuxi,Jiangsu214126,China)

The tectonic characteristics of the Wuxia fault belt in the northwestern margin of the Junggar Basin were systemically studied based on a large number of 3D seismic data. The study area has experienced five tectonic evolution periods, including the collision of foreland basin, thrust-nappe structure of foreland basin, thrusting fold of interior depression, uplift and denudation of interior depression, and cover setting. The forward modeling of the tectonic evolution of some typical profiles ever since the Late Carboniferous fit well with actual geological conditions. The results showed that the Wuxia fault belt transformed from the thrust with piggy back in foreland basin to thrusting with overstep in the interior depression. Experimental results also indicated that the formation and development of Wuxia fault belt were mainly controlled by the NW-SE oriented compressive stress field, which was mainly affected by some key factors such as plate activities, orogenic belts, and frontier basin-controlling faults.

tectonic characteristics; physical simulation experiment; Wuxia fault belt; northwestern margin of Junggar Basin

1001-6112(2017)04-0567-06

10.11781/sysydz201704567

2017-01-25；

2017-05-31。

李蕭(1987—)，碩士，助教，從事構造地質學研究。E-mail：136930394@qq.com。

勞海港(1978—)，博士，講師，從事盆地構造解析研究。E-mail:haigang1609@126.com。

國家自然科學基金項目(40772132)和中國石油大學勝利學院重點科研計劃項目(KY2015013)資助。

TE121 .2

A