蔡 超，李 楠，蘭學(xué)槿，秦艷麗，溫耀林，張園惠，王承淳，劉 鵬

(1.中國(guó)礦業(yè)大學(xué) 煤炭資源與安全開(kāi)采國(guó)家重點(diǎn)實(shí)驗(yàn)室，江蘇 徐州 221116；2.中國(guó)礦業(yè)大學(xué) 礦業(yè)工程學(xué)院，江蘇 徐州 221116；3.中國(guó)礦業(yè)大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院，江蘇 徐州 221116；4.中國(guó)礦業(yè)大學(xué) 安全工程學(xué)院，江蘇 徐州 221116；5.中國(guó)礦業(yè)大學(xué) 數(shù)學(xué)學(xué)院，江蘇 徐州 221116)

水力壓裂技術(shù)最早應(yīng)用于油氣及地?zé)衢_(kāi)發(fā)領(lǐng)域[1]。近年來(lái)，水力壓裂技術(shù)在煤礦堅(jiān)硬頂板控制、沖擊地壓和煤與瓦斯突出等煤巖動(dòng)力災(zāi)害防控、煤層瓦斯高效抽采等方面得到了廣泛應(yīng)用，對(duì)于煤礦安全高效開(kāi)采具有重要意義[2-4]。水力壓裂是通過(guò)向目標(biāo)煤巖層注入高壓流體，在水壓力作用下，煤巖體內(nèi)部產(chǎn)生壓裂裂縫，起到卸壓增透改性效果，從而達(dá)到煤巖動(dòng)力災(zāi)害防控和促進(jìn)瓦斯高效抽采的作用[5-7]。在水力壓裂過(guò)程中，研究裂縫的幾何形狀及裂縫的擴(kuò)展規(guī)律對(duì)壓裂工程師制定壓裂方案、評(píng)價(jià)壓裂效果和預(yù)測(cè)產(chǎn)氣量具有重要意義，同時(shí)也為現(xiàn)場(chǎng)壓裂鉆孔數(shù)量的調(diào)整、定向壓裂等提供依據(jù)。目前國(guó)內(nèi)外學(xué)者針對(duì)水力壓裂裂縫的擴(kuò)展，采用實(shí)驗(yàn)室實(shí)驗(yàn)、力學(xué)模型、數(shù)值模擬等方法進(jìn)行研究，然由于現(xiàn)場(chǎng)環(huán)境的復(fù)雜性，上述研究成果與真實(shí)的水力壓裂裂縫存在差異[8-10]。

微震監(jiān)測(cè)技術(shù)是指利用材料內(nèi)部失穩(wěn)破壞產(chǎn)生的震動(dòng)波信號(hào)來(lái)分析研究材料內(nèi)部結(jié)構(gòu)與穩(wěn)定性、破壞過(guò)程及機(jī)理的一種地球物理監(jiān)測(cè)手段[11，12]。1978年，Schuster C L.[13]將微震監(jiān)測(cè)技術(shù)引入石油領(lǐng)域，用于研究水力壓裂裂縫傳播的幾何形狀和激發(fā)范圍，并取得較好的效果。此后微震監(jiān)測(cè)技術(shù)被不斷的發(fā)展，產(chǎn)生一系列的理論成果，同時(shí)也逐漸商業(yè)化，產(chǎn)生較多的微震監(jiān)測(cè)公司[14]。美國(guó)的MicroSeismic (MSI)、南非的ISS、加拿大的ESG、道爾頓、中國(guó)的東方地球物理勘探等公司開(kāi)發(fā)出了具有自主知識(shí)產(chǎn)權(quán)的微震監(jiān)測(cè)系統(tǒng)軟件，能夠?qū)τ蜌馓锼毫蚜芽p的空間位置、長(zhǎng)度、高度、寬度、傾角等進(jìn)行定量計(jì)算，在油氣儲(chǔ)層改造壓裂監(jiān)測(cè)中得到了廣泛應(yīng)用，極大地推動(dòng)了非常規(guī)油氣田的勘探開(kāi)發(fā)進(jìn)程。但這些商用軟件大多數(shù)應(yīng)用于油氣開(kāi)采領(lǐng)域。與油氣開(kāi)采領(lǐng)域相比，煤礦井下水力壓裂的范圍和規(guī)模很小，產(chǎn)生的波形更微弱。利用這些軟件進(jìn)行微震數(shù)據(jù)處理效果不明顯，不能很好的拾取微弱波形的高精度到時(shí)，對(duì)裂縫定位效果較差。

