注采參數(shù)對松遼盆地干熱巖物理力學及波動特征的影響

崔翰博，唐巨鵬，姜昕彤，邱于曼

遼寧工程技術(shù)大學力學與工程學院，阜新123000

1 前言

增強型地熱系統(tǒng)（Enhanced Geothermal System，EGS）將開發(fā)干熱巖（Hot Dry Rock，HDR）作為研究目的，通過水力壓裂（Hydrofracturing method）的方式致裂儲層，利用水與HDR熱交換，提取熱量進行應用。HDR溫度是EGS熱采選址的首要因素，其決定了靶區(qū)的開采前景與運營成本；水作為攜熱介質(zhì)，水溫是熱采的重要控制方法之一；高溫遇水循環(huán)次數(shù)也與儲層壽命密切相關(guān)，因此開展巖樣溫度、水溫、高溫遇水循環(huán)次數(shù)與花崗巖物理力學、波動特征間關(guān)系的研究對EGS高質(zhì)量開采具有重要意義。

近年來，國內(nèi)外相關(guān)專家針對注采參數(shù)與EGS關(guān)系已進行了大量的數(shù)值模擬研究。Jing等（2014）、Vogt等（2012）、趙陽升等（2004）發(fā)現(xiàn)儲層初始溫度與EGS產(chǎn)能密切相關(guān)，溫度越高，HDR產(chǎn)熱效率越好。隨采熱進行，當儲層溫度降幅達到一定程度，產(chǎn)能降低。Fox等（2013）和翟海珍等（2017）指出產(chǎn)能降低，采熱效率明顯下降，此時需停止開采，待熱能恢復后再進行回采，在此期間HDR會經(jīng)歷多次高溫遇水冷卻循環(huán)過程。雷宏武（2014）和肖勇（2017）指出水溫與產(chǎn)熱效率、裂縫開度、儲層沉降存在一定關(guān)系。此外，部分專家學者對高溫條件下巖石的物理力學、波動特征進行了試驗研究。Jiang等（2018）進行了高溫后花崗巖波動特征測試試驗，指出250～300℃為裂縫加速擴展階段。Gautam等（2018）進行了高溫下花崗巖單軸抗壓試驗，指出300℃是花崗巖力學特征出現(xiàn)下降的起始溫度。秦嚴（2017）對自然冷卻后花崗巖進行了單軸壓縮實驗，指出600℃是花崗巖力學性質(zhì)產(chǎn)生突變的閾值。黃真萍等（2016）進行了高溫遇水冷卻后大理巖、石灰?guī)r物理力學特征試驗研究，并對損傷因子進行了簡單界定。郤保平等（2010）研究發(fā)現(xiàn)高溫遇水冷卻后花崗巖的力學特征會產(chǎn)生明顯下降。Li和Ju（2018）進行了多次熱循環(huán)后花崗巖力學、波動特征試驗，發(fā)現(xiàn)高溫循環(huán)后花崗巖力學特征損失、波速降低多發(fā)生在5個熱循環(huán)之內(nèi)。

基于以上研究背景發(fā)現(xiàn)，前人多以產(chǎn)熱效率為主要目的對EGS熱采進行了數(shù)值模擬，對高溫花崗巖的探究也多為自然冷卻后或?qū)崟r高溫條件下巖樣的物理力學特征。以HDR為研究背景，進行高溫遇水冷卻或高溫遇水循環(huán)冷卻對巖樣物理力學特征影響的研究極為少見。此外，注采、熱恢復過程中低溫巖體受高溫區(qū)域影響，含水量是逐漸變化的，進而引起力學、波動特征發(fā)生改變，考慮擱置過程能夠提高對恢復過程中壓裂區(qū)巖體力學、波動特征測量的準確性。另外，HDR處于深部地層，采用大量拉、壓試驗盡管能夠獲取巖樣的物理力學性質(zhì)，但在實際研究中并不現(xiàn)實。而通過超聲波測試的方法能夠很好的檢測出熱采靶區(qū)巖樣物理力學特征，且操作簡單，適于在EGS熱采中應用。基于以上因素，為探究注采參數(shù)對HDR物理力學和波動特征的作用，以松遼盆地HDR開采為研究背景，對不同注采參數(shù)下巖樣進行單軸抗壓試驗和縱、橫波波速測試試驗，研究巖樣溫度、水溫、高溫遇水循環(huán)次數(shù)對波動特征、物理力學特征的作用，并將其建立聯(lián)系，同時考慮擱置過程中HDR物理力學及波動特征變化規(guī)律。此研究有望為HDR力學特征無損檢測提供一種手段，同時為松遼盆地地熱田注采參數(shù)的選取提供參考。

