趙弘，吳婷婷，徐泉，趙毅鑫

1 中國(guó)石油大學(xué)(北京)機(jī)械與儲(chǔ)運(yùn)工程學(xué)院，北京 102249

2 中國(guó)石油大學(xué)(北京)重質(zhì)油國(guó)家重點(diǎn)實(shí)驗(yàn)室/新能源研究院，北京 102249

3 中國(guó)礦業(yè)大學(xué)(北京)資源與安全工程學(xué)院，北京 100083

*通信作者, ttabbmmjj@163.com

石油機(jī)械

3D打印類巖石材料的實(shí)驗(yàn)優(yōu)化研究

趙弘1*，吳婷婷1，徐泉2，趙毅鑫3

1 中國(guó)石油大學(xué)(北京)機(jī)械與儲(chǔ)運(yùn)工程學(xué)院，北京 102249

2 中國(guó)石油大學(xué)(北京)重質(zhì)油國(guó)家重點(diǎn)實(shí)驗(yàn)室/新能源研究院，北京 102249

3 中國(guó)礦業(yè)大學(xué)(北京)資源與安全工程學(xué)院，北京 100083

*通信作者, ttabbmmjj@163.com

3D打印被譽(yù)為工業(yè)制造2025的一項(xiàng)關(guān)鍵技術(shù)，已在高強(qiáng)度鋼材、高溫合金、塑料橡膠上取得不同程度的工業(yè)化應(yīng)用，也是當(dāng)前發(fā)展熱點(diǎn)。然而，3D打印在石油巖芯領(lǐng)域仍然存在較多技術(shù)難點(diǎn)，例如巖芯的高獲得成本及較低的重復(fù)一致性等均是制約其進(jìn)一步發(fā)展的關(guān)鍵因素。本文首次嘗試從3D打印材料組合出發(fā)，通過(guò)類材料的方法展開(kāi)研究。首先通過(guò)對(duì)幾種常用類巖石類材料及頁(yè)巖粉的性質(zhì)分析，成功選擇出了一種適用于3D打印的巖芯材料組合。結(jié)合響應(yīng)面法進(jìn)一步分析了不同配比對(duì)類巖石類材料滲透率、孔隙度和凝結(jié)時(shí)間的影響。最后結(jié)合中心組合設(shè)計(jì)(CCD)對(duì)模型進(jìn)行優(yōu)化，成功得到最佳的配方配比方案，并與真實(shí)數(shù)據(jù)進(jìn)行比較，結(jié)果顯示擬合效果較好。

類巖石類材料；3D打印技術(shù)；響應(yīng)面法；中心組合設(shè)計(jì)；人造巖芯

0 引言

3D打印(3D Printing，又稱三維打印)作為第三次工業(yè)革命的代表性技術(shù)之一，越來(lái)越多地受到人們關(guān)注[1]，已廣泛應(yīng)用于材料、機(jī)械、電子、軟件等領(lǐng)域[2]。在航空航天方面，美國(guó)Aero Met公司制備的F/A-18E/F戰(zhàn)機(jī)翼根吊環(huán)，其疲勞壽命是要求值的4倍，靜力加載到225%仍未被破壞[3]。在生物植入件方面，利用Ni、Ti金屬制備的載重植入體，其多孔和功能梯度結(jié)構(gòu)與人體具有良好的相容性，使用壽命可高達(dá)7年～12年[4-7]。在醫(yī)療方面，3D打印的個(gè)體化支架、器官或腫瘤模型等，其生物力學(xué)性能與正常狀態(tài)非常接近[8-11]。

3D打印技術(shù)能否有更廣泛的應(yīng)用領(lǐng)域取決于打印材料的發(fā)展是否能適應(yīng)打印技術(shù)的要求[12]。目前廣泛應(yīng)用的3D打印材料有光敏樹(shù)脂、工程塑料、金屬材料、橡膠類材料、陶瓷類材料等，除此之外，還有少量的人造骨粉材料、彩色石膏材料、細(xì)胞生物材料以及食品類材料等[13]。

如何準(zhǔn)確地表征與直觀地顯示固體內(nèi)部的復(fù)雜結(jié)構(gòu)是解決諸多工程實(shí)際問(wèn)題的基礎(chǔ)和關(guān)鍵。在石油、天然氣和礦業(yè)工程領(lǐng)域中，能否高效安全地開(kāi)采儲(chǔ)層中的石油、天然氣、煤層氣和礦物資源依賴于對(duì)地下儲(chǔ)層結(jié)構(gòu)及地應(yīng)力場(chǎng)的準(zhǔn)確認(rèn)知與科學(xué)描述[14-18]。

