李元松，吳春華，段鑫，李洋，董開源

1.武漢工程大學(xué)資源土木工程學(xué)院，湖北 武漢 430074；2.黃石市礦山安全檢測檢驗所，湖北黃石 435000

大志山注漿帷幕防治水工程大型抽水試驗

李元松1，吳春華2，段鑫1，李洋1，董開源2

1.武漢工程大學(xué)資源土木工程學(xué)院，湖北 武漢 430074；2.黃石市礦山安全檢測檢驗所，湖北黃石 435000

礦山帷幕注漿防治水工程具有典型隱蔽工程的特點，其施工質(zhì)量檢查鑒定和堵水效果評價一直沒有可依照和遵循的標(biāo)準(zhǔn)，現(xiàn)有的單因素評價或多因素混合評價都難以給出定量結(jié)論，大型抽水試驗是目前最為直接且有效的方法，但因工程復(fù)雜費用昂貴，很少使用.本文基于Dupuit穩(wěn)定井流基本原理，在充分調(diào)查、分析礦區(qū)生產(chǎn)現(xiàn)狀、水文地質(zhì)條件、帷幕勘察、設(shè)計、施工資料的基礎(chǔ)之上，設(shè)計并成功實施一大型抽水試驗，使用單井抽水多個觀測孔計算公式求解注漿后復(fù)合巖體的等效滲透系數(shù)，采用擬“大井法”對礦井深部的涌水量進行預(yù)測，其成果直接用于大志山帷幕注漿工程的質(zhì)量檢測與評價鑒定，取得較為滿意的效果.

巖溶礦床；地下水防治；帷幕注漿；穩(wěn)定井流；抽水試驗

0 引言

隨著礦產(chǎn)資源的開發(fā)逐漸向深部發(fā)展，帷幕注漿堵水方法越來越多受到廣泛的應(yīng)用.我國大多數(shù)富水巖溶礦床的開采，均或多或少使用了帷幕注漿技術(shù)進行地下涌水的封堵與阻隔［1－4］.然而，帷幕注漿工程具有典型的隱蔽工程特點，其工程施工質(zhì)量檢查和效果評價仍然沒有可依照與遵循的標(biāo)準(zhǔn).利用既有主井進行抽水試驗是目前檢驗帷幕注漿效果較為直接且有效的方法，根據(jù)抽水試驗取得的綜合成果，可采用地下水動力學(xué)法或涌水量曲線方程法預(yù)測礦坑涌水量、判定帷幕堵水率，并利用幕內(nèi)外水位觀測孔水位動態(tài)資料，預(yù)測地下水位降落漏斗擴展情況及地面巖溶塌陷可能的擴展范圍［5－7］.

本文根據(jù)大冶市大志山注漿帷幕防治水工程的特點，結(jié)合多個溶巖礦山帷幕堵水工程質(zhì)量檢測評定實踐，設(shè)計并實施大型主井抽水試驗.系統(tǒng)分析試驗工程背景、過程及成果，以供類似注漿工程質(zhì)量檢測鑒定時參考.

1 工程背景

1.1 礦山基本情況

大冶市大志山礦區(qū)為一典型的矽卡巖型巖溶富水礦床，礦體產(chǎn)于接觸帶外矽卡巖中，其產(chǎn)狀基本與接觸帶一致.礦區(qū)水文地質(zhì)條件復(fù)雜，曾進行過開采，因突水淹井而停產(chǎn).之前礦山完成基建工程：a．主井掘砌安裝工程總深501 m；b．風(fēng)井掘到－300 m標(biāo)高，總深350 m；c．－300 m中段已掘石門平巷125 m及沿脈平巷360 m，主井與風(fēng)井貫通；d．－350 m中段掘石門平巷135 m，沿脈巷180 m，水泵房及水倉形成21 312 m3/d的排水能力，已投入使用；e．－400 m中段已掘石門平巷135 m，沿脈平巷140 m，并有天井與－350 m中段聯(lián)通；f．－450 m中段掘石門平巷165 m，沿脈巷160 m，水泵房及水倉開掘安裝完成并投入使用，形成16 281 m3／d的排水能力；有天井與－400 m中段貫通.為開采利用銅資源，礦山?jīng)Q定恢復(fù)開采，經(jīng)防治水技術(shù)研究，決定實施帷幕注漿，并于2014年11月底全部竣工.

