黃 金， 于洪洋， 李 鋒， 張 宇， 張 坤， 王 俊， 李翠勤

(1.東北石油大學 化學化工學院 石油與天然氣化工省重點實驗室，黑龍江 大慶 163318；2.大慶師范學院 油田應用化學與技術省重點實驗室，黑龍江 大慶 163712)

隨著對能源需求的不斷增加以及常規(guī)油氣資源開發(fā)難度的日益增大，致密油藏和頁巖油藏等非常規(guī)油氣資源受到廣泛關注[1]。然而，這些非常規(guī)油藏油層物性差、滲透率低、縱向分布零散、非均質嚴重、黏土含量高，表面活性劑吸附損耗嚴重[2]。傳統(tǒng)驅油用的烷基苯磺酸鹽和石油磺酸鹽陰離子型表面活性劑由于具有耐鹽耐溫性差、吸附損耗大等缺點[3]，難以滿足非常規(guī)油氣資源的開采。

兩性離子型表面活性劑隨所在體系酸堿性的變化呈現不同的離子特性，不易受無機電解質的影響，不易形成疏水表面，相較于傳統(tǒng)表面活性劑具有更好的耐高溫和耐鹽性能，作為一種驅油助劑廣泛應用于高溫高鹽油藏地區(qū)[4]。目前國內對甜菜堿型兩性表面活性劑的研究十分踴躍，但由于壬基酚甜菜堿兩性表面活性劑和辛基酚甜菜堿型兩性表面活性劑等石油基甜菜堿型表面活性劑的原料對環(huán)境污染較大，生物基表面活性劑成為該領域研究的熱點之一[5-7]。Moslemizadeh等[8]以桑樹葉為原料合成了一種新型生物基表面活性劑，該表面活性劑作為驅油劑可提高原油采收率達16.4%。Bian等[9]以可再生油酸為基本原料制備了一種新型含有支鏈烷基團的生物基苯磺酸鹽表面活性劑，該表面活性劑具有較好的界面活性，可將水的表面張力降低至32.54 mN/m，將油/水界面張力降低至10-2mN/m。生物基表面活性劑因其可再生性、低毒性、易生物降解性及較高的乳化性能和界面活性受到廣泛關注，成為未來表面活性劑的主要發(fā)展方向之一[10]。

以筆者所在課題組合成的一類生物基兩性離子型表面活性劑(結構如圖1所示)為研究對象，通過實驗和分子模擬研究了該生物基兩性離子型表面活性劑的界面特性，分析其結構對油/水界面特性的影響以及對不同體系的適應性。重點采用分水體積法研究了該類生物基兩性離子型表面活性劑的乳化穩(wěn)定性動力學，考察了乳化條件和體系組成對模擬乳狀液穩(wěn)定性動力學的影響，確定了該生物基兩性離子型表面活性劑形成穩(wěn)定乳狀液的最佳體系和組成。

圖1 生物基兩性離子型表面活性劑的結構Fig.1 Chemical structure of bio-based zwitterionic surfactant

1 實驗部分

1.1 原料和試劑

氯化鎂、氯化鈣、氯化鈉和氫氧化鈉，均為分析純，科密歐化學試劑有限公司產品；部分水解聚丙烯酰胺(HPAM，黏均相對分子質量為2.5×107)，中國石油大慶煉化公司提供；生物基兩性離子型表面活性劑(CNBS)，實驗室自制[11]；脫水原油，大慶油田有限責任公司第四采油廠提供，20 ℃標準密度為0.8557 g/cm3，45 ℃運動黏度為20.208 mm2/s，凝點為22.75 ℃。

1.2 界面張力的測定

將氯化鈉、氯化鎂和氯化鈣按照質量比60∶3∶4配制一定礦化度的模擬注入水，然后采用模擬注入水配制不同濃度的系列生物基兩性離子型表面活性劑水溶液，采用TX-500型全自動界面張力儀測定系列生物基兩性離子型表面活性劑水溶液的界面張力，測試條件：溫度為45 ℃，轉速為5000 r/min。

1.3 乳化穩(wěn)定動力學參數的測定