在石油科學(xué)技術(shù)研究中，巖芯的勘探與研究作為一個(gè)重要的課題，一直以來(lái)備受關(guān)注。巖體介質(zhì)深埋于地下，其內(nèi)部有復(fù)雜的孔隙、裂隙結(jié)構(gòu)、由地應(yīng)力及開(kāi)采應(yīng)力引發(fā)的結(jié)構(gòu)演化與形變、流體滲流、應(yīng)力-滲流相互作用等物理力學(xué)過(guò)程，利用現(xiàn)有的方法和技術(shù)條件，難以準(zhǔn)確地獲知，定量表征難度更大[19]。傳統(tǒng)的巖石結(jié)構(gòu)面原位取樣或原位剪切試驗(yàn)有制作過(guò)程中結(jié)構(gòu)面表面形態(tài)容易破壞、上下盤(pán)吻合度差、取樣周期長(zhǎng)、成本高、不能進(jìn)行大量的結(jié)構(gòu)面直剪實(shí)驗(yàn)等缺點(diǎn)[20]，導(dǎo)致巖芯探測(cè)難度大、成本高、可靠性低等，亟待發(fā)展人造巖芯制造技術(shù)。研究發(fā)現(xiàn)，影響人造巖芯性能的物性參數(shù)主要有滲透率、孔隙度和孔道半徑等。

常規(guī)的優(yōu)化方法是通過(guò)改變單一變量，逐個(gè)評(píng)估確定材料配比，不僅實(shí)驗(yàn)工作量大，而且結(jié)果可靠性還差。響應(yīng)面法是數(shù)學(xué)方法和統(tǒng)計(jì)方法相結(jié)合的產(chǎn)物，以較少的實(shí)驗(yàn)次數(shù)對(duì)多個(gè)變量進(jìn)行建模分析，優(yōu)化相應(yīng)模型。中心組合設(shè)計(jì)(Central Composite Design,CCD)是國(guó)內(nèi)外近年來(lái)興起的一種基于統(tǒng)計(jì)技術(shù)的優(yōu)化方法，具有實(shí)驗(yàn)次數(shù)少，精度高，預(yù)測(cè)性好等優(yōu)點(diǎn)，已成功應(yīng)用于醫(yī)藥學(xué)等領(lǐng)域[21-23]。

本文根據(jù)3D打印技術(shù)對(duì)打印材料的要求[24]，對(duì)現(xiàn)有的類巖石材料進(jìn)行分析研究，從中選取出適合3D打印的類巖石材料，并基于不同配比對(duì)滲透率、孔隙度及凝結(jié)時(shí)間的影響進(jìn)行優(yōu)化。

1 實(shí)驗(yàn)設(shè)計(jì)

中心組合試驗(yàn)設(shè)計(jì)(CCD)，又稱二次回歸旋轉(zhuǎn)設(shè)計(jì)，應(yīng)用各影響因子的二次多項(xiàng)式來(lái)預(yù)測(cè)其對(duì)評(píng)價(jià)指標(biāo)的作用[25]，將傳統(tǒng)的插值方式與全因子或部分因子設(shè)計(jì)相結(jié)合，在全因子或部分因子設(shè)計(jì)的試驗(yàn)點(diǎn)之外，增加位于設(shè)計(jì)空間中心的試驗(yàn)點(diǎn)。在沿著通過(guò)試驗(yàn)中心且平行于各個(gè)坐標(biāo)軸的直線上，選取距離中心±δ距離的設(shè)計(jì)點(diǎn)[26]，分析他們對(duì)設(shè)計(jì)結(jié)果的影響，從而得出非線性的關(guān)系式(1)，用來(lái)表示輸入變量x與輸出變量y的關(guān)系。

其中n為試驗(yàn)次數(shù)，α0為常數(shù)項(xiàng)待定系數(shù)，αi為一次項(xiàng)待定系數(shù)，αij為二次項(xiàng)待定系數(shù)，F(xiàn)為全因子或部分因子設(shè)計(jì)的試驗(yàn)點(diǎn)數(shù)。若設(shè)計(jì)變量個(gè)數(shù)為n，每個(gè)設(shè)計(jì)變量分量取l，則試驗(yàn)點(diǎn)的個(gè)數(shù)為 ln。以n=3，l=3為例，中心組和試驗(yàn)點(diǎn)的分布形式如圖1所示。

1.1 實(shí)驗(yàn)材料

圖1 中心組合設(shè)計(jì)Fig. 1 Central Composite Design

3D打印技術(shù)對(duì)材料一般有如下要求：(1)有利于快速、精確地加工原型零件；(2)快速成型制件接近最終要求，盡量滿足對(duì)強(qiáng)度、剛度、耐潮濕性、熱穩(wěn)定性能等要求；(3)有利于后續(xù)的處理工藝[27-29]。

石膏以其來(lái)源廣、成本低等特性被廣泛應(yīng)用于3D打印工業(yè)生產(chǎn)中[30-32]，與其他材料相比石膏具有以下優(yōu)勢(shì)：

(1)表面光滑飽滿，粉末顆粒易于調(diào)整，打印質(zhì)量可以保證；

(2)顏色潔白，且便于上色，可實(shí)現(xiàn)彩色打?。?/p>

(3)質(zhì)地細(xì)膩，具有良好的裝飾性和加工性；

(4)制造成本低，加工的性價(jià)比高；

(5)安全環(huán)保，無(wú)毒害性。

普通水泥通常初凝時(shí)間較長(zhǎng)且早期的強(qiáng)度略低，但是速凝水泥不僅具有良好的耐腐蝕性、耐熱性、后期強(qiáng)度發(fā)展高且二次水化反應(yīng)慢等特點(diǎn)，還提高了初凝時(shí)間及早期的強(qiáng)度。

因此，本文選擇石膏粉和水泥作為打印材料。所選石膏粉和水泥的基本特征及適用范圍如表1和表2所示。