根據(jù)“安全第一，預(yù)防為主，綜合治理”的安全生產(chǎn)方針，判斷帷幕構(gòu)筑后井下開采是否安全，檢驗注漿帷幕的連續(xù)性和整體截流能力，明確是否具備安全驗收條件，為此設(shè)計并實施本次抽水試驗.

1.2 礦區(qū)水文地質(zhì)條件

礦區(qū)北部以碳酸鹽含水層同數(shù)條構(gòu)造斷裂以及坑道組成了一個大志山礦區(qū)導(dǎo)水和積水的網(wǎng)絡(luò)，給礦體開采造成了重大隱患，南部以巖漿巖體形成了隔水邊界.礦區(qū)水文地質(zhì)條件屬復(fù)雜類型巖溶充水礦床.

a．含水層.1）第四系殘坡積粘土夾碎石含水層；沖積、洪積粉質(zhì)粘土夾砂石孔隙含水層；2）下三疊系大冶灰?guī)r、嘉陵江灰?guī)r巖溶裂隙含水層；3）下二疊系茅口灰?guī)r溶洞裂隙含水層，為礦區(qū)內(nèi)主要含水層，受構(gòu)造斷裂帶的影響，發(fā)育深度達到-500 m，平均滲透系數(shù)2.886 9 m/d；4）下二疊系棲霞灰?guī)r中等富水溶洞裂隙含水層；5）石炭系黃龍、船山灰?guī)r中等富水溶洞裂隙含水層.

b．隔水層.1）石英閃長巖隔水層，分布在矽卡巖下盤；2）志留系砂頁巖隔水層.由于分布規(guī)模的局限性，僅具局部隔水作用.

c．構(gòu)造含水帶.1）斷裂含水帶，礦區(qū)內(nèi)斷裂構(gòu)造以NWW一組較為發(fā)育，其中較大斷裂F1斷層涌水量達25 L／s；F2斷層涌水量11.46～14.03 L／s；F4斷層起關(guān)鍵的導(dǎo)水作用，發(fā)育深度大，連通性好，將地表水體和淺層地下水緊密連通起來，向礦坑補給；F5推測斷層：自水南灣村西緣沿河谷北西向延展；F6推測斷層：放水中水位有明顯反映，為推測導(dǎo)水?dāng)鄬樱?）接觸含水帶，茅口灰?guī)r和大冶灰?guī)r與陽新侵入體之邊緣相石英閃長巖相接觸，形成礦區(qū)構(gòu)造接觸帶.富水性僅次于大理巖溶洞裂隙含水層，坑道中揭露涌水量1.0～10 L／s.

d．水文地質(zhì)邊界條件.礦區(qū)南側(cè)已被陽新閃長巖體侵吞，東起葉家莊，西至曹家灣方向分布，基本阻隔了南面的基巖地下水.帷幕線兩端勘察孔均進入該閃長巖隔水邊界；志留系頁巖分布于礦區(qū)北部，受構(gòu)造擠壓及風(fēng)化作用影響，淺部巖石風(fēng)化破碎，距礦坑1 km范圍內(nèi)基本沒有完整的隔水體存在，礦區(qū)北東部為敞開邊界.

e．地下水天然補給、徑流與排泄.天然條件下，分布在低山丘陵地帶的裸露灰?guī)r，為本區(qū)地下水的主要補給區(qū)，南北方向展布的水南灣河谷是地下水徑流區(qū)，北部的大冶湖盆低地是地下水排泄區(qū).

f．采區(qū)充水因素.1）礦坑充水主要來自二疊系茅口灰?guī)r溶洞裂隙含水層；2）接觸構(gòu)造的充水，由蝕變石英閃長巖、蝕變大理巖、矽卡巖組成，受構(gòu)造作用影響，裂隙發(fā)育，巖心破碎，地下水沿構(gòu)造裂隙流出；3）礦區(qū)斷層的存在，特別是斷層與礦體斜交，是引起本礦床充水主要因素之一；4）越流與地表水對礦坑的充水，沖洪積亞粘土夾砂礫石含水層，主要分布在礦區(qū)東部水南灣河床一帶.

g．帷幕軸線巖層特征.帷幕軸線上巖溶發(fā)育規(guī)律與巖石透水率見表1.

