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

1.中國(guó)地質(zhì)大學(xué)能源學(xué)院；2. 中國(guó)石油大學(xué)石油工程學(xué)院（青島）；3. 海洋水下設(shè)備試驗(yàn)與檢測(cè)技術(shù)國(guó)家工程實(shí)驗(yàn)室；4.中國(guó)石油大學(xué)石油工程學(xué)院（克拉瑪依）

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

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

1.中國(guó)地質(zhì)大學(xué)能源學(xué)院；2. 中國(guó)石油大學(xué)石油工程學(xué)院（青島）；3. 海洋水下設(shè)備試驗(yàn)與檢測(cè)技術(shù)國(guó)家工程實(shí)驗(yàn)室；4.中國(guó)石油大學(xué)石油工程學(xué)院（克拉瑪依）

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

蒙脫石；分子動(dòng)力學(xué)模擬；水化膨脹；頁(yè)巖抑制劑；甲酰胺；尿素

在石油工程技術(shù)發(fā)展過程中，如何維持井壁穩(wěn)定一直是困擾石油工程界的重大技術(shù)難題之一。井壁失穩(wěn)會(huì)造成井下復(fù)雜情況，導(dǎo)致井壁坍塌、儲(chǔ)層污染、影響后續(xù)作業(yè)實(shí)施、誘發(fā)各種嚴(yán)重鉆井事故等，甚至對(duì)經(jīng)濟(jì)財(cái)產(chǎn)和現(xiàn)場(chǎng)人員安全造成威脅［1］。因此，開展井壁穩(wěn)定性問題的研究對(duì)鉆井和油氣田開發(fā)都有十分重要的意義［2-3］。

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

1 模型的建立

Modeling

1.1 鈉蒙脫石

Na-montmorillonite

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

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

1.2 甲酰胺

Formamide

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

1.3 尿素

Urea

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

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

2 有機(jī)胺的抑制作用機(jī)理

Inhibition mechanism of organic amide

2.1 構(gòu)象

Conformation

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

（1）甲酰胺。向蒙脫石層間水分子中添加一個(gè)甲酰胺分子時(shí)的初始構(gòu)象和結(jié)構(gòu)優(yōu)化后的構(gòu)象分別如圖3、圖4所示。

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

從圖3可以看出，隨水分子數(shù)量增加，甲酰胺分子中原子所在平面與蒙脫石表面的傾角逐漸增大。蒙脫石層間只有單層水分子時(shí)，甲酰胺位于層水化層中，甲酰胺分子中原子所在平面與蒙脫石表面斜角較?。浑p層水分子時(shí)，甲酰胺位于2個(gè)水分子層中間，與蒙脫石表面的傾角增大，以分子中的氧原子靠近蒙脫石層間鈉離子，同一端的氫原子向蒙脫石表面靠近；當(dāng)為3層水分子時(shí)，甲酰胺分子近乎垂直于蒙脫石表面，并且分子中的氧原子及對(duì)應(yīng)距離最遠(yuǎn)的氫原子分別伸向蒙脫石表面。

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

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

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

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

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

2.2 氫鍵

Hydrogen bond

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

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

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

2.3 自擴(kuò)散系數(shù)

Self-diffusion coef fi cient

自擴(kuò)散系數(shù)，是指物體中某一點(diǎn)的溫度的擾動(dòng)傳遞到另一點(diǎn)的速率的量度?？杀硎緸?/p>

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

自擴(kuò)散系數(shù)Da也可以通過三維愛因斯坦方程計(jì)算得到，計(jì)算式如下

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

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

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

2.4 鈉離子的水化參數(shù)

Hydration parameter of sodion

徑向分布函數(shù)可以理解為系統(tǒng)的區(qū)域密度與平均密度的比，計(jì)算式為

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

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

根據(jù)離子配位數(shù)可進(jìn)一步計(jì)算出離子水化數(shù)和水化半徑

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

通過蒙脫石層間鈉離子Na+與水中氧O（w）的徑向分布函數(shù)，計(jì)算各水化層數(shù)條件下鈉離子的水化參數(shù)，見表2?？煽闯?，有機(jī)胺分子與層間水分子和蒙脫石表面形成的氫鍵，對(duì)鈉離子的水化起到了一定的抑制作用，使鈉離子周圍的水分子產(chǎn)生定向排列，各對(duì)應(yīng)水化層鈉離子的水化參數(shù)均有所降低。

