• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      HIFU治療中生物組織對(duì)聲學(xué)焦域及聲的非線性影響的研究

      2021-01-08 08:31:50徐遨璇王月兵鄭慧峰曹永剛
      聲學(xué)技術(shù) 2020年6期
      關(guān)鍵詞:衰減系數(shù)聲壓聲場(chǎng)

      徐遨璇,王月兵,鄭慧峰,曹永剛,陳 璐

      (中國計(jì)量大學(xué)計(jì)量測(cè)試工程學(xué)院,浙江杭州 310018)

      0 引 言

      近年來,隨著超聲診斷和治療的廣泛應(yīng)用,其安全性問題越來越受到人們的關(guān)注。高強(qiáng)度聚焦超聲(High Intensity Focused Ultrasound, HIFU)治療技術(shù)已應(yīng)用于多種實(shí)體性腫瘤的臨床治療中,其作為一種非侵入性的方法備受關(guān)注[1-2],而我國在 HIFU技術(shù)領(lǐng)域更是位于國際前列,在不斷的研究和臨床試驗(yàn)中,我國的HIFU技術(shù)越來越成熟[3]。HIFU起源于二十世紀(jì)四十年代,由Lynn等[4]提出了HIFU,并利用它熱凝固體內(nèi)組織。HIFU技術(shù)的治療過程是把水介質(zhì)作為媒介,將人體外的超聲波通過水體聚焦傳播到人體內(nèi)某一治療靶區(qū),期間需經(jīng)過生物組織層,然后在焦點(diǎn)區(qū)域形成瞬間高溫,從而殺滅腫瘤組織形成損傷灶,達(dá)到無創(chuàng)治療的目的[5]。在人體組織內(nèi)形成的損傷區(qū)域是HIFU治療的關(guān)鍵,而損傷區(qū)域受聲學(xué)焦域的位置、形狀、聲強(qiáng)和諧波成分共同影響。在HIFU治療過程中,超聲束需穿過皮膚、脂肪、肌肉等多層組織才能到達(dá)靶區(qū),超聲能量沉積在聲通道組織而發(fā)生衰減。由于生物組織的不均一性以及聲學(xué)特性的差異,致使組織的厚度及介質(zhì)內(nèi)的傳播聲速存在差異,這使得HIFU焦域的位置和形狀發(fā)生改變,從而影響靶區(qū)的能量分布和損傷灶的形成。同時(shí)生物組織影響HIFU焦域的強(qiáng)度變化,其非線性效應(yīng)會(huì)發(fā)生改變,對(duì)治療也會(huì)產(chǎn)生一定的影響。

      在生物損傷區(qū)域研究方面,國外Watkin等[6]進(jìn)行了HIFU在生物組織中形成損傷的實(shí)驗(yàn),研究表明,隨著聲強(qiáng)的提高,焦域位置前移且形狀由長橢球形逐漸變?yōu)椤膀蝌叫巍薄鴥?nèi)重慶醫(yī)科大學(xué)實(shí)驗(yàn)研究了體外HIFU不同頻率對(duì)損傷區(qū)域的影響,得到的結(jié)果與文獻(xiàn)[6]中的實(shí)驗(yàn)結(jié)果一致[7]。數(shù)值仿真和有限元仿真是預(yù)測(cè)HIFU治療中聲場(chǎng)分布的有效方法,眾多學(xué)者對(duì)HIFU治療過程進(jìn)行了仿真研究和臨床試驗(yàn)研究。在文獻(xiàn)[8]中,中國科學(xué)院的周玉祿和程建政通過數(shù)值仿真,研究了層狀生物媒質(zhì)中的聚焦聲場(chǎng),利用水、蓖麻油等介質(zhì)代替生物組織進(jìn)行實(shí)驗(yàn)驗(yàn)證。在文獻(xiàn)[9]中,浙江大學(xué)的解卓麗和鄭音飛為分析醫(yī)學(xué)超聲在非均勻組織中的分布特性,建立了超聲發(fā)射聲場(chǎng)的計(jì)算模型。在文獻(xiàn)[10]中,常詩卉和菅喜岐采用Westervelt方程的近似式,并結(jié)合Pennes生物熱傳導(dǎo)方程進(jìn)行數(shù)值仿真,以豬肝腫瘤為例,預(yù)測(cè)高強(qiáng)度聚焦超聲治療的溫度分布和確定治療劑量。在文獻(xiàn)[11]中,楊含和李發(fā)琪利用連續(xù)高強(qiáng)度聚焦超聲與脈沖高強(qiáng)度聚焦超聲輻照離體牛肝組織,進(jìn)而觀察靶區(qū)形態(tài),測(cè)算壞死體積。

      可見,前人或通過實(shí)驗(yàn)觀察現(xiàn)象得出結(jié)論,或進(jìn)行純粹的仿真,進(jìn)而分析得出結(jié)論。本文希望仿真和實(shí)驗(yàn)相輔相成,互相驗(yàn)證,并且自己測(cè)量組織聲參數(shù),貼合實(shí)際測(cè)量環(huán)境。同時(shí)前人多采用KZK方程作為出發(fā)點(diǎn),本文將利用Westervelt方程來建立聲學(xué)傳播模型。

      1 基本理論

      為探究生物組織對(duì) HIFU聲場(chǎng)的影響,HIFU探頭采用 61陣元的球殼點(diǎn)聚焦陣列換能器,探頭的示意圖如圖1所示。球殼陣由一個(gè)曲率半徑df=150 mm、橫向直徑L=160 mm的剛性球冠體和61個(gè)直徑d=16 mm的平面圓形陣元組成。

      圖1 HIFU探頭示意圖Fig.1 Schematic diagram of HIFU probe

      提高換能器的工作功率,產(chǎn)生高強(qiáng)度聚焦超聲,當(dāng)其輻射聲波達(dá)到有限振幅波范圍時(shí),在傳播過程中會(huì)產(chǎn)生非線性諧波,因此線性波動(dòng)方程不再適用于聲波傳播建模。目前多數(shù)非線性聲學(xué)研究工作采用 KZK方程作為出發(fā)點(diǎn),該方程以拋物近似為前提,求取的結(jié)果是有限孔徑聲源在無限大平面自由空間內(nèi)傳播時(shí)軸向和遠(yuǎn)場(chǎng)的普遍特征[12]。而本文將利用Westervelt方程來建立聲學(xué)傳播模型,該方程由線性波動(dòng)方程轉(zhuǎn)化而來,是累積非線性效應(yīng)在局部非線性效應(yīng)中占主導(dǎo)時(shí)全二階波動(dòng)方程的近似,具體表示如下[13]:

      式中:p為總聲壓;ρ0和c0分別為介質(zhì)的密度和傳播聲速;β=1+B/(2A)為非線性系數(shù),其中B/(2A)為非線性參數(shù);δ為聲擴(kuò)散率;t為傳播時(shí)間;?為拉普拉斯算符。

      2 仿真模型及仿真結(jié)果

      2.1 HIFU治療的仿真模型

      圖2為直徑200 mm,長230 mm的圓柱形仿真模型。其中換能器為曲率半徑R=150 mm的聚焦換能器,激勵(lì)信號(hào)為頻率f=1.0 MHz的正弦波。邊界采用完美匹配層(Perfectly Matched Layer, PML)作為吸收條件。為了精確解析整個(gè)聲學(xué)域的壓力梯度變化,使用二次單元進(jìn)行離散化處理,在進(jìn)行頻域仿真時(shí),對(duì)其采用大小為λ/6的細(xì)化網(wǎng)格,而時(shí)域仿真由于網(wǎng)格大小受解析諧波數(shù)量的影響,因此將網(wǎng)格大小定義為dx=c0/(6N0f0)(c0為基波聲速,N0為諧波數(shù)量,f0為基波頻率)。

      圖2 高強(qiáng)度聚焦超聲治療的仿真模型(單位:mm)Fig.2 Simulation model of HIFU treatment (unit: mm)

      2.2 仿真參數(shù)

      為達(dá)到理想的仿真效果,需要獲取真實(shí)的組織聲學(xué)參數(shù),因此利用水浸式脈沖透射法進(jìn)行豬肉組織的聲學(xué)參數(shù)測(cè)量,實(shí)驗(yàn)示意圖如圖3所示。在裝滿水的小水槽中,發(fā)射換能器和接收換能器相距Lw,面對(duì)面放置。當(dāng)發(fā)射換能器發(fā)射聲壓為p0的平面聲波時(shí),接收換能器接收到的聲壓為pi。

      圖3 聲學(xué)參數(shù)測(cè)量實(shí)驗(yàn)示意圖Fig.3 Experimental layout of acoustic parameter measurement

      在充滿除氣水的水槽中,未放入生物組織時(shí),設(shè)接收聲壓為p1,則p1的表達(dá)式為

      式中:αw為水的衰減系數(shù);Lw為換能器表面之間的距離。

      在除氣水中加入生物組織,設(shè)接收到的聲壓為p2,則p2的表達(dá)式為

      式中:αt為生物組織的聲衰減系數(shù),Lt為生物組織的厚度,Twt和Ttw分別為從水到生物組織和從生物組織到水的聲傳播透射系數(shù),其表達(dá)式為[14]

      式中:zt和zw為生物組織和水的聲阻抗。在本實(shí)驗(yàn)中,豬肉樣本切割得較為規(guī)整,形狀為長方體,密度通過聯(lián)系體積和質(zhì)量,多次取平均得到。生物組織的聲衰減系數(shù)的表達(dá)式為

      另外,為獲取樣本的聲速,假設(shè)在插入樣品的情況下,從發(fā)射換能器發(fā)射的聲波到接收換能器接收所需的時(shí)間為twt,則聲速和聲波傳播距離之間的關(guān)系表達(dá)式為

      式中:ct為生物組織的聲速;cw為水的聲速。因此,生物組的聲速ct為

      本實(shí)驗(yàn)中,測(cè)量多塊豬肉組織成分,通過讀取聲路徑傳播時(shí)間和接收聲壓信號(hào),將測(cè)量所得數(shù)據(jù)換算為聲速和聲衰減系數(shù),分別進(jìn)行線性擬合,得到各組織的平均聲速和聲衰減系數(shù)。其中發(fā)射換能器和接收換能器的工作頻率為2 MHz,忽略衍射效應(yīng),同時(shí)考慮到聲介質(zhì)的聲阻抗有差異,聲界面的反射效應(yīng)不可忽略。因此由式(6)計(jì)算出頻率為2 MHz下的聲波在各組織中的聲衰減系數(shù)(dB·cm-1)。然而生物組織的聲衰減系數(shù)與頻率并非完全線性關(guān)系,生物組織的超聲衰減系數(shù)與頻率的關(guān)系為

      式中:γ為生物組織在 1 MHz下的聲衰減系數(shù);f為測(cè)量時(shí)超聲頻率;n為頻率依賴指數(shù)。通過文獻(xiàn)整理得到皮膚的頻率依賴指數(shù)為 2,脂肪和肌肉的頻率依賴指數(shù)為 1.2[15]。上述測(cè)量結(jié)果是在頻率為2 MHz下測(cè)量,需要進(jìn)行換算,得到頻率為1 MHz下的聲衰減系數(shù),換算后的具體結(jié)果如表1所示。

      表1 仿真用的生物組織參數(shù)(實(shí)測(cè))Table 1 Biological tissue parameters (measured) for simulation

      2.3 有限元仿真結(jié)果

      圖4 焦域分布變化Fig.4 Change of focal region distribution

      為了能直觀地觀察聲學(xué)焦域整體的變化,并且考慮到三維組織體積過大,劃分細(xì)化網(wǎng)格,會(huì)占用大量的計(jì)算內(nèi)存,因此減少模擬計(jì)算量,此處將圖2作為二維模型進(jìn)行仿真計(jì)算。焦域聲場(chǎng)的仿真結(jié)果如圖4所示,顯示了焦域的形狀和位置的變化。圖4(a)為HIFU自由場(chǎng)的聲場(chǎng)分布,聲場(chǎng)焦域?yàn)闄E圓形。圖4(b)為HIFU透過厚度為72 mm的豬肉組織后的聲場(chǎng)分布,該仿真結(jié)果不考慮組織結(jié)構(gòu)上的不均一性,即各組織層的分界面為平整面,焦域的形狀基本不變。圖 4(c)同樣為 HIFU透過厚度為72 mm的豬肉組織后的聲場(chǎng)分布,該結(jié)果在圖4(b)的基礎(chǔ)上考慮組織結(jié)構(gòu)的不均一性,其組織的各層分界面為實(shí)際測(cè)量豬肉組織的分界面,焦域聲場(chǎng)發(fā)生略微的畸變,而且聲學(xué)焦域除了向換能器一側(cè)靠近,在焦平面方向也發(fā)生偏移。組織結(jié)構(gòu)的不均一性以及聲學(xué)特性的差異,導(dǎo)致HIFU透過組織后的聲學(xué)焦域相比自由場(chǎng)發(fā)生畸變,焦域位置也發(fā)生改變。同時(shí)組織具有較高的聲衰減能力,焦域的聲壓也發(fā)生了改變,具體仿真結(jié)果如表2所示,其中,聲壓值比為加入生物組織后的聲壓與未加入組織時(shí)自由場(chǎng)聲壓的比值,焦域?qū)挾仁前凑章晧鹤畲笾?6 dB處的寬度進(jìn)行計(jì)算。

      表2 焦域仿真結(jié)果Table 2 Results of focal field simulation

      由圖4和表2可知,由于生物組織聲阻抗的差異,導(dǎo)致HIFU焦域在聲軸上發(fā)生偏移;組織的結(jié)構(gòu)不均一性導(dǎo)致HIFU焦域在焦平面內(nèi)也會(huì)發(fā)生偏移;同時(shí)根據(jù)仿真結(jié)果可知,焦域的寬度基本不變,無散焦現(xiàn)象,即超聲束經(jīng)過生物組織不會(huì)發(fā)生聲束擴(kuò)散;此處由于采用的實(shí)際離體豬肉組織分層界面較為平整,因此結(jié)構(gòu)上的非均一性并沒有引起巨大的誤差,從仿真結(jié)果分析,兩者的聲壓值比只相差2.8%。

      由于高強(qiáng)度聚焦超聲的輻射聲波達(dá)到有限振幅波范圍,因此會(huì)產(chǎn)生非線性諧波。在進(jìn)行非線性仿真時(shí),通過上述二維焦域仿真結(jié)果發(fā)現(xiàn),平整組織層(理想)和非平整組織層(實(shí)際)聲壓衰減后的聲壓值比非常接近。因此此處仿真模型按照?qǐng)D 2所示,采用二維軸對(duì)稱進(jìn)行仿真,組織層為平整界面,仿真結(jié)果如圖5所示。

      圖5 高強(qiáng)度聚焦超聲的諧波分布Fig.5 Distribution of HIFU harmonics

      生物組織作為一種似流體組織,其結(jié)構(gòu)具有不均一性,因而超聲波在生物組織中傳播時(shí)引起的聲衰減機(jī)理是復(fù)雜的。在這種似流體介質(zhì)中,超聲衰減由超聲吸收和超聲散射兩部分組成。通過觀察圖5中的仿真結(jié)果發(fā)現(xiàn),由于生物組織具有較強(qiáng)的超聲衰減能力,HIFU進(jìn)入生物組織后,非線性效應(yīng)會(huì)發(fā)生較大程度的改變,高次諧波相對(duì)基波所占比例下降幅度更大。在HIFU治療過程中,HIFU的非線性效應(yīng)是極其重要的組成部分。因此,此處對(duì)HIFU的非線性效應(yīng)進(jìn)行了仿真研究。

      3 實(shí)驗(yàn)與結(jié)果

      3.1 實(shí)驗(yàn)系統(tǒng)搭建

      為了通過實(shí)驗(yàn)進(jìn)一步驗(yàn)證仿真結(jié)果的可靠性,搭建如圖6所示的實(shí)驗(yàn)測(cè)試系統(tǒng)。首先激勵(lì)球殼點(diǎn)聚焦陣列換能器工作,在水箱中形成聚焦聲場(chǎng)。通過夾具夾持光纖水聽器,利用三維運(yùn)動(dòng)機(jī)械結(jié)構(gòu)對(duì)HIFU聲場(chǎng)焦域進(jìn)行測(cè)量,輸出信號(hào)呈現(xiàn)在示波器上,最終將數(shù)據(jù)讀取并存儲(chǔ)在電腦中。實(shí)驗(yàn)器材:球殼點(diǎn)聚焦陣列換能器、AR功率放大器、信號(hào)源、示波器、光纖水聽器、光纖解調(diào)器、PLC機(jī)箱、上位機(jī)。設(shè)備型號(hào):示波器(RIGOL MSO5104)、信號(hào)源(RIGOL DG4062)、光纖解調(diào)器(FOH-880)、功率放大器(rf/microwave instrumentation model 800A3B)。

      3.2 HIFU聲學(xué)焦域測(cè)量

      圖6 實(shí)驗(yàn)系統(tǒng)Fig.6 Experimental system

      選取新鮮離體豬肉組織,均浸泡在0.9%的生理鹽水中。標(biāo)本的準(zhǔn)備:洗凈、切塊、浸泡。豬肉組織共切成120 mm×120 mm的4塊長方體,其厚度分別為72.0、70.5、58.0、76.0 mm,并分別標(biāo)記為1~4號(hào),復(fù)溫到20℃,放置于除氣水(水中氣體含量為0.14 mg·L-1)中,等待實(shí)驗(yàn)使用。各塊豬肉組織成分厚度如表3所示。

      表3 各塊離體豬肉組織成分Table 3 Components of isolated pork tissues

      在放入生物組織前,先通過3D定位系統(tǒng)將光纖水聽器移動(dòng)到HIFU自由場(chǎng)聲學(xué)焦域的焦點(diǎn)處,測(cè)量焦平面聲場(chǎng)分布。然后將準(zhǔn)備好的生物組織插入到球殼點(diǎn)聚焦陣列換能器的聲通道中。在對(duì)焦域聲場(chǎng)進(jìn)行測(cè)量時(shí),首先使生物組織能覆蓋HIFU探頭輻射的聲束,所切割的豬肉組織為長方體,長為120 mm,寬為120 mm,放置于距離換能器50 mm處,通過幾何分析,能夠滿足覆蓋全聲束的測(cè)量條件。其次通過3D定位系統(tǒng)微調(diào)運(yùn)動(dòng)機(jī)械結(jié)構(gòu),移動(dòng)光纖水聽器尋找焦點(diǎn)所在位置。最后,測(cè)量焦平面y軸和z軸上的聲壓分布。為精準(zhǔn)測(cè)量其聲壓分布,測(cè)量范圍為 10 mm×10 mm,采集間距為0.2 mm,測(cè)量結(jié)果如圖7所示。

      圖7 焦平面聲場(chǎng)分布Fig.7 Distribution of acoustical field in focal plane

      表4 焦域測(cè)量結(jié)果Table 4 Results of focal field measurement

      圖7分別為自由場(chǎng)和加入生物組織后焦平面的聲場(chǎng)分布對(duì)比。從圖7中可以看出:聲通道存在生物組織時(shí),會(huì)引起 HIFU焦域在焦平面內(nèi)發(fā)生偏移;超聲波束通過生物組織,焦域?qū)挾然静蛔儯瑳]有發(fā)生擴(kuò)散,沒有出現(xiàn)散焦現(xiàn)象。具體實(shí)驗(yàn)結(jié)果如表4所示。從表4中可以發(fā)現(xiàn),由于介質(zhì)聲阻抗的差異,除了導(dǎo)致HIFU的焦點(diǎn)在聲軸方向發(fā)生偏移,還產(chǎn)生了焦平面內(nèi)的偏移。從測(cè)量結(jié)果中發(fā)現(xiàn),平均焦域?qū)挾戎挥?.6 mm,焦平面方向的偏移量接近1 mm,因此焦平面內(nèi)的偏移對(duì)HIFU的治療影響非常大。同時(shí)由于生物組織的黏滯導(dǎo)致大量能量沉積在組織中,焦點(diǎn)處的聲壓幅值顯著下降,在測(cè)量的結(jié)果中,聲壓值比均下降至50%以下。

      由于生物組織會(huì)影響 HIFU在焦域的強(qiáng)度變化,HIFU非線性效應(yīng)也會(huì)改變,隨之HIFU治療效果也會(huì)產(chǎn)生變化。本文選取自由場(chǎng)焦域和放入1號(hào)豬肉組織后焦域的非線性效應(yīng)作為研究對(duì)象。函數(shù)發(fā)生器發(fā)射頻率為1 MHz,幅值為600 mV的正弦脈沖波,AR功率放大器的增益為65%,輻射聲功率為25.77 W。測(cè)量結(jié)果如表5所示。

      表5 諧波測(cè)量結(jié)果Table 5 Results of harmonic measurement

      從表5可以發(fā)現(xiàn),聲通道中的生物組織,極大程度減弱了焦域的聲強(qiáng),進(jìn)而影響非線性效應(yīng)。當(dāng)未加入生物組織時(shí),非線性明顯,諧波次數(shù)達(dá)到5 MHz。當(dāng)在聲通道加入生物組織,由于生物組織的黏滯,導(dǎo)致焦域的聲強(qiáng)衰減,非線性效應(yīng)不顯著,諧波發(fā)生衰減,且諧波次數(shù)越高,聲衰減值越大。在本研究中,當(dāng)加入生物組織后,五次及以上的諧波已經(jīng)不存在,且高次諧波相對(duì)基波所占比例下降幅度更大。

      4 結(jié) 論

      本文通過測(cè)量新鮮離體豬肉組織的聲學(xué)參數(shù)進(jìn)行有限元仿真,研究了HIFU治療過程中聲束通過生物組織后的聲學(xué)焦域變化。同時(shí)搭建測(cè)量系統(tǒng),利用離體組織進(jìn)行實(shí)驗(yàn)來驗(yàn)證仿真結(jié)論。

      仿真和實(shí)驗(yàn)研究表明:(1) 當(dāng)聲通道存在生物組織時(shí),HIFU焦域會(huì)發(fā)生偏移,焦平面內(nèi)的偏移量相較于焦域?qū)挾?,?duì)HIFU治療精確性影響較大;(2) 當(dāng) HIFU經(jīng)過生物組織,沒有散焦現(xiàn)象,即聲學(xué)焦域的寬度基本不變,這避免了因焦域?qū)挾茸兇?,在治療過程中損害正常組織的情況;(3) 焦域內(nèi)高次諧波所占比例下降,且諧波次數(shù)越高,衰減越明顯,而非線性效應(yīng)的變化會(huì)影響到治療效果,需要進(jìn)行適當(dāng)?shù)哪芰垦a(bǔ)償。

      HIFU作為一種新興的腫瘤微創(chuàng)外科治療技術(shù),具有獨(dú)特的優(yōu)勢(shì),但是HIFU技術(shù)在推廣的過程中,還有許多的問題等待解決如損傷灶的形成等。本文研究過程中,將生物組織進(jìn)行了簡(jiǎn)化,具體臨床情況相對(duì)更加復(fù)雜,因此后期需繼續(xù)完善仿真模型,使其更加貼切實(shí)際,為HIFU試驗(yàn)提供有利意見,對(duì)于HIFU通過各種復(fù)雜生物組織形成聲學(xué)焦域有待進(jìn)一步研究。

      猜你喜歡
      衰減系數(shù)聲壓聲場(chǎng)
      基于嘴唇處的聲壓數(shù)據(jù)確定人體聲道半徑
      基于BIM的鐵路車站聲場(chǎng)仿真分析研究
      車輛結(jié)構(gòu)噪聲傳遞特性及其峰值噪聲成因的分析
      汽車工程(2018年12期)2019-01-29 06:46:36
      探尋360°全聲場(chǎng)發(fā)聲門道
      復(fù)合材料孔隙率的超聲檢測(cè)衰減系數(shù)影響因素
      近岸及內(nèi)陸二類水體漫衰減系數(shù)的遙感反演研究進(jìn)展
      對(duì)《電磁波衰減系數(shù)特性分析》結(jié)果的猜想
      HT250材料超聲探傷中的衰減性探究
      基于GIS內(nèi)部放電聲壓特性進(jìn)行閃絡(luò)定位的研究
      板結(jié)構(gòu)-聲場(chǎng)耦合分析的FE-LSPIM/FE法
      拜泉县| 德江县| 桐柏县| 巴南区| 哈尔滨市| 鄄城县| 恩施市| 凤翔县| 滦南县| 于田县| 太白县| 宁南县| 巴中市| 井研县| 陇川县| 文山县| 布尔津县| 武山县| 任丘市| 和林格尔县| 金坛市| 安图县| 晋城| 西和县| 潮安县| 沅江市| 芜湖市| 广南县| 古交市| 桑日县| 上栗县| 平泉县| 清远市| 化德县| 东兰县| 旅游| 湘乡市| 肃宁县| 宜昌市| 岑巩县| 涞水县|