1.3 帷幕設(shè)計參數(shù)

a．帷幕型式：平面封閉式，剖面半封底式；

b．帷幕長度：1 616.8 m；

表1 帷幕軸線巖溶發(fā)育規(guī)律與巖石透水率Table1 Karst development law and rock permeable rate on grouting curtain axis

c．帷幕深度（帷幕底邊界標(biāo)高）：幕底深度沿軸線呈V字型，幕底最深-568.25 m.

d．帷幕高度：幕墻起始標(biāo)高為+18.46 m，幕底最深標(biāo)高-568.25 m，帷幕高度介于97.14～601.53 m之間；

e．帷幕厚度：10 m；

f．鉆孔孔徑：開孔孔徑Φ130～150 mm，終孔孔徑Φ91 mm；

g．鉆孔偏斜率：終孔最大偏斜率≤孔深的1.5%；

h．鉆孔間距：10 m；

i．注漿擴散半徑：7.1 m；

j．帷幕滲透系數(shù)：目標(biāo)值≤0.09 m／d.

1.4 帷幕施工

a．注漿材料：采用水泥尾砂漿、水泥粘土漿和水泥尾砂粘土漿.

b．漿液配制：漿液配制采用集中制漿站攪拌配制，用泵經(jīng)過輸漿管道送到各注漿孔的二次攪拌桶，再用注漿泵向鉆孔內(nèi)注漿.

c．注漿方式：采用孔口封閉，孔內(nèi)循環(huán)，自上而下分段注漿.為保證注漿質(zhì)量，加快施工進度，采用φ42鉆桿作為射漿管，射漿管的下置深度進入受注段底部距孔底0.5～5.0 m，注漿過程的觀測和記錄采用灌漿自動記錄儀.

d．分段高度：注漿段長原則上為30 m，在溶洞發(fā)育區(qū)，注漿段長5～10 m.在遇到大于1 m的溶洞時，停鉆注漿.

e．壓水試驗.試驗采用孔口封閉或孔內(nèi)膠塞止水封閉，單點式壓水，壓力表安裝在孔口，試驗壓力取灌漿壓力的80%.

f．注漿壓力：注漿終壓采用靜水壓力的1.5～2.0倍.

2 抽水試驗

2.1 試驗方案

本次抽水試驗采用既有主井抽水，礦山原進行過開采，形成一定的空區(qū)，原有采空區(qū)有一定的靜儲量.抽水試驗的目的是通過試驗求得礦井各階段的涌水量，抽水試驗設(shè)計要分出靜儲量與動涌水量，與相關(guān)抽水試驗規(guī)范的條件不完全一致.解決這一問題有兩種方法：控制循環(huán)時間法與控制水位法.考慮大志山礦井抽水降深大（－350 m），水中泥砂含量高，試驗條件困難，本次試驗采用控制水位法.即假定某一降深段單位時間的涌水量為常數(shù)q，控制井筒水位在一定范圍內(nèi)進行抽水、恢復(fù)水位循環(huán)，經(jīng)過若干循環(huán)后，達到每一循環(huán)的抽水量與循環(huán)時間趨于相同時，視為穩(wěn)定.

2.2 試驗步驟

抽水試驗步驟如下：

a．一次抽水至－350 m水平，停機，記錄時間，抽水量；

b．水位恢復(fù)至－300 m，開機，記錄時間，水位降至－350 m水平，停機，記錄時間，抽水量；

c．水位恢復(fù)至－300 m時，開機，記錄時間，水位降至－350 m水平，停機，記錄時間，抽水量；

d．如此循環(huán)，直到每一循環(huán)單位時間抽水量連續(xù)5次相差不到5%時，視為滲流穩(wěn)定，可停止抽水.

e．恢復(fù)水位至－252 m水平，按上述步驟進行第二降深段抽水試驗.

抽水試驗各時段的流量觀測，主井與觀測孔水位觀測同步進行，試驗期安排專職人員進行地表塌陷，主要建（構(gòu)）物的變形觀測與巡查.

2.3 儀器設(shè)備及試驗過程