1.2 不同配比對(duì)滲透率及孔隙度的影響

1.2.1 實(shí)驗(yàn)設(shè)計(jì)及操作

結(jié)合響應(yīng)面法，以石膏、水泥的質(zhì)量與水的體積為響應(yīng)變量，分別以A、B、C表示，巖石的滲透率、孔隙度為分析指標(biāo)，以R1、R2表示。當(dāng)石膏的質(zhì)量范圍為5 g～15 g，水泥質(zhì)量5 g～15 g，水的體積6 mL～12 mL時(shí)，試驗(yàn)分析水平及因素如表3所示。得到9組不同的配比方案，如表4所示。

實(shí)驗(yàn)操作如下：

(1)取樣

按照表4所確定的配比，用電子秤稱取各種材料。

(2)拌勻

將石膏、水泥依次倒入攪拌器中，然后邊攪拌邊倒入水，直到攪拌均勻。

(3)填料

將攪拌均勻的材料填入3D打印機(jī)的擠出裝置中，邊放入邊攪拌，防止類巖石材料凝結(jié)。

(4)打印

將裝好料的擠出裝置固定在3D打印機(jī)上，將材料從擠出裝置中擠壓出來(lái)，并在底盤(pán)上進(jìn)行加熱處理，加快類巖石材料的凝結(jié)速度，使之快速成型。

(5)清洗器具

打印完成后，迅速對(duì)實(shí)驗(yàn)所用器件進(jìn)行清理、養(yǎng)護(hù)，防止類巖石材料凝結(jié)在器具上，造成不必要的損壞。

1.2.2 實(shí)驗(yàn)結(jié)果及分析

試驗(yàn)所需檢測(cè)條件如表5所示，得到相應(yīng)滲透率及孔隙度值如表6所示，經(jīng)高壓滲透實(shí)驗(yàn)后的樣品如圖2所示。

實(shí)驗(yàn)中選取的石膏粉的密度為2.8 g/cm3～ 3.0 g/cm3，水泥的密度為3.0 g/cm3～3.15 g/cm3。為保證混合物表面的平整度，取石膏的質(zhì)量比例大于水泥的質(zhì)量比例。分別將水泥和石膏單獨(dú)與水混合進(jìn)行實(shí)驗(yàn)發(fā)現(xiàn)，水灰比為0.45時(shí)，混合物的收縮值最小[33]；膏水比為1:1.4時(shí)，混合物既有較好的流動(dòng)性，又能保證一定的耐壓強(qiáng)度[34]。根據(jù)這一結(jié)果，確定石膏、水泥、水按1.4:1:1.45的比例混合，水的質(zhì)量比例近似為35%，當(dāng)水的質(zhì)量比例小于35%時(shí)，無(wú)法保證混合物的流動(dòng)性。

用響應(yīng)面分析法進(jìn)行優(yōu)化分析[35]，由于孔隙度、滲透率系統(tǒng)存在彎曲，所以選取二階模型，更加接近于實(shí)際系統(tǒng)。分別以A、B、C代表石膏、水泥、水的質(zhì)量，以R1、R2表示孔隙度及滲透率，建立擬合方程：滲透率及孔隙度的實(shí)驗(yàn)值與擬合值對(duì)比情況如圖3所示，擬合方程中滲透率、孔隙度的方差分析如表7和表8所示。由圖表分析可得，擬合曲線可近似等價(jià)于實(shí)際系統(tǒng)的方程。

表1 實(shí)驗(yàn)材料-石膏粉Table 1 The experimental materials—gypsum powder

表2 實(shí)驗(yàn)材料-水泥Table 2 The experimental materials—cement

表3 實(shí)驗(yàn)水平及因素Table 3 Test levels and factors

表4 配比方案Table 4 The solution composition

表5 檢測(cè)條件Table 5 Measure conditions

表6 正交表Table 6 Orthogonal array

圖2 試樣表面視圖Fig. 2 The sample surface view

圖3 實(shí)驗(yàn)值與擬合值對(duì)比Fig. 3 Test value and fi tting value contrast

表7 滲透率的方差分析Table 7 Permeability analysis of variance

表8 孔隙度的方差分析Table 8 Analysis of variance of porosity

由表7和表8可知，對(duì)試樣的滲透率影響：石膏＞水＞水泥；對(duì)孔隙度的影響：水＞石膏＞水泥；綜合分析對(duì)滲透率與孔隙度的影響：石膏＞水＞水泥。類巖石類材料的滲透率、孔隙度與石膏、水泥、水的含量關(guān)系如圖4和圖5所示。該結(jié)論與試樣表面粗糙度分析相一致，證明擬合合理。

2 凝結(jié)時(shí)間優(yōu)化實(shí)驗(yàn)

2.1 實(shí)驗(yàn)操作及結(jié)果

選取表4所示9組配比方案進(jìn)行打印實(shí)驗(yàn)，得到各配比材料的初凝時(shí)間和終凝時(shí)間，如表9所示。

圖4 滲透率與石膏、水泥、水含量關(guān)系Fig. 4 The relationship between permeability and composition content of gypsum, cement and water

圖5 孔隙度與石膏、水泥、水含量關(guān)系Fig. 5 The relationship between porosity and composition content of gypsum, cement and water

表9 凝結(jié)時(shí)間Table 9 Setting time