按照表1的實驗方案配制含生物基兩性離子型表面活性劑的水溶液，將溶液與脫水原油置于45 ℃的恒溫水浴中預熱30 min；采用FA25型高速剪切分散乳化機以10000 r/min對含生物基兩性離子型表面活性劑溶液剪切1 min，然后與預熱后的脫水原油按一定比例混合，在45 ℃下以10000 r/min剪切乳化5 min，獲得模擬乳狀液。將模擬乳狀液置于100 mL具塞量筒中，45 ℃恒溫靜止并每隔5 min記錄析出水體積，連續(xù)記錄至90 min。設定φA為t時刻油相含水率(體積分數，下同)，根據乳狀液破壞過程(液滴絮凝，聚結，油、水分離) 及乳狀液破乳的一級速率方程(式(1))[12]，通過積分得到積分速率方程(式(2))，將乳狀液破乳一半(1/2φA0)的時間定義為半衰期T1/2(式(3))。根據實驗測得不同時間下乳狀液析出水的體積，結合式(2)和式(3)，得到不同體系和不同條件下形成乳狀液的乳化穩(wěn)定動力學參數：破乳速率常數hc、半衰期T1/2、90 min時油相的含水率φA90。

v=dφA/dt=hcφA

(1)

(2)

T1/2=ln2/hc

(3)

式中：v為乳狀液的破乳速率，%/min；φA為不同時間下油相的含水率(體積分數)，%；φA0為乳狀液中初始含水率(體積分數)，%；φA90為90 min時油相的含水率(體積分數)，%；hc為破乳速率常數，min-1；t為靜沉時間，min；T1/2為半衰期，min。

表1 乳化穩(wěn)定動力學研究的實驗方案Table 1 Experimental scheme for the research of emulsion-stabilizing kinetics

1.4 界面微觀特性分子模擬方法

生物基兩性離子型表面活性劑在油/水界面的微觀特性模擬采用全原子分子動力學模擬技術，全部模擬計算使用Gromacs2020程序[13]完成，力場選擇Gromos54a7力場，以十一烷代表油[14]。生物基表面活性劑分子使用Automated topology builder version 3.0進行結構優(yōu)化并獲得分子模擬用文件[15]。使用Gromacs命令構建5 nm×5 nm×9 nm的盒子(見圖2)，沿y軸劃分為2 nm(油)+5 nm(水)+2 nm(油)共3個區(qū)域；油區(qū)域填充入十一烷，水區(qū)域填充水獲得油+水體系(如圖2(a)所示)；水區(qū)域填充水和生物基兩性離子型表面活性劑獲得油+水+CNBS體系(如圖2(b)所示)。對油/水、油/水/CNBS、油/水/CNBS/NaCl和油/水/CNBS/NaCl/NaOH等4個體系進行模擬，各體系中十一烷分子數量為127個，除油/水體系外其他體系中CNBS分子數量為4個，油/水/CNBS/NaCl體系和油/水/CNBS/NaCl/NaOH體系中按照質量濃度32 g/L加入NaCl(根據界面活性實驗結果)，油/水/CNBS/NaCl/NaOH體系按照質量分數1.00%加入NaOH(根據界面活性實驗結果)。分子模擬全過程溫度為45 ℃，先進行能量最小化優(yōu)化，然后進行壓力為0.1 MPa的等溫等壓(NPT)系綜模擬，模擬總時長為300 ps，時間步長為2 fs。根據模擬后體系獲得的密度曲線，計算油/水界面厚度、總能量Etotal、油相能量Eoil和水相能量Ewater，根據式(4)計算油/水界面能Einterface，能量單位均為kJ/mol。

Einterface=Etotal-Eoil-Ewater

(4)

圖2 分子模擬體系盒子Fig.2 Box of molecular simulation system(a) Oil+Water system； (b) Oil+Water+CNBS system

2 結果與討論

2.1 乳化體系組成對生物基兩性離子型表面活性劑界面活性的影響