2014年度10月30日10：00開始坑內(nèi)試抽水，主井采用兩臺深水潛水泵，型號分別為HN QB1300－550／470和HNQB560－270／460.采用電測水位計觀測水位，TDS－100F1超聲波流量計計量.平均每天排水約12 000 m3，日均降深約2 m，至12月27日凌晨4點水位降至－350 m.隨后停機恢復(fù)水位至－300 m，27日凌晨5點開始正式循環(huán)抽水試驗.經(jīng)5個循環(huán)，單位時間的涌水量基本穩(wěn)定，折算日排水量約9 200 m3／d，隨后停機恢復(fù)水位至－252 m，至28日上午10點開始第二降深段抽水試驗，經(jīng)6個循環(huán)，單位時間的涌水量基本穩(wěn)定，至12年30日上午8時兩階段抽水試驗圓滿結(jié)束.

3 試驗結(jié)果分析

a．抽水量－循環(huán)－時間關(guān)系

將抽水試驗結(jié)果整理，繪制抽水量－循環(huán)－歷時曲線如圖1所示.

圖1 抽水量－時程曲線Fig.1 Water output-time history curve

從圖1中可以看出，兩降深段曲線相似，均從第二循環(huán)后近似于平行水平軸的直線，抽水量與循環(huán)時間趨于平穩(wěn)，符合設(shè)計穩(wěn)定流試驗要求.－350降深段循環(huán)5.5 h，抽水量約2 115 m3，折算成日排水量為9 257 m3/d；－300降深段循環(huán)時間4.05 h，抽水量為1 360 m3，折算成日排水量為8 047 m3／d.

b．觀測孔的降深

將抽水試驗過程中各觀測點水位觀測數(shù)據(jù)，整理繪制降深－歷時曲線如圖2所示.從圖2中可以看出，從10月30日開始至12月30日止，帷幕內(nèi)外側(cè)觀測孔水位變化特點不同，幕內(nèi)觀測孔（觀9，觀11）降深較大，與主井水位聯(lián)動性明顯.幕外各觀測孔（觀1，觀3，觀4，觀6，觀7，觀13）水位降深小，且與主井水位聯(lián)動性差.

c．Q－S，q－s曲線

對于穩(wěn)定流抽水試驗，除及時繪制出Q－t和s－t曲線外，尚需繪制出Q－s和q－s關(guān)系曲線（q為單位降深涌水量）.Q－s、q－s曲線可幫助了解抽水試驗進行是否正常.常見的Q－s和q－s曲線類型如圖3和圖4所示.圖3中：

1）曲線Ⅰ表示承壓井流；

2）曲線Ⅱ表示潛水或承壓轉(zhuǎn)無壓井流；

3）曲線Ⅲ表示從某一降深值起，涌水量隨降深的加大而增加很少；

4）曲線Ⅳ表示補給衰竭或水流受阻，隨s加大Q反而減少；

5）曲線Ⅴ表示涌水量隨降深增加而加速增大，異常.

由于各種原因，試驗只進行兩個降深，根據(jù)本次抽水試驗結(jié)果，結(jié)合先前－160 m的涌水記錄，繪制Q－S曲線與q－s曲線如圖5和圖6.從圖5可以看出，試驗Q－S比較符合圖3中曲線Ⅱ－潛水或承壓轉(zhuǎn)無壓井流的特征，這與事先設(shè)計的抽水試驗特征相符.對比圖4和圖6，可以看出本次試驗所得q－s曲線符合曲線Ⅱ的特征.

圖2 主井與部分觀測孔降深－歷時曲線Fig.2 Main shaft and the observation hole drawdown-time curve

圖3 常見抽水試驗的Q－S曲線類型Fig.3 Q－S curve in common type of pumping test

圖4 常見抽水試驗的q－s曲線類型Fig.4 q－s curve common type of pumping test

圖5 本次抽水試驗的Q-S曲線Fig.5 Q-S curve in this pumping test

圖6 本次抽水試驗的q-s曲線Fig.6 q-s curve in this pumping test

4 涌水量計算

目前，礦山帷幕注漿施工后分階段涌水量預(yù)測還沒有統(tǒng)一的理論計算方法，也無相應(yīng)規(guī)范可以參照.實際工程中往往以水文地質(zhì)學(xué)中關(guān)于均質(zhì)水體的滲透理論為基礎(chǔ)進行計算，然后根據(jù)工程經(jīng)驗進行適當(dāng)修正.常用的方法有曲線方程法、大井法、相似比擬法和數(shù)值模擬法.本次采用大井法進行深部兩階段的涌水量預(yù)測.

