陳益民，程文超，周建旭，楊 飛，孫洪亮

(1.中國電建集團(tuán)華東勘測設(shè)計研究院有限公司，浙江 杭州 311122； 2.河海大學(xué)水利水電學(xué)院，江蘇 南京 210098)

水電站輸水發(fā)電系統(tǒng)的穩(wěn)定運(yùn)行是保障電網(wǎng)安全可靠的重要前提。在水電站運(yùn)行實(shí)踐中，水力系統(tǒng)、機(jī)械系統(tǒng)以及電氣系統(tǒng)的特性相互耦合，直接影響水電站的安全穩(wěn)定運(yùn)行[1]。因此，應(yīng)針對水電站輸水發(fā)電系統(tǒng)開展全面的穩(wěn)定性分析，對于兩機(jī)或多機(jī)共水力單元的情況，甩荷或增荷機(jī)組對運(yùn)行機(jī)組的影響較為突出，應(yīng)著重進(jìn)行水力干擾穩(wěn)定性分析。

圖1 水電站輸水發(fā)電系統(tǒng)布置示意圖

在對水電站水力-機(jī)械系統(tǒng)進(jìn)行穩(wěn)定性分析時，常采用剛性模型或低階彈性模型[2-3]來描述水體的動態(tài)特性，并將描述系統(tǒng)動態(tài)特性的微分方程線性化，得到狀態(tài)方程組，通過數(shù)值求解以評估系統(tǒng)的穩(wěn)定性和機(jī)組的調(diào)節(jié)品質(zhì)[4-6]，或者推導(dǎo)出相應(yīng)的特征方程或傳遞函數(shù)，進(jìn)行調(diào)節(jié)控制研究[7-9]。針對水電站穩(wěn)定性，諸多學(xué)者開展了較為廣泛而深入的研究。楊建東等[10]研究了上下游雙調(diào)壓室系統(tǒng)穩(wěn)定域的變化規(guī)律，推導(dǎo)出穩(wěn)定域干涉點(diǎn)和共振點(diǎn)的解析公式。郭文成等[11]建立了水輪機(jī)調(diào)節(jié)系統(tǒng)非線性數(shù)學(xué)模型，將Hopf分岔理論應(yīng)用到變頂高尾水洞水電站水輪機(jī)調(diào)節(jié)系統(tǒng)的穩(wěn)定性研究。Yu等[12]研究了水輪機(jī)-調(diào)速器系統(tǒng)狀態(tài)方程的表述特征，利用矩陣變換方法，建立了一種求解狀態(tài)方程系數(shù)矩陣的新方法。Zhou等[13]針對單管單機(jī)水電站，考慮調(diào)速器不同調(diào)節(jié)方式，建立了相應(yīng)的狀態(tài)方程，分析了小負(fù)荷擾動下的動態(tài)過程。翟曉娟等[14]考慮水力、機(jī)械、電氣系統(tǒng)的耦合作用，建立了水電機(jī)組穩(wěn)定性非線性數(shù)學(xué)模型，分析了孤網(wǎng)運(yùn)行時水電機(jī)組的水力干擾穩(wěn)定性。Chen等[15]提出了一種多目標(biāo)綜合指標(biāo)，并使用該指標(biāo)對調(diào)速器的參數(shù)進(jìn)行優(yōu)化，為在尾水波動條件下的穩(wěn)定性分析提供了理論支持。Yu等[16]利用Jordan標(biāo)準(zhǔn)形式的狀態(tài)矩陣解耦，推導(dǎo)了無量綱轉(zhuǎn)輪轉(zhuǎn)速的顯式公式，采用直接求解方法，優(yōu)化選擇PID控制器參數(shù)。馬安婷等[17]基于水力系統(tǒng)振動特性分析方法，建立輸水發(fā)電系統(tǒng)總體矩陣，能夠準(zhǔn)確求解復(fù)雜輸水發(fā)電系統(tǒng)的運(yùn)行穩(wěn)定域。周建旭等[18]結(jié)合設(shè)置上下游雙調(diào)壓室的長引水式水電站輸水系統(tǒng)的水力特性，探討了上下游調(diào)壓室系統(tǒng)水力干擾穩(wěn)定性等特殊水力學(xué)問題。Huang等[19]對孤立電網(wǎng)下復(fù)雜水電站小波動過渡過程中的調(diào)速器參數(shù)敏感性進(jìn)行了研究。李玲等[20]構(gòu)建了超長引水隧洞水電站單機(jī)無窮大電網(wǎng)的仿真模型，研究了系統(tǒng)的小波動穩(wěn)定性。上述研究主要側(cè)重于小波動穩(wěn)定性分析，而水力干擾穩(wěn)定性研究相對較少。