2 試驗概況

2.1 地質(zhì)背景

松遼盆地地跨內(nèi)蒙古、黑龍江、吉林、遼寧四省，受地質(zhì)構(gòu)造運動影響，盆地基底斷裂發(fā)育明顯，火山活動頻繁，存在多個地熱異常點（章鳳奇等，2008）。地熱田熱源主要分為以下三種：地幔對流、巖漿體侵入、巖體內(nèi)部放射性元素衰變。在三者共同作用下，大地熱流值在51.5～90.0 mW/m2間，平均地熱梯度為5.7℃/100 m，HDR資源總量可達0.62×106EJ（李野，2017）。目標層上部主要由淺變質(zhì)巖系、花崗巖構(gòu)成；下部以深變質(zhì)巖系和片麻狀花崗巖為主。上覆泥巖蓋層，能夠保證熱能的貯存，且松花江水系位于地熱田附近。綜合以上地質(zhì)因素，松遼盆地地熱田具有良好地熱開發(fā)前景。

2.2 巖樣制備

選取松遼盆地北部為EGS熱采目標靶區(qū)進行取樣。巖樣表觀呈灰白色，密度介于2.71～2.99 g/cm3間。利用取芯機、切割機、磨平機依次對巖樣進行加工、打磨。巖樣尺寸為φ50 mm×100 mm，樣高、直徑誤差均小于3 mm），共進行15組試驗，每組3塊。

2.3 試驗方案

試驗方案設計：HDR溫度介于100～650℃，因此巖樣分別加熱至100℃、200℃、300℃、400℃、500℃、600℃，以下簡稱100～600℃（水溫20℃，高溫遇水循環(huán)次數(shù)1次）；松花江水溫介于1～25℃，考慮注水過程水溫上升，因此水溫分別選取1℃、20℃、40℃、60℃、80℃，以下簡稱1～80℃（巖樣溫度300℃，高溫遇水循環(huán)次數(shù)1次）；考慮循環(huán)注水過程和熱恢復作用（Fox，2013；翟海珍，2017），高溫遇水循環(huán)次數(shù)選取1、2、3、4、5次，以下簡稱循環(huán)1～5次（Li等，2018）（巖樣溫度300℃，水溫20℃）。

2.4 試驗設備

表1為主要試驗設備的詳細參數(shù)。

表1 試驗設備Table 1 Experiment equipments

2.5 試驗內(nèi)容

試驗內(nèi)容分為：高溫遇水循環(huán)冷卻試驗、縱、橫波波速測試試驗、單軸抗壓試驗。遵循《巖石物理力學試驗規(guī)程》（中華人民共和國國土資源部，2015），并根據(jù)試驗條件進行調(diào)整：

天然狀態(tài)巖樣性質(zhì)測定：天平測巖樣質(zhì)量；游標卡尺測樣高、直徑；超聲檢測儀測縱、橫波波速。

縱、橫波波速測試試驗：為避免巖樣上下表面凹凸對試驗影響，用凡士林將上、下面和探頭間耦合，然后將巖樣置于收、發(fā)傳感器間，發(fā)射聲波，采集數(shù)據(jù)（波速測量精度為0.01 m/s）。