據(jù)礦區(qū)水文地質(zhì)條件，礦區(qū)北邊為大理巖含水層分布區(qū)，南邊為閃長巖隔水巖體，可近似看作為一直線隔水邊界.第四系土層上部有一層粘土，為隔水層.粘土下部為洪積粉質(zhì)粘土夾砂石孔隙含水層，滲透系數(shù)采用承壓潛水非完整井公式，如圖7所示.

圖7 承壓潛水非完整井計算示意圖Fig.7 Computing sketch of partially penetrating well in phreatic aquifer under pressure

a．等效滲透系數(shù)計算

式（1）中：Q－涌水量（m3/d）；K為滲透系數(shù)（m/d）；Sw為主井的水位降深（m）；S1為觀測孔的水位降深（m）；Ha為靜止水位至有效帶下限的水柱高度（m）；r0為大井引用半徑（m）；r1為主井距觀測孔的距離（m）；β為大理巖地下水進入礦坑范圍與礦坑（大井）周長的比值.

計算參數(shù)的選擇與確定：

1）靜止地下水位標(biāo)高：取抽水試驗前靜止水位標(biāo)高18．46 m.

2）含水層厚度（Ha）：有效帶下限取鉆孔揭露的溶洞，綜合考慮?。?00 m，故含水層總厚度（Ha）為518．46 m；

3）大井引用半徑r0：大井引用半徑r0：分別?。?00 m，－350 m等水位線圈定的范圍折算出大井引用半徑r0分別為20 m和25 m.鉆孔剖面揭露的粉質(zhì)粘土最低點標(biāo)高為0．66 m；

4）β大理巖地下水進入礦坑范圍與礦坑大井周長的比值（近似直線隔水邊界）為0．5.

將上述參數(shù)與實測結(jié)果代入公式（1），計算－300 m水平和－350水平的等效滲透系數(shù)，其結(jié)果如表2所示.

表2 -300/-350水平等效滲透系數(shù)K值計算表Table 2 Equivalent permeability coefficient K value calculation at-300/-350 level

預(yù)測中段－400 m和－450 m的等效滲透系數(shù)按線性插值法近似計算，分別為0．086 m/d和0．085 m/d.

b．涌水量計算

為了說明帷幕注漿以后，幕內(nèi)礦坑排水量與帷幕注漿前比較其堵水效果、堵水率，采用大井法估算－400 m（標(biāo)高）和－450 m的涌水量，仍采用承壓－潛水非完整井涌水量計算公式：

式（2）中：

Ma為含水層頂板至有效帶下限的厚度（m）；

ha為動水位至有效帶下限的水柱高度（m）；

其余參數(shù)同式（1）.

根據(jù)礦山水文地質(zhì)確定的參數(shù)與實測水位降深代入式（2），計算結(jié)果如表3所示.

表3 地下水動力學(xué)法礦坑涌水量計算結(jié)果表Table 3 Calculation result of mine water inflow by groundwater dynamics method

根據(jù)表3計算結(jié)果，可算出礦坑水位降至相應(yīng)中段的堵水率，見表4.

表4 分中段堵水率表Table 4 Water plugging rate in typical level

從以上分析計算可知，帷幕堵水率均超過70%的設(shè)計目標(biāo)，最高達到75．5%，最低也達到72．1%，從而判定注漿效果明顯，達到設(shè)計要求.

5 結(jié)語

a．主井抽水試驗規(guī)模大，涉及環(huán)節(jié)多，成本高，管理復(fù)雜.本次試驗方案既遵循了相關(guān)規(guī)范規(guī)程，取得數(shù)據(jù)科學(xué)有效，又考慮了礦山的實際情況，因地制宜，切實可行，達到試驗?zāi)康?

b．抽水試驗與“大井法”計算結(jié)果表明，大志山礦－300、－350、－400和－450四個中段的堵水率分別為75．5%、72．3%、72．1%和72．7%.

致謝

感謝大冶市大志山礦業(yè)有限司提供項目經(jīng)費支持，同時感謝中南基礎(chǔ)工程公司的全力配合.

［1］高建軍，祝瑞勤，徐大寬．巖溶充水礦床帷幕注漿堵水技術(shù)研究［J］．水文地質(zhì)工程地質(zhì)，2007，34（5）：123－127.

