有機胺抑制蒙脫石水化機理的分子模擬

黃小娟徐加放丁廷稷李小迪，4李力文

1.中國地質大學能源學院；2. 中國石油大學石油工程學院（青島）；3. 海洋水下設備試驗與檢測技術國家工程實驗室；4.中國石油大學石油工程學院（克拉瑪依）

有機胺抑制蒙脫石水化機理的分子模擬

黃小娟1徐加放2，3丁廷稷2李小迪2，4李力文2

1.中國地質大學能源學院；2. 中國石油大學石油工程學院（青島）；3. 海洋水下設備試驗與檢測技術國家工程實驗室；4.中國石油大學石油工程學院（克拉瑪依）

井壁失穩一直是困擾石油鉆井的重大技術難題，泥頁巖水化是造成井壁失穩的最重要因素之一，作為最易吸水和泥頁巖主要組分的蒙脫石是引起井壁失穩的根本性因素。利用分子模擬的方法，建立了有機胺類抑制劑與鈉蒙脫石反應的分子結構模型，研究了甲酰胺和尿素等小分子有機胺類抑制劑對蒙脫石水化的抑制作用機理。胺類有機小分子抑制蒙脫石水化的能力與分子結構及官能團有關，甲酰胺和尿素均可與層間水分子及蒙脫石表面的氧原子形成氫鍵，抑制層間鈉離子的水化，限制層間水分子及鈉離子的自擴散。甲酰胺和尿素分子在不同濃度時對鈉蒙脫石體系的密度和力學性能影響不同，但有相同的變化規律，即存在一臨界濃度，低于此濃度時，胺類分子的加入可增大蒙脫石的密度，提高體系的力學性能，而高于此濃度時，體系密度減小，力學性能變差。甲酰胺和尿素都是在低濃度時能更好地發揮抑制作用，提高鈉蒙脫石的穩定性，綜合各參數對比分析，甲酰胺比尿素的整體作用效果好。

蒙脫石；分子動力學模擬；水化膨脹；頁巖抑制劑；甲酰胺；尿素

在石油工程技術發展過程中，如何維持井壁穩定一直是困擾石油工程界的重大技術難題之一。井壁失穩會造成井下復雜情況，導致井壁坍塌、儲層污染、影響后續作業實施、誘發各種嚴重鉆井事故等，甚至對經濟財產和現場人員安全造成威脅［1］。因此，開展井壁穩定性問題的研究對鉆井和油氣田開發都有十分重要的意義［2-3］。

井壁失穩是力學和物理化學問題，造成井壁失穩的原因主要有力學因素、工程施工因素和泥頁巖水化膨脹因素［4-6］。保持井壁穩定的最有效途徑之一是在鉆井液中添加頁巖抑制劑［7］，用來抑制泥頁巖中黏土礦物的水化膨脹及分散，防止井壁坍塌。目前油田常用的頁巖抑制劑有硅酸鹽、鉀鹽、銨鹽以及大分子聚合物等，其中胺基鉆井液是國內外近幾年發展起來的新型水基鉆井液，較適用于水敏性泥頁巖地層［8-10］。胺基頁巖抑制劑分子鏈中由于引入了胺基官能團，因此它的分子結構較為獨特。胺基頁巖抑制劑作用機理有：（1）電荷中和作用；（2）吸附作用；（3）形成半透膜。但當前對于有機胺頁巖抑制劑研究仍主要集中在宏觀實驗的基礎上，微觀分子層面上的研究十分有限。筆者從井壁穩定性角度出發，研究了甲酰胺和尿素作為頁巖抑制劑抑制泥頁巖水化的作用效果及作用機理，從理論上優選作用效果更好的抑制劑種類及相應的適宜濃度范圍，為頁巖抑制劑在現場的應用提供一定的理論指導，同時研究了泥頁巖水化及抑制劑作用下的力學參數變化，實現了力學和化學的耦合。

1 模型的建立

Modeling

1.1 鈉蒙脫石

Na-montmorillonite