不同試驗方案冷卻過程：不同巖樣溫度：利用電爐以30℃/min，分別加熱巖樣至設置溫度（100～600℃），為保證巖樣內(nèi)外受熱均勻，恒溫4 h（閆治國等，2006），然后將其置于堿骨料反應箱中遇水冷卻4 h（水溫20℃，試驗發(fā)現(xiàn)巖樣溫度不同，遇水冷卻4 h后均降至20℃左右，說明水與巖樣間不再傳熱，且基本飽水）；不同水溫：將巖樣加熱至300℃恒溫4 h后，遇水冷卻4 h（水溫分別保持在1～80℃）；高溫遇水冷卻循環(huán)次數(shù)：將300℃巖樣恒溫4 h，置于20℃水中降溫冷卻4 h，此過程循環(huán)1～5次。

擱置過程：將冷卻后巖樣置于干燥箱中（溫度20℃，只控溫不鼓風），擱置0、2、4、6、8 h，測量巖樣質(zhì)量和縱、橫波波速。

單軸抗壓試驗：擱置8 h，取出巖樣，在巖樣兩側(cè)中間對稱位置貼應變片，連接動態(tài)信號測試分析系統(tǒng)，然后將其置于試驗機上，以0.5 MPa/s加載速度進行單軸抗壓試驗，待試件破壞后，導出數(shù)據(jù)。

3 試驗結(jié)果及分析

3.1 不同注采參數(shù)下巖樣外觀特征

提高巖樣溫度、增加遇水循環(huán)次數(shù)，巖樣均由灰白向土黃色轉(zhuǎn)變。圖1為不同注采參數(shù)下巖樣外觀特征。巖樣溫度：100～200℃，巖樣與天然條件下相比變化微小，表面呈灰白色，局部存在黑斑；300～400℃，顏色變淺，黑斑略有減少；500～600℃，部分區(qū)域呈土黃色，局部存在灰斑，產(chǎn)生白色晶狀物。水溫：水溫變化，巖樣外觀特征改變不明顯。高溫遇水循環(huán)次數(shù)：循環(huán)次數(shù)增加，巖樣表觀黑斑減少，循環(huán)至5次，表面產(chǎn)生土黃色區(qū)域。由此可見，對巖樣外觀特征影響由大到小依次為巖樣溫度、高溫遇水循環(huán)次數(shù)、水溫，通過外觀特征的變化能夠判斷熱采區(qū)域所經(jīng)歷的注采條件。

3.2 擱置過程中巖樣力學及波動特征變化規(guī)律

注采、熱恢復過程中低溫區(qū)域受高溫巖體影響，含水量逐漸下降，進而使壓裂區(qū)巖體物理力學、波動特征產(chǎn)生一定改變，忽略此現(xiàn)象會使測量結(jié)果產(chǎn)生較大誤差。因此考慮擱置過程對HDR物理力學、波動特征的作用十分必要。

圖1 不同注采參數(shù)下巖樣外觀特征Fig.1 Granite superficial morphology under different injection-production parameters

3.2.1 擱置過程中不同溫度下巖樣物理力學及波動特征變化規(guī)律

試驗發(fā)現(xiàn)：擱置初期，巖樣溫度越高，質(zhì)量、縱橫波波速、彈性模量降幅越大。

由于擱置初期是巖樣物理力學及波動特征發(fā)生改變的主要時刻，因此本文以擱置0～2 h為例，以100℃、600℃的試驗結(jié)果作為研究對象展開分析，彈性模量、泊松比、相鄰時間段參數(shù)變化率公式如（1）～（3）。

式中：E為彈性模量/GPa， μ為泊松比，VP為縱波波速/（m/s），VS為橫波波速/（m/s），P為相鄰時間段參數(shù)變化率，X為參數(shù)（依次為m為質(zhì)量/g，VP，VS，E， μ），i為上一時間段/h，i′為下一時間段/h。注：變化率正為降幅，變化率負為增幅。

