，, 2， ，，

(1.同濟大學 地下建筑與工程系，上海 200092；2.紹興文理學院 土木工程學院，浙江 紹興 312000)

1 研究背景

節(jié)理巖體的力學和水力性質(zhì)主要受其內(nèi)部的節(jié)理影響，研究單個節(jié)理的力學和水力性質(zhì)是研究節(jié)理巖體的基礎(chǔ)。節(jié)理表面是粗糙不平的，節(jié)理面的粗糙特性對節(jié)理的剪切特性起到了主導作用[1-3]；上下節(jié)理面的粗糙程度及二者所組成的節(jié)理空腔更是影響節(jié)理滲流特性的主要因素，在節(jié)理剪切過程中，上下節(jié)理面在水平向發(fā)生剪切的同時在法向因為突起的存在而發(fā)生爬坡作用，2個方向的位移變化導致了節(jié)理空腔不斷地發(fā)生變化[4-5]。由于現(xiàn)有技術(shù)水平的限制，難以獲得節(jié)理剪切過程中的表面形貌數(shù)據(jù)及節(jié)理空腔數(shù)據(jù)。目前常用的解決方法是首先獲得初始狀態(tài)下節(jié)理表面形貌數(shù)據(jù)及節(jié)理空腔數(shù)據(jù)，在節(jié)理剪切過程中，將上下節(jié)理面數(shù)據(jù)根據(jù)剪切位移-法向位移進行錯位，在這一過程中上下節(jié)理面相交的區(qū)域被認為是剪切破壞區(qū)域或者是發(fā)生彈性變形區(qū)域予以刪除，最后獲得節(jié)理剪切過程中的節(jié)理面形貌數(shù)據(jù)及節(jié)理空腔數(shù)據(jù)，進而為相應(yīng)的研究提供模型基礎(chǔ)[4,6]。故節(jié)理面初始形貌數(shù)據(jù)及接觸狀態(tài)是進行節(jié)理剪切和滲流特性研究的基礎(chǔ)。

節(jié)理初始表面形貌數(shù)據(jù)的獲得在技術(shù)方面已經(jīng)比較成熟，目前常用的方法主要是使用三維激光掃描儀、光柵投影式形貌儀或者其他攝影測量設(shè)備來對單個節(jié)理面進行掃描[7-9]，進而獲得節(jié)理面形貌數(shù)據(jù)。節(jié)理空腔的測量技術(shù)主要有：①利用低熔點的金屬或者樹脂填充節(jié)理空腔，對填充物厚度分布進行測量進而獲得節(jié)理空腔數(shù)據(jù)[6,10-11],但是由于節(jié)理接觸的存在，在試驗過程中難以保證節(jié)理空腔被完全填充，而且此方法只能在試驗前或試驗后使用[12];②Re等[13]、Renshaw等[14]分別使用X射線和核磁共振的方法來測量節(jié)理空腔，但是這些方法只能定性地確定節(jié)理初始隙寬分布范圍，難以作為后續(xù)計算分析的基礎(chǔ);③Koyama等[15]、Tang等[16]應(yīng)用數(shù)值方法來計算節(jié)理初始空腔數(shù)據(jù)，即首先獲得上下節(jié)理面數(shù)據(jù)，然后將上節(jié)理面置于下節(jié)理面之上，通過在豎向移動、旋轉(zhuǎn)上節(jié)理面，使上下節(jié)理面有3個穩(wěn)定接觸點且上節(jié)理面重心位于3點所組成的三角形內(nèi)則計算完畢，此時的節(jié)理接觸狀態(tài)便認為是節(jié)理的初始接觸狀態(tài)，進而可以得到節(jié)理的初始空腔。但是此種方法的弊端在于首先在接觸狀態(tài)計算之前難以確定上下節(jié)理面在平面方向的相對位置，其次3點接觸的判斷標準與上下節(jié)理的實際接觸狀態(tài)存在偏差。

為了解決初始節(jié)理空腔計算存在的上述問題，利用自行研制的TJXW-3D 型便攜式巖石節(jié)理表面形貌測量儀分別對貼有標志點的上下節(jié)理面和上下節(jié)理面在初始接觸狀態(tài)下的整個試樣進行掃描，掃描結(jié)束后利用標志點將上下節(jié)理面分別對齊到其在初始接觸狀態(tài)下的坐標系中，進而可以獲得節(jié)理初始空腔數(shù)據(jù)。

2 試驗概況

2.1 試樣制備

通過劈裂花崗巖塊來獲得3組不同形貌的節(jié)理，節(jié)理平面尺寸為200 mm×100 mm；采用硅膠復制節(jié)理面的形貌，然后以此為模板澆筑石膏以制備節(jié)理試樣，分別命名為J1，J2，J3。Barton[2]提出了節(jié)理粗糙度系數(shù)JRC來評價節(jié)理表面粗糙程度，并給出JRC為0～20的10條典型的節(jié)理剖面線。Tse和Cruden[17]提出了JRC與剖面線高度均方根Z2的關(guān)系式,即：

