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

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

1 中國石油大學(xué)(北京)機械與儲運工程學(xué)院，北京 102249

2 中國石油大學(xué)(北京)重質(zhì)油國家重點實驗室/新能源研究院，北京 102249

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

*通信作者, ttabbmmjj@163.com

石油機械

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

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

1 中國石油大學(xué)(北京)機械與儲運工程學(xué)院，北京 102249

2 中國石油大學(xué)(北京)重質(zhì)油國家重點實驗室/新能源研究院，北京 102249

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

*通信作者, ttabbmmjj@163.com

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

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

0 引言

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

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

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

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

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

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

1 實驗設(shè)計

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

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

1.1 實驗材料

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

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

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

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

(2)顏色潔白，且便于上色，可實現(xiàn)彩色打印；

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

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

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

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

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

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

1.2.1 實驗設(shè)計及操作

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

實驗操作如下：

(1)取樣

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

(2)拌勻

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

(3)填料

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

(4)打印

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

(5)清洗器具

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

1.2.2 實驗結(jié)果及分析

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

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

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

表1 實驗材料-石膏粉Table 1 The experimental materials—gypsum powder

表2 實驗材料-水泥Table 2 The experimental materials—cement

表3 實驗水平及因素Table 3 Test levels and factors

表4 配比方案Table 4 The solution composition

表5 檢測條件Table 5 Measure conditions

表6 正交表Table 6 Orthogonal array

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

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

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

2.1 實驗操作及結(jié)果

選取表4所示9組配比方案進行打印實驗，得到各配比材料的初凝時間和終凝時間，如表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é)時間Table 9 Setting time

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

2.2 實驗結(jié)果分析

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

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

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

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

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

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

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

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

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

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

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

表10 頻率檢驗Table 10 Frequency check

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

圖9 材料A(石膏)初凝時間預(yù)測值與測定值圖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可知，當材料A(石膏)為5 g時，在水泥和水的取值范圍內(nèi)，等高線基本成水平分布，此時凝結(jié)時間隨水泥質(zhì)量變化不顯著，但隨著水的體積增加，水與水泥的交互作用顯著增加。當石膏取10 g，等高線的切線角度增大，且出現(xiàn)圓心，即石膏取10 g左右時凝結(jié)時間達到最佳狀態(tài)。當石膏取15 g時，水泥與水的交互作用下降，凝結(jié)時間基本隨水泥質(zhì)量的增加而增大。

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

(2)當材料B(水泥)在區(qū)間范圍內(nèi)取5 g、10 g、15 g時，材料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)，當材料B(水泥)的值一定時，材料A(石膏)與材料C(水)的質(zhì)量呈線性相關(guān)關(guān)系，且相關(guān)性大于水泥與水，即石膏的水溶性更好。

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

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

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

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

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

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

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

圖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作為空白組與最佳配比組做比較。空白組的拉壓比為1.27，最佳配比組的拉壓比為1.52。可見優(yōu)化后的材料凝結(jié)時間更短，抗壓強度更高，拉壓比更大，如圖16和圖17所示。

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

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

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

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

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

圖17 最佳配比組與空白組抗壓強度對比Fig17 The optimal proportion compared with blank group compressive strength

3 結(jié)論

通過以上實驗分析可得結(jié)論如下：

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

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

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

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

[1] 王忠宏, 李揚帆, 張曼茵. 中國3D打印產(chǎn)業(yè)的現(xiàn)狀及發(fā)展思路[J]. 經(jīng)濟縱橫, 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] 費群星, 張雁, 譚永生, 等．激光近凈成形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] 周新桂, 操成杰, 袁嘉音. 儲層構(gòu)造裂縫定量預(yù)測與油氣滲流規(guī)律研究現(xiàn)狀和進展[J] .地球科學(xué)進展, 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] 侯泉林, 李會軍, 范俊佳, 等. 構(gòu)造煤結(jié)構(gòu)與煤層氣賦存研究進展[J]. 中國科學(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)力場的可視化方法[J]. 科學(xué)通報, 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)友, 杜時貴, 黃曼, 等. 一種直剪試驗結(jié)構(gòu)面模型的制作方法:中國, 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è)計法優(yōu)化3′,5′-二辛酰基-氟苷藥質(zhì)體制備工藝[J]. 藥學(xué)學(xué)報, 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]. 中國經(jīng)濟信息, 2012, (2): 18-19. [QIU B S. The power of the third industrial revolution[J].China Economic Information, 2012, (2): 18-19.]

[25] 翁德凱, 程寓, 夏玲玲. 基于結(jié)合面的立式加工中心進給系統(tǒng)的動態(tài)特性分析[J]. 機械設(shè)計與制造, 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)面方法在多種實際優(yōu)化問題中的應(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] 張永莉, 霍書娟, 高建平, 等. 硫酸鈣復(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打印一個時代的開啟[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] 張巨松, 楊奇, 李宗陽. 水灰比/礦物外加劑對水泥的收縮變形實驗[J]. 沈陽建筑大學(xué)學(xué)報(自然科學(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] 張萬勝. 石膏模具生產(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)面方法的改進及其對工程優(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] 趙開斌. 關(guān)于學(xué)生化殘差的一點注記[J]. 淮北煤師院學(xué)報, 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)可靠性計算方法—響應(yīng)面法[J]. 上海交通大學(xué)學(xué)報, 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打印噴頭泥漿流動規(guī)律分析[J]. 機械設(shè)計與制造工程, 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

國家自然科學(xué)基金面上項目(51575528)資助

趙弘, 吳婷婷, 徐泉, 趙毅鑫. 3D打印類巖石材料的實驗優(yōu)化研究. 石油科學(xué)通報, 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