圖2a為相鄰時間段巖樣質(zhì)量變化率。擱置0～2 h，100、600℃巖樣質(zhì)量降幅依次為0.066%、0.135%，提高溫度，質(zhì)量降幅增加，600℃巖樣質(zhì)量降幅約為100℃的2.05倍。圖2b,c為相鄰時間段巖樣縱、橫波波速變化率。擱置0～2 h，100℃、600℃縱波波速降幅依次為0.48%、3.46%，橫波波速降幅依次為0.52%、3.45%，提高溫度，縱、橫波波速降幅增加，600℃巖樣縱、橫波波速降幅依次為100℃的7.19和6.63倍。圖2d為相鄰時間段巖樣彈性模量變化率。彈性模量與縱、橫波波速的變化規(guī)律近似。擱置0～2 h，100℃、600℃彈性模量降幅依次為1.08%、6.94%，溫度越高，彈性模量降幅越多，600℃巖樣彈性模量降幅為100℃的6.43倍；圖2e為相鄰時間段巖樣泊松比變化率。能夠看出隨巖樣溫度和擱置時間變化，泊松比變化率呈無規(guī)則波動，因此擱置時間、巖樣溫度對泊松比的影響有待進一步研究。

圖2 相鄰時間段不同溫度下巖樣物理力學及波動特征變化規(guī)律Fig.2 The variationsof physico-mechanical and wave characteristics of rock samplesatdifferenttemperatures in adjacenttime periods

通過試驗結(jié)果看出，擱置初期，不同溫度下巖樣質(zhì)量與力學、波動特征變化規(guī)律呈一定關(guān)聯(lián)性。提高溫度，巖樣內(nèi)部裂縫增多，遇水后吸水量隨之增加（林睦曾，1991）。隨著擱置時間增長，巖樣內(nèi)部的水產(chǎn)生逃逸。溫度越高，水逸出量越多，相同時間段巖樣質(zhì)量降幅越大。而聲波在巖石中的波速遠低于水中的傳播速率，所以縱、橫波波速、彈性模量降幅增大。由此可見，巖樣內(nèi)部含水率的損失是導致其力學、波動特征產(chǎn)生變化的重要因素之一。擱置后期，水大幅逸出，巖樣物理力學及波動特征的變化受試驗條件影響較大，試驗器材的精度、測量間隔的誤差、耦合劑涂抹的均勻程度、巖樣上、下面的粗糙度等均會對試驗結(jié)果產(chǎn)生一定影響。此外，在今后試驗中應縮短擱置初期（0～2 h）的時間間隔，能夠有效提高對熱恢復時間預測的準確性。

3.2.2 擱置過程中不同水溫下巖樣物理力學及波動特征變化規(guī)律

試驗結(jié)果表明：擱置初期，隨水溫升高，縱、橫波波速、彈性模量降幅先增大后減小。

圖3 相鄰時間段不同水溫下巖樣物理力學及波動特征變化規(guī)律Fig.3 The variations of physico-mechanical and wave characteristics of rock samples at different water temperatures in adjacent time periods

圖3 a為相鄰時間段巖樣質(zhì)量變化率。擱置0～2 h，1～80℃ 巖樣質(zhì)量降幅依次為 0.088%、0.086%、0.069%、0.066%、0.069%，1℃時質(zhì)量降幅約為80℃的1.28倍。圖3b，c為相鄰時間段巖樣縱、橫波波速變化率。擱置0～2 h，1～80℃，縱波波速分別降幅0.73%、1.04%、1.40%、1.19%、0.89%，橫波波速降幅0.93%、1.04%、1.50%、1.34%、0.90%，提高溫度，縱、橫波波速降幅均呈先增大后減小的趨勢。圖3d為相鄰時間段巖樣彈性模量變化率。擱置0～2 h，1～80℃彈性模量分別降幅1.82%、2.16%、3.02%、2.67%、1.89%，與縱、橫波波速的變化規(guī)律近似。圖3e為相鄰時間段巖樣泊松比變化率。能夠看出隨水溫和擱置時間變化，泊松比變化率呈上下波動的趨勢，因此水溫對泊松比變化規(guī)律的影響有待進一步討論。