固定總礦化度為5 g/L，配制不同生物基兩性離子型表面活性劑濃度的溶液，按照1.2節(jié)中方法檢測界面張力并繪制曲線，見圖3(a)中曲線(1)。由圖3(a)中曲線(1)可知，單一的生物基兩性離子型表面活性劑體系，油/水界面張力未達到超低界面張力(1×10-3mN/m)。當體系中生物基兩性離子型表面活性劑質量濃度低于1.5 g/L時，隨其質量濃度增大，油/水界面張力下降；當表面活性劑質量濃度達到1.5 g/L時，繼續(xù)增大表面活性劑濃度，界面張力變化很小。這主要是因為，較低濃度的兩親性表面活性劑分子在油/水界面上排布比較疏松；隨著濃度增大，吸附在油/水界面的表面活性劑分子排布比較緊密；當表面活性劑質量濃度接近1.5 g/L時，表面活性劑在油/水界面達到飽和吸附。且通過圖3(a)中曲線(1)得到該生物基兩性離子型表面活性劑在該體系中的臨界膠束質量濃度為1.29 g/L，此濃度下的油/水界面張力為3.28×10-2mN/m。與課題組合成的單一陰離子型生物基表面活性劑(3.68 g/L)[16]相比，由于該生物基兩性離子型表面活性劑雙親水基電荷的中和作用，使其在較低的質量濃度下(1.29 g/L)就能在油/水界面達到飽和吸附。

固定生物基兩性離子型表面活性劑的質量濃度為1.29 g/L，配制不同礦化度的生物基兩性離子型表面活性劑溶液，按照1.2節(jié)中測定界面張力并繪制曲線，見圖3(a)中曲線(2)。由圖3(a)中曲線(2)可知：隨著體系礦化度增大，生物基兩性離子型表面活性劑的界面活性增加；但當體系礦化度達到32 g/L時，繼續(xù)增大體系的礦化度，界面張力變化較小。這主要是由于該生物基表面活性劑為兩性離子型表面活性劑，分子中含有陽離子親水基團和陰離子親水基團；隨著體系礦化度的增大，體系中無機鹽在水中電離的陽離子和陰離子濃度增大，這些離子與表面活性劑分子中的離子相互作用，使得表面活性劑分子中親水基團之間的電荷斥力作用減弱，表面活性劑分子在油/水界面排列緊密，界面張力降低，界面活性增大；但當體系中的礦化度增大到一定值時，表面活性劑親水基團附近與其所帶電荷相反的無機離子濃度增大到一定值時，電荷之間的作用達到平衡，繼續(xù)增大體系礦化度對表面活性劑分子在油/水界面上的吸附性能影響較小，界面張力幾乎不變化[17]。因此，針對該類生物基兩性離子型表面活性劑所形成的油/水乳狀液體系，較適宜的礦化度為32 g/L。

固定生物基兩性離子型表面活性劑質量濃度為1.29 g/L、礦化度為32 g/L，NaOH質量分數對含有生物基兩性離子型表面活性劑溶液的油/水界面張力的影響見圖3(b)中曲線(1)。由圖3(b)中曲線(1)可知，隨著體系中NaOH質量分數增大，界面張力逐漸降低。當NaOH質量分數增加至1.0%時，界面張力為9.96×10-3mN/m，達到了超低界面張力(10-3mN/m)；繼續(xù)增大NaOH質量分數，界面張力繼續(xù)降低，在NaOH質量分數增大到2.5%時，界面張力降為9.60×10-4mN/m。這表明堿的加入有利于降低生物基兩性離子型表面活性劑體系的油/水界面張力，這是由于一方面堿與油相中的脂肪酸反應形成表面活性物質，與表面活性劑分子在油/水界面形成混合界面膜[18]；另一方面堿電離出的離子與生物基兩性離子型表面活性劑的親水基反應，降低表面活性劑分子間的斥力，使得表面活性劑分子在油/水界面排列更加緊密，從而降低界面張力。但當堿的用量過大時，會對地層造成傷害。因此，針對該類生物基兩性離子型表面活性劑所形成的油/水乳狀液體系，較適宜的NaOH質量分數為1.0%。