JRC=32.2+32.47lgZ2；

(1)

(2)

式中：xi，zi為剖面線坐標；L為剖面線長度；N為采樣點數(shù)。計算得到的3組節(jié)理的JRC分別為：JRC1=11.5，JRC2=13.2，JRC3=15.3。

2.2 試驗儀器

圖1 形貌測量儀Fig.1 Apparatus of surface topography measurement

節(jié)理表面形貌的測量是本試驗的主要內(nèi)容，本試驗采用的是基于3D光柵投影測量原理的TJXW-3D型便攜式巖石節(jié)理表面形貌測量儀，該測量系統(tǒng)具有精度高、單面掃描速度快、易于操作、可實現(xiàn)多個單幅測量數(shù)據(jù)的自動拼接等特點。系統(tǒng)主要包括機頭、三腳架、計算機和標定板，見圖1。機頭利用螺栓固定在三腳架頂部，方便拆卸，機頭主要由光柵投影裝置和2個CCD攝像機組成。測量時機頭朝下向被測物體投射測量光柵，機頭與被測物體距離在50 cm左右為最佳。系統(tǒng)單次測量精度為0.02 mm，其他具體參數(shù)參見文獻[18]。

壓縮試驗使用巖石節(jié)理全剪切-滲流耦合試驗系統(tǒng)，該試驗設(shè)備由主機加載框架、液壓系統(tǒng)、伺服控制系統(tǒng)、計算機控制及處理系統(tǒng)4大部分組成。該設(shè)備能夠?qū)崿F(xiàn)PC控制全自動加載，同步記錄各項變形、位移、荷載數(shù)據(jù)。法向位移采用位移傳感器測量，位移傳感器測試精度為0.001 mm。

2.3 試驗方法

對于一個節(jié)理試樣，首先能夠分別獲得上下節(jié)理面的形貌數(shù)據(jù)和上下節(jié)理面處于初始閉合狀態(tài)下的相對位置。如果能夠?qū)⑸舷鹿?jié)理面的數(shù)據(jù)分別對齊到上下節(jié)理面處于初始閉合的狀態(tài)，便可以獲得節(jié)理在初始閉合狀態(tài)下的空腔數(shù)據(jù)。TJXW-3D型便攜式巖石節(jié)理表面形貌測量儀自帶的測量軟件可以識別特定的標志點為特征點，并對具有共同特征點的數(shù)據(jù)按照特征點之間的距離不變原則進行拼接。首先要對側(cè)面貼有標志點的處于初始閉合狀態(tài)的節(jié)理試樣進行掃描，在此之后，分別對上下節(jié)理面進行掃描，側(cè)面的標志點同時也被掃描記錄，軟件會自動根據(jù)貼在試樣側(cè)面的標志點將上下節(jié)理面分別對齊到初始偶合接觸狀態(tài)下的掃描數(shù)據(jù)。這樣，初始偶合狀態(tài)下節(jié)理空腔可以根據(jù)式(3)計算得到。

zaperture(xi,yi)=zupper(xi,yi)-zlower(xi,yi)。

(3)

式中：zaperture(xi,yi)為節(jié)理空腔數(shù)據(jù);zupper(xi,yi)，zlower(xi,yi)分別為上、下節(jié)理面形貌數(shù)據(jù)。

在獲得節(jié)理形貌數(shù)據(jù)的過程中，每個節(jié)理面要在不同方位進行多次掃描才可以拼接得到一個完整的節(jié)理面數(shù)據(jù)，因為只有同時出現(xiàn)在2個CCD鏡頭中的對象才會被掃描成功。

節(jié)理閉合試驗在巖石節(jié)理全剪切-滲流耦合試驗系統(tǒng)上進行，采用4支位移傳感器采集節(jié)理法向閉合變形，并以其平均值來表征節(jié)理法向閉合變形。法向荷載為80 kN(4 MPa)，加載速率為1 kN/s。

3 初始接觸狀態(tài)下節(jié)理空腔分布

利用本文第2節(jié)介紹的方法分別對試樣J1，J2，J3的上下節(jié)理面的形貌數(shù)據(jù)和節(jié)理空腔結(jié)果進行了計算，計算過程見圖2，計算結(jié)果見圖3。

圖2 節(jié)理空腔計算過程Fig.2 Computation of void space of joint

圖3 節(jié)理空腔計算結(jié)果Fig.3 Computation result of joint’s void space

圖3(a)給出的各節(jié)理上下表面形貌數(shù)據(jù)可以看到，節(jié)理J1表面具有較大的一階起伏度，J2表面則明顯以二階粗糙度為主，J3節(jié)理表面一階起伏度和二階粗糙度都可以清楚地分辨出來。