通過試驗結(jié)果看出，擱置初期水溫是影響巖樣力學、波動特征的主要因素。聲波在40～45℃的水中傳播速度最快（蘭朝鳳，2012）。擱置2 h巖樣溫度恢復至室溫，因此水溫40℃的巖樣所受影響最大，力學、波動特征降幅最大。擱置4 h后，5組巖樣均接近室溫，力學及波動特征的變化規(guī)律主要受巖樣內(nèi)部含水率和試驗條件的影響。綜上所述，考慮擱置過程能夠提高對不同水溫條件下注采、熱恢復過程測量的精確度。

3.3 注采參數(shù)對巖樣波動特征的影響

試驗得到：提高巖樣溫度、增加高溫遇水循環(huán)次數(shù)，波速降幅增大；水溫升高，波速降幅減小。

圖4 不同注采參數(shù)下巖樣縱、橫波波速變化規(guī)律Fig.4 Changes of longitudinal and transverse wave velocity of rock samples under different injection-production parameters

參數(shù)相對值如公式（4）。

式中：Q為參數(shù)相對值，R為天然狀態(tài)下參數(shù)（依次為VP，VS，σ為峰值強度（/MPa），E， μ），Rn為高溫后參數(shù)。

波速一定程度上能夠體現(xiàn)熱采過程中壓裂層巖體內(nèi)部結(jié)構(gòu)及裂縫擴展特征。圖4（a-c）為不同注采參數(shù)下巖樣縱、橫波波速變化規(guī)律。考察巖樣溫度對波速的影響可知，100～200℃，巖樣溫度的改變對波速影響較??；300～400℃，波速受巖樣溫度影響明顯，產(chǎn)生較大降幅，結(jié)果與朱振南等（2018）結(jié)論相似。400℃時降幅分別達到30.21%、40.02%；500～600℃，降幅趨勢減弱，經(jīng)歷600℃高溫，巖樣縱、橫波波速降幅依次達到天然狀態(tài)的53.44%、58.02%。由此可見，靶區(qū)巖層溫度越高，對波速降幅影響越大。這是因為巖樣內(nèi)顆粒成分不同，注水冷卻后顆粒間產(chǎn)生差異性膨脹（即熱應力），黏結(jié)力下降，促使拉應力區(qū)結(jié)構(gòu)破碎，原本致密的巖體內(nèi)部產(chǎn)生裂縫，熱應力作用下，裂縫延伸、貫通，形成裂縫區(qū)域，導致聲波傳遞受到影響，波速有較大降幅（方新宇等，2016；Zhang et al.，2016；朱振南等，2018；崔翰博等，2019）。溫度高于300℃，巖樣內(nèi)部結(jié)晶水爆裂轉(zhuǎn)化為自由水，產(chǎn)生大量裂縫，導致外界冷水注入，高溫巖體與冷水間產(chǎn)生復雜反應，促使巖石發(fā)生熱破裂（朱振南等，2018），因此300～400℃間巖樣波速降幅明顯。

考察水溫對波速的影響發(fā)現(xiàn)，其與巖樣溫度不同時，波速間變化規(guī)律相反。水溫1℃，縱、橫波波速分別降低20.55%、18.65%；40～60℃時波速損失相差較小；80℃時縱、橫波波速分別損失12.89%、10.33%。低溫注水，使水巖間溫差增大，提高了降溫速率，使巖樣所受熱沖擊作用增強，產(chǎn)生了更多裂縫，促使波速明顯降低。

考察高溫遇水循環(huán)次數(shù)對縱、橫波波速影響得到：循環(huán)1次，縱、橫波波速分別下降19.55%、16.44%；循環(huán)3次，縱、橫波波速出現(xiàn)明顯下降；循環(huán)至5次，縱、橫波波速損失分別達到33.61%、33.63%，隨循環(huán)次數(shù)增加，波速降幅增大。這是由于經(jīng)過多次遇水循環(huán)過程巖樣會產(chǎn)生軟化現(xiàn)象（趙陽升等，2004；康健，2008；郤保平等，2010），其結(jié)果必然使巖石內(nèi)顆粒間粘結(jié)力減弱，破壞了內(nèi)部骨架結(jié)構(gòu)。此外，熱巖循環(huán)遇水，巖樣遭受多次熱沖擊，導致內(nèi)部產(chǎn)生大量裂縫，也會促使波速降低。綜合以上因素可知，對波速變化規(guī)律影響由大到小的注采參數(shù)依次為靶區(qū)溫度、高溫遇水循環(huán)次數(shù)、注水溫度。