白鶴灘水電站單機(jī)容量達(dá)百萬千瓦，機(jī)組運(yùn)行水頭為163.9～243.1 m，具有單機(jī)容量巨大、運(yùn)行水頭變幅大的特點(diǎn)。該水電站輸水發(fā)電系統(tǒng)采用單管單機(jī)引水、兩機(jī)共用尾水調(diào)壓室和尾水隧洞的布置方式，其中尾水系統(tǒng)采用“三段變坡”的形式，包括有壓尾水隧洞、坡度較大的反坡銜接隧洞和平坡尾水隧洞，而平坡尾水隧洞中會發(fā)生明流、滿流，以及明滿流現(xiàn)象。因此，本文基于有壓輸水系統(tǒng)特征線法，引入明滿流尾水隧洞瞬變流分析的特征隱式格式，以便準(zhǔn)確描述尾水系統(tǒng)中可能的復(fù)雜流態(tài)，同時考慮調(diào)速器頻率調(diào)節(jié)模式，建立輸水發(fā)電系統(tǒng)穩(wěn)定性分析和調(diào)控模型，并結(jié)合白鶴灘水電站典型水力單元開展水力干擾穩(wěn)定性分析，闡明受擾機(jī)組的穩(wěn)定運(yùn)行及其影響因素，揭示尾水系統(tǒng)復(fù)雜流態(tài)對系統(tǒng)穩(wěn)定運(yùn)行的影響，為制定可靠有效的調(diào)節(jié)控制策略提供技術(shù)支撐。

1 輸水發(fā)電系統(tǒng)穩(wěn)定性分析模型

1.1 輸水系統(tǒng)瞬變流分析的數(shù)學(xué)模型

白鶴灘水電站輸水發(fā)電系統(tǒng)布置見圖1，其尾水系統(tǒng)的反坡銜接隧洞和平坡尾水隧洞在過渡過程中可能出現(xiàn)無壓流、有壓流，以及明滿流現(xiàn)象。圖1中1—1斷面為反坡銜接隧洞進(jìn)口斷面，2—2斷面為反坡銜接隧洞出口斷面，亦為明滿流隧洞段串聯(lián)節(jié)點(diǎn)。

輸水發(fā)電系統(tǒng)的有壓輸水部分，即從水庫進(jìn)水口至反坡銜接隧洞進(jìn)口(斷面1—1)，采用有壓輸水系統(tǒng)特征線法來描述其瞬態(tài)特性，可求解出計算斷面的測壓管水頭和流量，相應(yīng)的特征方程為

C+：HP=Cp-SpQP

(1)

C-：HP=Cm+SmQP

(2)

式中：H為斷面瞬時測壓管水頭；Q為斷面瞬時流量；Cp、Sp、Cm、Sm均為有壓輸水系統(tǒng)特征方程的已知系數(shù)；下標(biāo)P對應(yīng)的參數(shù)為計算節(jié)點(diǎn)未知量。

輸水發(fā)電系統(tǒng)的明滿流隧洞部分(運(yùn)行尾水位低于尾水隧洞出口洞頂或在尾水隧洞出口洞頂附近)，即從始終為有壓流的反坡銜接隧洞進(jìn)口斷面1—1至平坡尾水隧洞出口，采用明渠水流特征隱式格式來描述其動態(tài)特性，可求解出計算斷面的瞬時水深和流量，相應(yīng)的特征方程為

a1jΔhj-1+b1jΔQj-1+c1jΔhj+d1jΔQj=e1j

(3)

a2jΔhj+b2jΔQj+c2jΔhj+1+d2jΔQj+1=e2j

(4)

式中：Δh、ΔQ分別為斷面水深和流量的增量形式；aij、bij、cij、dij均為系數(shù)(i=1,2)；下標(biāo)j為計算斷面編號。