分別以A、B、C代表石膏、水泥和水的質(zhì)量，用響應(yīng)面分析法[35]對(duì)實(shí)驗(yàn)值進(jìn)行二次函數(shù)的回歸擬合，建立初凝時(shí)間Rt1、終凝時(shí)間Rt2的回歸模型，其表達(dá)式分別為：

2.2 實(shí)驗(yàn)結(jié)果分析

2.2.1 材料與初凝時(shí)間關(guān)系分析

(1)將凝結(jié)時(shí)間的實(shí)驗(yàn)值與擬合值進(jìn)行比較，如圖6所示，雖然存在差別，但大致圍繞在預(yù)測(cè)值左右，可以認(rèn)為擬合方程合理。

圖6 凝結(jié)時(shí)間的實(shí)驗(yàn)值與擬合值Fig. 6 Setting time of the experimental value and fi tting value

(2)如圖7所示，以材料A(石膏)輸出值的學(xué)生化殘差為橫坐標(biāo)[36]，以其百分比為縱坐標(biāo)，其輸出值(初凝時(shí)間)圍繞在直線上下無(wú)規(guī)律分布，可近似看成是一條直線。輸出值在各個(gè)區(qū)間內(nèi)的布情況如表10所示，由此可見(jiàn)，雖然擬合存在一定的誤差，但材料A(石膏)的凝結(jié)時(shí)間可近似看作正態(tài)分布曲線。

(3)當(dāng)以凝結(jié)時(shí)間的預(yù)測(cè)值為橫坐標(biāo)，學(xué)生化殘差為縱坐標(biāo)時(shí)，所描繪的點(diǎn)圍繞殘差值等于0的直線上下隨機(jī)散布，且都分布在(-2，2)之間，沒(méi)有任何系統(tǒng)趨勢(shì)。認(rèn)為回歸直線與原觀測(cè)值的擬合情況良好，材料A(石膏)的質(zhì)量與初凝時(shí)間有較為明顯的線性關(guān)系。初凝時(shí)間預(yù)測(cè)值與學(xué)生化殘差圖如圖8所示。

(4)以材料A(石膏)初凝時(shí)間的測(cè)定值為橫坐標(biāo)，預(yù)測(cè)值為縱坐標(biāo)做圖。剔除無(wú)效點(diǎn)后，預(yù)測(cè)值在測(cè)定值的上下附近分布，可近似擬合為一條直線，凝結(jié)時(shí)間在2 s～40 s之間，如圖9所示。

(5)以實(shí)驗(yàn)次數(shù)為橫坐標(biāo)，學(xué)生化殘差為縱坐標(biāo)，剔除不合理的點(diǎn)之后，發(fā)現(xiàn)殘差值無(wú)規(guī)律的分布在(-2，2)之間，由此可以認(rèn)為預(yù)測(cè)值的擬合方程基本符合實(shí)際系統(tǒng)，學(xué)生化殘差的運(yùn)行結(jié)果如圖10所示。

綜上所述，材料A(石膏)初凝時(shí)間的實(shí)驗(yàn)值與預(yù)測(cè)值在合理誤差范圍內(nèi)，其質(zhì)量比例與初凝時(shí)間基本符合線性相關(guān)關(guān)系。水泥實(shí)驗(yàn)分析得到與圖6～圖10相近的結(jié)果，說(shuō)明，水泥的質(zhì)量、水的體積與初凝時(shí)間也符合線性相關(guān)關(guān)系。

2.2.2 材料之間相互作用關(guān)系

(1)當(dāng)材料A(石膏)在區(qū)間范圍內(nèi)分別取5 g、10 g、15 g時(shí)，材料B(水泥)與材料C(水)之間的相互影響如圖11所示。

圖7 材料A(石膏)初凝時(shí)間殘差與百分比圖Fig. 7 The residual (gypsum) initial setting time vs percentage fi gure of material A (gypsum)

表10 頻率檢驗(yàn)Table 10 Frequency check

圖8 材料A(石膏)初凝時(shí)間的殘差與預(yù)測(cè)值Fig. 8 The predictive values vs residuals of initial setting time of material A (gypsum)

圖9 材料A(石膏)初凝時(shí)間預(yù)測(cè)值與測(cè)定值圖Fig. 9 The measured value and prediction of initial setting time of material A (gypsum)

圖10 材料A(石膏)殘差Fig. 10 The residuals of material A (gypsum)

由圖11可知，當(dāng)材料A(石膏)為5 g時(shí)，在水泥和水的取值范圍內(nèi)，等高線基本成水平分布，此時(shí)凝結(jié)時(shí)間隨水泥質(zhì)量變化不顯著，但隨著水的體積增加，水與水泥的交互作用顯著增加。當(dāng)石膏取10 g，等高線的切線角度增大，且出現(xiàn)圓心，即石膏取10 g左右時(shí)凝結(jié)時(shí)間達(dá)到最佳狀態(tài)。當(dāng)石膏取15 g時(shí)，水泥與水的交互作用下降，凝結(jié)時(shí)間基本隨水泥質(zhì)量的增加而增大。