3.4 注采參數(shù)對巖樣力學特征的影響

熱采過程中HDR力學特征是逐漸改變的，進而影響到采熱效率。選擇合適的注采參數(shù)能夠使巖石力學特征向有利于熱采的方向發(fā)展。

3.4.1 不同注采參數(shù)下巖樣受壓后破壞形態(tài)

試驗發(fā)現(xiàn)：巖樣溫度300℃、500℃、水溫低于60℃、遇水循環(huán)至5次，巖樣破壞形態(tài)會發(fā)生改變。

圖5為不同注采參數(shù)下巖樣受壓后破壞形態(tài)。巖樣溫度：100～200℃，巖樣呈軸向拉剪破壞；300～400℃，為斜向剪切破裂；500～600℃，則近似呈橫向錐形破裂。水溫：1～60℃，巖樣破壞形態(tài)為斜向剪切破裂；80℃，呈軸向拉剪破裂。高溫遇水循環(huán)次數(shù)：循環(huán)1～4次，巖樣破壞形態(tài)為斜向剪切破裂；循環(huán)5次，趨近于橫向錐形破裂。通過巖樣受壓后的破壞形態(tài)，可為采熱過程中熱儲區(qū)變化規(guī)律的研究提供一定參考。

3.4.2 不同注采參數(shù)下巖樣單軸受壓下的應力—應變曲線

圖6（a-c）為不同注采參數(shù)下巖樣應力—應變曲線。注采參數(shù)不同巖樣應力—應變曲線均呈現(xiàn)4個過程，分別為（1）孔隙壓密；（2）彈性變形；（3）屈服-破壞；（4）卸荷（郤保平等，2010）。

圖5 不同注采參數(shù)下巖樣受壓后破壞形態(tài)Fig.5 Destruction morphology of rock samples under different injection and production parameters

圖6 不同注采參數(shù)下巖樣應力—應變曲線Fig.6 Stress-strain curves of rock samples under different injection-production parameters

考察巖樣溫度對應力—應變曲線影響可知：提高溫度，巖樣孔隙率增大，縱向應變增加，彈性過程斜率降低；100～300℃對屈服—破壞過程影響較小，破壞時失穩(wěn)過程為突發(fā)失穩(wěn)；400～600℃，屈服—破壞過程延長，卸荷過程斷面滑移增大，縱向應變增加。失穩(wěn)過程由400℃時的突發(fā)失穩(wěn)過渡為500℃時的準突發(fā)失穩(wěn)。

考察水溫對應力—應變曲線影響得到：水溫改變，應力—應變曲線無明顯變化。

考察高溫遇水循環(huán)次數(shù)對應力—應變曲線影響得到：增加循環(huán)次數(shù)與提高巖樣溫度對巖樣應力—應變曲線變化規(guī)律影響相似。循環(huán)次數(shù)增加，彈性過程斜率降低，經(jīng)歷5次循環(huán)后屈服—破壞過程巖樣的失穩(wěn)過程由突發(fā)失穩(wěn)過渡為準突發(fā)失穩(wěn)。

由此可見，巖樣溫度和高溫遇水循環(huán)次數(shù)改變會使巖樣應力—應變曲線發(fā)生較大改變，水溫對其影響較小。

3.4.3 不同注采參數(shù)下巖樣峰值強度變化規(guī)律

試驗表明：提高巖樣溫度、增加高溫遇水循環(huán)次數(shù)、低溫注水均會使巖樣峰值強度產(chǎn)生較大損失。