運用Materials Studio中的Crystal Builder模塊建立蒙脫石模型，空間群選擇C2/m，對稱型選擇L2PC；單晶胞參數a=0.523 nm，b=0.906 nm，c值可隨層間水分子含量及層間陽離子的類型和數量不同而出現變化［11］；根據從文獻中查到的晶胞相關參數的實驗值［12］建立晶胞單元，考慮到晶體層間距和體系帶電量及范德華力等客觀條件的限制，超晶胞由4a×2b×1c個單位晶胞組成；根據Marry V［13］等關于鈉蒙脫石的研究中所建模型的電荷設置，賦予晶胞中各原子相應的電荷；理想的蒙脫石結構中晶格取代規則［14］，在鋁氧八面體層中有4個Al3+被Mg2+取代，硅氧四面體層有2個Si4+被Al3+取代，取代位隨機分布。模擬應用三維周期性邊界條件。最后在超晶胞中充填6個Na+以補償晶格取代產生的負電荷，形成鈉蒙脫石多晶胞模型，如圖1，其中白色為氫，紅色為氧，綠色為鎂，黃色為硅，紫色為鋁，中間球狀為鈉。

圖1 鈉蒙脫石超晶胞模型Fig. 1 Super-cell model of Na-montmorillonite

1.2 甲酰胺

Formamide

用MS軟件建立甲酰胺分子模型，如圖2（a）所示。模型中甲酰胺分子的幾何結構為C O鍵長 RCO=1.214?，C H 鍵 長 RCH=1.085?，N—C鍵 長RCH=1.356 ?，N—H鍵長RCH=1.044 ?。N—C—O鍵角為 119.481°，N—C—H 鍵角為 120.745°，C—N—H鍵角為 119.420°。

1.3 尿素

Urea

通過MS軟件構建尿素模型，模型中尿素分子的幾何結構為C O鍵長RCO=1.212?，N—C鍵長RNC=1.349 ?，N—H 鍵長 RNH=1.041 ?。N—C—O 鍵角為 119.916°，N—C—N 鍵角為 120.167°，H—N—H鍵角為10.123°，C—N—H鍵角為118.99°，如圖2。對構建得到的尿素分子模型進行結構優化，并將優化后不同數量的尿素分子與水分子同時添加到鈉蒙脫石層間，進行體系結構優化后進行1 000 ps的分子動力學模擬，主要分不同水化層添加一個尿素分子、不同水化層添加不同濃度的尿素以及不同水化層相同質量濃度尿素時的對比分析。

圖2 甲酰胺與尿素分子模型Fig. 2 Molecular model of formamide and urea

2 有機胺的抑制作用機理

Inhibition mechanism of organic amide

2.1 構象

Conformation

隨著蒙脫石層間水分子數量的增加，蒙脫石的層間距逐漸增大，水分子也逐漸開始分層，蒙脫石的層間水分子通?？煞譃?層［9］。有機分子受層間距的限制，其存在狀態會隨層間距的大小而改變。

（1）甲酰胺。向蒙脫石層間水分子中添加一個甲酰胺分子時的初始構象和結構優化后的構象分別如圖3、圖4所示。

圖3 含不同水化層鈉蒙脫石添加1個甲酰胺分子構象Fig. 3 Conformation of adding one formamide molecule in Namontmorillonite with different hydration layers

從圖3可以看出，隨水分子數量增加，甲酰胺分子中原子所在平面與蒙脫石表面的傾角逐漸增大。蒙脫石層間只有單層水分子時，甲酰胺位于層水化層中，甲酰胺分子中原子所在平面與蒙脫石表面斜角較??；雙層水分子時，甲酰胺位于2個水分子層中間，與蒙脫石表面的傾角增大，以分子中的氧原子靠近蒙脫石層間鈉離子，同一端的氫原子向蒙脫石表面靠近；當為3層水分子時，甲酰胺分子近乎垂直于蒙脫石表面，并且分子中的氧原子及對應距離最遠的氫原子分別伸向蒙脫石表面。

圖4 1個甲酰胺在2層水分子的鈉蒙脫石中的動力學過程構象Fig. 4 Conformation of dynamic process of one formamide molecule in Na-montmorillonite with two layers’ hydrone

如圖4所示，隨著時間的推移［t=0時刻即圖3（b）］，鈉蒙脫石層間鈉離子始終處于兩層水分子中間，水分子由明顯分層擴散至分層界限不再明顯，部分水分子進一步靠近蒙脫石層表面，甲酰胺則向最鄰近蒙脫石表面擴散，懸掛在硅氧四面體片的六氧元環的上方，以羰基中的氫靠近蒙脫石表面，其擴散受六氧元環作用而受到限制。