綜上所述，當(dāng)材料A(石膏)的質(zhì)量取值一定時(shí)，材料B(水泥)的質(zhì)量與材料C(水)的體積呈線性相關(guān)關(guān)系，且這種交互作用在石膏質(zhì)量的取值范圍內(nèi)先減小后增大。隨著材料A(石膏)質(zhì)量的改變，混合物的初凝時(shí)間也隨之變化，且變化趨勢(shì)呈二次分布，變化率先減小后增大，當(dāng)石膏質(zhì)量取10 g左右時(shí)，變化率最小，達(dá)到最佳凝結(jié)時(shí)間。

(2)當(dāng)材料B(水泥)在區(qū)間范圍內(nèi)取5 g、10 g、15 g時(shí)，材料A(石膏)與材料C(水)之間的相互影響如圖12所示。

圖11 水泥與水等高線圖(a)石膏5 g；(b)石膏10 g；(c)石膏15 gFig. 11 Cement and water contour map(a) gypsum 5 g; (b) gypsum 10 g; (c) gypsum 15 g

圖12 石膏與水等高線圖(a)水泥5 g; (b)水泥10 g; (c)水泥15 gFig. 12 Gypsum contour map with water(a) cement 5 g; (b) cement 10 g; (c) cement 15 g

結(jié)果發(fā)現(xiàn)，當(dāng)材料B(水泥)的值一定時(shí)，材料A(石膏)與材料C(水)的質(zhì)量呈線性相關(guān)關(guān)系，且相關(guān)性大于水泥與水，即石膏的水溶性更好。

隨著材料B(水泥)質(zhì)量的改變，混合物的初凝時(shí)間也隨之變化，且變化趨勢(shì)呈二次分布，當(dāng)材料B(水泥)的質(zhì)量在5 g左右，石膏與水的交互作用最小，初凝時(shí)間達(dá)到最佳。

(3)當(dāng)材料C(水)在區(qū)間范圍內(nèi)取不同值時(shí)，材料A(石膏)與材料B(水泥)之間的相互影響，初凝時(shí)間的三維響應(yīng)曲線如圖13所示，縱坐標(biāo)為初凝時(shí)間。

從圖中可以看出，當(dāng)水取6 mL或12 mL時(shí)，石膏和水泥的質(zhì)量變化呈二次函數(shù)關(guān)系，交互作用較為顯著，當(dāng)水為9 mL時(shí)，石膏和水泥的交互作用明顯降低，且在材料A(石膏)為9 g，材料B(水泥)為5 g時(shí)出現(xiàn)最低點(diǎn)，即凝結(jié)時(shí)間最快。因此在混合物的配比中，加入的水量較少或較多時(shí)，將出現(xiàn)難以混合均勻或凝結(jié)時(shí)間較長(zhǎng)等問(wèn)題，這將影響打印件的表面粗糙度。

當(dāng)水含量小于20%時(shí)，混合物難以均勻攪拌，凝結(jié)產(chǎn)物的硬度和表面粗糙度都難以得到保障。水含量大于50%時(shí)，將大大增加混合物凝結(jié)時(shí)間，無(wú)法在打印過(guò)程完成材料的迅速凝結(jié)，無(wú)法達(dá)到打印材料的硬度標(biāo)準(zhǔn)，甚至引起坍塌等現(xiàn)象的發(fā)生。

綜上所述，根據(jù)中心組合設(shè)計(jì)，利用二次函數(shù)的擬合方程[37-39]，得到21種優(yōu)化方案，如表11所示。其中A為石膏質(zhì)量，單位是g，B為水泥質(zhì)量，單位是g，C為水的體積，單位mL，R1、R2分別為初凝時(shí)間和終凝時(shí)間，單位為s。

在噴頭直徑為2 mm，出料速度為30 mm/s±10 mm/s的條件下，為滿足x、y軸方向收縮量為0.1 mm～0.2 mm，z軸方向的收縮量不大于0.2 mm的打印要求[40]，材料的初凝時(shí)間應(yīng)在1 s～3 s之間，終凝時(shí)間在3 s～5 s之間。這樣既不至于堵塞打印噴頭，又能保證打印件不坍塌堆疊。據(jù)此選擇初凝時(shí)間較短，終凝時(shí)間合理，且期望值也較小的6組優(yōu)化設(shè)計(jì)進(jìn)行進(jìn)一步的實(shí)驗(yàn)驗(yàn)證。根據(jù)這6組的實(shí)驗(yàn)配比，在同等條件下重新實(shí)驗(yàn)，驗(yàn)證結(jié)果如表12所示，初凝時(shí)間、終凝時(shí)間的預(yù)測(cè)值與實(shí)驗(yàn)值的擬合柱形圖如圖14和圖15所示。

從表12及圖14和圖15可以看出，這6組優(yōu)化參數(shù)所對(duì)應(yīng)的初凝時(shí)間和終凝時(shí)間與實(shí)驗(yàn)結(jié)果非常接近，相對(duì)誤差都在4%左右，說(shuō)明所建立的模型精確、可靠。其中第6組實(shí)驗(yàn)(石膏9.35 g、水泥5 g、水6.12 mL)初凝時(shí)間最短，期望最小，雖然終凝時(shí)間較前5項(xiàng)長(zhǎng)一點(diǎn)，但對(duì)打印效果并無(wú)明顯影響，因此最終選擇由該比例配制而成的類巖石材料為3D打印機(jī)的打印材料。