圖7（a-c）為不同注采參數(shù)下巖樣峰值強度變化規(guī)律?？疾鞄r樣溫度對巖樣峰值強度影響得到：100～200℃，溫度對峰值強度影響較?。惶岣邷囟?，峰值強度降低，300～400℃峰值強度產(chǎn)生明顯降低，降幅為36.60%；經(jīng)歷600℃高溫，與天然狀態(tài)相比，峰值強度損失達到66.56%，巖樣溫度升高，峰值強度損失隨之增大，所得結(jié)果與郤保平等（2010）得到的規(guī)律相似。

考察水溫對巖樣峰值強度影響得到：1～80℃時，巖樣峰值強度分別為127.57 MPa、129.86 MPa、133.66 MPa、 132.19 MPa、 135.15 MPa， 水 溫 升高，巖樣峰值強度隨之上升。低溫注水會加速巖層破碎產(chǎn)生新裂縫。

考察高溫遇水循環(huán)次數(shù)對巖樣峰值強度影響得到：循環(huán)次數(shù)增加，巖樣峰值強度降低，循環(huán)5次，強度損失達到天然狀態(tài)的34.22%。

由此可見，巖樣峰值強度與縱、橫波波速變化規(guī)律相似，為達到巖層破碎增透儲層的目的，可選取高溫儲層作為靶區(qū)鉆井，并選取低溫注水、多次冷水循環(huán)的方法。

3.4.4 不同注采參數(shù)下巖樣彈性模量變化規(guī)律

試驗發(fā)現(xiàn)：巖樣溫度300～400℃、循環(huán)至3次，彈性模量降幅最為明顯；注水溫度對彈性模量影響較小。

彈性模量決定了HDR對水力壓裂及循環(huán)過程中對彈性變形的抵抗能力。圖8（a-c）為不同注采參數(shù)下巖樣彈性模量變化規(guī)律?？疾鞄r樣溫度變化對彈性模量影響得到：300～400℃，巖樣彈性模量明顯下降，降幅為64.04%；經(jīng)歷600℃高溫，彈性模量降幅為79.48%，由此可見，經(jīng)歷高溫后巖樣抵抗彈性變形的能力明顯降低。

考察水溫對巖樣彈性模量影響得到：水溫升高，巖樣彈性模量增大，0～40℃，彈性模量呈階梯形增長，40～60℃水溫對彈性模量影響不大，60～80℃彈性模量又有所增加。1℃與80℃彈性模量相差為13.95%。因此采用低溫注水能夠使HDR變形增大，促進儲層裂縫的發(fā)育。

圖7 不同注采參數(shù)下巖樣峰值強度變化規(guī)律Fig.7 Changes of compressive strength of rock samples under different injection-production parameters

考察高溫遇水循環(huán)次數(shù)對巖樣彈性模量影響得到：循環(huán)次數(shù)增加，彈性模量逐漸下降。循環(huán)次數(shù)2～3次，彈性模量降幅最大，為15%左右，經(jīng)歷5次高溫循環(huán)，彈性模量降幅達到56%。由此可知，高溫多次遇水循環(huán)后目標熱儲層會產(chǎn)生再次變形，應對其進行實時監(jiān)測，保證采熱效率。

圖8 不同注采參數(shù)下巖樣彈性模量變化規(guī)律Fig.8 Changes of elastic modulus of rock samples under different injection-production parameters

3.4.5 不同注采參數(shù)下巖樣泊松比變化規(guī)律

泊松比為水力壓裂及循環(huán)過程中巖體橫、縱向應變之比的絕對值。圖9（a-c）為不同注采參數(shù)下巖樣泊松比變化規(guī)律。試驗結(jié)果表明：溫度100～200℃對巖樣泊松比的影響可以忽略，300℃時泊松比逐漸下降，600℃時泊松比降幅達到21.35%。因此經(jīng)歷高溫后巖體縱向應變變化率明顯高于橫向應變。水溫改變對巖樣泊松比影響較小。高溫遇水循環(huán)次數(shù)增加，巖樣泊松比逐漸上升，循環(huán)3次時泊松比上升趨勢明顯，循環(huán)5次，比天然狀態(tài)下巖樣泊松比高0.08%。由此可見，增加注水循環(huán)次數(shù)與提高巖樣溫度對泊松比變化規(guī)律的影響效果相反。

