馬武剛
1 中國(guó)地震局地震研究所,武漢市洪山側(cè)路40號(hào),430071 2 中國(guó)地震局地震大地測(cè)量重點(diǎn)實(shí)驗(yàn)室,武漢市洪山側(cè)路40號(hào),430071 3 湖北省地震局,武漢市洪山側(cè)路48號(hào),430071
VP型寬頻帶垂直擺傾斜儀是用來(lái)記錄由天體相互作用力引起的地傾斜固體潮變化的地震觀(guān)測(cè)設(shè)備,其基本原理是利用垂直于地面的懸掛擺系來(lái)拾取地傾斜變化,由高精度電容式微位移傳感器將這種變化轉(zhuǎn)變?yōu)槲⑷蹼娦盘?hào),通過(guò)放大、整形、濾波等環(huán)節(jié),由數(shù)采記錄、整理,還原出實(shí)際的地傾斜變化狀況的過(guò)程。研究表明[1-2],秒采樣gPhone重力儀可以準(zhǔn)確記錄到地震波的同震響應(yīng)信號(hào),且記錄的P波信號(hào)也能真實(shí)反映震源破裂過(guò)程。VP寬頻帶傾斜儀自2012年投入觀(guān)測(cè)以來(lái),將形變觀(guān)測(cè)儀器采樣率從min提升到s,采樣率的大幅提升使其能記錄到更多的高頻信息。與地震計(jì)、秒采樣gPhone重力儀等類(lèi)似,該儀器也能夠觀(guān)測(cè)到地震的同震波、地脈動(dòng)等豐富的地球物理信息[3]。
本文利用河北地區(qū)8套VP型寬頻帶垂直擺傾斜儀的連續(xù)觀(guān)測(cè)資料,從映震特征、映震能力等方面討論傾斜儀(儀器參數(shù)修改前后)的映震效能,并對(duì)其時(shí)頻響應(yīng)特征進(jìn)行分析,為VP型寬頻帶垂直擺傾斜儀連續(xù)觀(guān)測(cè)數(shù)據(jù)在同震信號(hào)分析中的應(yīng)用和儀器的參數(shù)修改提供科學(xué)參考。
河北地區(qū)共有12個(gè)定點(diǎn)形變臺(tái)站,其中8個(gè)臺(tái)站布設(shè)有VP型寬頻帶垂直擺傾斜儀,包括懷來(lái)臺(tái)、張家口臺(tái)、易縣臺(tái)、承德臺(tái)、赤城臺(tái)、豐寧臺(tái)、寬城臺(tái)及崇禮臺(tái),主要分布在黑山寺與狼山斷裂帶、北東向蔚縣山前大斷裂帶、桑干河斷裂帶、東西向狼窩溝-崇禮-赤城斷裂帶和懷安-宣化斷裂帶、紫荊關(guān)斷裂帶上。除崇禮地震臺(tái)外,VP型寬頻帶垂直擺傾斜儀的儀器硬件技術(shù)參數(shù)均一致。2015年,武漢地震科學(xué)儀器研究院在對(duì)崇禮地震臺(tái)VP型寬頻帶垂直擺傾斜儀進(jìn)行設(shè)備質(zhì)量自查后,為解決調(diào)零頻繁和頻帶范圍高頻端不足的問(wèn)題,將儀器的輸出信號(hào)放大倍數(shù)和低通濾波電路參數(shù)進(jìn)行調(diào)整,其中輸出信號(hào)放大倍數(shù)降低至原來(lái)的1/8,低通濾波電路濾波常數(shù)降低至原來(lái)的1/2,之后儀器于2020-09入網(wǎng)產(chǎn)出觀(guān)測(cè)數(shù)據(jù),具體觀(guān)測(cè)效果正在實(shí)際觀(guān)測(cè)中進(jìn)行檢驗(yàn)。
自觀(guān)測(cè)以來(lái),除張家口臺(tái)儀器NS分量由于儀器原因觀(guān)測(cè)曲線(xiàn)存在干擾、赤城臺(tái)儀器NS分量存在不明原因畸變外,其他各臺(tái)儀器觀(guān)測(cè)資料均穩(wěn)定、可靠且連續(xù)性好。為方便對(duì)各臺(tái)觀(guān)測(cè)資料進(jìn)行比對(duì),在本文后續(xù)分析過(guò)程中,選取各臺(tái)站VP型寬頻帶垂直擺傾斜儀EW分量2020-09至今的秒采樣數(shù)據(jù)進(jìn)行處理和分析。
VP型寬頻帶垂直擺傾斜儀觀(guān)測(cè)數(shù)據(jù)通常包括潮汐和非潮汐2個(gè)部分,其中非潮汐部分包含有地震波信號(hào)。考慮到固體潮信號(hào)非常強(qiáng),在對(duì)地震波信號(hào)進(jìn)行分析時(shí),為更好地識(shí)別地震波信息,需要利用小波分解去除固體潮及長(zhǎng)趨勢(shì)等影響。此外,原始觀(guān)測(cè)資料中還會(huì)有突跳、間斷等錯(cuò)誤記錄和儀器自身的漂移,這些影響也應(yīng)一并去除,最后才能得到地震波信號(hào)占主要成分的非潮汐觀(guān)測(cè)序列。圖1為記錄到地震的VP型寬頻帶垂直擺傾斜儀原始數(shù)據(jù)和經(jīng)過(guò)預(yù)處理后的觀(guān)測(cè)序列,可以看出,去除固體潮及長(zhǎng)趨勢(shì)等信號(hào)的影響之后,地震波信息更為突顯。
設(shè)信號(hào)為x(j),j=1,2,…,n,則信號(hào)的功率譜密度(PSD)為:
(1)
式中,Δt為采樣間隔,N為采樣點(diǎn)總數(shù)。一定頻率范圍內(nèi)的平均功率譜密度為:
(2)
式中,k1和k2分別為所選取頻率范圍的上限和下限,M為所選取頻率范圍內(nèi)的采樣點(diǎn)數(shù)。
3.1.1 近震映震特征
據(jù)中國(guó)地震臺(tái)網(wǎng)測(cè)定,北京時(shí)間2021-04-16 16:06河北灤州(39.75°N,118.72°E)發(fā)生M4.3地震,震源深度9 km。河北省內(nèi)8套VP型寬頻帶垂直擺傾斜儀均記錄到本次地震,震中距為96~344 km,其中離震中最近的是寬城臺(tái),最遠(yuǎn)的是張家口臺(tái)。
圖2(a)為各臺(tái)站經(jīng)預(yù)處理之后的震后觀(guān)測(cè)序列,考慮到近震同震持續(xù)時(shí)間較短,這里只選取震后0.5 h的數(shù)據(jù),由于大部分臺(tái)站處理結(jié)果一致,本文僅給出部分臺(tái)站結(jié)果;圖2(b)為承德臺(tái)同址觀(guān)測(cè)的寬頻帶垂直擺傾斜儀和寬頻帶地震儀的震后觀(guān)測(cè)序列曲線(xiàn)對(duì)比??梢钥闯觯瑥牡卣鸢l(fā)生時(shí)刻16:06開(kāi)始,各臺(tái)觀(guān)測(cè)數(shù)據(jù)陸續(xù)出現(xiàn)小幅突跳,1 min左右出現(xiàn)最大脈沖,隨后逐漸減小并衰減至震前水平。各臺(tái)觀(guān)測(cè)曲線(xiàn)在地震前后趨勢(shì)一致,未出現(xiàn)漂移及大的階躍現(xiàn)象。除崇禮臺(tái)外,相較于寬頻帶地震儀,其他各臺(tái)儀器均未記錄到較完整的地震波信息,無(wú)法獲得P波、S波等震相,而崇禮臺(tái)記錄的同震信號(hào)則較為完整。這是由于崇禮臺(tái)低通濾波電路參數(shù)有改變,濾波常數(shù)降低,高頻成分就會(huì)大量保留,對(duì)于近震而言,地震波幅度不會(huì)衰減太多。但值得注意的是,崇禮臺(tái)震中距并非最遠(yuǎn),但該臺(tái)儀器在地震發(fā)生2 min之后才出現(xiàn)同震響應(yīng),這一現(xiàn)象除了與地震信號(hào)傳輸介質(zhì)相關(guān)外,其他成因還需作進(jìn)一步探尋。
圖2 2021-04-16河北寬頻帶傾斜儀和寬頻帶地震儀觀(guān)測(cè)序列Fig.2 Observation sequence of Hebei broadband tiltmeter and broadband seismograph on April 16, 2021
3.1.2 遠(yuǎn)震映震特征
據(jù)中國(guó)地震臺(tái)網(wǎng)測(cè)定,北京時(shí)間2021-01-15 02:28蘇拉威西(西里伯斯)島(3.04°S,118.99°E)發(fā)生M6.2地震,震源深度10 km。圖3(a)為各臺(tái)站經(jīng)預(yù)處理的震后1 h觀(guān)測(cè)序列,震中距為4 727~4 926 km;圖3(b) 為承德臺(tái)同址觀(guān)測(cè)的寬頻帶垂直擺傾斜儀和寬頻帶地震儀的震后觀(guān)測(cè)序列對(duì)比曲線(xiàn)。河北省內(nèi)8套VP型寬頻帶垂直擺傾斜儀均明顯地記錄到了本次地震的同震響應(yīng),本文只給出具有代表性的結(jié)果。由圖3可以看出,幾套儀器的映震特征較為一致,表現(xiàn)為地震波初至?xí)r曲線(xiàn)小幅震蕩,面波到達(dá)時(shí)振幅達(dá)到最大值,之后振幅逐漸減小、衰減至震前水平;且P波與S波的起始時(shí)刻非常清晰,S波振幅大于P波振幅,其中崇禮臺(tái)記錄到的地震波信號(hào)中P波和S波信號(hào)成分占比更大、能量更強(qiáng)。由圖3(b)可以看出,承德臺(tái)寬頻帶垂直擺傾斜儀EW分量與寬頻帶地震儀E分量記錄到的遠(yuǎn)震地震波波形基本吻合,兩者具有較高的一致性。
圖3 2021-01-15河北寬頻帶傾斜儀和寬頻帶地震儀觀(guān)測(cè)序列Fig.3 Observation sequence of Hebei broadband tiltmeter and broadband seismograph on January 15, 2021
本文擬從震級(jí)、震源深度及震中距3個(gè)方面對(duì)河北地區(qū)VP型寬頻帶垂直擺傾斜儀的映震能力進(jìn)行分析。對(duì)于VP型寬頻帶垂直擺傾斜儀來(lái)說(shuō),遠(yuǎn)震的地震波信息能夠記錄得比較完整,因此本文共選取了6個(gè)遠(yuǎn)震事件作為研究對(duì)象。
3.2.1 震級(jí)
為分析不同震級(jí)對(duì)河北地區(qū)VP型寬頻帶垂直擺傾斜儀映震能力的影響,本文選取發(fā)震地點(diǎn)相近的2個(gè)不同震級(jí)的地震事件,基本參數(shù)見(jiàn)表1。2次地震震中位置最遠(yuǎn)相距約39 km,河北地區(qū)8套VP型寬頻帶垂直擺傾斜儀的震中距分別為10 236~10 496 km和10 222~10 481 km,遠(yuǎn)超2次地震之間的距離,可將其視為同一震源。河北地區(qū)8套VP型寬頻帶垂直擺傾斜儀對(duì)這2次地震的同震響應(yīng)結(jié)果見(jiàn)表2及圖4,由于篇幅限制,本文僅展示部分結(jié)果??梢钥闯?,相比于新西蘭克馬德克群島M5.5地震,各儀器記錄的新西蘭克馬德克群島M7.8地震的地震波振幅要大得多,而對(duì)于新西蘭克馬德克群島M5.5地震而言,張家口臺(tái)和豐寧臺(tái)儀器由于儀器自身或環(huán)境背景噪聲略大,因此記錄到的地震波信號(hào)較不明顯;崇禮臺(tái)由于調(diào)整了儀器的電路參數(shù),導(dǎo)致儀器記錄曲線(xiàn)噪聲較大,在震級(jí)較小時(shí)地震波信號(hào)基本被噪聲湮沒(méi)。
表1 地震基本參數(shù)
表2 2次地震同震響應(yīng)參數(shù)
圖4 河北寬頻帶傾斜儀觀(guān)測(cè)序列Fig.4 Observation sequence of the broadband tiltmeter in Hebei
綜上所述,地震震級(jí)對(duì)VP型寬頻帶垂直擺傾斜儀映震能力的影響是顯而易見(jiàn)的,震級(jí)越大,儀器記錄到的地震波信號(hào)越強(qiáng)。
3.2.2 震源深度
為分析震源深度對(duì)河北地區(qū)VP型寬頻帶垂直擺傾斜儀映震能力的影響,選取發(fā)震地點(diǎn)相近、震級(jí)相同的2個(gè)地震,基本參數(shù)見(jiàn)表3。2021-01-08 08:28與2021-01-06 14:20新西蘭克馬德克群島海域發(fā)生2次M6.1地震,震源深度分別為200 km和10 km,河北地區(qū)8套VP型寬頻帶垂直擺傾斜儀的震中距分別為10 156~10 414 km和10 175~10 436 km,遠(yuǎn)超2次地震之間的距離(156 km),可視為發(fā)震地點(diǎn)相同。河北地區(qū)8套VP型寬頻帶垂直擺傾斜儀對(duì)這2次地震的同震響應(yīng)的部分結(jié)果見(jiàn)表4及圖5??梢钥闯?,除崇禮臺(tái)外,不管是從最大振幅還是同震持續(xù)時(shí)間上來(lái)看,其他7套VP型寬頻帶垂直擺傾斜儀均對(duì)淺源地震的映震效能更加突出;而崇禮臺(tái)VP型寬頻帶垂直擺傾斜儀則恰好相反,對(duì)震源深度較深的地震有更好的映震表現(xiàn),這一現(xiàn)象可能是由于儀器的頻帶寬度和輸出信號(hào)放大倍數(shù)改變?cè)斐傻摹?/p>
表3 地震基本參數(shù)
表4 2次地震同震響應(yīng)參數(shù)
圖5 河北寬頻帶傾斜儀觀(guān)測(cè)序列Fig.5 Observation sequence of the broadband tiltmeter in Hebei
3.2.3 震中距
為分析不同震中距的地震對(duì)河北地區(qū)VP型寬頻帶垂直擺傾斜儀映震能力的影響,選取2個(gè)震級(jí)和震源深度相同但發(fā)震地點(diǎn)不同的地震,基本參數(shù)見(jiàn)表5。其中,2020-12-21 01:23日本本州東岸近海的M6.3地震距離河北較近,震中距為2 040~2 334 km;2021-02-16 08:49瓦努阿圖群島的M6.3地震較遠(yuǎn),震中距為8 241~8 495 km。河北地區(qū)8套VP型寬頻帶垂直擺傾斜儀對(duì)這2次地震的部分同震響應(yīng)結(jié)果見(jiàn)表6及圖6。可以看出,就同震最大振幅而言,8套寬頻帶垂直擺傾斜儀對(duì)震中距最近的日本本州東岸近海M6.3地震的波幅均大于瓦努阿圖群島的M6.3地震;就同震持續(xù)時(shí)間而言,8套VP型寬頻帶垂直擺傾斜儀表現(xiàn)一致:同震持續(xù)時(shí)間與震中距呈現(xiàn)出負(fù)相關(guān),即對(duì)震中距較遠(yuǎn)的瓦努阿圖群島M6.3地震的同震持續(xù)時(shí)間長(zhǎng),對(duì)震中距較近的日本本州東岸近海M6.3地震的同震持續(xù)時(shí)間短。這是由于,隨著震中距增加,P波和S波的周期也明顯增大,同震持續(xù)時(shí)間也越長(zhǎng),記錄的震相也越完整。
表6 2次地震同震響應(yīng)參數(shù)
圖6 河北寬頻帶傾斜儀觀(guān)測(cè)序列Fig.6 Observation sequence of the broadband tiltmeter in Hebei
為分析VP型寬頻帶垂直擺傾斜儀記錄地震波信號(hào)的頻譜特征,以2021-04-16 16:06河北灤州M4.3地震(震中距96~344 km)及2021-01-15 02:28蘇拉威西(西里伯斯)島M6.2地震(震中距4 727~4 926 km)為例,選擇懷來(lái)、赤城、崇禮3個(gè)臺(tái)站記錄的震后1 h數(shù)據(jù)進(jìn)行分析。
圖7為懷來(lái)、赤城、崇禮3個(gè)臺(tái)站VP型寬頻帶垂直擺傾斜儀河北灤州M4.3地震震后1 h數(shù)據(jù)的波形圖、時(shí)頻譜和功率譜。由于采樣率為1 s,根據(jù)采樣定理,寬頻帶垂直擺傾斜儀只能記錄到0.5 Hz以?xún)?nèi)的信號(hào)和能量。由懷來(lái)臺(tái)和赤城臺(tái)時(shí)頻譜結(jié)果可以看出,地震發(fā)生時(shí),信號(hào)能量在0~0.5 Hz之間均有增強(qiáng),特別是0~0.2 Hz之間,能量增強(qiáng)尤為明顯,功率譜圖也顯示出同樣結(jié)果;而崇禮臺(tái)時(shí)頻譜和功率譜圖結(jié)果與前兩臺(tái)完全不同:信號(hào)能量集中在較高頻段的0.2~0.5 Hz之間。這主要是由于廠(chǎng)家對(duì)崇禮臺(tái)VP型寬頻帶垂直擺傾斜儀的低通濾波電路參數(shù)進(jìn)行調(diào)整,使得儀器頻帶寬度的高頻端得到擴(kuò)展,所以崇禮臺(tái)能記錄到更多的高頻信號(hào),故信號(hào)能量在高頻端較為集中。
圖7 河北灤州M4.3地震震后1 h波形、時(shí)頻譜和功率譜Fig.7 1-hour waveform, time spectrum and power spectrum after Luanzhou M4.3 earthquake in Hebei province
圖8為懷來(lái)、赤城、崇禮3個(gè)臺(tái)站VP型寬頻帶垂直擺傾斜儀蘇拉威西(西里伯斯)島M6.3地震震后1 h數(shù)據(jù)的波形、時(shí)頻譜和功率譜。由時(shí)頻譜結(jié)果可以看出,P波到時(shí),信號(hào)能量開(kāi)始有增強(qiáng)跡象,頻域較寬,能量集中在0.1~0.3 Hz之間;S波到時(shí),信號(hào)能量有明顯增強(qiáng),能量集中頻段較P波稍低,集中在0~0.2 Hz之間;面波在P波后十幾min到達(dá),強(qiáng)度遠(yuǎn)高于其他信號(hào),能量在整個(gè)頻域內(nèi)最強(qiáng),頻域隨后收窄,集中在0~0.1 Hz之間,而后逐漸減弱。懷來(lái)、赤城的功率譜圖也顯示出本次地震VP型寬頻帶垂直擺傾斜儀所記錄到的地震波能量集中在0~0.1 Hz之間;崇禮臺(tái)功率譜圖結(jié)果與前兩臺(tái)稍有不同,信號(hào)能量集中的頻帶更寬,在0.1~0.3 Hz之間。
圖8 蘇拉威西(西里伯斯)島M6.2地震震后1 h波形、時(shí)頻譜和功率譜Fig.8 1-hour waveform, time spectrum and power spectrum after the M6.2 earthquake in Sulawesi (Celebes)
本文利用河北地區(qū)8套VP型寬頻帶垂直擺傾斜儀的連續(xù)觀(guān)測(cè)資料,從映震特征、映震能力等方面討論河北地區(qū)VP型寬頻帶垂直擺傾斜儀的映震效能,并對(duì)其時(shí)頻響應(yīng)特征進(jìn)行分析,結(jié)果表明:
1)VP型寬頻帶垂直擺傾斜儀能夠記錄到較為完整的遠(yuǎn)震地震波信息,從波形上能夠看到清晰的P波與S波起始時(shí)刻;對(duì)于近震而言,除崇禮臺(tái)外,其他各儀器均未記錄到較完整的地震波信息,只能夠記錄到面波,無(wú)法獲得P波、S波等震相,而崇禮臺(tái)儀器由于低通濾波電路參數(shù)改變,使得記錄的同震信號(hào)較為完整,但同震響應(yīng)存在延時(shí)現(xiàn)象,其原因還需進(jìn)一步探尋。
2)從6個(gè)震例的同震響應(yīng)結(jié)果可以初步推斷,在映震能力方面張家口臺(tái)VP型寬頻帶垂直擺傾斜儀的映震能力最優(yōu),赤城和易縣臺(tái)儀器映震能力相對(duì)較弱,更全面的結(jié)果還需更多的震例來(lái)驗(yàn)證。此外,從震級(jí)、震源深度和震中距3個(gè)方面對(duì)河北地區(qū)VP型寬頻帶垂直擺傾斜儀映震能力進(jìn)行分析,結(jié)果表明:
①地震震級(jí)和震中距對(duì)VP型寬頻帶垂直擺傾斜儀映震能力的影響較為明顯,震級(jí)或震中距越大,儀器記錄到的地震波信號(hào)越強(qiáng);
②除崇禮臺(tái)外,不管是從最大振幅還是同震持續(xù)時(shí)間上來(lái)看,其他臺(tái)站7套VP型寬頻帶垂直擺傾斜儀均對(duì)淺源地震的映震效能更加突出,而崇禮臺(tái)VP型寬頻帶垂直擺傾斜儀則恰好相反;
③對(duì)于同一地震,各臺(tái)站的震中距遠(yuǎn)近與VP型寬頻帶垂直擺傾斜儀的映震能力沒(méi)有表現(xiàn)出明顯的相關(guān)性。震中距遠(yuǎn),映震能力卻不一定差。這可能跟臺(tái)站所處環(huán)境和地下介質(zhì)密度差異有關(guān):由于地下介質(zhì)密度不均,因此即使震級(jí)和震源深度相同,地震波傳播方向不同,也會(huì)造成地表振動(dòng)強(qiáng)度不一致。
3)VP型寬頻帶垂直擺傾斜儀記錄地震波信號(hào)的頻譜特征分析結(jié)果顯示,該儀器能夠記錄到較為完整的面波信號(hào),但由于受到采樣率的限制,體波(P波和S波)信號(hào)記錄不全,無(wú)法記錄到0.5 Hz以上的信號(hào)。但目前研究結(jié)果顯示,地震計(jì)記錄的P波的優(yōu)勢(shì)頻帶在0.5 Hz以?xún)?nèi)[2],因此認(rèn)為,VP型寬頻帶垂直擺傾斜儀能夠記錄到大部分的P波信號(hào)。
4)崇禮臺(tái)VP型寬頻帶垂直擺傾斜儀對(duì)淺源地震的映震效能產(chǎn)生了與其他儀器相反的結(jié)論,并且對(duì)震級(jí)較小地震的映震能力明顯降低。產(chǎn)生這一現(xiàn)象的原因可能是由于儀器廠(chǎng)家在2018年對(duì)儀器參數(shù)進(jìn)行了調(diào)整,將儀器電路放大倍數(shù)降低至原來(lái)的1/8,低通濾波電路時(shí)間常數(shù)降低至原來(lái)的1/2。其主要目的是通過(guò)降低放大倍數(shù)擴(kuò)大儀器的量程,在相同零漂的情況下減少調(diào)零次數(shù);通過(guò)減小濾波電路常數(shù)使得高頻截止頻率變大,以達(dá)到對(duì)儀器頻帶寬度的要求。這樣的調(diào)整是在儀器的機(jī)械部分不變動(dòng)的前提下進(jìn)行的,所以對(duì)儀器的安裝沒(méi)有影響。儀器參數(shù)調(diào)整之后,高頻端信號(hào)響應(yīng)能力與其他各臺(tái)儀器區(qū)別明顯,能夠記錄到更多的高頻信號(hào),但同時(shí)也引入了更多的噪聲信號(hào)。而對(duì)于長(zhǎng)周期形變儀器而言,信噪比是其重要的參數(shù)之一,儀器參數(shù)的調(diào)整雖然滿(mǎn)足了部分測(cè)試要求,但引入的大量的噪聲信號(hào)使儀器記錄曲線(xiàn)變得不再平滑,對(duì)于幅度較小的同震信號(hào)或異常的高頻信號(hào)而言,其提取將會(huì)變得更加困難。因此認(rèn)為,發(fā)揮儀器的優(yōu)勢(shì)觀(guān)測(cè)頻段才是儀器設(shè)計(jì)者的初衷,過(guò)多地跨頻段使用儀器,效果值得商榷。
致謝:廖成旺研究員、呂永清研究員、趙瑩博士、馬廣慶高級(jí)工程師、焦成麗高級(jí)工程師為本文提供專(zhuān)業(yè)支持,中國(guó)地震局地殼形變臺(tái)網(wǎng)中心、河北省地震局承德中心地震臺(tái)、易縣地震臺(tái)等為本文提供數(shù)據(jù)支持,在此一并表示感謝!