圖13 石膏與水泥三維響應(yīng)曲面圖(a)水6 mL；(b)水9 mL；(c)水12 mLFig. 13 Gypsum and cement 3D response surface fi gure(a) water 6 mL;(b) water 9 mL;(c) water 12 mL

以石膏10 g、水泥10 g、水6 mL作為空白組與最佳配比組做比較?？瞻捉M的拉壓比為1.27，最佳配比組的拉壓比為1.52?？梢?jiàn)優(yōu)化后的材料凝結(jié)時(shí)間更短，抗壓強(qiáng)度更高，拉壓比更大，如圖16和圖17所示。

表11 優(yōu)化方案列表Table 11 Optimization scheme list

表12 驗(yàn)證實(shí)驗(yàn)條件及結(jié)果Table 12 Validation test conditions and results

圖14 初凝時(shí)間擬合情況柱形圖Fig. 14 Initial setting time fi tting bar charts

圖15 終凝時(shí)間擬合柱形圖Fig15 The fi nal setting time fi tting bar chart

圖16 最佳配比組與空白組凝結(jié)時(shí)間對(duì)比Fig. 16 The optimal proportion compared with blank group setting time

圖17 最佳配比組與空白組抗壓強(qiáng)度對(duì)比Fig17 The optimal proportion compared with blank group compressive strength

3 結(jié)論

通過(guò)以上實(shí)驗(yàn)分析可得結(jié)論如下：

(1)適應(yīng)3D打印的類巖石類材料為高強(qiáng)硬度特白石膏粉-B和FDB-401速凝型水泥。

(2)隨著石膏質(zhì)量增加、水泥質(zhì)量減少、水質(zhì)量減少，滲透率呈上升趨勢(shì)；隨著石膏質(zhì)量增加、水泥質(zhì)量增加、水質(zhì)量增加孔隙度呈上升趨勢(shì)。

(3)各因素中石膏對(duì)滲透率和孔隙度的影響最大，水泥的影響最小。

(4)當(dāng)石膏、水泥、水按照9.35 g、5 g、6.12 mL配比混合時(shí)，材料凝結(jié)時(shí)間短，抗壓強(qiáng)度高，打印效果最符合實(shí)驗(yàn)要求。

[1] 王忠宏, 李揚(yáng)帆, 張曼茵. 中國(guó)3D打印產(chǎn)業(yè)的現(xiàn)狀及發(fā)展思路[J]. 經(jīng)濟(jì)縱橫, 2013, (1): 90-93. [WANG Z H, LI Y F, ZHANG M Y,et al. China’s present situation and development of 3D printing industry[J]. Economic Review, 2013, (1): 90-93.]

[2] WALTERS P, DAVIES K. 3D printing for artists: Research and creative practice[J]. Journal of the Norwegian Print Association, 2010,(1):12-15.

[3] VANDENBROUCKE B. Selective laser melting of biocompatible metals for rapid manufacturing of medical parts[J]. Rapid Prototyping Journal, 2007, 13: 196-203.

[4] VAMSI KRISHNA B, XUE W C, BOSE S, et al. Engineered porous metals for implants[J]. JOM, 2008 (5) : 45-48．

[5] VAMSI KRISHNA B, BOSE S, BANDYOPADHYAY A. Low stiffness porous Ti structures for load-bearing implants[J]. Acta Biomaterialia, 2007, (3): 997-1 006.

[6] ZHANF B. Preliminary study on some properties of Co-Cr dental alloy formed by selective laser melting technique[J]. Journal of Wuhan University of Technology (Mater Sci Ed), 2012, 27: 665-668.

[7] 費(fèi)群星, 張雁, 譚永生, 等．激光近凈成形Ni-Cu-Sn合金[J]．稀有金屬材料與工程, 2007, (36): 2 052-2 056．[FEI Q X, ZHANG Y, TAN Y S, et al. Laser near net forming of Ni-Cu-Sn alloy [J]. Rare Metal Materials and Engineering, 2007, (36): 2 052-2 056.]

[8] FISCHBACH C, CHEN R, MASTSUMOTO T, et al. Engineering tumors with 3D scaffolds [J]. Nat Methods, 2007, 4(10); 855-860.

[9] DING C, QIAO Z, JIANG W, et al. Regeneration of a goat femoral head using a tissue-speci fi c, biphasic scaffolds fabricated with CAD/CAM technology [J]. Biomaterials, 2013, 34(28): 6 706-6 716.

[10] CIOCCA L, DONATI D, FANTINI M, et al. CAD-CAM-generated hydroxyapatite scaffolds to replace the mandibular condyle in sheep:Preliminary results [J]. J Biometer Apply, 2013, 28(2): 207-218.

[11] LEE J Y, CHOL B, WU B, et al. Customized biomimetic scaffolds created by indirect three-dimensional printing for tissue engineering[J]. Bio fabrication, 2013 ,5(4): 045 003-045 011.

[12] TIM O. 7 things you should know about 3D printing [EB/OL]. [2013-05-22]. http://net. educause. Edu/ir/library/pdf/ELI7086. pdf.

[13] 李滌塵, 田小永, 王永信, 等．增材制造技術(shù)的發(fā)展[J]. 電加工與模具, 2012, (S1): 20-22. [LI D C, TIAN X Y, WANG Y X, et al.Add material manufacturing technology development [J]. Electro Machining & Mold, 2012, (S1): 20-22.]