（2）尿素。尿素分子在蒙脫石層間的初始構象和結構優化后的構象分別如圖5、圖6所示。從圖5中可以看出，單層水分子時，尿素位于層中間位置，呈平鋪狀，尿素原子所在平面與蒙脫石表面夾角很??；雙層和三層水分子，尿素主要位于其中的單層水分子中，并且靠近蒙脫石表面，尿素原子所在平面與蒙脫石表面有一定的傾角，但角度小于甲酰胺。圖6表明，整個過程水分子不斷向四周擴散，兩層水分子分層逐漸模糊；層間鈉離子主要在層間運動，有向蒙脫石表面擴散的趨勢，但仍主要位于兩水分子層中間；尿素分子由中間位置向一側擴散，最后靠近蒙脫石層表面，將周圍的水分子擠開。由俯視圖［圖6（d）］可以看出，尿素以其中一個NH2基吸附在蒙脫石六氧元環上方，與K+鑲嵌機理類似（NH2基團半徑與六氧元環半徑相近）。

圖5 含不同水化層鈉蒙脫石添加1個尿素分子構象Fig. 5 Conformation of adding one urea molecule in Namontmorillonite with different hydration layers

圖6 1個尿素分子在2層水分子的鈉蒙脫石中的動力學過程構象Fig. 6 Conformation of dynamic process of one formamide molecule in Na-montmorillonite with two layers’ hydrone

2.2 氫鍵

Hydrogen bond

與無機鹽抑制蒙脫石水化機理不同［15-17］，有機胺通過氫鍵與蒙脫石表面和層間水分子作用，從而抑制蒙脫石的水化。由圖7可以看出，有機胺中的O和NH2基中的H可以分別與蒙脫石硅氧四面體中的O和層間水分子的H之間形成氫鍵，同時影響了層間水分子的排列。

圖7 蒙脫石層間氫鍵示意圖Fig. 7 Schematic hydrogen bond between montmorillonite layers

通過計算，未加入有機胺抑制劑時，蒙脫石在1、2、3層水化條件下對應的氫鍵個數為8、31、65；一個甲酰胺分子存在時，對應的1、2、3層水分子蒙脫石層間的氫鍵個數為16，47，83；一個尿素分子為17、42、78，有機胺的加入增加了蒙脫石層間的氫鍵數量，增強了晶層間的吸附。

2.3 自擴散系數

Self-diffusion coef fi cient

自擴散系數，是指物體中某一點的溫度的擾動傳遞到另一點的速率的量度?？杀硎緸?/p>

式中，λ為熱導率，W/（m·K）；ρ為物體的密度，g/cm3；Cp為物體的定壓比熱容，J/（kg·K）。

自擴散系數Da也可以通過三維愛因斯坦方程計算得到，計算式如下

自擴散系數為均方位移曲線平衡點斜率的1/6。計算加入甲酰胺和尿素前后蒙脫石層間鈉離子和水分子的自擴散系數，如表1。

從表1可以看出，與未加有機胺抑制劑時相比，加入甲酰胺和尿素后，鈉離子和水分子對應相同水化層時的自擴散系數都有減小，表明有機胺對蒙脫石層間鈉離子和水分子的自擴散有一定的限制作用。通過對比可知，甲酰胺在蒙脫石初始水化時的抑制能力較強，而尿素更適合抑制蒙脫石后期水化。

表1 有機胺加入前后蒙脫石層間水分子和鈉離子的自擴散系數Table 1 Self-diffusion coeffcient of hydrogen and sodion between montmorillonite layers before and after the addition of organic amide

2.4 鈉離子的水化參數

Hydration parameter of sodion

徑向分布函數可以理解為系統的區域密度與平均密度的比，計算式為

式中，gαβ(r)為β粒子的分布機率；nβ為半徑為r+dr中圍繞α粒子的β粒子個數；ρβ為粒子β的密度；r為α粒子與β粒子之間的距離。

可通過VS軟件編程計算出Na+的水分子配位數，對應含有一、二、三個水化層的鈉蒙脫石中鈉離子配位數分別為9.47、7.67、6.07?？梢钥闯?，隨著層間水分子個數的增多，鈉離子配位數呈現下降趨勢，這是因為離子配位數與體系的溫度、濃度相關，低溫、高濃度有利于形成高配位數。

