Physical characteristics of cores used in this study.
The permeability of tight gas reservoirs is usually lower than 1 md. When the external fluids from drilling and completion processes invade such reservoirs, formation damage occurs and causes serious damage to oil and gas production. Fluorocarbon surfactants are most often recommended for removing such damage because they have extremely low surface tension, which means that they can change the reservoir wettability from water wet to gas or oil wet. However, they are not normally applied in the field because they are not cost-effective. Besides, some environmental concerns also restrict their application. In this work, we studied the effects of an oligomeric organosilicon surfactant (OSSF) on wettability modification, surface tension reduction, invasion of different fluids, and fluid flow back. It was found that the amount of spontaneous imbibition and remaining water could be reduced by the surfactant as a result of surface tension reduction and wettability alteration. Compared to the distilled water, the concentration of 0.20 wt% OSSF could decrease water saturation of cores by about 4%. At a flow-back pressure of 0.06 and 0.03 MPa after 20 PV displacement, permeability recovery could increase from 8 to 7–93% and 86%, respectively. We also found that the mechanism of OSSF includes the physical obstruction effect, surface tension reduction of external fluids, and wettability alteration of the reservoir generated. Meanwhile, quantum chemical calculations indicated that adsorbent layer of polydimethylsiloxane could decrease the affinity and adhesion of CH4 and H2O on the pore surface.
- low-permeability reservoir
- silicone hydrophobic polymer
- spontaneous imbibition
- antiwater blocking
- water saturation
- gas permeability improvement
Water blocking can damage the low permeability of a normally tight gas reservoir due to the increase of water saturation and the reduction of gas phase permeability in the process of drilling, completion, and stimulation. Zhong et al.  showed that the permeability damage rate in a gas reservoir could reach 70–90%, and the gas well production could decrease by more than 70% when water-blocking damage occurred. Therefore, it is important to develop an antiwater blocking agent with high efficiency and low cost for the development of a superior low-permeability gas reservoir.
At present, most researchers consider water blocking as caused by capillary thermodynamics and dynamics . The commonly used antiwater blocking agents include lower alcohol content agents (especially methanol), alcohol ethers, silyl ethers, hydrocarbon surfactants , and fluorocarbon surfactants, and “… the alcohols react with formation water in the reservoir to form low boiling point azeotrope, which is helpful for gasification and flow back of the injected water” . Nasr et al. indicated that the glycol ether and polyethylene glycol monobutyl ether could remove the water blocking damage and improve the permeability of oil fields in Arab countries . Zhang et al. evaluated the effects of methanol, ethanol, and ethylene glycol on alleviating the water blocking damage of the low-permeability sandstone gas reservoir. The results showed that methanol had the most favorable performance, followed by ethanol and ethylene glycol, respectively . “Bai et al. compared the anti-water blocking effects of methanol and petroleum sulfonate on a low permeability gas reservoir and found methanol to have better performance…” Surfactants can reduce the surface tension and change the reservoir wettability from water wetting to gas or oil wetting . Bang et al. reduced the water blocking damage and increased the fluid flow in the fracture reservoirs with an alcohol solution containing a fluorocarbon surfactant. They also reduced the water blocking damage of the condensate in gas reservoirs by applying the fluorocarbon surfactant. The results showed that the hydrated silanol groups could adsorb water from pore surfaces via covalent bonding with the alkenyl groups of fluorocarbon surfactants. This made fluorine-containing alkyl directionally arranged, changing the reservoir wettability from water wetting to oil wetting and finally increasing the gas phase permeability [8, 9, 10]. Li et al. developed an antiwater blocking agent containing perfluoroalkyl side chains by using a stepwise emulsion polymerization method, which showed low surface tension and interfacial tension. The surfactant can adsorb on the surface of the pore via chemical adsorption and modify the surface to preferential gas wetting . Li et al. used two nonionic fluorocarbon surfactant methanol solutions for water phase displacement experiments and found that it could improve the permeability . Liu et al. synthesized a cationic fluoride Gemini surfactant and recorded ultralow surface tension; moreover, its solution could remarkably reduce the damage of water blocking on low-permeability formations . Liu et al. found that a wettability variation from water wetting to gas wetting was achieved by adding 0.1 wt% fluoride. The core flow tests indicated that both the flow back rate and gas relative permeability were significantly improved . Biosurfactants such as
Due to the high cost and environment concerns of fluorocarbon surfactants, an oligomeric organosilicon surfactant (OSSF) contained functional groups such as the silicon hydroxyl group. The silicon oxygen chain and silicon methyl were prepared by polycondensation reaction . Gas preferential wettability could be achieved by oriented adsorption of silicon methyl via hydrogen bonding or chemical condensation. It could also decrease the dynamic surface tension of fluids owing to faster interfacial adsorption of small molecules during diffusion adsorption. Besides, it possesses super-strong surface activity, and a very small amount can decrease the surface tension to less than 25 mN/m, while the cost is less than 1/10 of fluorocarbon surfactants.
In this study, the influence of OSSF on the aqueous phase trapping damage was evaluated through the measurement of spontaneous imbibition distilled water saturation, water blocking damage rate, permeability recovery, and remaining water saturation. The mechanisms were explained using aggregates blocking, directional adsorption, wettability alteration, and quantum chemistry calculations.
2. Experimental approach
In this study, OSSF and ABSN were used to alleviate the aqueous phase trapping damage to low-permeability sand formation. Both surfactants are water soluble and transparent liquids and were made in our laboratory. The relative molecular mass of an OSSF is 2000–3000. It includes a hydrophobic group of permethylated siloxane, hydrophilic groups of sulfonic acid groups, and hydroxyls. The polydimethylsilioxane chain enables the surfactant to have water-repelling characteristics and low surface tension. These properties enable its adsorption on the rock via multiple points of attachment and give the surfactant a long-term effectiveness. ABSN is quaternary ammonium salt cationic surfactant, and the hydrophobic tail is dodecyl benzene. Figure 1(a) and (b) shows the chemical structure of dodecyl benzene sulfonate triethanolamine (ABSN) and oligomeric organosilicon surfactant (OSSF), respectively. The artificial and reservoir cores used in the same experiment have the same composition and similar porosity and permeability.
The cores used in this study include artificial cores and reservoir cores. They have same compositions and similar porosity and permeability. The artificial cores (obtained from Haian Petroleum Scientific Research Instrument Co., Ltd. in Nantong, China) are mainly composed of quartz sand, which does not contain any water sensitive substance. Reservoir cores were drilled at the Yingxi exploratory area in the Qinghai oilfields. The physical characteristics of these cores are listed in Table 1.
|Type||L/(mm)||D/(mm)||m/(g)||Ф/(%)||K/(10−3 μm2)||Standard error|
Triply distilled water was used in all experiments. About 4% of sodium bentonite mud was added with 0.2% Na2CO3 and then stirred for 2 days. It was used as the basic mud. Sodium bentonite was obtained from Weifang Hua Bentonite Group Co., Ltd. Na2CO3 is a commercially pure reagent.
2.2 Experimental method
2.2.1 Determination of contact angle, wettability, and surface energy
Cores were immersed in NaOH solution with a pH of 9 at 150°C for 16 h. Each core was then taken out and dried at 150°C for 4 h. After that, it was cooled down to room temperature. The contact angles among the cores, distilled water, and ethylene glycol were measured with a JC200D3 contact angle analyzer. When the contact angle of water is less than 75°, it is termed “water wetting.” When the contact angle of water is greater than 110°, it is considered “oil wetting” . When the contact angle is between 75° and 110°, it is considered “intermediate wetting.” This is a condition of gas preferential wetting. The surface energy of cores was calculated using the Owens-Wendt formula . Aging in this paper means rolling under the high temperature of 150°C for 16 h.
2.2.2 Determination of surface tension
The surface tension measurements of NaOH solution with a pH of 9 and basic mud filtrates before and after aging (rolling in the temperature of 150°C oven for 16 h) were performed with an interfacial tensiometer sigma 701 (KSV, Finland) using the dynamic Wilhelmy plate method.
2.2.3 Spontaneous imbibition experiments
Figure 2 shows the experiment’s core inside the core holder, and a certain confinement pressure was added. The distilled water was added into the metering tube followed by the OSSF solution. The initial contact volume of the liquid with cores was determined based on the different flow rates. The increased water saturation, defined as spontaneous imbibition volume divided by the core’s pore volume, was recorded during the experiment . Core permeability change rates at different times were defined as the water blocking damage rate. This rate can be calculated :
2.2.4 Gas driven flow-back experiment
The cores were gas driven and saturated with flow-back spontaneous imbibition. They were used to weigh the core mass at different flow-back stages and to calculate the residual water saturation. The gas permeability recovery was calculated as follows [21, 22]:
3. Results and discussion
3.1 Wettability properties
This section presents an evaluation of the wettability alteration on both artificial cores and reservoir cores using different surfactants in alkaline solutions. Both the distilled water and ethylene glycol experiments were carried out under high temperatures. As shown in Table 2 (Figure 3), water and ethylene glycol spreading on the artificial cores were treated by 0.2 and 0.4 wt% ABSN; thus, the contact angles were too low to measure. The contact angles of 0.2 wt% OSSF treated cores were 110.12° and 27.54°, respectively. Considering the microscopic anisotropy of the cores and the measurement error, we believe that the real contact angles should not vary dramatically from one another. The surface energy of the cores decreased to 20–25 mJ/m2, which indicated that OSSF could greatly change the wettability from hydrophilic to highly hydrophobic.
|Concentration/wt%||Contact angle/°||Surface energy/(mJ/m2)||Standard error|
|Distilled water||Ethylene glycol||Distilled water||Ethylene glycol||Surface energy|
|0.00||Droplet infiltration||Droplet infiltration||—||—||—|
|0.20% ABSN||Droplet infiltration||Droplet infiltration||—||—||—|
|0.40% ABSN||Droplet infiltration||Droplet infiltration||—||—||—|
As shown in Table 3 (Figure 4), the 0.2 and 0.4 wt% ABSN-treated cores displayed lower water contact angle, and the ethylene glycol droplet completely spread on the surface. However, the contact angle of water on the OSS-modified cores still reached 110°, and the surface energy decrease was between 26 and 28 mJ/m2. These results indicated that water did not spread on the surface of the pore due to the formation of low energy adsorption film treated by OSSF. Since ABSN does not affect wettability alteration, its surface energy is not evaluated here.
|Concentration/wt%||Contact angle/ (°)||Surface energy/(mJ/m2)||Standard error|
|Distilled water||Ethylene glycol||Distilled water||Ethylene glycol||Surface energy|
|0.20% ABSN||8.75||Droplet spreading||—||—||—|
|0.40% ABSN||6.42||Droplet spreading||—||—||—|
3.2 Interfacial properties
To evaluate the surface activities of OSSF and ABSN in a high-temperature basic environment, the surface tension of solutions and fluid filtrates was measured after aging at 150°C for 16 h. In Figures 5 and 6, both 0.2–0.4 wt% ABSN and OSSF could decrease the surface tension under high temperature. Meanwhile, OSSF was better at reducing surface tension, meaning that it could better reduce filtrate adsorption and aqueous phase trapping damage in the rock.
3.3 Evaluation of antiwater blocking properties of OSSF
3.3.1 Spontaneous imbibition property of cores
As shown in Figure 7, water saturation of cores gradually increased with time. Obvious spontaneous imbibition and diffusion stages could also be seen. The increase in water saturation was dramatic at the spontaneous imbibition stage, but water saturation was compromised at the diffusion stage. Effective pores can imbibe the liquid by capillary force, whereas liquid accessing to unconnected pores only relies on diffusion from connected pores, which requires more time. Water saturation decreased from about 71–68% with a 0.20 wt% OSSF in the spontaneous imbibition stage. The spontaneous imbibition time then changed from about 33 to 28 min.
Figure 8 shows the cores treated by 0.2% OSSF solution. Water saturation decreased to 56% and the saturation time decreased to 68 min in the spontaneous imbibition stage. These results indicate that the spontaneous imbibition rate and the amount of external liquid in a reservoir can easily be reduced.
3.3.2 Water blocking damage rate of cores
As shown in Figure 9, the water blocking damage rate of distilled water and OSSF solutions increased sharply at the beginning. However, after 20 min, it slowly started to increase. Compared with distilled water, OSSF solution could achieve an additional damage rate reduction of 15%. The results suggest that the water blocking damage rate of the reservoir increased remarkably when the external liquids invaded into reservoir at the beginning. After reaching a certain water saturation, the water blocking damage did not increase significantly.
3.3.3 Permeability recovery of cores
As shown in Figure 10, the permeability recovery of the cores gradually increases with flow-back PV. Meanwhile, the recovery rate gradually decreases. When gas pressure was 0.06 MPa and the flow-back PV is 20, the cores’ permeability recovery changes from 7 to 77%. It is evident that the higher the displacing pressure, the higher the permeability recovery. A permeability recovery of 0.20% OSSF solution increased from about 8 to 93%; then, with the distilled water, the permeability recovery was 77%. Both increases were at 0.06 MPa and 20 PV. At 0.03 MPa and 20 PV, permeability recovery of 0.20% OSSF solution increased from about 7 to 86%, with the distilled water 67% by contrast. These results illustrate that the increase of flow-back volume and gas pressure can improve the permeability recovery and reduce the water blocking damage of cores. In addition, OSSF is beneficial to the permeability recovery of the reservoirs damaged by external fluids.
3.3.4 Retained water saturation of cores
Accumulation of external fluids in gas reservoir is the principle factor responsible for water blocking damage [23, 24]. As shown in Figure 11, the remaining water saturation decreased with PV, but the decreasing rate reduced gradually. Retained water saturation of 0.20% OSSF solution decreased from about 83 to 30%; moreover, with the distilled water, the water saturation decreased from 90 to 34% by contrast. The results reveal that increasing flow-back time can reduce the water saturation and the water blocking damage of the reservoir. Besides, OSSF is also beneficial to the flow back of invasive external fluid in a reservoir.
3.4 Mechanism analysis
3.4.1 Physical blocking
Figure 12 shows that aggregates can be formed when the concentration exceeds its critical micelle concentration (CMC). Adjusting the OSSF concentrations to match reservoir needs with varied pore and throat sizes is necessary when the aggregates invade reservoir with external fluids. This causes friction to form on the surface in which pores or throats and fluid come in contact with each other.
When the aggregates move into the pores or throats of reservoirs, as shown in Figure 13(b), the interfacial tensions on the two contact surfaces are even. When the aggregates flow from the pores to the reservoir throats (Figure 13(c)), friction can decrease the flow rate. When aggregates move from the throats to the pores (Figure 13(d)), the interfacial tension between aggregates and external fluid becomes smaller. In the process of production, the fluid begins to flow back, the aggregates move back from the throats to the pores, and then, the differences in the interfacial tensions between two contact surfaces could become the driving force for flowing, which could accelerate these fluids to flow back.
3.4.2 Wettability alteration
When the alkaline solution invades into the reservoir, the structural groups of the silicon-oxygen bond (Si-O-Si) on the surface of the silicate minerals can produce silanol groups (Si-OH) by hydrolysis. The density can reach up to 6–7 nm per square. These silanol groups can adsorb the water molecules owing to the hydrogen bonding and the Van Der Waals force upon the water-pore interface. Within both forces, hydrogen bonding plays an important role. It is needed to overcome the adhesion force between rock and water, thereby decreasing the remaining water in the formation . The basic structural unit of the OSSF is polydimethylsiloxane. The silicon-oxygen chain is the polar part of the OSSF composition, and the methyl group represents the nonpolar part.
Under high temperature and catalyst conditions, polarization occurs in the silicon oxygen main chain. This leads to multipoint adsorption of silanol groups (Si-OH) and silicon-oxygen bond (Si-O-Si) through chemical bonding and hydrogen bonding. Meanwhile, the methyl groups directionally rotate and sequentially arrange on the rock surface. This process is shown in Figure 14 . These directionally adsorbed methyl groups reduce the surface energy of the rock and lead to the hydrophobization on the rock surface, changing the composition of the rock surface. Meanwhile, the capillary force direction of fluids and the physical properties of the gas reservoir change .
As shown in Figure 15, the cores treated by OSSF have a higher weight and atom Si number as well as a lower weight and atom O number. This is because the cores are composed of silicon and oxygen, with a chemical structural much like that of OSSF. Figure 16 shows how the OSSF adsorption on the reservoir core surface can significantly increase the weight and number of atom Si and Ca but slightly decrease that of atom O. Since we use carbonate cores, their surface components change greatly as a result of OSSF adsorption. Both artificial and reservoir core surfaces show more Si atoms. The outermost areas of the core are composed of Si-CH3, and the cores have become hydrophobic.
3.4.3 Adsorption energy
The weak interaction among the molecules can be calculated accurately by the computational technique based on quantum chemistry. In this paper, (a) silicon-oxygen tetrahedron and (b) dimethyl siloxane were applied to simulate the surface of the untreated cores and OSSF treated cores, respectively. The oxygen atoms at the edge were saturated with hydrogen atoms.
Gaussian 09 W software (Gaussian Inc., Wallingford, CT, USA) is used as the simulation tool, and 6-31G is selected as the basic set. Correction is performed using a basis set superposition error (BSSE). In the adsorption system, the distance between the C atoms of methane and the O atoms of water is the equilibrium distance
As shown in Figure 17 and Table 4, the distances of the dimethyl siloxane models between CH4 and H2O are larger or similar than that of silicon-oxygen tetrahedron. The adsorption energy of silicon-oxygen tetrahedron and the dimethyl siloxane models to CH4 is 11.21 and 0.46 kJ/mol, respectively, which indicates that the adsorption of methane on the core surfaces is nonhydrogen-bonding physical adsorption. The adsorption energy of silicon-oxygen tetrahedron and the dimethyl siloxane models to H2O is 81.78 and 52.18 kJ/mol, respectively, suggesting the adsorption of water is hydrogen-bonding adsorption. Since the adsorption binding energy of CH4 and H2O on the silicon-oxygen tetrahedron model is higher than that on the dimethyl siloxane model, the stability of CH4 and H2O adsorbed on the pore surface of cores treated by OSSF is lower than that on the untreated cores. At the same flow-back pressure, CH4 and H2O adsorbed on the surface of treated cores are often more easily desorbed. Under the same temperature and pressure conditions, after the adsorption of OSSF, the affinity and adhesion of the rock surface to CH4 and H2O are reduced, leading to the reduction of the shear stress of CH4 and H2O on the pore surfaces of the reservoir, when gas or external fluids flow in the reservoir . This can increase the fluidity of CH4 and H2O and finally improve the permeability of CH4 and H2O, thereby decreasing water saturation of the reservoir.
|Adsorption model||Adsorbate||re/ (nm)||Ee/ (kJ·mol−1)|
|Silicon-oxygen tetrahedron (a)||CH4||4.33||−11.21|
1. An ultralow surface tension could be obtained using OSSF. A 0.20 wt% OSSF solution can alter cores from water wetting to preferential gas wetting.
2. OSSF can remarkably decrease the water saturation and spontaneous imbibition rate of cores, resulting in a reduction in the cores’ water blocking damage rate. In this case, the effect of cores’ wettability alteration to water saturation and the spontaneous imbibition rate is higher.
3. In the gas driven flow-back experiment, OSSF improved core permeability recovery. The remaining water saturation was also found to be low.
4. The main mechanism for improving core permeability includes the physical obstruction effect, surface tension reduction, and adsorption energy reduction, as well as the capillary force, adhesion, and shear stress of methane and liquid on rock surfaces.
This work was supported by the National Major Science and Technology Project of China (Nos. 2016ZX05051, 2016ZX05052, and 2016ZX05030) and Research on Qaidam’s Excellent Drilling and Completion Technology (2016E-0109).
Conflicts of interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
The supervisor, Jie Zhang, conducted the experimental work and analyzed the results of the research. All authors contributed to the analysis of results. Jie Zhang wrote the manuscript. Xu-Yang Yao contributed to experiment’s design and analysis of the results. Bao-Jun Bai contributed to analysis and discussion of the results, reviewing, and editing of the manuscript.
Zhong XY, Huang L, Wang LH. Research progress on water blockage effect in low permeability gas reservoirs. Special Oil & Gas Reserves. 2008; 15:12-23
He CZ, Hua MQ. Quantitative study of water blocking mechanism. Drilling Fluid & Completion Fluid. 2000; 3:4-7
Pham VP, Jo YW, Oh JS, Kim SM, Park JW, Kim SH, et al. Effect of plasma-nitric acid treatment on the electrical conductivity of flexible transparent conductive films. Japanese Journal of Applied Physics. 2013; 52:075102
Liu M, Wu Q. Research progress and prospect of water blocking agent. Guangzhou. Chemistry & Industry. 2013; 41:32-42
Nasr-EI-Din HA, Lynn JD, AI-Dossary KA. Formation damage caused by a water blockage chemical: Prevention through operator supported test programs. SPE. 2002; 33:297-325
Zhang WB, Li CF, Shi XW. Application of alcoholic acid acidification in Qinghai oil field. Oilfield Chem. 2009; 26:24-28
Bai FL. Water phase trapping damage and the study of relieving measures in gas reservoir. Petrochemical Industry Application. 2010; 29:14-17
Bang V, Pope G, Sharma M, Baran JR Jr. Development of a successful chemical treatment for gas wells with water and condensate blocking damage. In: SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 4–7 October 2009. Paper SPE 124977. 2009
Bang VSS, Yuan CW, Pope GA, Baran JR Jr, Skildum J. Improving productivity of hydraulically fractured gas condensate wells by chemical treatment. In: OSSFhore Technology Conference, Houston, Texas, 5–8 May 2008. Paper OTC 19599. 2008
Bang V, Pope G, Mukul M, Baran JR Jr, Ahmadi M. A new solution to restore productivity of gas Wells with condensate and water blocks. In: SPE Annual Technical Conference and Exhibition, Denver, Colorado, 21–24 September 2008. SPE 116711. 2008
Li YY, Jiang GC, Xuan Y. Synthesis and performance evaluation of drilling fluid anti-blockage agent of low porosity and low permeability reservoir. Drilling Fluid and Completion Fluid. 2014; 31:9-12
Li XQ, Wang YC, Wang YQ. Experimental study on the effect of mixed surfactants on water blockage. Chemistry & Bioengineering. 2013:85-89
Liu Y, Guo L, Bi K, Ren W, Liu L, Li YG, et al. Synthesis of cationic gemini surfactant and its application in water blocking. Fine and Specialty Chemicals. 2011; 19:8-10
Liu XF, Kang YL, Luo PY, You LJ, Tang Y, Kong L. Wettability modification by fluoride and its application in aqueous phase trapping damage removal in tight sands to ne reservoirs. Journal of Petroleum Science and Engineering. 2015; 133:201-207
Zhang XQ. Researching on the use of new type surfactant to weaken water blocking and jamin effect. Master thesis. Qingdao, China: Ocean Univ. China; 2013
Yin D, Luo P, Zhang J. Synthesis of Oligomeric silicone surfactant and its interfacial properties. Applied Sciences. 2019; 9(3):497
Jiang GC, Wei YJ, Zhang M, Luo SJ, Yang Z, Pang JT. Review of evaluation methods to gas wettability. Science Technology & Engineering. 2012; 13:5562-5565
Owens DK, Wendt RC. Estimation of the surface free energy of polymers. Journal of Applied Polymer Science. 1969; 13:1741-1747
Dong B, Lan L, Cheng ZH, Xia HY, Niu J. Study on capillary self-absorption characteristics of dense gas reservoir rocks. Drilling and Production Technology. 2012; 35:34-37
Lai NJ, Ye ZB, Liu XJ, Yang JJ, Zhang JF. Study on water blocking damage in low permeability tight sandstone gas reservoir. Natural Gas Industry. 2005; 25:125-127
Stand. Adm. China. Determination Method of Relative Permeability of Two-Phase Fluids in Rocks (GB/T 28912–2012). Beijing, China: China Stand. Press; 2012
Guo TS, Li TW, Li HB. Laboratory study of water locking damage under imbibition in volcanic gas reservoirs. Science Technology & Review. 2011; 29:62-66
Bennion DB, Thomas FB, Ma T. Formation damage processes reducing productivity of low permeability gas reservoirs. In: SPE Rocky Mountain Regional/Low Permeability Reservoirs Symposium and Exhibition, Denver, Colorado, 12–15 March 2000. Paper SPE 60325. 2000
Bennion DB, Bietz RF, Thomas FB. Reductions in the productivity of oil and gas reservoirs due to aqueous phase trapping. Petroleum Society of Canada. 1994; 33:45-54
Wu JF, Wu DQ. Surface ionization and complexation of mineral-water interface. Advances in Earth Science. 2000; 1:90-96
Cao Y, Li HL. Adsorption behavior of polymeric surfactants at solid-liquid interface. Acta Physico-Chimica Sinica. 1999; 10:895-899
Koretsky CM, Sverjensky DA, Sahai N. A model of surface site types on oxide and silicate minerals based on crystal chemistry. Implications for site types and densities, multi-site adsorption, surface infrared spectroscopy, and dissolution kinetics. American Journal of Science. 1998; 298:349-438
Castro EAS, Gargano R, Martins JBL. Theoretical study of benzene interaction on kaolinite. Computer-Aided Materials Design. 2012; 112:2828-2831
Taghvaei E, Moosavi A, Nouri-Borujerdi A. Superhydrophobic surfaces with a dual-layer micro-and nanoparticle coating for drag reduction. Energy. 2017; 12:1-10