GAO Jianjun，ZHU Ruiqin，XU Dakuan．Re－discussion on plugging technique with curtain grouting in karst water－filled ore deposits［J］．Hydrogeology＆EngineeringGeology，2007，34（5）：123－127．（in Chinese）

［2］祝世平，王伏春，曾夏生．大紅山礦帷幕注漿治水工程及其評價［J］．金屬礦山，2007（9）：79－83.

ZHUShiping，WANGFuchun，ZENGXiasheng．Groundwater blockage by grouting curtain in dahongshanm ine and its evaluation［J］．Metal Mine，2007（9）：79－83.（in Chinese）

［3］韓貴雷，于同超，劉殿鳳，等．礦山帷幕注漿方案研究及堵水效果綜合分析［J］．礦業(yè)研究與開發(fā)，2010，30（3）：95－98.

HAN Guilei，YU Tongchao，LIU Dianfeng，et al．Study on curtain grouting scheme and comprehensive analysis on its water plugging effection in Zhongguan ironmine［J］．Mining R＆D，2010，30（3）：95－98.（in Chinese）

［4］陳靜，王軍．巖溶礦區(qū)帷幕注漿截流效果的抽水試驗法評價［J］．采礦技術(shù)，2008，8（6）：53－55.

CHEN Jing，WANG Jun．The effect evaluation of water prevention after curtain grouting engineering in karst mine area using pumping test［J］．Mining Technology，2008，8（6）：53－55.（in Chinese）

［5］康鳳新，魏東，張新文，等．大型抽水試驗的水文地質(zhì)意義［J］．水文地質(zhì)工程地質(zhì)，2005（5）：27－30.

KANG Fengxin，WEI Dong，ZHANG Xinwen，et al．Hydrogeological functions of a large pumping test［J］．Hydrogeology＆Engineering Geology，2005（5）：27－30.（in Chinese）

［6］薛禹群，吳吉春．地下水動力學(xué)［M］．3版.北京：地質(zhì)出版社，2010．

XUE Y Q，WU J C．Groundwater dynamics［M］．3rd ed. Beijing：The Geological Publishing House，2010．（in Chinese）

［7］雷國輝，葛國昌，雷國剛，等．承壓完整井非穩(wěn)定流抽水試驗的觀測時間要求［J］．水文地質(zhì)工程地質(zhì)，2009，36（1）：40－42.

LEI G H，GE G C，LEI G G，et al．Requirement of survey time of unsteady flow pumping test of complete bearing［J］．Hydrogeology＆Engineering Geology，2009，36（1）：40－42．（in Chinese）

Steady－flow pumping test of groundwater blockage by grouting curtain in Dazhishan mine

LI Yuan－song1，WU Chun－h(huán)ai2，DUAN Xin1，LI Yang1，DONG Kai－yuang2
1．School of Resource and Civil Engineering，Wuhan Institute of Technology，Wuhan 430074，China；2．Huangshi Institute of Mine Safety Inspection，Wuhan 435000，China

The grouting curtain，with the typical characteristics of underground engineering，has no specific standard to inspect its construction quality and evaluate the effect of grouting，and the quantitative conclusion is difficult to be given by single factor or multiple factors combination，so the large pumping test，which is rarely be implemented in practice because of its expensive and complicated characteristics，could be the most direct and effective method．Based on the basic principle of Dupuit’s steady well flow，a large pumping test was designed and completed after the sufficient investigation and analysis of production status，hydrogeological conditions and grouting curtain document of investigation，design and construction about Dazhishan mine．In the test，the equivalent permeability coefficient of composite rock body after grouting curtain was solved by using single well pumping-multiple observation hole calculation formula，and the amount of water in deep was forecasted using“virtual large diameter well”method，which was used in quality evaluation and appraisal in grouting curtain water plugging engineering in Dazhishan mine，and the results are satisfactory．

Karst water－filled deposit；groundwater blockage；grouting curtain；steady－flow；pumping test

TU 443

A

10.3969/j.issn.1674-2869.2015.08.002

1674－2869（2015）08－0006－07

本文編輯：龔曉寧

2015－04－30

李元松（1964－），男，湖北應(yīng)城人，博士，教授.研究方向：巖土工程領(lǐng)域的教學(xué)與研究.