根據離子配位數可進一步計算出離子水化數和水化半徑

式中，nhyd為離子水化數；h為離子水化因子，其中，Na+水化因子為0.69；rhyd為離子水化半徑，nm；v為1個水分子的體積，2.99×10-29m3；reff為離子有效半徑，定義為徑向分布函數的第一峰位與水的有效半徑之差，其中水分子的有效半徑為0.138 nm。

通過蒙脫石層間鈉離子Na+與水中氧O（w）的徑向分布函數，計算各水化層數條件下鈉離子的水化參數，見表2。可看出，有機胺分子與層間水分子和蒙脫石表面形成的氫鍵，對鈉離子的水化起到了一定的抑制作用，使鈉離子周圍的水分子產生定向排列，各對應水化層鈉離子的水化參數均有所降低。

2.5 力學參數

Mechanical parameter

提取有機胺加入前后蒙脫石的彈性力學參數，見表3。

表2 加入有機胺前后蒙脫石層間鈉離子的水化參數Table 2 Hydration parameter of sodion between montmorillonite layers before and after the addition of organic amide

表3 鈉蒙脫石體系彈性力學參數Table 3 Elasticity parameter of Na-montmorillonite system

從表3可以看出，隨著層間水化層數的增加，蒙脫石的體積模量、切變模量和楊氏模量逐漸減小，壓縮系數逐漸增大，鈉蒙脫石的穩定性逐漸變差。但與未加有機胺時相比，有機胺（甲酰胺和尿素）的加入明顯提高了彈性模量，降低了壓縮系數，提高了蒙脫石體系的穩定性。

3 有機胺濃度的影響

Effect of organic amide concentration

3.1 水化參數

Hydration parameter

由圖8可以看出，相同甲酰胺濃度時，層間鈉離子的水化參數（配位數、水化數和水化半徑）隨層間距增加而逐漸減小；隨著甲酰胺濃度的增加，一個水化層時鈉離子的水化參數先減小、后增大，2個和3個水化層時水化參數整體趨于減小。表明隨甲酰胺的濃度增加，對鈉離子的水化抑制增強。

圖8 蒙脫石層間鈉離子的水化參數Fig. 8 Hydration parameter of sodion between montmorillonite layers

3.2 自擴散系數分析

Analysis on self-diffusion coef fi cient

圖9為蒙脫石層間粒子的自擴散系數。由圖9可以看出，相同甲酰胺濃度時，蒙脫石層間鈉離子、水分子和甲酰胺的自擴散系數隨著水化層數的增加而逐漸增大；隨著甲酰胺濃度的增加，三者的自擴散系數均表現出先減小、后增大的趨勢。

圖9 蒙脫石層間粒子的自擴散系數Fig. 9 Self-diffusion coeffcient of particles between montmorillonite layers

3.3 力學參數分析

Analysis on mechanical parameter

圖10為甲酰胺加入前后蒙脫石體系的力學參數變化。從圖10中可以看出，同濃度甲酰胺條件下，隨蒙脫石水化層數的增加，其體積模量、剪切模量和楊氏模量均逐漸減小，壓縮系數則逐漸增大；隨甲酰胺濃度的增加，體積模量、剪切模量和楊氏模量先增大后減小，壓縮系數為先減小后增大。

圖10 蒙脫石的力學參數Fig. 10 Mechanical parameter of montmorillonite

綜上可以看出，甲酰胺對蒙脫石水化的抑制作用存在一最佳濃度，即質量分數為5.21%～13.02%，而并非越大越好。研究表明，尿素具有同樣的作用效果和作用規律，經計算尿素的最佳作用濃度為6.94% ～15.63%。

4 結論

Conclusions

（1）有機胺主要通過與蒙脫石層間水分子和硅氧四面體表面的氧形成氫鍵來達到抑制蒙脫石水化的作用，同時限制了層間鈉離子和水分子的自擴散運動，降低了鈉離子的水化性。

（2）有機胺的加入可以提高蒙脫石水化體系的力學性能，最佳作用濃度甲酰胺質量分數為5.21% ～13.02%，尿素質量分數為6.94%～5.63%。

［1］ 黃熠.南海高溫高壓勘探鉆井技術現狀及展望［J］.石油鉆采工藝，2016，38（6）：737-745.HUANG Yi. Drilling technology for HTHP exploration in South China Sea and its prospect［J］. Oil Drilling &Produc-tion Technology, 2016, 38(6): 737-745.