基于特征方程(3)和(4)，結(jié)合尾水閘門井節(jié)點(diǎn)的流量連續(xù)條件、水頭平衡條件和閘門井水位波動方程，可建立尾水閘門井邊界的數(shù)學(xué)模型；同時，考慮反坡銜接隧洞進(jìn)口斷面、明滿流隧洞串聯(lián)節(jié)點(diǎn)、明滿流隧洞出口的流量連續(xù)條件和水頭平衡條件，建立相應(yīng)邊界的數(shù)學(xué)模型。

有壓輸水和明滿流模擬段的交界面(即始終為有壓流的反坡銜接隧洞進(jìn)口斷面1—1)，滿足流量連續(xù)條件和水頭平衡條件(即有壓模擬段出口斷面流量等于明滿流模擬段進(jìn)口流量，以及有壓模擬段出口測壓管水頭等于明滿流模擬段進(jìn)口水深加上洞底高程)，據(jù)此可建立融合有壓輸水系統(tǒng)特征線法和明滿流隧洞特征隱式格式的輸水系統(tǒng)瞬變流分析模型。在瞬變流計算分析過程中，若明滿流模擬段某計算斷面的瞬時水深近似等于斷面洞高，則該計算斷面為明滿流分界面。

1.2 水輪發(fā)電機(jī)組的動態(tài)分析模型

白鶴灘水電站輸水發(fā)電系統(tǒng)的穩(wěn)定性分析模型應(yīng)能夠準(zhǔn)確描述百萬千瓦機(jī)組的綜合特性及其對系統(tǒng)運(yùn)行穩(wěn)定性的影響。白鶴灘水電站左岸輸水發(fā)電系統(tǒng)的水輪機(jī)型號為D545A-F15，結(jié)合機(jī)組模型特性曲線，可得到反映水輪機(jī)全特性的離散數(shù)據(jù)，包括流量特性Q11=f1(n11,a)和效率特性η=f2(n11,a)，其中Q11、n11、η和a分別為機(jī)組單位流量、單位轉(zhuǎn)速、效率和導(dǎo)葉開度，f1和f2表示不同的函數(shù)關(guān)系，見圖2。

圖2 水輪機(jī)D545A-F15綜合特性曲線

白鶴灘水電站右岸輸水發(fā)電系統(tǒng)的水輪機(jī)型號為A1181a，其全特性曲線見圖3，包括流量全特性曲線Q11=f1(n11,a)和力矩全特性曲線T11=f3(n11,a)，其中T11為機(jī)組單位力矩。

描述機(jī)組節(jié)點(diǎn)的方程有：

H=(Cp1-Cm2)-(Sp1+Sm2)Q

(5)

P=γQHη

(6)

(7)

(8)

(9)

式中：H、Q、P、n分別為機(jī)組運(yùn)行水頭、流量、出力、轉(zhuǎn)速；γ為水的容重；D1為轉(zhuǎn)輪直徑；Cp1、Sp1、Cm2和Sm2均為基于有壓輸水系統(tǒng)特征方程的已知系數(shù)。

結(jié)合式(5)～式(9)，以及水輪機(jī)模型綜合特性曲線Q11=f1(n11,a)和η=f2(n11,a)(左岸機(jī)組)，或者Q11=f1(n11,a)和T11=f3(n11,a)(右岸機(jī)組)，則可求解得到任一瞬時水輪機(jī)節(jié)點(diǎn)的相關(guān)參數(shù)，包括蝸殼進(jìn)口測壓管水頭和流量、尾水管進(jìn)口測壓管水頭和流量、機(jī)組瞬時出力或力矩等。

描述擾動機(jī)組動態(tài)特性的狀態(tài)方程為

(10)

式中：φ、p分別為機(jī)組相對轉(zhuǎn)速變化量和相對出力變化量；t為時間；x為外擾動；Sg為負(fù)荷自調(diào)節(jié)系數(shù)；Tm為機(jī)組啟動時間常數(shù);P0、n0、GD2分別為機(jī)組初始出力、初始轉(zhuǎn)速和轉(zhuǎn)動慣量?；谑?10)和機(jī)組的綜合特性，可得到受擾機(jī)組的瞬時轉(zhuǎn)速。

1.3 調(diào)速器的狀態(tài)方程