2.5 力學(xué)參數(shù)

Mechanical parameter

提取有機(jī)胺加入前后蒙脫石的彈性力學(xué)參數(shù)，見表3。

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

表3 鈉蒙脫石體系彈性力學(xué)參數(shù)Table 3 Elasticity parameter of Na-montmorillonite system

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

3 有機(jī)胺濃度的影響

Effect of organic amide concentration

3.1 水化參數(shù)

Hydration parameter

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

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

3.2 自擴(kuò)散系數(shù)分析

Analysis on self-diffusion coef fi cient

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

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

3.3 力學(xué)參數(shù)分析

Analysis on mechanical parameter

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

圖10 蒙脫石的力學(xué)參數(shù)Fig. 10 Mechanical parameter of montmorillonite

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

4 結(jié)論

Conclusions

（1）有機(jī)胺主要通過與蒙脫石層間水分子和硅氧四面體表面的氧形成氫鍵來達(dá)到抑制蒙脫石水化的作用，同時(shí)限制了層間鈉離子和水分子的自擴(kuò)散運(yùn)動(dòng)，降低了鈉離子的水化性。

（2）有機(jī)胺的加入可以提高蒙脫石水化體系的力學(xué)性能，最佳作用濃度甲酰胺質(zhì)量分?jǐn)?shù)為5.21% ～13.02%，尿素質(zhì)量分?jǐn)?shù)為6.94%～5.63%。

［1］ 黃熠.南海高溫高壓勘探鉆井技術(shù)現(xiàn)狀及展望［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］ 閆傳梁. 硬脆性泥頁(yè)巖井壁失穩(wěn)機(jī)理研究［D］. 北京：中國(guó)石油大學(xué)（北京），2014.YAN Chuanliang. Borehole instability mechanism of hardbrittle shale［D］. Beijing: China University of Petroleum （Beijing）, 2014.

［3］ 鐘漢毅. 聚胺強(qiáng)抑制劑研制及其作用機(jī)理研究［D］.青島：中國(guó)石油大學(xué)（華東），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］ 黃桃，蒲曉林，羅霄，于浩，徐生江 . 基于流固耦合理論的泥頁(yè)巖井壁穩(wěn)定性分析［J］. 石油天然氣學(xué)報(bào)，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］ 趙寶祥，陳江華，徐一龍，鄭浩鵬，殷水平.北部灣盆地開發(fā)井表層鉆井模式［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ī)鹽抑制蒙脫石水化機(jī)理［J］. 石油學(xué)報(bào)，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］ 王進(jìn). 蒙脫石層間結(jié)構(gòu)的分子力學(xué)和分子動(dòng)力學(xué)模擬研究［D］. 太原：太原理工大學(xué)，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］ 王進(jìn)，王軍霞，曾凡桂，吳秀玲. 蒙脫石晶體結(jié)構(gòu)構(gòu)型及其 XRD、IR 的分子模擬［J］. 礦物學(xué)報(bào)，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］ 袁俊生，包捷. 鉀、鈉、氯離子水化現(xiàn)象的分子動(dòng)力學(xué)模擬［J］. 計(jì)算機(jī)與應(yīng)用化學(xué)，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ī)胺抑制蒙脫石水化機(jī)理的分子模擬［J］.石油鉆采工藝，2017，39（4）：442-448.

TE254.4< class="emphasis_bold">文獻(xiàn)標(biāo)識(shí)碼；

；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.

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃973項(xiàng)目“海洋深水油氣安全高效鉆完井基礎(chǔ)研究”（編號(hào)：2015CB251200）；油氣重大專項(xiàng)“大型油氣田及煤層氣開發(fā)”（編號(hào)：2016ZX05044003-001）。

黃小娟（1973-），1995年畢業(yè)于江漢石油學(xué)院，現(xiàn)主要研究方向?yàn)橛蜌馓镩_發(fā)與工程應(yīng)用，在讀博士生。通訊地址：（100083）北京海淀區(qū)學(xué)院路29號(hào)能源學(xué)院。E-mail：huangxj@cugb.edu.cn

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