[14] 周新桂, 操成杰, 袁嘉音. 儲(chǔ)層構(gòu)造裂縫定量預(yù)測(cè)與油氣滲流規(guī)律研究現(xiàn)狀和進(jìn)展[J] .地球科學(xué)進(jìn)展, 2003, 18: 398–404. [ZHOU X G, CAO C J, YUAN J Y. Quantitative forecasting reservoir structural fractures and oil and gas seepage law research status and progress[J]. Advance in Earth Sciences, 2003, 18: 398–404.]

[15] 侯泉林, 李會(huì)軍, 范俊佳, 等. 構(gòu)造煤結(jié)構(gòu)與煤層氣賦存研究進(jìn)展[J]. 中國(guó)科學(xué):地球科學(xué), 2012,42: 1 487–1 495. [HOU Q L, LI H J, FAN J, et al. Tectonic coal structure and occurrence of coal bed methane (CBM) is reviewed[J]. Scientia Sonica(Terrae), 2012, 42:1 487–1 495.]

[16] ZHANG L K, WANG Z L, QU Z H, et al. Physical simulation experiment of gas migration in sandstone porous media[C]. Acta Geol Sin, 2007, 81: 539–544.

[17] MAZUMDER S, WOLF K. Differential swelling and permeability change of coal in response to CO2injection for ECBM[C]. Int J Coal Geol, 2008,74:123–138.

[18] SONG H Z, ZENG H R, SUN J X, et al. Methods of server tectonic fracture prediction and its application[C]. Seismol Geol, 1999, 21:205–213.

[19] 鞠楊, 謝和平, 鄭澤民, 等. 基于3D打印技術(shù)的巖體復(fù)雜結(jié)構(gòu)與應(yīng)力場(chǎng)的可視化方法[J]. 科學(xué)通報(bào), 2014, 59: 3 109–3 119. [JU Y,XIE H P, ZHENG Z M, et al. Based on 3D printing visualization method of complex structure and stress fi eld of rock mass[J]. Chinese Science Bulletin, 2014, 59: 3 109–3 119.]

[20] 羅戰(zhàn)友, 杜時(shí)貴, 黃曼, 等. 一種直剪試驗(yàn)結(jié)構(gòu)面模型的制作方法:中國(guó), CN200910154781.9 [P]. 2010-06-02. [LUO Z Y, DU S G,HUANG M, et al. A kind of direct shear test method of structural plane model: China CN200910154781.9 [P]. 2010-06-02.]

[21] ISMAIL A, SOUTANE S, GHOUL M. Optimization of the enzymatic synthesis of butyl glucoside using response surface methodology[J].Bioethanol Prog, 1998, 14(6): 874-878.

[22] OOIJKAASL P, WILKINSON E C, TRAMPER J. Medium optimization for spore production of coniothyrium minivans using statistically based experimental designs [J]. Biotechnology and Bioengineering, 1999, 64(1): 92–100.

[23] 張志榮, 王建新, 呂劍, 等. 中心組合設(shè)計(jì)法優(yōu)化3′,5′-二辛?；?氟苷藥質(zhì)體制備工藝[J]. 藥學(xué)學(xué)報(bào), 2001, 36(6): 456-461.[ZHANG Z R, WANG J X, LV J, et al. Central composite design method to optimize 3, 5-dactyl acyl fl uorine nucleoside drugs plasmid preparation technology[J]. Acta Pharmaceutical Sonica, 2001, 36(6): 456-461.].

[24] 丘北上. “第三次工業(yè)革命”的力量[J]. 中國(guó)經(jīng)濟(jì)信息, 2012, (2): 18-19. [QIU B S. The power of the third industrial revolution[J].China Economic Information, 2012, (2): 18-19.]

[25] 翁德凱, 程寓, 夏玲玲. 基于結(jié)合面的立式加工中心進(jìn)給系統(tǒng)的動(dòng)態(tài)特性分析[J]. 機(jī)械設(shè)計(jì)與制造, 2012(3): 130-132. [WENG D K, CHENG Y, XIA L L. The dynamic characteristics analysis of feeding system of a vertical machining center based on conjoint interfaces[J]. Machinery Design & Manufacture, 2012(3): 130-132.]

[26] 馬寶勝. 響應(yīng)面方法在多種實(shí)際優(yōu)化問(wèn)題中的應(yīng)用[D]. 北京: 北京工業(yè)大學(xué), 2007. [MA B S. Application of response surface method to several practical optimization problems[D]. Beijing: Beijing University of Technology, 2007.]

[27] BERDT J F, ANDERSON T C. Three-dimensional printing materials system: US Patent, No. 6416850[P]. 2002.

[28] SACHS E M, CYMA M J, CARADONNA M A, et al. Jetting layers of powder the formation of fi ne powder beds thereby: US Patent,No.6596224[P].2003.

[29] SHEN J. Material system for use in three-dimensional printing: US Patent, No.7049363[P].2006.

[30] 劉力源,類巖石材料的配制方法與加工制作. [EB/OL]. http://www.doc88.com/p-991529215414.html. [LIU L Y. Class preparation method of rock material and processing production. [EB/OL]. http://www.doc88.com/p-991529215414.html.]