圖3(b)給出了黑色區(qū)域代表節(jié)理接觸區(qū)域。J1上下節(jié)理接觸區(qū)域較小，J2上下節(jié)理面接觸區(qū)域零星分布于整個節(jié)理空腔，J3的節(jié)理接觸則主要分布于邊緣位置，且接觸面積為J3>J2>J1。

從圖3(c)給出的各節(jié)理隙寬概率分布也可以得到同樣的結(jié)論，J3節(jié)理的初始接觸面積最大，在3.3%左右。根據(jù)計算得到的節(jié)理空腔結(jié)果，統(tǒng)計了不同隙寬的概率分布，可見各組節(jié)理空腔分布遵循高斯分布。從圖3(c)中也可以看出節(jié)理初始平均隙寬J1(0.41 mm)>J2(0.35 mm)>J3(0.31 mm)，這與Hakami等[12]的研究結(jié)果基本吻合。Hakami等[12]使用樹脂材料填充節(jié)理空腔得到了耦合花崗巖的在低法向應(yīng)力(0.45 MPa)作用下的空腔分布，結(jié)果表明節(jié)理隙寬分布同樣均遵循高斯分布，且統(tǒng)計的初始接觸面積均在5%以下，平均隙寬為0.36 mm。

4 節(jié)理空腔數(shù)據(jù)的應(yīng)用

Greenwood等[19]將粗糙表面的微凸體視為球形，基于Hertz理論和概率統(tǒng)計理論提出了光滑平面與粗糙平面接觸的變形模型，以下簡稱Greenwood模型。Brown和Scholz[20]在Greenwood模型的基礎(chǔ)上，利用節(jié)理的組合形貌提出了節(jié)理閉合模型——Brown模型，即：

(4)

(5)

式中：σn為法向應(yīng)力;ηc為節(jié)理組合形貌的峰點密度；βc為節(jié)理組合形貌的峰點平均半徑；〈ψ〉為切向應(yīng)力相關(guān)因子的平均值，通常〈ψ〉≈1；E′為有效彈性模量；E，ν分別為彈性模量和泊松比；zc為隨機變量，表示組合形貌某峰頂?shù)絽⒖济娴木嚯x；d為上下節(jié)理面參考面之間的距離;δ為節(jié)理法向閉合變形。

計算節(jié)理的組合形貌是Brown模型的基礎(chǔ)，圖4給出了節(jié)理空腔與節(jié)理組合形貌之間的關(guān)系。在這里，利用本文的節(jié)理空腔數(shù)據(jù)來計算得到節(jié)理的組合形貌，進而可以方便地應(yīng)用Brown模型，并與實際的壓縮試驗結(jié)果進行對比。

圖4 節(jié)理組合形貌與節(jié)理空腔的關(guān)系Fig.4 Relationship between composite topography and void space

由圖5可見，在同一法向應(yīng)力下，Brown模型預測結(jié)果總是比法向位移實際試驗結(jié)果大。經(jīng)分析，主要原因在于:雖然Brown模型是在Greenwood模型上進行改進的，但其基本假設(shè)相同，即認為突起具有完全相同的突起半徑且互不影響;在實際情況下，突起在被壓縮的過程中會牽動突起周邊的節(jié)理面,受力面積不僅限于各突起(實際受力的不只是Greenwood所假設(shè)的突起)，而且各突起互相之間具有約束作用。所以Brown模型預測的節(jié)理變形會比試驗結(jié)果要大。

圖5 Brown模型預測結(jié)果與試驗結(jié)果的比較Fig.5 Comprison between Brown model’s result and test result

5 結(jié) 論

本文介紹了確定節(jié)理在初始接觸狀態(tài)下空腔分布的新方法，初步探索了3個不同粗糙度節(jié)理試樣的初始隙寬分布規(guī)律以及節(jié)理初始空腔數(shù)據(jù)在Brown法向閉合模型中的應(yīng)用。所介紹的方法在節(jié)理剪切-滲流的研究中，尤其是為節(jié)理剪切滲流的數(shù)值研究奠定了模型基礎(chǔ)。

(1)利用自行研制的TJXW-3D型便攜式巖石節(jié)理表面形貌測量儀分別對上下節(jié)理面和上下節(jié)理面在初始接觸狀態(tài)下的整個試樣進行掃描，掃描結(jié)束后利用標志點將上下節(jié)理面分別對齊到其在初始接觸狀態(tài)下的坐標系中，進而可以獲得節(jié)理初始空腔數(shù)據(jù)。

(2)利用本文介紹的節(jié)理初始空腔計算方法對3個不同粗糙度的節(jié)理進行計算，得到了各節(jié)理初始空腔數(shù)據(jù)，并統(tǒng)計了各節(jié)理隙寬的分布規(guī)律。各節(jié)理隙寬分布符合高斯分布，平均隙寬為0.31～0.41 mm，與Hakami等的研究所得結(jié)論基本一致。