固定表面活性劑質量濃度為1.29 g/L、礦化度為32 g/L、NaOH質量分數為1.0%，HPAM質量濃度對含有生物基兩性離子型表面活性劑溶液的油/水界面張力的影響見圖3(b)中曲線(2)。由圖3(b)中曲線(2)可知：生物基兩性離子型表面活性劑溶液中引入HPAM，油/水界面張力變化較?。划擧PAM質量濃度小于500 mg/L時，隨著HPAM濃度的增大，界面張力變化較??；當HPAM質量濃度達到500 mg/L時，繼續(xù)增大HPAM濃度，界面張力增加幅度較大。這是由于加入體系中的HPAM，會與表面活性劑在油/水界面上發(fā)生競爭吸附，導致油/水界面上吸附的有效表面活性劑分子濃度降低，界面張力增大[19]。

圖3 乳化體系組成對生物基兩性離子型表面活性劑界面活性的影響Fig.3 Effect of the composition of emulsion system on interfacial activity of bio-based zwitterionic surfactant(a)—(1) Effect of CNBS mass concentration on interfacial tension (γ), condition: ρ(Salinity)=5 g/L；(a)—(2) Effect of salinity on interfacial tension (γ), condition: ρ(CNBS)=1.29 g/L;(b)—(1) Effect of NaOH mass fraction on interfacial tension (γ), conditions: ρ(Salinity)=32 g/L；ρ(CNBS)=1.29 g/L;(b)—(2) Effect of HPAM mass concentration on interfacial tension (γ), conditions: ρ(Salinity)=32 g/L；ρ(CNBS)=1.29 g/L； w(NaOH)=1.0%

2.2 體系總含水量對生物基兩性離子型表面活性劑乳化穩(wěn)定動力學的影響

根據2.1節(jié)中界面張力的研究結果，當生物基兩性離子型表面活性劑質量濃度為1.29 g/L、NaOH質量分數為2.5%、HPAM質量濃度為500 mg/L、礦化度為32 g/L時界面張力最低，以此條件配制溶液，以不同總含水量進行油/水混合，按照1.3節(jié)中方法研究體系總含水體積分數對生物基兩性離子型表面活性劑乳化穩(wěn)定動力學性能的影響，結果見圖4和表2所示。

φA—The water volume fraction in the oil phase圖4 體系總含水體積分數對模擬乳狀液乳化穩(wěn)定性的影響Fig.4 Effect of the bulk volume fraction of water contenton stability of the simulated emulsionConditions: ρ(Salinity)=32 g/L； ρ(CNBS)=1.29 g/L；w(NaOH)=2.5%； ρ(HPAM)=0.5 g/L

表2 體系總含水體積分數對模擬乳狀液乳化穩(wěn)定動力學參數的影響Table 2 Effect of the bulk volume fraction ofwater content on parameters regarding theemulsion-stabilizing kinetics of the simulated emulsion

由圖4可知，隨著靜沉時間延長，不同含水量的乳狀液體系中油相含水率φA均降低，乳狀液最終因不穩(wěn)定而破乳。這是由于乳狀液是熱力學不穩(wěn)定體系，隨著靜沉時間的延長，乳狀液中的微油滴(微水滴)發(fā)生運動碰撞，界面膜破裂，油滴(水滴)聚并，導致乳狀液穩(wěn)定性下降。總含水體積分數為60%的乳狀液油相含水率下降幅度較小，表明當體系總含水率為60%時，該生物基兩性離子型表面活性劑的油/水界面膜較穩(wěn)定。由表2可知,隨著總含水量增大，破乳速率常數hc先增大后減小，半衰期T1/2和90 min時油相含水率φA90先升高后降低。破乳速率常數hc越小、半衰期T1/2和90 min時油相含水率φA90越大，對應的乳狀液越穩(wěn)定。當乳化體系總含水率為60%時，破乳速率常數hc最小，半衰期T1/2和90 min時油相含水率φA90達到最高值，分別為256.72 min和47.23%。這進一步表明，當總含水率等于60%時，形成的乳化體系穩(wěn)定性較高，且乳狀液以W/O型為主。另外，隨著體系總含水率的增大，半衰期T1/2降低，體系穩(wěn)定性下降，但90 min時油相含水率φA90大于0，這表明該生物基兩性離子型表面活性劑由于其原料飽和腰果酚分子中長鏈親油烷基以及在反應過程中表面活性劑分子中引入的2個丁烷基團增加了該表面活性劑的親油性，體系中該表面活性劑的添加易于形成W/O型乳狀液。體系總含水率影響結果初步證實該生物基兩性離子型表面活性劑在總含水體積分數為60%時具有良好的乳化穩(wěn)定性。

2.3 生物基兩性離子型表面活性劑濃度對其乳化穩(wěn)定動力學的影響

固定總礦化度為32 g/L和體系總含水體積分數為60%時，分析該生物基兩性離子型表面活性劑質量濃度對模擬乳狀液穩(wěn)定動力學的影響，結果如圖5和表3所示。

φA—The water volume fraction in the oil phase圖5 表面活性劑CNBS濃度對模擬乳狀液乳化穩(wěn)定性的影響Fig.5 Effect of concentration of the surfactant CNBSon stability of the simulated emulsionConditions: ρ(Salinity)=32 g/L；φ(Total water)=60%

表3 表面活性劑濃度對模擬乳狀液乳化穩(wěn)定動力學參數的影響Table 3 Effect of surfactant concentration on parametersregarding the emulsion-stabilizing kinetics ofthe simulated emulsion

由圖5可知:隨著靜沉時間的延長，不同質量濃度生物基表面活性劑的模擬乳狀液油相含水率φA下降幅度開始較大，隨后降幅較??；在靜沉初期，隨生物基兩性離子型表面活性劑質量濃度增大，模擬乳狀液油相含水率φA增大；當靜沉時間為90 min，生物基兩性離子型表面活性劑質量濃度為1.5 g/L時，模擬乳狀液油相含水率φA最大。這是由于當表面活性劑質量濃度較低時，吸附在油/水界面上的表面活性劑分子較少，分子處于平躺或斜躺狀態(tài)，形成的界面膜穩(wěn)定性較差；當表面活性劑濃度大于臨界膠束濃度時，其在油/水界面達到飽和吸附，形成穩(wěn)定的界面膜。當表面活性劑濃度大于臨界濃度時，一方面，表面活性劑在水相易于形成膠束，促進水相乳化，不利于W/O型乳狀液的形成；另一方面，隨著表面活性劑濃度增大，表面活性劑在溶液中形成膠束，模擬乳狀液體系的黏度增大，一定程度上阻止表面活性劑分子進一步在油/水界面的定向排列[20]。由表3的動力學參數可知：當生物基兩性離子型表面活性劑質量濃度為1.5 g/L時，模擬乳狀液破乳速率常數最小，為0.0010 min-1；半衰期T1/2和90 min時油相含水率φA90達到最大值，分別為693.15 min和48.52%，油/水體系主要形成W/O型乳狀液。在該實驗測量濃度范圍內，φA90大于0，說明油/水體系主要形成W/O型乳狀液。當生物基兩性離子型表面活性劑質量濃度低于1.5 g/L時，隨著表面活性劑質量濃度的增大，模擬乳狀液破乳速率常數減小，半衰期T1/2和油相含水率φA90均升高；當表面活性劑質量濃度高于1.5 g/L，隨著表面活性劑質量濃度的增大，模擬乳狀液破乳速率常數增大，半衰期T1/2和油相含水率φA90均降低，這表明，表面活性劑濃度是影響乳狀液穩(wěn)定性的主要因素。

2.4 NaOH質量分數對生物基兩性離子型表面活性劑乳化穩(wěn)定動力學的影響

固定礦化度為32 g/L、總含水體積分數為60%、生物基兩性離子型表面活性劑質量濃度為1.5 g/L，考察NaOH質量分數對生物基兩性離子型表面活性劑形成的二元模擬乳狀液穩(wěn)定動力學的影響，結果如圖6和表4所示。

φA—The water volume fraction in the oil phase圖6 NaOH質量分數對模擬乳狀液乳化穩(wěn)定性的影響Fig.6 Effect of NaOH mass fraction on emulsifyingstability of the simulated emulsionConditions: ρ(Salinity)=32 g/L；φ(Total water)=60%； ρ(CNBS)=1.5 g/L

表4 NaOH質量分數對乳狀液乳化穩(wěn)定動力學參數的影響Table 4 Effect of NaOH mass fraction on parametersregarding the emulsion-stabilizing kinetics ofthe simulated emulsion