[31] 張永莉, 霍書(shū)娟, 高建平, 等. 硫酸鈣復(fù)合生物材料[J]. 高分子材料科學(xué)與工程, 2006,(1):57-58. [ZHANG Y L, HUO S J, GAO J P,et al. Calcium sulfate composite biomaterials[J]. Polymer Materials Science and Engineering, 2006, (1): 57-58.]

[32] 陳潔. 3D打印一個(gè)時(shí)代的開(kāi)啟[J]. 高科技與產(chǎn)業(yè)化, 2013,(4):37-38. [CHEN J. 3D printing the opening of an era[J]. High-Technology & Industrialization, 2013, (4): 37-38.]

[33] 張巨松, 楊奇, 李宗陽(yáng). 水灰比/礦物外加劑對(duì)水泥的收縮變形實(shí)驗(yàn)[J]. 沈陽(yáng)建筑大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 28(03): 521-526. [ZHANG J S, YANG Q, LI Z Y. Water /cement ratio mineral admixture on cement shrinkage deformation experiment[J]. Journal of Shenyang Jianzhu University (Natural Science), 2012, 28(03): 521-526].

[34] 張萬(wàn)勝. 石膏模具生產(chǎn)中的質(zhì)量控制[J]. 江蘇陶瓷, 2004, 37(01): 32-37. [ZHANG W S. Quality control in gypsum mold production[J]. Jiangsu Ceramics, 2004, 37(01): 32-37.]

[35] 隋允康, 宇慧平. 響應(yīng)面方法的改進(jìn)及其對(duì)工程優(yōu)化的應(yīng)用[M]. 北京: 科學(xué)出版社, 2011. [SUI Y K, YU H P. Improvement of response surface method and its application to engineering optimization[M]. Beijing: Science Press, 2011]

[36] 趙開(kāi)斌. 關(guān)于學(xué)生化殘差的一點(diǎn)注記[J]. 淮北煤師院學(xué)報(bào), 2002,(4):10-13. [ZHAO K B. A little note about students the residuals[J].Journal of Huaibei coal teacher’s college, 2002, (4): 10-13.]

[37] KAMAS I, MCMAHON C A. A response surface method based on weighted regression for structural reliability analysis[J]. Probabilistic Engineering Mechanics, 2005, 20(1): 11-17.

[38] GAVIN H P, YAU S C. High-order limit state functions in the response surface method for structural reliability analysis[J]. Structural Safety, 2008, 30(2): 162-179.

[39] 朱靜, 郭軍, 陸鑫森. 一種新的結(jié)構(gòu)可靠性計(jì)算方法—響應(yīng)面法[J]. 上海交通大學(xué)學(xué)報(bào), 1995, 29(2): 26-31. [ZHU J, GUO J, LU X S. A new method of structural reliability calculation of the response surface method[J]. Journal of Shanghai Jiao tong University, 1995,29(2): 26-31.]

[40] 王康, 黃筱調(diào), 陳勇, 等. 噴墨式水泥3D打印噴頭泥漿流動(dòng)規(guī)律分析[J]. 機(jī)械設(shè)計(jì)與制造工程, 2015, 45(05): 28-31. [WANG K,HUANG X T, CHENG Y, et al. Analysis of mud fl ow pattern of inkjet cement 3D printing sprinkler [J]. Machine Design and Manufacturing Engineering, 2015, 45(05): 28-31].

Abstract3D printing is regarded as a key technology in “made in China 2025” and has been applied in high strength steel, high temperature alloys, plastic rubber etc. But 3D printing still has not been used in petroleum industry because of many technical challenges. For example, arti fi cial core manufacturing faces the obstacles of high production cost and low sample repeatability.Here a potential technological solutions for artificial core manufacturing were discussed. The characteristics of several rock materials and shale powders were analyzed, and a core material candidate for 3D printing was selected out. The in fl uence of components ratios of the materials on rock permeability, porosity and setting time was studied by the response surface method.The central composite design (CCD) was used to optimize the model and a material with best component ratio was gained. The simulation results fi t well with the experiment.

Keywordsrock-like material; 3D printing; response surface method; central composite design; arti fi cial core

(編輯 馬桂霞)

Optimization of rock-like materials for 3D printing

ZHAO Hong1, WU Tingting1, XU Quan2, ZHAO Yixin3
1 College of Mechanical and Transportation Engineering, China University of Petroleum-Beijing, Beijing 102249, China
2 State Key Laboratory of Heavy Oil/New Energy Research Institute, China University of Petroleum-Beijing, Beijing 102249,China
3 College of Resources & Safety Engineering, China University of Mining & Technology, Beijing 100083, China

2017-07-13

國(guó)家自然科學(xué)基金面上項(xiàng)目(51575528)資助

趙弘, 吳婷婷, 徐泉, 趙毅鑫. 3D打印類巖石材料的實(shí)驗(yàn)優(yōu)化研究. 石油科學(xué)通報(bào), 2017, 03: 399-412

ZHAO Hong, WU Tingting, XU Quan, ZHAO Yixin. Optimization of rock-like materials for 3D printing. Petroleum Science Bulletin,2017, 03: 399-412. doi: 10.3969/j.issn.2096-1693.2017.03.037

10.3969/j.issn.2096-1693.2017.03.037