考慮到白鶴灘水電站在電網(wǎng)中的重要性和機(jī)組的運(yùn)行穩(wěn)定性要求，基于調(diào)速器的頻率調(diào)節(jié)模式，進(jìn)行系統(tǒng)的調(diào)節(jié)控制分析。在頻率調(diào)節(jié)模式下，調(diào)速器采用PID調(diào)節(jié)規(guī)律，將導(dǎo)葉開度作為反饋值，實(shí)時跟蹤機(jī)組實(shí)際功率。將調(diào)速器頻率調(diào)節(jié)模式的狀態(tài)方程[13]耦合機(jī)組的狀態(tài)方程(10)以及輸入擾動，即可確定機(jī)組和調(diào)速器的狀態(tài)變量值，包括機(jī)組導(dǎo)葉相對開度等。

1.4 融合輸水系統(tǒng)復(fù)雜瞬變流分析的輸水發(fā)電系統(tǒng)調(diào)控模型

融合輸水系統(tǒng)明滿流等復(fù)雜瞬變流分析的數(shù)學(xué)模型、考慮機(jī)組綜合特性的水輪發(fā)電機(jī)組動態(tài)分析模型和調(diào)速器頻率調(diào)節(jié)模式的狀態(tài)方程，構(gòu)建形成白鶴灘水電站百萬千瓦機(jī)組及輸水發(fā)電系統(tǒng)的調(diào)節(jié)控制模型。

表2 左岸1號水力單元水力干擾過渡過程調(diào)節(jié)品質(zhì)

基于所建立的白鶴灘水電站輸水發(fā)電系統(tǒng)調(diào)節(jié)控制模型，進(jìn)行水力-機(jī)械系統(tǒng)水力干擾穩(wěn)定性分析，能夠準(zhǔn)確反映擾動機(jī)組和輸水系統(tǒng)瞬態(tài)參數(shù)的較大幅度變化，以及受擾機(jī)組參數(shù)的動態(tài)變化過程和調(diào)節(jié)品質(zhì)，并直觀揭示輸水發(fā)電系統(tǒng)在水力耦合下的水力干擾特性，特別是受擾機(jī)組的動態(tài)特性。

2 輸水發(fā)電系統(tǒng)水力干擾穩(wěn)定性

2.1 白鶴灘水電站參數(shù)

白鶴灘水電站地下廠房采用首部開發(fā)方式，左、右岸輸水發(fā)電系統(tǒng)的布置格局相同，各布置8臺混流式水輪發(fā)電機(jī)組，分為4個水力單元，引水系統(tǒng)采用單洞單機(jī)布置方式，尾水系統(tǒng)采用兩機(jī)合一洞的布置方式，典型水力單元布置如圖1所示。電站左岸水輪機(jī)型號為D545A-F15，單機(jī)額定出力為1 015 MW，額定流量為545.49 m3/s，額定水頭為202 m，額定轉(zhuǎn)速為111.1 r/min，轉(zhuǎn)動慣量為36萬 t·m2；左岸1號水力單元引水道長度l1為480.50 m，尾水支洞長度l2為259.75 m，尾水隧洞總長度l3為1 695.80 m。右岸水輪機(jī)型號為A1181a，單機(jī)額定出力為1 015 MW，額定流量為538.80 為m3/s，額定水頭為202 m，額定轉(zhuǎn)速為107.1 r/min，轉(zhuǎn)動慣量為37萬t·m2；右岸8號水力單元引水道長度l1為486.12 m，尾水支洞長度l2為240.15 m，尾水隧洞總長度l3為1 769.87 m。左岸1號和右岸8號水力單元尾水隧洞出口頂高程均為592.0 m，尾水調(diào)壓室有效面積均為1 590.43 m2。

2.2 水力干擾分析工況

結(jié)合白鶴灘水電站典型水力單元，開展水力干擾過渡過程計算分析，評估發(fā)生水力干擾時運(yùn)行機(jī)組的運(yùn)行穩(wěn)定性以及調(diào)節(jié)品質(zhì)。已知調(diào)速器調(diào)節(jié)參數(shù)的整定值為：永態(tài)轉(zhuǎn)差系數(shù)bp=0.04，相應(yīng)比例常數(shù)KB=2.5、積分常數(shù)KI=0.385 1/s、微分常數(shù)KD=2.5 s。