由圖6可知：隨著NaOH的加入，乳狀液體系的φA增加幅度較小；且隨著靜沉時間的延長，NaOH在二元模擬乳狀液體系的φA變化較小，實驗所測范圍內模擬乳狀液的φA范圍為53%～60%。這是由于當生物基兩性離子型表面活性劑質量濃度為1.5 g/L，形成的模擬乳狀液比較穩(wěn)定(見表3)。當體系中加入NaOH后，一方面NaOH通過電荷之間的相互作用，使溶液中的表面活性劑離子與溶液中無機離子達到平衡，同時影響生物基兩性離子型表面活性劑在水相分布能力和油/水界面的分布能力，使油/水界面張力降低；另一方面NaOH與原油中的酸性物質反應形成的表面活性物質也會吸附到油滴表面，與吸附在油/水界面上的生物基兩性離子型表面活性劑發(fā)生協(xié)同吸附[21]，增加界面膜上吸附的分子密度，提高了界面膜的穩(wěn)定性，多方面的協(xié)同作用增加了乳狀液的穩(wěn)定性。但是由于表面活性劑分子在油/水界面已達到飽和吸附，形成的表面活性物質在油/水界面有效吸附量較小，且受原油中酸性物質含量的影響，導致NaOH質量分數對生物基兩性離子型表面活性劑形成的二元模擬乳化體系的乳化穩(wěn)定動力學性能影響較小。從表4也可以看出，隨著NaOH質量分數增大，φA90和破乳速率常數hc變化均較小，且當NaOH質量分數在0～1.49%范圍內，φA90大于0，體系仍主要是W/O型乳狀液。然而，半衰期T1/2變化很大，當NaOH質量分數為0.40%時，半衰期T1/2為866.43 min；當NaOH質量分數為1.00%時，半衰期T1/2達到最大，為3465.74 min。這表明，在二元體系中，NaOH與生物基兩性離子型表面活性劑協(xié)同形成較穩(wěn)定的模擬乳狀液。但是NaOH質量分數過高會傷害地層，因此，針對該研究體系，NaOH適宜的質量分數為1.00%。

2.5 部分水解聚丙烯酰胺濃度對生物基兩性離子型表面活性劑乳化穩(wěn)定動力學的影響

固定礦化度為32 g/L、總含水體積分數為60%、表面活性劑質量濃度為1.5 g/L、堿質量分數為1.00%，考察HPAM質量濃度對模擬三元乳狀液穩(wěn)定動力學的影響，結果見圖7和表5所示。

由圖7可以看出，隨著靜沉時間的延長，φA逐漸下降，其中HPAM質量濃度為2000 mg/L的體系降低速率最快。這表明，隨著HPAM的加入，生物基兩性離子型表面活性劑形成的乳狀液乳化穩(wěn)定性下降。從表5也可以看出，隨著HPAM質量濃度增大，破乳速率常數hc增大，半衰期T1/2減小，這是由于HPAM除影響乳化體系的油、水分離過程外，還會影響乳狀液的形成。這是因為，在乳狀液制備過程中，隨著HPAM質量濃度增大，水相溶液黏度逐漸增大，會消耗一定量的機械能，導致促進乳狀液生成的能量減少。另外，隨著HPAM質量濃度增大，表面活性劑分子向油/水界面擴散的速率降低，溶液界面張力增大(圖3)；同時模擬乳狀液體系黏度的增大，導致供乳狀液生成的能量和乳化液滴運動速率和碰撞機率降低[22]。當HPAM質量濃度為900 mg/L時，φA90最大，這是由于當HPAM質量濃度小于900 mg/L時，提高HPAM質量濃度可以促進乳狀液生成和穩(wěn)定；當HPAM質量濃度大于900 mg/L時，提高HPAM質量濃度會減弱乳狀液的生成。因此，采用強堿三元復合驅油體系提高原油采收率時，結合驅油體系的黏度和界面張力，針對該生物基兩性離子型表面活性劑形成的強堿三元復合驅油體系，選擇適宜的HPAM質量濃度為900 mg/L。

φA—The water volume fraction in the oil phase圖7 聚合物質量濃度對模擬乳狀液乳化穩(wěn)定性的影響Fig.7 Effect of polymer mass concentration on stability ofthe simulated emulsionConditions: ρ(Salinity)=32 g/L；φ(Total water)=60%；ρ(CNBS)=1.5 g/L； w(NaOH)=1.00%