［2］ 閆傳梁. 硬脆性泥頁巖井壁失穩機理研究［D］. 北京：中國石油大學（北京），2014.YAN Chuanliang. Borehole instability mechanism of hardbrittle shale［D］. Beijing: China University of Petroleum （Beijing）, 2014.

［3］ 鐘漢毅. 聚胺強抑制劑研制及其作用機理研究［D］.青島：中國石油大學（華東），2012.ZHONG Hanyi. Development and mechanism study on high performance polyamine inhibitor in waterbased drilling fuid［D］. Qingdao: China University of Petroleum （East China）, 2012

［4］ BARTEZKO A, WESSLING S, DAHLT PEIJ. Challenges for wellbore stability workfow automation［J］. Oil Gas European Magazine, 2014, 37（2）: 65-71.

［5］ 黃桃，蒲曉林，羅霄，于浩，徐生江 . 基于流固耦合理論的泥頁巖井壁穩定性分析［J］. 石油天然氣學報，2014，36 （10）：117-126.HUANG Tao, PU Xiaolin, LUO Xiao, YU Hao, XU Shengjiang. Borehole stability analysis of shale wells based on the theory of fuid-solid coupling［J］. Journal of Oil and Gas Technology, 2014, 36 （10）: 117-126.

［6］ 趙寶祥，陳江華，徐一龍，鄭浩鵬，殷水平.北部灣盆地開發井表層鉆井模式［J］.石油鉆采工藝，2015，37（6）：23-26.ZHAO Baoxiang, CHEN Jianghua , XU Yilong, ZHENG Haopeng, YIN Shuiping .The pattern of surface drilling in development wells of Beibu Gulf Basin［J］. Oil Drilling& Production Technology, 2015, 37（6）: 23-26.

［7］ Lü Kaihe, ZHONG Hanyi, REN Guanlong, LIU Y.Properties evaluation and application of organic amine inhibitor on the properties of drilling fluids［J］. Open Petroleum Engineering Journal, 2014, 7（1）: 50-54.

［8］ ABASS H, SHEBATALHAMD A. Wellbore instability of Shale Formation, Zuluf feld, Saudi Arabia［R］. SPE 106345, 2006.

［9］ 徐加放，孫澤寧，劉洪軍，付顯威，孫中富，李小迪. 分子模擬無機鹽抑制蒙脫石水化機理［J］. 石油學報，2014，35（2）：377-384.XU Jiafang, SUN Zening, LIU Hongjun, FU Xianwei, SUN Zhongfu, LI Xiaodi. Molecular simulation for inorganic salts inhibition mechanism onmontmorillonitehydration［J］.Acta Petrolei Sinica, 2014, 35（2）: 377-384.

［10］ GHOLIZADEH-DOONECHALY N, TAHMASBI K,DAVANI E. Development of high-performance waterbased mud formulation based on amine derivatives［R］.SPE 121228, 2009.

［11］ 王進. 蒙脫石層間結構的分子力學和分子動力學模擬研究［D］. 太原：太原理工大學，2005.WANG Jin. Molecular mechanics and dynamic simulation study of Montmorillonite interlayer structure［D］.Taiyuan: Taiyuan Science & Technology University, 2005.

［12］ BOEK E S. Molecular dynamics simulations of interlayer structure and mobility in hydrated Li-, Na-and K-montmorilloniteclays［J］. Molecular Physics, 2014,112（9/10）: 1472-1483.

［13］ MARRYV, TURQP, CARTAILLER T, LEVESQUED.Microscopic simulation of structure and dynamics of water and counterions in a monohydrated montmorillonite［J］. The Journal of Chemical Physics,2002, 117（7）: 3454-3463.

［14］ YOUNG D A, SMITH D E. Simulations of clay mineral swelling and hydration: dependence upon interlayer ion size and charge［J］. The Journal of Physical Chemistry B, 2000, 104（39）: 9163-9170.

［15］ 王進，王軍霞，曾凡桂，吳秀玲. 蒙脫石晶體結構構型及其 XRD、IR 的分子模擬［J］. 礦物學報，2011，31（1）：133-138.WANG Jin, WANG Junxia, ZENG Fangui, WU Xiuling. Molecular simulations of crystal structurecon formation, X-Ray diffraction and infra-red spectrum in montmorillonites［J］. Acta Mineralogica Sinica, 2011,31（1）: 133-138.