因此，開(kāi)發(fā)一款適用于煤礦領(lǐng)域煤巖體水力壓裂裂縫微震監(jiān)測(cè)軟件已成為當(dāng)務(wù)之急。本文針對(duì)煤巖體水力壓裂微震信號(hào)高精度到時(shí)拾取、裂縫定位監(jiān)測(cè)難題，開(kāi)發(fā)了一款煤巖體水力壓裂裂縫微震監(jiān)測(cè)軟件，并通過(guò)開(kāi)展原煤真三軸水力壓裂實(shí)驗(yàn)，采集壓裂全過(guò)程數(shù)據(jù)，進(jìn)行到時(shí)拾取與震源定位，驗(yàn)證該軟件的可靠性及準(zhǔn)確性。

1 煤巖水力壓裂誘發(fā)信號(hào)特征及處理

1.1 煤巖水力壓裂誘發(fā)信號(hào)特征

煤巖體水力壓裂在較短時(shí)間內(nèi)會(huì)產(chǎn)生大量的波形信號(hào)，這些波形信號(hào)與傳統(tǒng)的礦震信號(hào)在幅值、頻率、信噪比等方面存在很大的差異，傳統(tǒng)的礦山微震監(jiān)測(cè)不能適用于煤巖體水力壓裂微震監(jiān)測(cè)。與致密砂巖、頁(yè)巖等致密巖層相比，煤層較為松軟、節(jié)理裂隙發(fā)育、非均質(zhì)性較強(qiáng)，煤巖體水力裂縫擴(kuò)展過(guò)程及水力壓裂誘發(fā)的微震信號(hào)具有獨(dú)特的特征。大量的研究表明，煤層水力壓裂誘發(fā)的信號(hào)具有低信噪比、高頻率、低幅值、能量低的特點(diǎn)[15，16]。

1.2 聲發(fā)射信號(hào)處理

聲發(fā)射波形的初至到時(shí)是震源定位必不可少的數(shù)據(jù)，其精度直接影響震源定位的準(zhǔn)確性[17]。目前波形到時(shí)拾取算法眾多，使用最為廣泛的是長(zhǎng)短時(shí)窗法(STA/LTA)。該方法是一種類似于信噪比的衡量方式，其中STA對(duì)時(shí)間序列振幅的快速變化比較敏感，而LTA則反映的是背景噪聲的信息[18]。當(dāng)兩者的比值大于預(yù)設(shè)的閾值時(shí)，判定為有效聲發(fā)射波形并記錄到時(shí)。但該方法在信號(hào)信噪比較低時(shí)，拾取誤差較大。此外，Akaike 信息準(zhǔn)則法(AIC)，也使用比較廣泛。AIC法是基于聲發(fā)射/微震信號(hào)是非平穩(wěn)的概念，利用自回歸過(guò)程，求取時(shí)窗內(nèi)聲發(fā)射/微震波形的AIC值[19]。AIC的求解過(guò)程包括三步：

1)通過(guò)式(1)計(jì)算背景噪音信號(hào)的閾值T：

T=m±3σ

(1)

式中，m、σ分別表示所選取噪音信號(hào)片段的平均值與標(biāo)準(zhǔn)差。

2)選擇AIC計(jì)算窗口，將微震/聲發(fā)射波形中，第一次超過(guò)背景噪音閾值T的點(diǎn)記為參考點(diǎn)，以參考點(diǎn)為起始點(diǎn)，向前、向后分別選取一定量的采樣點(diǎn)，構(gòu)成時(shí)間窗口。需要注意的是，該時(shí)間窗口內(nèi)需同時(shí)包含有效信號(hào)與背景噪音信號(hào)。

3)根據(jù)式(2)計(jì)算AIC(k)函數(shù)的最小值。

AIC(k)=klog(var{x(1，k)})+(N-k-1)

log(var{x(k+1，N)})

(2)

式中，k是輸入微震/聲發(fā)射波形的樣本數(shù)據(jù)；var{x}為方差函數(shù)。