圖9 不同注采參數(shù)下巖樣泊松比變化規(guī)律Fig.9 Changes of Poisson’s ratio of rock samples under different injection-production parameters

3.5 不同注采參數(shù)下巖樣力學及波動特征關(guān)系

通過所得結(jié)論能夠看出，熱采過程中不同注采參數(shù)下巖體力學及波動特征間存在關(guān)聯(lián)性，利用波動特征能夠很好的推斷熱采區(qū)巖體的力學特征。本文基于所得結(jié)論，對巖樣力學—波動特征曲線進行擬合。圖10（a-c）為不同注采參數(shù)下巖樣縱、橫波波速—峰值強度擬合曲線。

對巖樣縱、橫波波速—峰值強度進行回歸分析，得到下列方程：

巖樣溫度：

圖10 不同注采參數(shù)下巖樣縱、橫波波速-峰值強度擬合曲線Fig.10 Fitting curve of longitudinal and transverse wave velocity-compressive strength of rock samples under different injection-production parameters

根據(jù)式（5～10）看出，巖樣溫度不同，峰值強度和縱、橫波波速間存在較好關(guān)聯(lián)性，可決系數(shù)R2大于0.95；高溫遇水循環(huán)次數(shù)不同，關(guān)聯(lián)性次之，可決系數(shù)R2大于0.75；水溫不同，關(guān)聯(lián)性則較弱。這是因為多種因素會對試驗結(jié)果造成影響，既有巖體本身存在的差異，也有人為操作、試驗設備引起的誤差，各因素使試驗數(shù)據(jù)產(chǎn)生較大離散性。而巖樣溫度、高溫遇水循環(huán)次數(shù)是影響熱采區(qū)巖體力學及波動特征的主要因素，參數(shù)發(fā)生改變，巖體力學及波動特征差值較大，變化規(guī)律較為顯著；注水溫度對采熱過程中巖體力學及波動特征影響相對較弱，受客觀因素影響較大，但仍可體現(xiàn)出一定關(guān)聯(lián)性。因此通過波動特征對熱采區(qū)HDR峰值強度進行預估是一種可行的試驗方法。

4 結(jié)論

本文進行了不同注采參數(shù)下高溫遇水后巖樣物理力學及波動特征試驗研究，所得結(jié)論如下。

注采參數(shù)不同，擱置時間增加，巖樣物理力學及波動特征均產(chǎn)生一定的降幅。這主要與注水溫度、巖體內(nèi)部含水量、試驗條件有關(guān)。擱置初期，注水溫度相同，隨擱置時間增加，巖樣溫度越高，水流損失越多，縱、橫波波速、彈性模量降幅增大；巖樣溫度相同，注水溫度升高，縱、橫波波速、彈性模量受水溫影響較大，擱置初期降幅呈先增大后減小的趨勢。擱置后期受試驗條件影響較大，應考慮其對試驗結(jié)果產(chǎn)生的誤差?；谝陨弦蛩兀紤]擱置時間的作用，提高了對EGS注采、熱恢復過程中力學及波動特征預估的精確度。

提高巖樣溫度、增加高溫遇水循環(huán)次數(shù)，巖樣縱、橫波波速、峰值強度、彈性模量均逐漸下降。提升注水溫度對巖樣力學及波動特征變化規(guī)律的影響與其相反。由此可見，對熱采區(qū)巖體物理力學及波動特征影響由大到小的注采因素依次為靶區(qū)溫度、注水循環(huán)次數(shù)、注水溫度?；诖?，選取高溫儲層作為靶區(qū)鉆井，并采用低溫注水、多次冷水循環(huán)的方法，能夠達到巖層破碎增透儲層的目的。

影響力學及波動特征關(guān)聯(lián)性的注采參數(shù)依次為靶區(qū)溫度、注水循環(huán)次數(shù)、注水溫度，可以利用波動特征對熱采區(qū)巖體力學特征進行預測，節(jié)約實際工程中消耗的成本。