［16］ 袁俊生，包捷. 鉀、鈉、氯離子水化現象的分子動力學模擬［J］. 計算機與應用化學，2009，26（10）：1295-1299.YUAN Junsheng, BAO Jie. Molecular dynamics simulation of K+, Na+and Cl-hydration［J］. Computers and Applied Chemistry, 2009, 26（10）: 1295-1299.

［17］ LüG S, CHANGP H, XINGX B, JIANG W T, JEAN J S,LI Z H. Investigation of intercalation of diphenhydramine into the interlayer of smectite by XRD, FTIR, TG-DTG analyses and molecular simulation［J］. Arabian Journal of Chemistry, 2015, 70（1）: 1-8.

（修改稿收到日期 2017-06-27）

〔編輯 薛改珍〕

Molecular simulation on the inhibition mechanism of organic amine to montmorillonite hydration

HUANG Xiaojuan1, XU Jiafang2,3, DING Tingji2, LI Xiaodi2,4, LI Liwen2
1. School of Energy Resources, China University of Geosciences，Beijing 100083, China; 2. School of Petroleum Engineering, China Petroleum University, Qingdao 266580, Shandong, China; 3. National Engineering Laboratory for Subsea Equipment Testing and Detection Technology, Qingdao 266580, Shandong, China; 4. School of Petroleum Engineering, China Petroleum University, Karamay 834000, Xinjiang, China

Wellbore instability has always been the major diffculty in the process of well drilling, and mud shale hydration is one of the most important factors leading to wellbore instability. Montmorillonite, as the composition with the highest water absorption and the main component in mud shale, is the radical factor of wellbore instability. In this paper, the molecular structure model for the reaction between organic amine inhibitors and Na-montmorillonite was established by means of molecular simulation to study the inhibition mechanisms of small-molecule organic amine inhibitors (e.g. formamide and urea) to montmorillonite hydration. It is shown that the inhibition capacity of small-molecule organic amine mainly depends on its molecular structure and functional group. Formamide andurea can form hydrogen bond with the interlayer hydrone and the oxygen on the surface of montmorillonite to inhibit the hydration of interlayer sodion and limit the self-diffusion of interlayer hydrone and sodion. When the concentrations of formamide and urea change,their effects on the density and mechanical performance of Na-montmorillonite are different, but they follow the same change laws.There is a certain critical concentration, below which the addition of amine can increase the density of montmorillonite, increasing the mechanical performance of this system, and vice versa. And at low concentration, they both can present good inhibition and increase the stability of montmorillonite. Based on the parameter analysis and comparison, formamide is, on the whole more effective than urea.

montmorillonite; molecular dynamics simulation; hydrate expansion; shale inhibitor; formamide; urea

黃小娟，徐加放，丁廷稷，李小迪，李力文.有機胺抑制蒙脫石水化機理的分子模擬［J］.石油鉆采工藝，2017，39（4）：442-448.

TE254.4< class="emphasis_bold">文獻標識碼；

；A

1000 – 7393（ 2017 ） 04 – 0442 – 07

10.13639/j.odpt.2017.04.009

:HUANG Xiaojuan, XU Jiafang, DING Tingji, LI Xiaodi, LI Liwen. Molecular simulation on the inhibition mechanism of organic amine to montmorillonite hydration［J］. Oil Drilling & Production Technology, 2017, 39(4): 442-448.

國家重點基礎研究發展計劃973項目“海洋深水油氣安全高效鉆完井基礎研究”（編號：2015CB251200）；油氣重大專項“大型油氣田及煤層氣開發”（編號：2016ZX05044003-001）。

黃小娟（1973-），1995年畢業于江漢石油學院，現主要研究方向為油氣田開發與工程應用，在讀博士生。通訊地址：（100083）北京海淀區學院路29號能源學院。E-mail：huangxj@cugb.edu.cn

徐加放（1973-），1998年畢業于中國石油大學（華東），現主要從事水合物開發與防治、井壁穩定、深水鉆井液及油氣層保護、鉆井廢棄物處理與資源化利用等教學科研工作。通訊地址：（266580）山東省青島市黃島區長江西路66號石油工程學院。E-mail：xjiafang@upc.edu.cn