(3)將本文計算得到的節(jié)理初始空腔數(shù)據(jù)應(yīng)用到節(jié)理法向閉合Brown模型中，Brown模型預測結(jié)果比實際試驗結(jié)果相對較大，經(jīng)分析主要原因在于Brown模型無法考慮各突起之間的相互作用。

[1] PATTON F D. Multiple Modes of Shear Failure in Rock[C]∥International Society for Rock Mechanics and Rock Engineering，1st ISRM Congress，Lisbon，Portugal，September 25-October 1，1966：509-513.

[2] BARTON N. Review of a New Shear-strength Criterion for Rock Joints[J]. Engineering Geology，1973，7(4)：287-332.

[3] 唐志成，夏才初，宋英龍，等.Grasselli節(jié)理峰值抗剪強度公式再探[J]. 巖石力學與工程學報，2012，31(2) ：356-364.

[4] XIA C C，GUI Y，WANG W，etal. Numerical Method for Estimating Void Spaces of Rock Joints and the Evolution of Void Spaces under Different Contact States[J]. Journal of Geophysics and Engineering，2014，11(6):065004.doi:10.1088/1742-2132/11/6/065004.

[5] GUI Y，XIA C C，DING W Q，etal. A New Method for 3D Modeling of Joint Surface Degradation and Void Space Evolution under Normal and Shear Loads[J]. Rock Mechanics and Rock Engineering, 2017，50(10):2827-2836.

[6] PARK J W, SONG J J. Numerical Method for the Determination of Contact Areas of a Rock Joint under Normal and Shear Loads[J]. International Journal of Rock Mechanics and Mining Sciences，2013，58:8-22

[7]GRASSELLI G, WIRTH J, EGGER P. Quantitative Three-dimensional Description of a Rough Surface and Parameter Evolution with Shearing[J]. International Journal of Rock Mechanics and Mining Sciences，2002，39(6):789-800.

[8]BELEM T，SOULEY M，HOMAND E F. Method for Quantication of Wear of Sheared Joint Walls Based on Surface Morphology[J]. Rock Mechanics and Rock Engineering，2009，42(6):883-910.

[9] FATHI A，MORADIAN Z，RIVARD P，etal. Geometric Effect of Asperities on Shear Mechanism of Rock Joints[J]. Rock mechanics and Rock Engineering，2016，49(3):801-820.

[10] GENTIER S, BILLAUX D,VAN VLIET L. Laboratory Testing of the Voids of a Fracture[J]. Rock Mechanics and Rock Engineering，1989，22(2):149-157.

[11] NEMOTO K，WATANABE N，HIRANO N，etal. Direct Measurement of Contact Area and Stress Dependence of Anisotropic Flow Through Rock Fracture with Heterogeneous Aperture Distribution[J]. Earth and Planetary Science Letters，2009，281(1/2):81-87.

[12] HAKAMI E,LARSSON E.Aperture Measurements and Flow Experiments on a Single Natural Fracture[J].International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts，1996，33(4):395-404.

[13] RE F, SCAVIA C. Determination of Contact Areas in Rock Joints by X-ray Computer Tomography[J]. International Journal of Rock Mechanics and Mining Sciences，1999，36(7):883-890.

[14] RENSHAW C E，DADAKIS J S，BROWN S R. Measuring Fracture Apertures: A Comparison of Methods[J]. Geophysical Research Letters，2000，27(2):289-292.

[15] KOYAMA T，F(xiàn)ARDIN N，JING L，etal. Numerical Simulation of Shear-induced Flow Anisotropy and Scale-dependent Aperture and Transmissivity Evolution of Rock Fracture Replicas[J]. International Journal of Rock Mechanics and Mining Sciences，2006，43(1):89-106.

[16] TANG Z C, LIU Q S, XIA C C，etal. Mechanical Model for Predicting Closure Behavior of Rock Joints Under Normal Stress[J]. Rock Mechanics and Rock Engineering，2014，47(6):2287-2298.

[17] TSE R, CRUDEN D M. Estimating Joint Roughness Coefficients[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts，1979，16(5):303-307.

[18] 夏才初，王 偉，丁增志.TJXW-3D型便攜式巖石三維表面形貌儀的研制[J].巖石力學與工程學報，2008，27(7)：1505-1512.

[19]GREENWOOD J A，TRIPP J H. The Contact of Two Nominally Flat Rough Surfaces[J]. Proceedings of the Institution of Mechanical Engineers，1970，185：625-633.

[20] BROWN S R，SCHOLZ C H. Closure of Random Elastic Surfaces in Contact[J]. Journal of Geophysical Research，1985，90(B7)：5531-5545.