表1給出了白鶴灘水電站左岸和右岸控制輸水發(fā)電系統(tǒng)水力干擾穩(wěn)定性分析的典型工況，其中工況D1和D2水庫水位一列分別對應(yīng)左岸1號和右岸8號水力單元的數(shù)值。

表1 水力干擾工況

2.3 受擾機(jī)組穩(wěn)定性和調(diào)節(jié)品質(zhì)分析

結(jié)合白鶴灘水電站左岸1號和右岸8號水力單元，以及相應(yīng)的機(jī)組特性和導(dǎo)葉啟閉規(guī)律，結(jié)合表1給出的水力干擾穩(wěn)定性分析工況，進(jìn)行受擾機(jī)組的運(yùn)行穩(wěn)定性分析，調(diào)節(jié)品質(zhì)分析見表2和表3，表中給出了對應(yīng)轉(zhuǎn)速偏差為±0.2%和±0.4%時的調(diào)節(jié)時間。典型工況頻率調(diào)節(jié)模式下受擾機(jī)組的轉(zhuǎn)速和機(jī)組出力的動態(tài)過程見圖4和圖5。

圖4和圖5表明：依據(jù)整定的調(diào)速器參數(shù)、電網(wǎng)負(fù)荷自調(diào)節(jié)系數(shù)等相關(guān)參數(shù)，左岸1號和右岸8號水力單元受擾機(jī)組轉(zhuǎn)速是衰減的，水力干擾過渡過程是穩(wěn)定的；受擾機(jī)組的出力發(fā)生明顯的擺動，最大出力均未超出額定出力的105%，且超出力歷時較短；受擾機(jī)組轉(zhuǎn)速的衰減度均小于90%。隨著機(jī)組運(yùn)行水頭的增大，受擾機(jī)組的調(diào)節(jié)品質(zhì)明顯改善，主要表現(xiàn)為調(diào)節(jié)時間縮短，振蕩次數(shù)減小，最大偏差減小和衰減度增大，即甩荷工況D3相較于D1，增荷工況D4相較于D2，機(jī)組運(yùn)行水頭較大，受擾機(jī)組轉(zhuǎn)速的波動明顯減小，衰減明顯加快。因工況D1運(yùn)行水頭明顯較低，機(jī)組轉(zhuǎn)速和出力動態(tài)過程振蕩幅度較大，穩(wěn)定性較差，同時工況D1對應(yīng)運(yùn)行開度為額定開度，受擾過程中受調(diào)速器接力器行程限制的影響，故工況D1和D3的水力干擾過渡過程存在較大差別。

表3 右岸8號水力單元水力干擾過渡過程調(diào)節(jié)品質(zhì)分析

圖4 左岸1號水力單元受擾機(jī)組轉(zhuǎn)速、出力動態(tài)過程

圖5 右岸8號水力單元受擾機(jī)組轉(zhuǎn)速、出力動態(tài)過程

2.4 尾水隧洞典型流態(tài)對系統(tǒng)穩(wěn)定性的影響

左岸1號和右岸8號水力單元尾水隧洞洞頂高程為592.00 m，工況D1和D2尾水位為597.42 m，工況D3和D4尾水位為581.50 m，因此工況D1和D2尾水系統(tǒng)全線為有壓流流態(tài)，而工況D3和D4在反坡銜接隧洞段出現(xiàn)明滿交替流流態(tài)，高程較高的平坡尾水隧洞段出現(xiàn)明流流態(tài)。圖6和圖7分別給出左岸1號水力單元各工況尾水調(diào)壓室水位和左/右岸典型水力單元明滿流分界面動態(tài)過程線，以分析尾水隧洞典型流態(tài)對系統(tǒng)穩(wěn)定性的影響。

圖6 尾水調(diào)壓室水位動態(tài)過程線

圖7 明滿流分界面動態(tài)過程線

圖6表明：機(jī)組額定出力運(yùn)行時，工況D1和D2運(yùn)行水頭相對較低，運(yùn)行流量較大，調(diào)壓室水位波動幅度較大，且增荷工況衰減較快；因工況D1和D2尾水隧洞全線有壓，調(diào)壓室水位波動周期較長，而工況D3和D4尾水隧洞下游段為明滿流段或明流段，有壓段長度縮短，調(diào)壓室水位波動周期較短，同時運(yùn)行水頭相對較高，調(diào)壓室水位穩(wěn)定時間明顯縮短，有利于系統(tǒng)的運(yùn)行穩(wěn)定性。圖7表明：工況D3和D4尾水系統(tǒng)中出現(xiàn)明滿流流態(tài)，對應(yīng)的明滿流分界面位于反坡銜接隧洞段，即在斷面1—1和斷面2—2之間移動，其下游側(cè)始終為無壓流，并且很快趨于穩(wěn)定。