在微震/聲發(fā)射波形記錄中，到時(shí)以前被認(rèn)為是背景噪聲信號(hào)。噪聲信號(hào)與有效微震波形信號(hào)具有不一樣的統(tǒng)計(jì)特性，在最小二乘法的意義下，兩者擬合效果最差，對(duì)應(yīng)的AIC值最小，該點(diǎn)是噪音信號(hào)和有效波形信號(hào)最佳的分界點(diǎn)，即到時(shí)點(diǎn)。

AIC法拾取波形到時(shí)較為準(zhǔn)確，但波形信噪比較低時(shí)，會(huì)計(jì)算出多個(gè)局部最小值，導(dǎo)致波形到時(shí)拾取的結(jié)果有明顯錯(cuò)誤。因此軟件在后續(xù)編程時(shí)，先采用STA/LTA法拾取微震/聲發(fā)射波形的大致到時(shí)點(diǎn)，在通過(guò)AIC法，獲取波形精確到時(shí)。此外，為克服AIC存在的問(wèn)題，增加到時(shí)驗(yàn)證。從AIC法拾取到時(shí)點(diǎn)開(kāi)始，在后續(xù)40個(gè)采樣點(diǎn)內(nèi)，聲發(fā)射波形的幅值至少存在4次大于背景噪聲信號(hào)閾值T的情況，且在AIC法拾取到時(shí)點(diǎn)位置后100個(gè)采樣點(diǎn)，聲發(fā)射波形幅值的符號(hào)至少變化5次。當(dāng)上述條件滿足時(shí)，認(rèn)為AIC法拾取的到時(shí)不存在較大偏差。

在已知波形高精度到時(shí)的基礎(chǔ)上，方可進(jìn)行震源定位。震源定位是指利用傳感器接收的震動(dòng)波形信號(hào)，通過(guò)定位方法確定震源的發(fā)震時(shí)刻和空間位置[20，21]。定位過(guò)程可以用式(3)來(lái)描述。

式中，ti為微震傳感器接收到信號(hào)的時(shí)間；t0為震源的起震時(shí)刻；(x，y，z)為震源坐標(biāo)，傳感器坐標(biāo)為(xi，yi，zi)；v為P波波速。

當(dāng)P波波速已知，式(3)中存在四個(gè)未知數(shù)，故震源定位需要不少于四個(gè)傳感器采集的波形信息。為了增加方程的可解性，通常使用較多的微震傳感器構(gòu)建超定方程組。然方程組求解過(guò)程中，可能會(huì)出現(xiàn)不收斂、數(shù)值不穩(wěn)定、甚至無(wú)解等問(wèn)題。單純形定位算法可以很好地避免這類問(wèn)題。單純形定位算法是一種多維搜索尋優(yōu)的迭代方法，具有迭代收斂和迭代過(guò)程快速穩(wěn)健的優(yōu)點(diǎn)。單純形定位算法的實(shí)施需要經(jīng)過(guò)四個(gè)步驟：①首先在三維空間中選取4個(gè)點(diǎn)構(gòu)成初始單純形；②其次，選取合適的目標(biāo)函數(shù)，將初始單純形的四個(gè)頂點(diǎn)作為破裂源，計(jì)算四種情況下目標(biāo)函數(shù)的最值；③根據(jù)拉伸、收縮、映射等基本變換準(zhǔn)則，不斷構(gòu)造新的單純形；④當(dāng)目標(biāo)函數(shù)滿足迭代終止條件時(shí)，此時(shí)使得目標(biāo)函數(shù)取得最優(yōu)解的點(diǎn)即為破裂點(diǎn)的最佳位置。因此，軟件采用單純形定位算法，求解破裂源的空間坐標(biāo)。為了能清晰地刻畫水力壓裂裂縫的空間形態(tài)，增加了定位事件點(diǎn)的篩選。當(dāng)聲發(fā)射傳感器采集的波形大于6個(gè)且傳感器采集的到時(shí)與計(jì)算到時(shí)之間的標(biāo)準(zhǔn)差在3μs之內(nèi)方可參與定位。

2 微震監(jiān)測(cè)軟件設(shè)計(jì)與開(kāi)發(fā)

2.1 軟件開(kāi)發(fā)的原則