表5 HPAM質量濃度對三元模擬乳狀液乳化穩(wěn)定參數的影響Table 5 Effect of HPAM mass concentration on stabilityparameters regarding the ternary simulated emulsion

2.6 生物基兩性離子型表面活性劑的界面微觀特性

生物基兩性離子型表面活性劑在不同體系中的界面活性以及乳化穩(wěn)定動力學研究結果表明，該活性劑在油/水界面的分布是影響界面張力和乳液穩(wěn)定性的關鍵因素。為了進一步研究其在油/水界面的微觀特性，對油/水、油/水/CNBS、油/水/CNBS/NaCl和油/水/CNBS/NaCl/NaOH 4種體系形成的油/水界面的微觀特性進行分子模擬及模擬結果分析。研究油/水/表面活性劑體系的界面時，按照“90%～90%”原則進行界面厚度計算，計算油相和水相的密度分布見圖8(a)，其中油相密度的90%與水相密度90%之間的距離為界面厚度[23]?；诖朔椒ńy(tǒng)計所模擬的體系油相和水相密度沿Z軸分布情況并計算界面厚度，見圖8。根據式(4)計算油/水界面能，結果見表6。

由圖8和表6可知:純油/水體系的油/水界面膜厚度為0.3795 nm；而加入生物基兩性離子型表面活性劑后，油/水/CNBS體系的油/水界面厚度增加到0.5948 nm，且界面能由-2241.14 kJ/mol降至-2379.91 kJ/mol。表明，隨著生物基兩性離子型表面活性劑的加入，形成的模擬乳狀液的穩(wěn)定性增強。向油/水/CNBS體系中加入NaCl，油/水/CNBS/NaCl體系的油/水界面厚度由0.5948 nm降至0.5739 nm，油/水/CNBS/NaCl體系的界面能由-2379.91 kJ/mol降至-2464.24 kJ/mol(見表6)。向油/水/CNBS/NaCl體系中加入NaOH，油/水/CNBS/NaCl/NaOH體系的油/水界面厚度和界面能進一步下降，分別降至0.5422 nm和-2517.86 kJ/mol。這進一步表明，隨著NaCl和NaOH的加入，油/水/CNBS/NaCl體系和油/水/CNBS/NaCl/NaOH體系越穩(wěn)定，形成的乳狀液穩(wěn)定性增強。這是因為，在水相中加入NaCl或NaOH后，離子濃度發(fā)生了變化，導致界面層被壓縮，油/水界面膜厚度變小，但壓縮后界面層的界面能更低，界面的穩(wěn)定性得到增強。

圖8 模擬體系中油分子和水分子沿z軸密度分布曲線Fig.8 z-Axis density distribution of water and oil molecules in the simulation system(a) Oil/Water system； (b) Oil/Water/CNBS system； (c) Oil/Water/CNBS/NaCl system； (d) Oil/Water/CNBS/NaCl/NaOH system

表6 模擬體系界面能計算結果Table 6 Calculation results of the interfacial energy in the simulation system

3 結 論

(1)通過測定界面張力和分子動力學模擬技術研究了生物基兩性離子型表面活性劑在油/水界面的性質和分布，界面張力測定結果表明，在礦化度為32 g/L、NaOH質量分數為2.5%的條件下，該生物基兩性離子型表面活性劑的質量濃度為1.29 g/L的溶液的油/水界面張力最低可達9.60×10-4mN/m。分子模擬結果表明，在油/水體系中加入該生物基兩性離子型表面活性劑可將油/水界面膜厚度由0.3795 nm提高至0.5948 nm。

(2)NaCl和NaOH使油/水界面膜厚度和界面能均降低。

(3)乳化穩(wěn)定動力學研究結果表明，該生物基兩性離子型表面活性劑形成的模擬乳狀液以W/O型為主，表面活性劑質量濃度是影響乳狀液穩(wěn)定動力學的關鍵因素，當表面活性劑濃度接近其臨界膠束濃度時，形成的模擬乳狀液穩(wěn)定動力學參數半衰期和油相含水率均最大，破裂速率常數最小，此時形成的乳狀液最穩(wěn)定。