2.5 機(jī)組型號對系統(tǒng)穩(wěn)定性的影響

白鶴灘水電站左岸1號水力單元的機(jī)組型號為D545A-F15，右岸8號水力單元的機(jī)組型號為A1181a。通過對比不同型號機(jī)組在相同工況下的機(jī)組轉(zhuǎn)速和反坡銜接隧洞出口斷面2—2測壓管水頭動態(tài)過程，分析不同機(jī)組型號時系統(tǒng)的水力干擾穩(wěn)定性。圖8和圖9分別給出了甩荷工況D1和D3下機(jī)組轉(zhuǎn)速變化相對量φ和反坡銜接隧洞出口(斷面2—2)測壓管水頭H2變化過程線。

圖8 工況D1下受擾機(jī)組轉(zhuǎn)速及斷面2—2測壓管水頭動態(tài)過程

圖9 工況D3下受擾機(jī)組轉(zhuǎn)速及斷面2—2測壓管水頭動態(tài)過程

圖8和圖9表明：在機(jī)組甩荷條件下，工況D1轉(zhuǎn)速變化相對值的動態(tài)變化趨勢基本一致，對應(yīng)的機(jī)組調(diào)節(jié)品質(zhì)評價指標(biāo)接近(表2、表3)，而工況D3轉(zhuǎn)速變化相對值衰減過程存在較明顯的差別，即不同機(jī)組型號下，第1個調(diào)節(jié)周期內(nèi)機(jī)組最大轉(zhuǎn)速不同和振蕩衰減趨勢明顯不一致，在調(diào)節(jié)過程中相關(guān)參數(shù)動態(tài)過程亦受不同的尾水隧洞明滿流過渡段和明流段的影響；左岸1號和右岸8號水力單元輸水系統(tǒng)特別是尾水系統(tǒng)布置存在一定的差異，其中8號水力單元尾水隧洞較長，表現(xiàn)出兩個水力單元銜接隧洞出口斷面2—2測壓管水頭的變化趨勢基本一致，其中工況D1銜接隧洞出口斷面為有壓狀態(tài)，呈周期性振蕩衰減，右岸8號水力單元對應(yīng)的振蕩周期略長，工況D3銜接隧洞出口斷面為無壓狀態(tài)，測壓管水頭快速衰減并趨于穩(wěn)定。

3 結(jié) 論

a.依據(jù)整定的調(diào)速器參數(shù)、電網(wǎng)負(fù)荷自調(diào)節(jié)系數(shù)等相關(guān)參數(shù)，考慮尾水系統(tǒng)多流態(tài)特性，白鶴灘水電站典型輸水發(fā)電系統(tǒng)水力干擾過渡過程是穩(wěn)定的，隨著機(jī)組運(yùn)行水頭的逐漸增大，受擾機(jī)組的調(diào)節(jié)品質(zhì)逐步改善。

b.機(jī)組運(yùn)行水頭相對較低且尾水系統(tǒng)全線有壓時，尾水系統(tǒng)銜接隧洞出口斷面測壓管水頭呈周期性振蕩衰減，而運(yùn)行水頭較高且尾水隧洞下游段為明滿流段或明流段時，尾水系統(tǒng)銜接隧洞出口斷面流態(tài)為明流，測壓管水頭呈快速趨于穩(wěn)定的動態(tài)過程。

c.在機(jī)組甩荷條件下，左岸1號和右岸8號水力單元不同機(jī)組型號對應(yīng)的轉(zhuǎn)速、尾水系統(tǒng)銜接隧洞出口斷面測壓管水頭的動態(tài)變化趨勢基本一致，對應(yīng)的機(jī)組調(diào)節(jié)品質(zhì)評價指標(biāo)接近，表明兩種機(jī)組型號對應(yīng)水力干擾穩(wěn)定性和受擾機(jī)組調(diào)節(jié)品質(zhì)接近。