軟件開(kāi)發(fā)過(guò)程時(shí)不應(yīng)違背客觀自然規(guī)律，同時(shí)采用簡(jiǎn)潔的計(jì)算機(jī)編程語(yǔ)言，保證軟件在運(yùn)行過(guò)程中的流暢性，給用戶帶來(lái)良好的體驗(yàn)。此外，還應(yīng)充分考慮煤巖體水力壓裂誘發(fā)信號(hào)的特征。由于煤礦生產(chǎn)環(huán)境的復(fù)雜性，微震數(shù)據(jù)通常包含大量的噪音信號(hào)。如何從監(jiān)測(cè)數(shù)據(jù)中準(zhǔn)確地識(shí)別水力壓裂誘發(fā)的微弱波形信號(hào)，是該軟件開(kāi)發(fā)需要解決的關(guān)鍵問(wèn)題?？煽啃允俏⒄鸨O(jiān)測(cè)系統(tǒng)正常工作的基礎(chǔ)。

2.2 軟件總體架構(gòu)設(shè)計(jì)

煤巖體水力壓裂微震監(jiān)測(cè)軟件通過(guò)調(diào)用儲(chǔ)存在微震采集儀上的原始波形，進(jìn)行數(shù)據(jù)分析與處理。微震信號(hào)的處理需經(jīng)過(guò)數(shù)據(jù)調(diào)用，濾波、波速模型建立、參數(shù)設(shè)置、P/S波高精度到時(shí)拾取、震源定位幾個(gè)過(guò)程。根據(jù)微震事件處理過(guò)程，軟件在開(kāi)發(fā)過(guò)程中，把需實(shí)現(xiàn)的功能模塊化。數(shù)據(jù)接口層由以太網(wǎng)電纜直接連接微震采集儀，軟件對(duì)原始波形數(shù)據(jù)掃描和讀取，并將原始波形傳遞至功能驅(qū)動(dòng)層，實(shí)現(xiàn)波形數(shù)據(jù)的進(jìn)一步處理分析。功能驅(qū)動(dòng)層需要實(shí)現(xiàn)微震監(jiān)測(cè)軟件的具體功能，它將用戶操作命令用程序代碼實(shí)現(xiàn)。運(yùn)行控制層是溝通用戶界面層與功能控制層的橋梁，它根據(jù)用戶的實(shí)際需求執(zhí)行相應(yīng)的功能模塊。用戶層界面是呈現(xiàn)在用戶電腦桌面上的界面窗口，用戶通過(guò)點(diǎn)擊界面功能按鈕，下發(fā)相應(yīng)的指令，實(shí)現(xiàn)對(duì)應(yīng)功能，如圖1所示。

2.3 軟件開(kāi)發(fā)

采用python語(yǔ)言與Visual Studio Code編輯器作為開(kāi)發(fā)工具，在已搭建軟件框架的基礎(chǔ)上，將上述波形到時(shí)拾取、震源定位算法及其他功能程序化，設(shè)計(jì)開(kāi)發(fā)了煤巖體水力壓裂裂縫微震監(jiān)測(cè)軟件。軟件由參數(shù)配置、到時(shí)拾取、震源定位、三維可視化等模塊組成。

軟件主界面左上方為菜單欄，由項(xiàng)目、幫助組成。項(xiàng)目菜單欄中包含新建項(xiàng)目與導(dǎo)入項(xiàng)目。打開(kāi)本軟件，導(dǎo)入從微震采集儀下載的數(shù)據(jù)文件，可獲得傳感器坐標(biāo)及原始波形信息，此時(shí)，軟件自動(dòng)顯示項(xiàng)目的基本信息，包括項(xiàng)目名、創(chuàng)始人、創(chuàng)建時(shí)間、項(xiàng)目開(kāi)始時(shí)間、項(xiàng)目結(jié)束時(shí)間、備注。此外也可通過(guò)軟件主界面左側(cè)的項(xiàng)目管理，手動(dòng)填寫項(xiàng)目的基本信息。軟件主界面右側(cè)顯示傳感器的信息與參數(shù)配置。點(diǎn)擊參數(shù)配置按鈕，可彈出參數(shù)配置界面。參數(shù)配置菜單欄包含AIC時(shí)窗長(zhǎng)度、檢驗(yàn)時(shí)窗長(zhǎng)度、到時(shí)驗(yàn)證時(shí)窗長(zhǎng)度、檢波器參與個(gè)數(shù)、長(zhǎng)短時(shí)窗、速度模型參數(shù)設(shè)置功能。點(diǎn)擊參數(shù)配置，可手動(dòng)輸入層號(hào)、層厚、P波速度、S波速度等參數(shù)。參數(shù)設(shè)置完畢后，點(diǎn)擊軟件主界面左側(cè)的數(shù)據(jù)處理可實(shí)現(xiàn)微弱波形高精度到時(shí)拾取。在已知微震傳感器坐標(biāo)的基礎(chǔ)上，點(diǎn)擊震源定位按鈕可獲取震源的三維空間坐標(biāo)，并將其單獨(dú)顯示在軟件界面右側(cè)。

圖1 軟件總體架構(gòu)

3 煤巖水力壓裂裂縫監(jiān)測(cè)軟件驗(yàn)證

3.1 實(shí)驗(yàn)系統(tǒng)及實(shí)驗(yàn)方案

采用真三軸水力壓裂系統(tǒng)，如圖2所示，開(kāi)展煤樣水力壓裂實(shí)驗(yàn)，采集水力壓裂誘發(fā)的信號(hào)，并通過(guò)煤巖水力壓裂裂縫監(jiān)測(cè)軟件進(jìn)行波形到時(shí)拾取與高精度定位。

圖2 實(shí)驗(yàn)系統(tǒng)

實(shí)驗(yàn)采用的煤樣選自山東協(xié)莊煤礦，采用風(fēng)鉆和切割機(jī)將煤樣制備成長(zhǎng)寬高均為150mm的試樣。利用空心鉆在原煤試樣上表面中心位置處，垂直鉆取深度90mm、直徑14mm的壓裂孔，成孔后采用植筋膠將研制的封孔器密封在壓裂孔中，封孔器插入壓裂孔中的深度為60mm，剩下30mm的裸孔段作為壓裂段。煤樣表面共布設(shè)12個(gè)R15傳感器，如圖3所示，圖中黑色圓點(diǎn)表示聲發(fā)射傳感器。

圖3 傳感器布置及部分坐標(biāo)

在進(jìn)行水力壓裂實(shí)驗(yàn)之前，為了驗(yàn)證軟件的可靠性，首先利用已知震源空間坐標(biāo)的標(biāo)定實(shí)驗(yàn)進(jìn)行分析，即在原煤試樣表面上的已知位置處開(kāi)展斷鉛實(shí)驗(yàn)，采集斷鉛產(chǎn)生的聲發(fā)射波形，利用軟件拾取各聲發(fā)射傳感器位置處到時(shí)信息，并通過(guò)軟件進(jìn)行震源定位。震源定位的結(jié)果及誤差見(jiàn)表1。定位誤差的計(jì)算公式如式(4)。

式中，(x，y，z)為震源坐標(biāo)；傳感器坐標(biāo)為(xi，yi，zi)；N為傳感器的數(shù)量；煤樣的P波波速為2.0km/s。

由表1可知，在所有斷鉛位置，通過(guò)軟件求解的震源坐標(biāo)與實(shí)際斷鉛的誤差皆在2mm之內(nèi)，這證明了本軟件波形到時(shí)拾取與破裂源定位算法的準(zhǔn)確性及可靠性。

表1 破裂源定位的結(jié)果及誤差

3.2 到時(shí)拾取及震源定位驗(yàn)證

根據(jù)0.02MPa/s的加載速率進(jìn)行圍壓加載，待三個(gè)圍壓達(dá)到預(yù)設(shè)值后保持恒載。恒載一段時(shí)間后，按照100mL/min的注水流量進(jìn)行注水壓裂，直到煤樣破壞。將煤樣水力壓裂全過(guò)程采集的數(shù)據(jù)導(dǎo)入至本軟件中，進(jìn)行到時(shí)拾取及震源定位驗(yàn)證。

煤樣水力壓裂全過(guò)程曲線，如圖4所示。采用美國(guó)物理聲學(xué)軟件進(jìn)行聲發(fā)射波形到時(shí)拾取如圖5所示，聲發(fā)射波形到時(shí)存在明顯偏差，圖中紅色豎線表示聲發(fā)射波形的到時(shí)線。由于煤樣水力壓裂誘發(fā)的聲發(fā)射波形比較微弱，信噪比低，P波起跳不明顯，故單純采用長(zhǎng)短時(shí)窗法不能獲取精確到時(shí)。為驗(yàn)證本軟件到時(shí)拾取的可靠性，將煤樣水力壓裂誘發(fā)的波形導(dǎo)入本軟件，并進(jìn)行到時(shí)拾取，如圖6所示，該軟件到時(shí)拾取較為精確，圖中綠色豎線表示聲發(fā)射波形的到時(shí)線。

圖4 煤樣水力壓裂全過(guò)程曲線

圖5 美國(guó)物理聲學(xué)軟件聲發(fā)射波形到時(shí)

圖6 煤巖水力壓裂裂縫監(jiān)測(cè)軟件聲發(fā)射到時(shí)

煤樣水力壓裂試驗(yàn)結(jié)束后，試樣的破裂形態(tài)如圖7所示。美國(guó)物理聲學(xué)軟件定位的結(jié)果如圖8所示，圖8中有890個(gè)聲發(fā)射事件參與定位，定位事件點(diǎn)散亂的分布在空間中，無(wú)法辨別水力壓裂裂縫的空間分布形態(tài)。煤巖水力壓裂裂縫監(jiān)測(cè)軟件定位的結(jié)果如圖9所示，該軟件定位的過(guò)程中，增加了定位事件點(diǎn)的篩選，共有327個(gè)聲發(fā)射事件參與定位。由于部分聲發(fā)射事件的波形數(shù)不足6個(gè)或傳感器采集的到時(shí)與計(jì)算到時(shí)之間的標(biāo)準(zhǔn)差大于3μs，故超過(guò)一半的定位點(diǎn)沒(méi)有顯現(xiàn)。水力壓裂裂縫呈現(xiàn)平行最大主應(yīng)力(σ1)的空間分布特征，且裂縫在擴(kuò)展時(shí)存在明顯的偏轉(zhuǎn)現(xiàn)象，如圖9所示，圖中，黑色圓點(diǎn)為定位點(diǎn)在XZ面的投影，綠色圓點(diǎn)為定位點(diǎn)在YZ面的投影，藍(lán)色圓點(diǎn)為定位點(diǎn)在XY面的投影。將圖8與圖9兩者對(duì)比，可清晰凸顯出該軟件的可靠性與先進(jìn)性。

圖7 煤樣水力壓裂破壞實(shí)物

圖8 美國(guó)物理聲學(xué)軟件定位的結(jié)果

圖9 煤巖水力壓裂裂縫監(jiān)測(cè)軟件定位的結(jié)果

4 結(jié) 論

1)研究提出了STA/LTA與AIC聯(lián)合到時(shí)拾取算法，并增加了波形到時(shí)驗(yàn)證的功能，克服了AIC法在信噪比較低情況下，存在較多個(gè)極小值的問(wèn)題。此外，在破裂源定位方面，添加了定位事件點(diǎn)篩選功能，突出水力壓裂裂縫的空間分布特征。

2)利用python言語(yǔ)與Visual Studio Code，開(kāi)發(fā)了一款煤巖體水力壓裂裂縫微震監(jiān)測(cè)軟件。軟件從波形到時(shí)拾取與震源定位算法出發(fā)，基于微震數(shù)據(jù)處理過(guò)程，搭建了軟件總體框架；在此基礎(chǔ)上，采用功能模塊化設(shè)計(jì)，實(shí)現(xiàn)了項(xiàng)目管理、參數(shù)配置、數(shù)據(jù)調(diào)用、數(shù)據(jù)分析處理等功能。

3)開(kāi)展原煤真三軸水力壓裂實(shí)驗(yàn)，將聲發(fā)射采集的數(shù)據(jù)作為軟件輸入數(shù)據(jù)，進(jìn)行聲發(fā)射波形到時(shí)拾取與破裂源定位。測(cè)試結(jié)果表明：本軟件可精確拾取煤巖水力壓裂誘發(fā)微弱波形的到時(shí)，且震源定位的結(jié)果可以很好刻畫煤巖水力壓裂主裂縫的空間分布形態(tài)，這為煤礦井下煤巖體水力壓裂裂縫監(jiān)測(cè)及壓裂效果評(píng)價(jià)提供了技術(shù)手段。