[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Recent Trends in the Application of Photocatalytic Membranes in Removal of Emerging Organic Contaminants in Wastewater
Previous Article in Journal
A Two-Level Facility Layout Design Method with the Consideration of High-Risk Facilities in Chemical Industries
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microbially Induced Calcite Precipitation (MICP) Improved Drilling Fluid Optimization for Gravel Stratum

1
College of Civil Engineering and Transportation, Hohai University, Nanjing 210098, China
2
China Railway Major Bridge Reconnaissance & Design Institute Co., Ltd., Wuhan 430101, China
3
Department of Civil Engineering, Suzhou University of Science and Technology, Suzhou 215011, China
*
Author to whom correspondence should be addressed.
Processes 2025, 13(1), 162; https://doi.org/10.3390/pr13010162
Submission received: 18 December 2024 / Revised: 5 January 2025 / Accepted: 7 January 2025 / Published: 9 January 2025
(This article belongs to the Section Chemical Processes and Systems)
Figure 1
<p>Sample preparation.</p> ">
Figure 2
<p>Correlation analysis between various factors and responses. scatter plots: (<b>a</b>) the interactions between B and PV; (<b>b</b>) the interactions between B and D; (<b>c</b>) the interactions between B and YP. (A: Bentonite; B: biological solution; C: barite; PV: plastic viscosity; D: density; YP: yield point. The symbol ‘-’ indicates a negative correlation between variables; “Run” refers to the number of experiments conducted).</p> ">
Figure 3
<p>Three-dimensional response surface plots for various factors: (<b>a</b>) A and B, (<b>b</b>) A and C, (<b>c</b>) B and C; (<b>d</b>) diagnostic chart of predicted and actual values of plastic viscosity; (<b>e</b>) the normal graph of residuals; (<b>f</b>) residual and running chart.</p> ">
Figure 4
<p>Three-dimensional response surface plots of various factors: (<b>a</b>) A and B, (<b>b</b>) A and C, (<b>c</b>) B and C; (<b>d</b>) diagnostic chart of predicted and actual densities; (<b>e</b>) the normal graph of residuals; (<b>f</b>) residual and running chart.</p> ">
Figure 5
<p>Three-dimensional response surface plots of various factors: (<b>a</b>) A and B, (<b>b</b>) A and C, (<b>c</b>) B and C; (<b>d</b>) diagnosis chart of predicted and actual values of yield point; (<b>e</b>) the normal graph of residuals; (<b>f</b>) residual and running chart.</p> ">
Figure 6
<p>EDS spectrum of biological mud (<b>a</b>–<b>d</b>).</p> ">
Figure 7
<p>Microscopic morphology diagram of biological mud: (<b>a</b>) Calcium carbonate distribution; (<b>b</b>) ribbon-like connection; (<b>c</b>) calcium carbonate distribution on the fibrous connector.</p> ">
Versions Notes

Abstract

:
During the exploration of the gravel stratum, incidents such as wellbore leakage, stuck drilling, and unstable wellbore walls frequently occur. These issues lead to diminished drilling efficiency and prolonged construction timelines, ultimately adversely affecting the core recovery rate, resulting in a significant waste of manpower and material resources. To address the issue of hole collapse during drilling, the microbially induced calcite carbonate precipitation (MICP) technique was employed to enhance the properties of bentonite mud drilling fluids. This study analyzed the effects of three factors, i.e., bentonite, biological solution, and barite powder, on the bentonite mud bio-cementation effectiveness through an orthogonal experiment and response surface methodology (RSM). The biological mechanism was examined using scanning electron microscopy (SEM). The experimental results indicated that optimal formulation was achieved when the mass fraction of bentonite was 13.96%, the biological solution comprised 0.6% xanthan gum and 0.4% carboxymethyl cellulose, and the mass fraction of barite was 25%. This research explores the application potential of MICP in enhancing the rheological properties of bentonite mud drilling fluids, which provides new insights and technical references for optimizing their performance.

1. Introduction

The gravel stratum is composed of gravel particles of different shapes and sizes [1]. Its structure is characterized by large pores and minimal cohesion. During drilling exploration in gravel strata, several challenges frequently arise, including falling blocks, collapsing holes, and the formation of sand bridges, which may impede progress and necessitate additional precautions [2]. In addition, vibration during the drilling process may cause wellbore wall collapse, thereby leading to buried drilling problems and incurring economic losses. To mitigate borehole collapse during drilling, mud circulation is often employed to ensure wellbore stability and prevent issues such as borehole collapse and stuck drilling [3]. Drilling mud plays a vital role in suspending drilling cuttings, cooling drill bits, maintaining stable wellbore walls, and preventing wellbore collapse throughout the drilling process. However, conventional bentonite mud, due to its low viscosity and density, as well as its poor rheological properties, fails to meet the requirements for drilling in gravel strata [4,5]. Therefore, effectively enhancing the rheological properties of drilling fluids is an urgent challenge that must be addressed. Currently, improvements in drilling issues are primarily focused on three facets: the enhancement of drilling equipment, the refinement of drilling techniques, and the optimization of drilling fluids [6,7,8]. The optimization of drilling fluids has generally improved their quality and performance by introducing various additives into conventional drilling fluids. Keren et al. [9] demonstrated that the addition of bentonite significantly impacts the rheological properties of drilling fluids. In a paper by Sun and Zhao [10], various scholars utilized foam, polymers, polymer suspensions, and high-density slurries as regulators in laboratory experiments targeting gravel stratum to determine the optimal mud ratios. Au et al. [4] studied the influence of bentonite content on the fundamental properties of mud and concluded that the bentonite content determines the yield stress magnitude of the kaolin–bentonite composite suspensions. However, this approach is neither economical nor environmentally friendly [11]. Experimental investigations on bentonite application in gravelly soil drilling have revealed an optimal range of mud dosages, which have been successfully implemented in practical engineering projects [10,12]. Akpan et al. [13] found that the addition of xanthan gum significantly enhances the viscosity and film-forming properties of mud, thus effectively reducing filtration loss [14]. Sun et al. [10] further analyzed the modification effect of carboxymethyl cellulose sodium (CMC) on mud. Nonetheless, current research predominantly emphasizes enhancing the cementing effect of drilling fluids through various additives, which often result in difficult preparation, environmental pollution, and non-biodegradability [15]. Rheological characteristics are fundamental properties of drilling mud that reflect its flow behavior and film-forming potential in various soil types [16]. The rheological properties of drilling fluids refer to their flow and deformation characteristics under external forces and are typically described by rheological parameters such as rheological curves, plastic viscosity, and yield point [17]. Good rheological properties effectively transport drilling debris, enhance mechanical drilling efficiency [18], and maintain wellbore integrity [19]. A common practice in drilling operations in gravel formations is to use a dispersion-blocking system for drilling muds, which is achieved by adding specific particulate materials as blocking agents to the base drilling fluid. However, this practice often has a significant impact on the environment because of the difficulty in degrading the components of these special particulate materials, which leads to serious environmental pollution issues, and the application costs are relatively high. To address this situation, we propose to introduce MICP technology into dispersion-based drilling fluid systems to improve the wall-consolidation performance of drilling fluids.
The MICP technique essentially utilizes specific urease-producing bacteria found in nature for decomposing urea to generate carbonate ions that react with metal ions in the environment to form carbonate precipitation [20,21,22]. It effectively reduces greenhouse gas emissions and improves environmental carbon sequestration efficiency [23,24,25,26,27]. Currently, Bacillus subtilis is commonly employed as the experimental strain due to its clear reaction mechanism, controllable reaction conditions, robust growth activity, and favorable environmental adaptability [22,26,27,28]. Chu et al. [29] proposed utilizing the MICP technique to pre-generate calcium carbonate particles, also referred to as “biogrout”, to reinforce coarse sand. Their experimental results demonstrated that this technique effectively reduces the permeability coefficient of coarse sand. Li [30] and Shi et al. [31] found that xanthan gum can delay the bacterial growth rate without affecting bacterial activity. Saracho et al. [32] highlighted the beneficial role of carboxymethyl cellulose in enhancing CaCO3 generation when the MICP technique is used. Choi. et al. [33] showed that the addition of CMC increases microbial calcite precipitation. Through experiments, Jia et al. [34] verified that incorporating montmorillonite effectively enhances CaCO3 precipitation compared to traditional MICP approaches [35]. The content of solid substances and the viscosity of biological mud directly influence the effectiveness of mud wall protection [36]. Drilling fluids with higher viscosity can significantly enhance bonding and wall protection in fractured and unstable strata [15,37]. The MICP technique facilitates the generation of calcite crystals, which possess excellent cementation properties. These crystals play a crucial role in forming calcium carbonate within loose and fractured bottom strata, thereby stabilizing the wellbore wall [37]. MICP technology offers a novel approach to addressing wellbore instability in gravel stratum. It could enable the creation of a microbial drilling fluid system that induces calcium carbonate precipitation within gravel stratum during drilling. This would enhance the mechanical strength of the wellbore, improve the integrity of the formation, and provide a new technical solution for effectively stabilizing wellbores in gravel stratum.
This study explored the potential of implementing the MICP technique into a dispersed drilling fluid system to prepare microbial drilling fluid, aiming to improve the wall-consolidation performance of drilling muds. The response surface methodology (RSM) was employed to analyze data and optimize the drilling fluid formulation [38]. The center composite design (CCD) method was utilized to analyze the effects of bentonite, biological solution, and barite powder on the rheological properties of drilling mud, including plastic viscosity (PV) (mPa·s ), density (g/cm3), and yield point (YP) (Pa). Finally, a quadratic regression equation was developed to fit each factor, leading to the determination of the optimal formulation based on RSM.

2. Materials and Methods

2.1. Test Materials

Table 1 describes the mixing ratios of the culture medium compositions that are utilized as primary experimental materials and the additive admixtures.. All experimental reagents utilized in this study, purchased from Sinopharm Chemical Reagent Co., Ltd., Shanghai, were of analytical grade. Polydimethylsiloxane lotion was employed as an additive to reduce the amount of foam generated during mixing. CMC and xanthan gum are chemical additives that are commonly used in drilling. Chen et al. [23,27] demonstrated that the addition of CMC increased calcium carbonate production in microbial clays. Saracho [28] and Xiao et al. [31] found that xanthan gum hinders bacterial growth without affecting activity.
In this study, Sporosarcina pasteurii, which was acquired from the China General Microbiological Culture Collection Center (CGMCC), was utilized as the bacterial culture. The composition of the bacterial activation culture medium was prepared based on the liquid culture medium provided by the CGMCC. The culture medium was sterilized with high-pressure steam at 121 °C for 20 min and then placed on the ultra-clean workbench until cooled [21]. The primary bacterial strain was inoculated into the culture medium and then shaken and cultured at 180 r/min and 30 °C for 24 h [23,31]. The bacteria concentration was expressed by measuring the absorbance of the bacterial solution (turbidimetric method). This principle was mainly based on the fact that the bacterial concentration was directly proportional to the turbidity of the bacterial solution. Therefore, it was also directly proportional to the absorbance. In this paper, the absorbance (OD600) of an ultraviolet visible spectrophotometer (model uv-Molecules, INESA Company, Shanghai, China) at a wavelength of 600 nm was used to determine the bacteria concentration. Following cultivation, the absorbance value of the bacterial solution at 600 nm (OD600) was measured to be 1.396.

2.2. Preparation of Biological Solution

To investigate the combined effect of XG and CMC on bacterial growth activity during MICP, different proportions of XG and CMC were introduced during bacterial culture. The 200 mL culture medium solution was subjected to high-pressure steam sterilization, after which XG and CMC were added in specific proportions. Mixed solutions of 0%:0%, 0.4%:0.2%, 0.6%:0.4%, 0.8%:0.6%, and 0.4%:0.6% were prepared in accordance with the mass of the culture medium. After mixing, 1 mL of bacterial solution was added to every 200 mL of culture medium, which was then placed in a constant-temperature oscillation box for cultivation. The unit urease activity and density OD600 of the bacteria were measured.
Bacterial activity was assessed using the conductivity method, and the average conductivity change per minute (ms·min−1) was recorded using the conductivity meter [28,31,39]. According to Whitfin’s empirical function, a conductivity change of 1 ms · min−1 corresponded to 11 mM of urea hydrolysis per minute [40]. Considering that the bacterial solution was diluted 10 times during the mixing process with urea, the conductivity value was multiplied by 10 to determine the amount of urea hydrolyzed per minute in the tested bacterial solution, which served as a measure of urease activity [40]. During the bacterial culture, the OD600 value was measured at regular intervals. When it reached approximately 1.500, the bacterial culture was considered to be complete. Then, the activity was measured after a cultivation duration of 9 h. After the completion of the bacterial culture, the OD600 was measured using a visible-light spectrophotometer to characterize the bacterial content. Due to the colorless and transparent nature of xanthan gum and CMC solutions, it was experimentally determined that the addition of these solutions had minimal effect on the determination of OD600.
The patterns of change in activity and OD600 were consistent. The highest-unit urease activity and OD600 values were observed at the end of the 9 h cultivation period in the culture medium without added admixtures, yielding values of 15.86 U/mL (1 U = 1 μM of urea hydrolyzed/min) and 1.696, respectively. The remaining four groups exhibited a slight decrease in activity compared to OD600 as the content increased, but the decrease was not significant, and the activity remained above 11.5 U/mL. When the cultivation time was shorter (4 h), the bacterial concentration in the culture medium without any added substances was the lowest at 0.418. It was speculated that this is due to xanthan gum acting similarly to nutrients in the culture medium, without adversely affecting the bacteria in the short term. In the other four experimental groups, the OD600 value increased as the mixture content in the medium increased, while the activity value remained relatively stable (Figure 1).
The results for the unit urease activity and OD600 value indicated that as the amounts of xanthan gum and CMC added increased, there was a slight effect on bacterial growth; however, the difference in unit urease activity was not significant. Thus, the addition of mixed powder slowed the bacterial growth rate and prolonged the bacterial culture time, but it had minimal effect on the bacterial growth activity. Simultaneously, it effectively increased the viscosity of the bacterial solution, which created favorable conditions for mud production. Therefore, 0.6% XG and 0.4% CMC were added during the cultivation process.
Urea and CaCl2 were added in moderate proportions to the cultured Bacillus subtilis solution; the urea provided a nitrogen source for biological growth and CaCl2 supplied calcium for the MICP process. The biological solution with added CMC and XG was utilized as one of the variables, while bentonite and barite powder were used as the other two variables to ensure the basic properties of the slurry and to enhance its viscosity and density. The required biological solution for the experiment was obtained, and rheological parameter testing was conducted according to the experimental plan derived from the response model.

2.3. Response Surface Experimental Design

For multivariate analysis, RSM provides a research approach that employs statistical methods to study parameter interactions [41]. RSM was initially utilized for experimental response modeling and has since been widely adopted for numerical experimental modeling [42]. Essentially, it encompasses a collection of statistical and mathematical techniques used to establish empirical models [43]. Central Composite Design (CCD) is a type of response surface design that can optimize multiple response variables, generate mathematical models, and determine the optimal solution for drilling fluid under various conditions [40,44].
To ensure the fundamental performance of the mud and enhance its viscosity and density, specific amounts of bentonite and barite powder were incorporated into the drilling fluid. Independent variables (bentonite (g), biological solution (mL), and barite powder (g)) were investigated at five levels with three replications in the center domain. Each independent variable was coded at five levels, as detailed in Table 2. A three-factor, five-level response surface design was conducted using a CCD model [38]. The experimental design of this model includes 8 factorial points, 8 axial points ± 2 distance from the center, and 5 repeated center points, which can be utilized to describe the relationship between the response (i.e., PV, density, and YP) and variables (bentonite, biological solution, and barite powder) [38,43].

2.4. Density and Rheological Parameter Determination

Relative density was measured using a density meter with a volume of 140 m3. The biological mud was mixed using the HTD-3S digital high-speed mixer at a speed of 800 rpm. Following the preparation of the mud, it was transferred to a mud cup, covered, and placed on a bracket. The vernier was adjusted to achieve a horizontal lever, and the left indicator displayed the relative density of the mud [17]. The following formula can also be used for calculation:
ρ x = m 3 m 1 m 2 m 1
where ρ x is the relative density; m 1 is the mass of any fixed container; m 2 is the mass of the container filled with clean water; m 3 is the mass of the container filled with drilling mud.
The rheological parameters were primarily measured using a ZNN-D6 rotational viscometer. In this experiment, all results were conducted at consistent ambient temperature and humidity conditions [45,46]. The rheological parameters (i.e., PV and YP) were calculated based on the Bingham plastic fluid model [19]. The Bingham model represents the relationship between shear stress and shear rate:
σ = σ y + η B γ ˙
where σ is the shear stress of the fluid at different shear rates; σ y is the resistance to initial flow and represents the stress required to initiate the fluid movement, i.e., the yield point; η B is the plastic viscosity of the fluid; and γ ˙ is the shear rate of the fluid [19].
The Bingham plasticity model is a widely employed rheological model that can aid in determining the plastic viscosity and yield point of various drilling fluids. This model can predict the minimum force required to initiate fluid flow [17]. The PV is calculated with a rotational viscometer based on Equation (3), while the YP is determined using Equation (4).
P V = φ 600 φ 300 m P a s
Y P = 0.511 φ 300 P V P a

3. Results and Discussion

3.1. Response Surface Test Results

Table 3 summarizes the operational combinations and corresponding results of the 36 conducted experiments. To minimize the impact of unexpected variability in the observed response, all experimental tests in this study were conducted based on a random arrangement. The statistical significance of the polynomial model was ascertained by evaluating the coefficient of determination (R2) and the F-statistic, which represents the ratio of the mean square to the residual mean square of the model. This analysis aided in evaluating the fitting accuracy of the generated polynomial model. Additionally, the significance of the model terms was evaluated using the P-value (probability value) at a 95% confidence level [38]. To determine the interrelationships between the variables, multivariate regression analysis was conducted, and the polynomial model was obtained using the least squares method. The regression coefficients were determined through analysis of variance (ANOVA) [41].
Figure 2 presents the correlation analysis heatmap among various factors based on experimental results. From this graph, we observe a strong positive correlation between plastic viscosity (PV) and the biological solution. This indicates that as the amount of biological solution increases, plastic viscosity also rises correspondingly. Conversely, plastic viscosity shows a negative correlation with barite content; that is, higher barite content results in lower plastic viscosity. Additionally, the correlation between plastic viscosity and bentonite content is weak, suggesting that bentonite has a limited direct impact on plastic viscosity. There is a strong positive correlation between density values and barite content, which is expected, as barite is a weighting agent that is commonly utilized to increase the density of drilling fluids. Furthermore, while density values also demonstrate some correlation with bentonite content, this relationship is not as pronounced as the one with barite. Notably, the addition of biological solutions significantly reduces density values, indicating that as biological solution content increases, the density of the drilling fluid decreases. The yield point (YP) is influenced by the combination of bentonite and biological solution content, with the latter having a more pronounced effect. This highlights the biological solution’s importance in adjusting the yield point of drilling fluids. To further illustrate these relationships, scatter plots showing the interactions between plastic viscosity, yield point, density, and biological solution content are presented in Figure 2, using biological mud content as a reference. These scatter plots offer concrete visual evidence to reinforce the findings of the correlation analysis discussed above.

3.2. Statistical Analysis

3.2.1. Plastic Viscosity Analysis

Plastic viscosity can be utilized to evaluate drilling fluid’s ability to suspend unstable cuttings, which represents a component of the flow resistance caused by mechanical friction, specifically referring to the shear stress exceeding the yield point [17]. For general drilling fluids, the plastic viscosity primarily depends on the solid concentration within the drilling fluid [15]. In the samples with higher solid content, the increased frictional force between particles leads to enhanced plastic viscosity. Moreover, within the same volume, smaller particle sizes result in larger surface areas, thereby increasing frictional force and consequently yielding higher PV values.
The experimental results for the first response value of PV were used to fit the model and obtain the best-fitting result through variance analysis. The optimal fitting equation is essentially quadratic [38]. The second-order polynomial mathematical model representing the relationship between the independent variable and the response variable PV (Y1) is presented in Equation (5):
Y1 = 49.03 + 1.85A + 36.43B + 0.58C + 0.77AB + 0.18AC − 0.21BC − 0.69A2 − 8.91B2
where A represents the bentonite content; B represents the biological solution content; and C represents the barite content. Among them, positive values of coefficients A, B, C, AB, and AC indicated that linear effects increase the PV, while coefficients BC, A2, and B2 exerted negative effects on the increase in PV.
The orthogonal analysis results for the PV response surface quadratic model are demonstrated in Table 4. The model demonstrated good accuracy at a 95% confidence level with an F-statistic of 675.99, indicating the significance and a minimal 0.01% probability of such a high F-statistic resulting from interference. Based on the F-statistic, the influence of the three factors on PV followed the order B > A > C. This indicates that the biological solution has a significantly greater impact compared to the contents of bentonite and barite. The analysis yielded an acceptable AP ratio of 103.7526, exceeding the threshold value of four.
To evaluate the accuracy of the model, a diagnostic chart comparing predicted and actual values was utilized, which confirmed the acceptable accuracy of the model [41]. In cases where the fitting degree of the model is insufficient, continued optimization of the response surface could yield adverse consequences [44]. The appropriate distribution of data points near the line represents the model’s ability to predict response values [43]. Figure 3d illustrates the diagnostic graph comparing the predicted and actual PV values. The experimental design points were distributed diagonally, which indicated that the actual PV values align with the predicted values of the model. According to Table 4, the R2 and Adj.R2 values were similar, with a difference of less than 0.2, which indicated that the generated mathematical model can effectively reflect the relationship between mud viscosity and various factors and can predict the optimal viscosity.
Figure 3e and Figure 4f demonstrated that the residuals obtained from the operating results conform to the basic assumption of the error term, with an R2 of the residuals fitted at 0.95684. It is evident that the error is within an acceptable range, and there is no violation of the assumptions of independence or constant variables [36]. In addition, the interaction between various influencing factors and PV is illustrated in Figure 3, where the addition of each material ultimately increases the solid content. Figure 3a–c depict the relationship between plastic viscosity (Y1) and bentonite (A), biological solution (B), and barite powder (C). Figure 3a,c show that as the containment of the biological solution increased, the PV value gradually increased and exhibited a rapid increase in the range of 0 mL to 200 mL. The additional amount significantly influenced the PV value. However, in the range of 200 mL to 400 mL, while the PV value continued to increase, the growth rate gradually slowed down, and the curvature of the surface decreased. Figure 3 illustrates that bentonite positively affects the increase in PV value, but the effect is not as significant compared to biological solutions. The barite content had a minimal effect on the PV value.

3.2.2. Density Analysis

The fitting equation derived from the variance analysis on density is presented in Equation (6):
Y2 = 1.1714 + 0.0096A − 0.0240B + 0.0439C − 0.0021AB − 0.0004AC −
0.0018BC + 0.0004A2 − 0.0001B2 + 0.0015C2
where positive values of coefficients A, C, BC, and C2 indicate that linear effects increase the density, while coefficients B, AB, AC, A2, and B2 exert negative effects on the increase in density.
The orthogonal analysis results of the response surface quadratic model are displayed in Table 5. The F-statistic of the model was 408.09, indicating significant results. The p-values of A, B, and C were all less than 0.0001, designating them as important model terms. The importance of the three factors on the response results was as follows: C > B > A. Barite C had a significant impact on density, while biological solution and density had a relatively minor effect. The R2 value for this response was 0.9930, which was close to 1. The Adj.R2 value was 0.9905, which was in good agreement with the R2 value. The AP ratio was 81.7644, exceeding 4, which was acceptable. A ratio greater than four was expected and indicated model differentiation, which suggested that the model can guide the design space defined by CCD [38].
Figure 4 illustrates the interaction between various influencing factors and the density of drilling fluid. Figure 4a–c show the relationships between bentonite (A), biological solution (B), barite powder (C), and density (Y2). The density of drilling fluid was positively correlated with the contents of bentonite and barite. The addition of bentonite resulted in a slight increase in density, though the effect was not significant. Conversely, the addition of barite effectively increased the density. There was a negative correlation between biological solution and density value, and the addition of biological solution significantly decreased the density of drilling fluid. An increase in biological solution reduces the density; however, it may slightly increase the density at lower concentrations. This result aligned with expectations, as barite possessed a very high density (typically 4.3 × 103 kg/m3) and functioned as a density enhancer in water-based drilling fluids [47].
Figure 4d presents a diagnostic chart of predicted and actual densities, where points above the diagonal indicate overestimation, while points below the diagonal indicate underestimation [44,48]. The actual response values were relatively consistent with the predicted values, which indicates that the model provides a good estimate. The predicted values closely aligned with the actual experimental value, and the model yielded a high degree of fit, which can be utilized for analysis and prediction in this experiment. Figure 4e,f reveal that the R2 value of the residual fitting is 0.98115, and the error term follows a normal distribution.

3.2.3. The Yield Point Analysis

The yield point is an important component of the flow resistance of drilling fluid, which is also referred to as yield value. It reflects drilling fluid’s ability to form structures during laminar flow and indicates the drilling fluid’s capacity to carry and suspend drilling cuttings [17]. These forces arise from the negative and positive charges on the surface of clay particles, electrolytes, and polymers within the drilling fluid [9,19].
The experimental results for the third response value, i.e., YP, were used to fit the model and obtain the best-fitting result through analysis of variance. The second-order polynomial mathematical model representing the relationship between independent variables A, B, and C and the response variable YP (Y3) is expressed in Equation (7):
Y3 = 59.6259 + 5.6167A + 44.2132B − 1.3064C + 2.8236AB + 0.04338AC +
0.4196BC − 0.5138A2 − 12.1055B2 − 0.0110C2
where positive values of coefficients A, B, AB, AC, and BC indicate that the linear effect increases the yield point, while coefficients C, A2, B2, and C2 decrease the yield point.
The orthogonal analysis results of the response surface quadratic model of YP are displayed in Table 6. The F-statistic of the model was 654.37, which indicates that the model is significant at a 95% confidence level, and the probability of such a large “model F-statistic” was only 0.01% due to interference [43]. According to the F-statistic, the order of influence of the three factors on plastic viscosity was B > A > C, which showed that the biological solution exerts a much greater influence than barite and bentonite contents. Figure 5d shows a diagnostic graph comparing the predicted and actual values of YP. The experimental design points were distributed diagonally, indicating that the actual YP values align with the predicted values of the model [41]. Table 6 demonstrated that the R2 and Adj.R2 values are similar, with a difference of less than 0.2. Figure 5e,f indicate that the mathematical model is reliable, and the residual distribution exhibits normality. The R2 value of the residual fitting was 0.97109, which revealed that the error is not significant.
Figure 5a–c show that for both bentonite and biological solutions, YP increased as the amount of content increased. When the content of the biological solution was low, the interaction effect of bentonite on YP was relatively minor. As the content of the biological solution increased, the effect of bentonite on the YP value gradually became significant. An increase in barite content led to a decrease in YP value; however, the magnitude of the change was relatively small. When the bentonite content was maintained at a constant level, the YP value initially increased rapidly along with the content of the biological solution before the growth rate slowed down; overall, however, it underwent an increasing trend.

3.3. Analysis of Optimization Results

3.3.1. Model Validation

Given that the primary properties of drilling fluids include density and rheological properties, separate optimization analyses were conducted for density and rheological properties. Under the premise of meeting quality requirements while ensuring economic and environmental protection, the optimal solution was selected to determine the optimal concentration of each material, so as to provide the best microbial drilling fluid formulation suitable for drilling applications in loose formations. The optimal ratio of bentonite mud under the conditions of maximum viscosity, maximum density, and minimum yield point was obtained with the aid of the optimization function in the Design-Expert 12.0.3.0 software. Specifically, the optimal content was as follows: 50 g of bentonite, 200 mL of biological solution, and 120 g of barite. At this point, the predicted viscosity, density, and yield point of mud were 40.174 mPa·s, 1.223 g/cm3, and 46.584 Pa, respectively. The results obtained from three sets of repeated tests are presented in Table 7. The actual test results matched well with the predicted values derived from the response surface-based model, validating the accuracy of the model’s predictions. This confirmed that the results for plastic viscosity, density, and yield point predicted by the response surface are correct and that the response model results are of practical reference value. The b-value of the dynamic–plastic ratio roughly reflect the strength of mud shear dilution. The b-value obtained through testing met the requirements with the standards. [49].

3.3.2. Characterization of Biological Mud

The mud cake was obtained by subjecting the optimally formulated drilling fluid to a gas pressure of 0.69 MPa for 30 min at a recorded filtration loss of 5 mL in accordance with the API standard. The resulting mud cake was dried using the Critical Point Dryer, Tousimis Autosamdri-815 Series, Rockville, MD, USA, and the samples were subsequently gold-sputtered using the Ion Sputterer, Hitachi High-Technologies Corporation, Tokyo, Japan.. Finally, the samples were subjected to scanning electron microscope (SEM) testing with the SU8010 Field Emission Scanning Electron Microscope, Hitachi High-Technologies Corporation, Tokyo, Japan to observe their microstructures.
Figure 6 displays the energy-dispersive spectroscopy (EDS) spectrum, with a small image in the upper-right corner displaying the SEM image of the test material. Figure 6a,c depict samples incorporated with biological solution, while Figure 6b,d represent ordinary mud samples without added biological solution. In Figure 6a,c, characteristic peaks corresponding to calcium (Ca), sulfur (S), and barium (Ba) are evident. Notably, there is a significant increase in the peaks for elemental calcium compared to those of the untreated samples shown in Figure 6b,d. This suggests a marked rise in calcium content in the mud samples treated with the biological solution, potentially due to calcium being released from the solution or due to its reaction with the mud cake samples, aligning with the experimental expectations. In contrast, the changes in other elements aside from calcium were relatively minor. Specifically, in Figure 6b,d, the atomic ratio of sulfur (S) to barium (Ba) is 1, consistent with the atomic ratio of the molecular formula for barium sulfate, the primary component of the selected barite powder. This indicates that the large granular material with a smooth surface seen in the sample image is the barite component.
Additionally, experimental observations revealed that the addition of the biological solution caused the slurry to undergo volumetric expansion during high-speed mixing, leading to a decrease in density. Given that colloidal substances tend to expand, it is supposed that this phenomenon resulted from the introduction of the biological media and colloidal materials.
Figure 7 illustrates the mud cake samples obtained from the previously described experiments. Figure 7a reveals that the generated calcium carbonate crystals are uniformly distributed on the sample surface, with particle sizes ranging from approximately 1 to 5 μm. Under the influence of microbial solidification, abundant calcium carbonate crystals can be recognized on the surfaces and interfaces of solid particles and fibers, leading to a compact structure [29,34]. The small particle sizes of calcium carbonate solid satisfy the requirement that a smaller solid phase in drilling fluid yields a larger specific surface area, thereby increasing friction force and YP. The better the rheological properties of drilling fluid, the more favorable it is for engineering applications. It is speculated that the irregular shape of calcium carbonate crystals is due to the recombination effect on the calcium carbonate crystals produced by biological reactions during high-speed stirring. Distinct fibrous connections can be observed in Figure 7b,c, where Figure 7b exhibits a ribbon-like connection that aligns with the typical microscopic appearance of xanthan colloid. The locally enlarged image in Figure 7c shows a number of small solid particles distributed along the fiber connectors formed by the colloid, presumed to be calcium carbonate precipitates. The figure illustrates a large number of clustered calcium carbonate precipitate crystals at both ends of the connecting colloid, with additional crystal distribution on the connecting body. This indicates that the addition of colloid components provides more nucleation sites in the mud, resulting in uniform precipitation distribution and increased YP. The interaction between colloids and small solid particles forms a combination in the fluid, creating large solid clusters that contribute to an increase in the PV of the drilling fluid. This finding corroborates the conclusions reported by Mandala et al. and Bayas, who state that “there is a synergistic effect between the two biopolymers, leading to the formation of stable associations that increase PV” [15,36,50,51].

3.3.3. Biodegradability and Heavy Metal Detection Results

The drilling fluid prepared based on the optimized process formulation has been proven to meet the requirements of national standards regarding physical, chemical, and microbiological indicators. It does not contain heavy metals and falls within the permissible safety range established by national regulations. The drilling fluid underwent tests for biodegradability, and heavy metal elements (cadmium, chromium, lead, and arsenic) were assessed separately.
According to the Technical Specification for Environmental Performance Evaluation of Drilling Fluids, a DXY-3 biological toxicity tester was employed to monitor the BOD5 (Biochemical Oxygen Demand, an important indicator of organic matter pollution in water that is indirectly represented by the dissolved oxygen consumed by microbial metabolism) and COD (Chemical Oxygen Demand, a chemical method for measuring the amount of reducing substances that need to be oxidized in water samples) of the biological mud drilling fluid. The results, presented in Table 8, are in accordance with the “Identification Standard for Hazardous Wastes” and the “Technical Specification for Environmental Performance Evaluation of Drilling Fluid”.

3.3.4. Analysis of Cost Advantages

The raw materials for microbial drilling fluid are easily sourced and readily available. The unit cost of producing 1 kg of microbial drilling fluid can be calculated as Table 9:
We have selected other types of drilling fluids that are top sellers in the market and calculated their unit prices, as shown in the Table 10:
From the table above, it is evident that the unit cost of microbial drilling fluid is significantly lower than that of other drilling fluids on the market. This shows that microbial drilling fluid has a clear cost advantage. In practice, microbial drilling fluid is particularly suitable for gravel formations with larger pores, which are more sensitive to pressure changes. However, before this fluid can be utilized, attention must be given to the cultivation of microbial fluids. Currently, the expansion of primary bacteria still requires laboratory work. Some researchers are focusing on the induced domestication of microbial strains to identify those better suited for deep-well environments. With ongoing advancements in culture technology, it is expected that bacterial fluids will soon be available directly on the market, avoiding the need for laboratory processes. This would broaden the potential applications of biological mud significantly.

4. Discussion

This study evaluated the optimal preparation formulation for biological mud using the CCD and explored the mechanism by which bio-cementation improves drilling fluid performance at the microscopic level. Analyses of the experimental results yielded the following conclusions:
  • The analyses of the rheological and density test results of drilling fluids containing different amounts of biological solution indicated that the PV of biological drilling fluids increased with the addition of biological solutions, exhibiting a growth rate that initially increased and then decreased. The density increased monotonically with the addition of barite, while the increased rate of YP exhibited a trend of first increasing and then decreasing as the biological solution content increased. Interactions among various factors were evident, and the effect of biological solutions on density was manifested as an increase in response to low content and a decrease in response to high content, which was significantly influenced by the interaction. The influence of bentonite and barite on YP exhibited a trend of first increasing and then decreasing; as such, they could generally offset each other. However, as amount of biological solution increased, the influence of bentonite on the YP gradually became significant.
  • Based on the response surface optimization results, the formulation containing bentonite, biological solution, and barite in quantities of 50 g, 200 mL, and 120 g, respectively, was optimal, and the rheological parameters of drilling fluid under this formulation were obtained through experimental measurements. The plastic viscosity, density, and yield point were 40.174 mPa·s, 1.223 g/cm3, and 46.584 Pa, respectively, all of which met the regulatory requirements and expected standards.
  • The experimental results indicate that the addition of colloidal substances does not affect the growth activity of Bacillus subtilis. The microstructure of the mud cake samples obtained from the experiment was analyzed using SEM. The SEM and EDS analyses demonstrated that the addition of colloidal substances does not hinder the formation of solid calcium carbonate precipitates and facilitates the formation of fibrous and ribbon connections within the mud, improving the plastic viscosity of the mud. The generated small calcium carbonate solids contributed to the increased yield stress of the drilling fluid. After inspection, the mud was found to comply with national environmental protection requirements, and its heavy metal content and biodegradability met the standard requirements. The unit cost of the mud was deemed to be lower, indicating its cost advantage in terms of pricing.

Author Contributions

Conceptualization, R.P., Z.S., Y.C., X.S. and Y.H.; data curation, R.P., X.Z. and Y.H.; methodology, R.P., Z.S., Y.C. and X.S.; software, R.P., X.S., X.Z. and Y.H.; validation, Z.S., Y.C., X.S. and Y.H.; formal analysis, R.P. and Y.C.; resources, X.Z. and Y.H.; writing—original draft preparation, R.P. and Z.S.; writing—review and editing, R.P., Z.S., Y.C., X.S., X.Z. and Y.H.; supervision, Y.C., X.S., X.Z. and Y.H.; project administration, Y.C., X.S. and X.Z.; funding acquisition, Z.S., Y.C. and X.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Natural Science Foundation of China (Grant No. 51879090 and Grant No. 52179101), the National Natural Science Foundation of China (Grant No. 41831282), Jiangsu Province Natural Science Foundation (BK20220975).

Data Availability Statement

The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors, Zhou Shu, Xiaobing Sha and Xinquan Zhang, are employees of China Railway Major Bridge Reconnaissance & Design Institute Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ANOVAAnalysis of variance.
APIAmerican Petroleum Institute.
CCDCentral Composite Design.
CGMCCChina General Microbiological Culture Collection Center.
CMCCarboxymethyl cellulose sodium.
EDSEnergy-dispersive spectroscopy.
MICPMicrobially induced calcite carbonate precipitation.
OD600The absorbance value of the bacterial solution at 600 nm.
PVPlastic viscosity.
P-valueProbability value.
R2Coefficient of determination.
RSMResponse surface methodology.
SEMScanning electron microscopy.
XGXanthan gum.
YPYield point.

References

  1. Yu, L.; Han, Z.; He, J.; Wang, W. Experimental Study on Sand Control Gravel Particle Gradation Optimization in Mud Hydrate Reservoir. Energy Rep. 2022, 8, 9201–9210. [Google Scholar] [CrossRef]
  2. Sehly, K.; Chiew, H.-L.; Li, H.; Song, A.; Leong, Y.-K.; Huang, W. Stability and Ageing Behaviour and the Formulation of Potassium-Based Drilling Muds. Appl. Clay Sci. 2015, 104, 309–317. [Google Scholar] [CrossRef]
  3. Yeu, W.J.; Katende, A.; Sagala, F.; Ismail, I. Improving Hole Cleaning Using Low Density Polyethylene Beads at Different Mud Circulation Rates in Different Hole Angles. J. Nat. Gas Sci. Eng. 2019, 61, 333–343. [Google Scholar] [CrossRef]
  4. Au, P.-I.; Leong, Y.-K. Rheological and Zeta Potential Behaviour of Kaolin and Bentonite Composite Slurries. Colloids Surf. Physicochem. Eng. Asp. 2013, 436, 530–541. [Google Scholar] [CrossRef]
  5. Lee, C.E.; Chandra, S.; Leong, Y.-K. Structural Recovery Behaviour of Kaolin, Bentonite and K-Montmorillonite Slurries. Powder Technol. 2012, 223, 105–109. [Google Scholar] [CrossRef]
  6. Stilwell, W.B. Drilling Practices and Equipment in Use at Wairakei. Geothermics 1970, 2, 714–724. [Google Scholar] [CrossRef]
  7. Kim, D.W.; Lee, Y.S.; Park, M.S.; Chu, C.N. Tool Life Improvement by Peck Drilling and Thrust Force Monitoring during Deep-Micro-Hole Drilling of Steel. Int. J. Mach. Tools Manuf. 2009, 49, 246–255. [Google Scholar] [CrossRef]
  8. Kong, L.; Cao, S.; Chin, J.-H.; Si, Y.; Miao, F.; Li, Y. Vibration Suppression of Drilling Tool System during Deep-Hole Drilling Process Using Independence Mode Space Control. Int. J. Mach. Tools Manuf. 2020, 151, 103525. [Google Scholar] [CrossRef]
  9. Keren, R. Rheology of Mixed Kaolinite-Montmorillonite Suspensions. Soil Sci. Soc. Am. J. 1989, 53, 725–730. [Google Scholar] [CrossRef]
  10. Sun, Y.; Zhao, D. Research and Experimental Application of New Slurry Proportioning for Slag Improvement of EPB Shield Crossing Sand and Gravel Layer. Coatings 2022, 12, 1961. [Google Scholar] [CrossRef]
  11. Lagaly, G. Principles of Flow of Kaolin and Bentonite Dispersions. Appl. Clay Sci. 1989, 4, 105–123. [Google Scholar] [CrossRef]
  12. Budach, C.; Thewes, M. Application Ranges of EPB Shields in Coarse Ground Based on Laboratory Research. Tunn. Undergr. Space Technol. 2015, 50, 296–304. [Google Scholar] [CrossRef]
  13. Akpan, E.U.; Enyi, G.C.; Nasr, G.G. Enhancing the Performance of Xanthan Gum in Water-Based Mud Systems Using an Environmentally Friendly Biopolymer. J. Pet. Explor. Prod. Technol. 2020, 10, 1933–1948. [Google Scholar] [CrossRef]
  14. Khan, R.; Kuru, E.; Trembla, B.; Saasen, A. An Investigation of Formation Damage Characteristics of Xanthan Gum Solutions Used for Drilling, Drill-in, Spacer Fluids, and Coiled Tubing Applications; OnePetro: Richardson, TX, USA, 2003. [Google Scholar]
  15. Baba Hamed, S.; Belhadri, M. Rheological Properties of Biopolymers Drilling Fluids. J. Pet. Sci. Eng. 2009, 67, 84–90. [Google Scholar] [CrossRef]
  16. Yao, L.; Anne Naeth, M.; Jobson, A. Soil Microbial Response to Waste Potassium Silicate Drilling Fluid. J. Environ. Sci. 2015, 29, 189–198. [Google Scholar] [CrossRef] [PubMed]
  17. Agwu, O.E.; Akpabio, J.U.; Ekpenyong, M.E.; Inyang, U.G.; Asuquo, D.E.; Eyoh, I.J.; Adeoye, O.S. A Critical Review of Drilling Mud Rheological Models. J. Pet. Sci. Eng. 2021, 203, 108659. [Google Scholar] [CrossRef]
  18. Crow, W.; Carey, J.W.; Gasda, S.; Brian Williams, D.; Celia, M. Wellbore Integrity Analysis of a Natural CO2 Producer. Int. J. Greenh. Gas Control 2010, 4, 186–197. [Google Scholar] [CrossRef]
  19. Liu, N.; Zhang, D.; Gao, H.; Hu, Y.; Duan, L. Real-Time Measurement of Drilling Fluid Rheological Properties: A Review. Sensors 2021, 21, 3592. [Google Scholar] [CrossRef] [PubMed]
  20. Yu, T.; Souli, H.; Pechaud, Y.; Fleureau, J.-M. Review on Engineering Properties of MICP-Treated Soils. Geomech. Eng. 2021, 27, 13–30. [Google Scholar] [CrossRef]
  21. Liu, H.; Xiao, P.; Xiao, Y.; Wang, J.; Chen, Y.; Chu, J. Dynamic behaviors of MICP-treated calcareous sand in cyclic tests. Chin. J. Geotech. Eng. 2018, 40, 38–45. [Google Scholar] [CrossRef]
  22. Liu, B.; Zhu, C.; Tang, C.-S.; Xie, Y.-H.; Yin, L.-Y.; Cheng, Q.; Shi, B. Bio-Remediation of Desiccation Cracking in Clayey Soils through Microbially Induced Calcite precipitatMicrobially Induced Carbonate Precipitation (MICP) for Soil Strengtheniion (MICP). Eng. Geol. 2020, 264, 105389. [Google Scholar] [CrossRef]
  23. Chen, L.; Song, Y.; Fang, H.; Feng, Q.; Lai, C.; Song, X. Systematic Optimization of a Novel, Cost-Effective Fermentation Medium of Sporosarcina pasteurii for Microbially Induced Calcite Precipitation (MICP). Constr. Build. Mater. 2022, 348, 128632. [Google Scholar] [CrossRef]
  24. Comadran-Casas, C.; Brüggemann, N.; Jorat, M.E. Greenhouse Gas Fluxes of Microbial-Induced Calcite Precipitation at Varying Urea-to-Calcium Concentrations. Eur. J. Soil Sci. 2024, 75, e13516. [Google Scholar] [CrossRef]
  25. Fouladi, A.S.; Arulrajah, A.; Chu, J.; Horpibulsuk, S. Application of Microbially Induced Calcite Precipitation (MICP) Technology in Construction Materials: A Comprehensive Review of Waste Stream Contributions. Constr. Build. Mater. 2023, 388, 131546. [Google Scholar] [CrossRef]
  26. Deng, X.; Li, Y.; Liu, H.; Zhao, Y.; Yang, Y.; Xu, X.; Cheng, X.; Wit, B. de Examining Energy Consumption and Carbon Emissions of Microbial Induced Carbonate Precipitation Using the Life Cycle Assessment Method. Sustainability 2021, 13, 4856. [Google Scholar] [CrossRef]
  27. Chen, L.; Song, Y.; Huang, J.; Lai, C.; Jiao, H.; Fang, H.; Zhu, J.; Song, X. Critical Review of Solidification of Sandy Soil by Microbially Induced Carbonate Precipitation (MICP). Crystals 2021, 11, 1439. [Google Scholar] [CrossRef]
  28. Fu, T.; Saracho, A.C.; Haigh, S.K. Microbially Induced Carbonate Precipitation (MICP) for Soil Strengthening: A Comprehensive Review. Biogeotechnics 2023, 1, 100002. [Google Scholar] [CrossRef]
  29. Chu, J.; Ivanov, V.; Naeimi, M.; Stabnikov, V.; Liu, H.-L. Optimization of Calcium-Based Bioclogging and Biocementation of Sand. Acta Geotech. 2014, 9, 277–285. [Google Scholar] [CrossRef]
  30. Li, F.; Li, X.; Ye, L.; Liu, X.; Zhu, J.; Yang, S.; Yan, Y.; Xu, L.; Yan, J. A Genetically Engineered Composite Biofilm for Microbial Induced Calcium Carbonate Precipitation by Synergic Effect of Urease, Protein Adhesive and Xanthan Gum. J. Environ. Chem. Eng. 2022, 10, 108431. [Google Scholar] [CrossRef]
  31. Shi, J.; Xiao, Y.; Fu, G.; Feng, C.; Hu, J.; Haegeman, W.; Liu, H. Calcareous Silt Earthen Construction Using Biopolymer Reinforcement. J. Build. Eng. 2023, 72, 106571. [Google Scholar] [CrossRef]
  32. Saracho, A.C.; Lucherini, L.; Hirsch, M.; Peter, H.M.; Terzis, D.; Amstad, E.; Laloui, L. Controlling the Calcium Carbonate Microstructure of Engineered Living Building Materials. J. Mater. Chem. A 2021, 9, 24438–24451. [Google Scholar] [CrossRef] [PubMed]
  33. Han, Z.; Cheng, X.; Ma, Q. An Experimental Study on Dynamic Response for MICP Strengthening Liquefiable Sands. Earthq. Eng. Eng. Vib. 2016, 15, 673–679. [Google Scholar] [CrossRef]
  34. Tang, G.; Jia, C.; Wang, G.; Yu, P.; Zhang, H. Role of Na-Montmorillonite on Microbially Induced Calcium Carbonate Precipitation. Molecules 2021, 26, 6211. [Google Scholar] [CrossRef] [PubMed]
  35. Duan, Y.; Niu, L.; Li, B.; He, Y.; Xu, X.; Yu, C.; Wang, Z.; Xiao, C.; Zheng, C. Montmorillonite-Coupled Microbially Induced Carbonate Precipitation (MICP) Enhanced Contaminant Removal and Carbon Capture in Cyanide Tailings. J. Environ. Chem. Eng. 2024, 12, 113498. [Google Scholar] [CrossRef]
  36. Lim, S.N.; Khalil, M.; Jan, B.M.; Ali, B.S. Lightweight Biopolymer Drilling Fluid for Underbalanced Drilling: An Optimization Study. J. Pet. Sci. Eng. 2015, 129, 178–188. [Google Scholar] [CrossRef]
  37. Phillips, A.J.; Troyer, E.; Hiebert, R.; Kirkland, C.; Gerlach, R.; Cunningham, A.B.; Spangler, L.; Kirksey, J.; Rowe, W.; Esposito, R. Enhancing Wellbore Cement Integrity with Microbially Induced Calcite Precipitation (MICP): A Field Scale Demonstration. J. Pet. Sci. Eng. 2018, 171, 1141–1148. [Google Scholar] [CrossRef]
  38. Saikia, T.; Mahto, V. Experimental Investigations and Optimizations of Rheological Behavior of Drilling Fluids Using RSM and CCD for Gas Hydrate-Bearing Formation. Arab. J. Sci. Eng. 2018, 43, 6541–6554. [Google Scholar] [CrossRef]
  39. Wang, K.; Wu, S.; Chu, J. Mitigation of Soil Liquefaction Using Microbial Technology: An Overview. Biogeotechnics 2023, 1, 100005. [Google Scholar] [CrossRef]
  40. Onal Okyay, T.; Frigi Rodrigues, D. Optimized Carbonate Micro-Particle Production by Sporosarcina pasteurii Using Response Surface Methodology. Ecol. Eng. 2014, 62, 168–174. [Google Scholar] [CrossRef]
  41. Salehnezhad, L.; Heydari, A.; Fattahi, M. Experimental Investigation and Rheological Behaviors of Water-Based Drilling Mud Contained Starch-ZnO Nanofluids through Response Surface Methodology. J. Mol. Liq. 2019, 276, 417–430. [Google Scholar] [CrossRef]
  42. Tartibu, L.K.; Sun, B.; Kaunda, M.A.E. Multi-Objective Optimization of the Stack of a Thermoacoustic Engine Using GAMS. Appl. Soft Comput. 2015, 28, 30–43. [Google Scholar] [CrossRef]
  43. Wang, X.; Zhang, Y.; Zhao, W.; Wang, Z.; Wang, Z.; Wang, Y. Research on Optimizing Performance of New Slurries for EPBS Soil Conditioning Based on Response Surface Method. Constr. Build. Mater. 2023, 375, 130818. [Google Scholar] [CrossRef]
  44. Pan, Z.; Cao, S. Optimization of Culture Medium to Improve Bio-Cementation Effect Based on Response Surface Method. Sci. Rep. 2024, 14, 8752. [Google Scholar] [CrossRef]
  45. Xie, B.; Ting, L.; Zhang, Y.; Liu, C. Rheological Properties of Bentonite-Free Water-Based Drilling Fluids with Novel Polymer Viscosifier. J. Pet. Sci. Eng. 2018, 164, 302–310. [Google Scholar] [CrossRef]
  46. Chu, Q.; Lin, L. Synthesis and Properties of an Improved Agent with Restricted Viscosity and Shearing Strength in Water-Based Drilling Fluid. J. Pet. Sci. Eng. 2019, 173, 1254–1263. [Google Scholar] [CrossRef]
  47. Yap, J.; Leong, Y.-K.; Liu, J. Structural Recovery Behavior of Barite-Loaded Bentonite Drilling Muds. J. Pet. Sci. Eng. 2011, 78, 552–558. [Google Scholar] [CrossRef]
  48. Yao, D.; Wu, J.; Wang, G.; Wang, P.; Zheng, J.-J.; Yan, J.; Xu, L.; Yan, Y. Effect of Wool Fiber Addition on the Reinforcement of Loose Sands by Microbially Induced Carbonate Precipitation (MICP): Mechanical Property and Underlying Mechanism. Acta Geotech. 2021, 16, 1401–1416. [Google Scholar] [CrossRef]
  49. Ali, I.; Ahmad, M.; Ridha, S.; Iferobia, C.C.; Lashari, N. Enhancing Drilling Mud Performance through CMITS-Modified Formulations: Rheological Insights and Performance Optimization. RSC Adv. 2023, 13, 32904–32917. [Google Scholar] [CrossRef]
  50. George, A.; Sanjay, M.R.; Srisuk, R.; Parameswaranpillai, J.; Siengchin, S. A Comprehensive Review on Chemical Properties and Applications of Biopolymers and Their Composites. Int. J. Biol. Macromol. 2020, 154, 329–338. [Google Scholar] [CrossRef]
  51. Benzerara, M.; Guihéneuf, S.; Belouettar, R.; Perrot, A. Combined and synergic effect of algerian natural fibres and biopolymers on the reinforcement of extruded raw earth. Constr. Build. Mater. 2021, 289. [Google Scholar] [CrossRef]
Figure 1. Sample preparation.
Figure 1. Sample preparation.
Processes 13 00162 g001
Figure 2. Correlation analysis between various factors and responses. scatter plots: (a) the interactions between B and PV; (b) the interactions between B and D; (c) the interactions between B and YP. (A: Bentonite; B: biological solution; C: barite; PV: plastic viscosity; D: density; YP: yield point. The symbol ‘-’ indicates a negative correlation between variables; “Run” refers to the number of experiments conducted).
Figure 2. Correlation analysis between various factors and responses. scatter plots: (a) the interactions between B and PV; (b) the interactions between B and D; (c) the interactions between B and YP. (A: Bentonite; B: biological solution; C: barite; PV: plastic viscosity; D: density; YP: yield point. The symbol ‘-’ indicates a negative correlation between variables; “Run” refers to the number of experiments conducted).
Processes 13 00162 g002
Figure 3. Three-dimensional response surface plots for various factors: (a) A and B, (b) A and C, (c) B and C; (d) diagnostic chart of predicted and actual values of plastic viscosity; (e) the normal graph of residuals; (f) residual and running chart.
Figure 3. Three-dimensional response surface plots for various factors: (a) A and B, (b) A and C, (c) B and C; (d) diagnostic chart of predicted and actual values of plastic viscosity; (e) the normal graph of residuals; (f) residual and running chart.
Processes 13 00162 g003
Figure 4. Three-dimensional response surface plots of various factors: (a) A and B, (b) A and C, (c) B and C; (d) diagnostic chart of predicted and actual densities; (e) the normal graph of residuals; (f) residual and running chart.
Figure 4. Three-dimensional response surface plots of various factors: (a) A and B, (b) A and C, (c) B and C; (d) diagnostic chart of predicted and actual densities; (e) the normal graph of residuals; (f) residual and running chart.
Processes 13 00162 g004
Figure 5. Three-dimensional response surface plots of various factors: (a) A and B, (b) A and C, (c) B and C; (d) diagnosis chart of predicted and actual values of yield point; (e) the normal graph of residuals; (f) residual and running chart.
Figure 5. Three-dimensional response surface plots of various factors: (a) A and B, (b) A and C, (c) B and C; (d) diagnosis chart of predicted and actual values of yield point; (e) the normal graph of residuals; (f) residual and running chart.
Processes 13 00162 g005
Figure 6. EDS spectrum of biological mud (ad).
Figure 6. EDS spectrum of biological mud (ad).
Processes 13 00162 g006
Figure 7. Microscopic morphology diagram of biological mud: (a) Calcium carbonate distribution; (b) ribbon-like connection; (c) calcium carbonate distribution on the fibrous connector.
Figure 7. Microscopic morphology diagram of biological mud: (a) Calcium carbonate distribution; (b) ribbon-like connection; (c) calcium carbonate distribution on the fibrous connector.
Processes 13 00162 g007
Table 1. Medium materials and experimental raw materials.
Table 1. Medium materials and experimental raw materials.
Raw MaterialsContents
Composition of culture medium (100 mL)Deionized water93 mL
0.5 g/L MnSO4·H2O solution2 mL
Yeast extract2 g
NH4Cl1 g
1.2 g/L NiCl2·6H2O solution2 mL
1 M/L NaOH solution3 mL
Other experimental reagentsBentonite-
Barite-
Carboxymethylcellulose Sodium (CMC)-
Xanthan gum (XG)-
Polydimethylsiloxane lotion-
Table 2. Experimental range and levels of independent variables.
Table 2. Experimental range and levels of independent variables.
VariablesFactor Level
NumberNameMinimumCoded LowMeanCoded HighMaximum
ABentonite16.3630.0047.3070.0083.64
BBiological solution0.000.00238.07500.00670.45
CBarite powder0.007591.72125.00142.05
Table 3. Operational combinations of independent variables in the CCD experiment.
Table 3. Operational combinations of independent variables in the CCD experiment.
Group NumberInfluencing FactorsResponse Index
ABCPVDensityYP
130.005.504.0052.001.1968.406
290.008.006.005.001.280.175
370.005.508.0064.001.1972.349
410.008.006.0046.001.0161.905
550.008.006.0047.001.0365.115
650.008.006.0068.001.2185.638
730.0010.504.0048.001.0962.905
850.003.006.004.001.192.425
930.005.508.0063.001.1174.888
1050.008.002.004.001.202.864
1150.008.006.0049.001.1757.655
1270.0010.508.0049.001.1757.427
1350.008.006.0050.001.1857.497
1450.0013.006.0025.001.2123.915
1550.008.006.0048.001.1760.934
1630.005.504.0024.001.1328.436
1750.008.006.0067.001.1288.137
1830.0010.504.0050.001.1760.457
1930.0010.508.0049.001.1660.822
2050.008.006.0048.001.1860.583
2130.0010.508.0047.001.1760.395
2270.0010.508.001.001.180.176
2350.008.006.0049.001.1760.891
2470.005.504.0051.001.1860.328
2550.008.006.0084.001.13101.412
2650.008.006.0045.001.0158.665
2730.005.508.0058.001.1772.301
2850.008.006.0081.001.1588.478
2950.008.006.0028.001.1539.165
3070.0010.504.0029.001.2436.915
3170.0010.504.0049.001.2059.978
3250.008.006.005.001.241.765
3350.008.0010.0045.001.1649.286
3470.005.508.0043.001.0053.785
3570.005.504.0050.001.2159.595
3650.008.006.0051.001.2558.957
Table 4. Analysis of variance results of response value Y1 (PV).
Table 4. Analysis of variance results of response value Y1 (PV).
ResponseF-StatisticR2pAdj.R2APSDPRESS
PV (Y1)675.990.9957<0.00010.9943103.75261.56153.65
Table 5. Analysis of variance results of response value Y2 (density).
Table 5. Analysis of variance results of response value Y2 (density).
ResponseF-StatisticR2pAdj.R2APSDPRESS
Density (Y2)408.090.9930<0.00010.990581.76440.00640.0020
Table 6. Analysis of variance results of response value Y3 (YP).
Table 6. Analysis of variance results of response value Y3 (YP).
ResponseF-StatisticR2pAdj.R2APSDPRESS
YP (Y3)654.370.9956<0.00010.994197.63781.99283.51
Table 7. Numerical and actual optimization results of response surface.
Table 7. Numerical and actual optimization results of response surface.
Factor Y1 (mPa·s)Y2 (g/cm3)Y3 (Pa)
Estimate40.1741.22346.584
Experimental value 1401.2348.545
Experimental value 2401.2350.078
Experimental value 3411.2349.056
Table 8. Results of biodegradability testing and heavy metal content testing of drilling fluid.
Table 8. Results of biodegradability testing and heavy metal content testing of drilling fluid.
ParametersTest Content (mg/L)Biodegradable BOD5/CODStandard
CODCR97580.426>5%
BOD54156
Element NameTest result (mg/kg)Standard value (mg/kg)
Pb0.9500<1000
Cr2.0600<1000
AS1.9100<75
Cd0.0151<20
Hg0.0019<25
Table 9. Cost of microbial drilling fluid.
Table 9. Cost of microbial drilling fluid.
Raw MaterialsPriceCost per Kilogram (CNY)
Bentonite0.5 CNY/kg1 × 10−3
Biological solutionComposition of culture medium12 yuan/L2.4
Carboxymethylcellulose sodium (CMC)20 CNY/kg1.6 × 10−3
Xanthan gum (XG)19 CNY/kg2.28 × 10−2
Barite0.44 CNY/kg1.056 × 10−1
Final cost2.531 CNY/kg
Table 10. Cost of other types of drilling fluids.
Table 10. Cost of other types of drilling fluids.
TypePrice per Kilogram
Lianyungang Chemical Drilling Fluid A7.4 CNY/kg
Shandong Compound Drilling Fluid B4.8 CNY/kg
Shandong Chemical Fiber Drilling Mud C8.8 CNY/kg
Xi’an Chemical Mud D3.92 CNY/kg
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pan, R.; Shu, Z.; Chen, Y.; Sha, X.; Zhang, X.; Han, Y. Microbially Induced Calcite Precipitation (MICP) Improved Drilling Fluid Optimization for Gravel Stratum. Processes 2025, 13, 162. https://doi.org/10.3390/pr13010162

AMA Style

Pan R, Shu Z, Chen Y, Sha X, Zhang X, Han Y. Microbially Induced Calcite Precipitation (MICP) Improved Drilling Fluid Optimization for Gravel Stratum. Processes. 2025; 13(1):162. https://doi.org/10.3390/pr13010162

Chicago/Turabian Style

Pan, Rui, Zhou Shu, Yumin Chen, Xiaobing Sha, Xinquan Zhang, and Yi Han. 2025. "Microbially Induced Calcite Precipitation (MICP) Improved Drilling Fluid Optimization for Gravel Stratum" Processes 13, no. 1: 162. https://doi.org/10.3390/pr13010162

APA Style

Pan, R., Shu, Z., Chen, Y., Sha, X., Zhang, X., & Han, Y. (2025). Microbially Induced Calcite Precipitation (MICP) Improved Drilling Fluid Optimization for Gravel Stratum. Processes, 13(1), 162. https://doi.org/10.3390/pr13010162

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop