国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

水-鋼界面凸面相控陣聲場(chǎng)仿真

2017-06-05 14:15:42郭偉燦錢盛杰繆存堅(jiān)
無(wú)損檢測(cè) 2017年5期
關(guān)鍵詞:柵瓣凸面管狀

郭偉燦,錢盛杰,繆存堅(jiān)

(1. 浙江省特種設(shè)備檢驗(yàn)研究院,杭州 310020;2.寧波市特種設(shè)備檢驗(yàn)研究院,寧波 315048)

水-鋼界面凸面相控陣聲場(chǎng)仿真

郭偉燦1,錢盛杰2,繆存堅(jiān)1

(1. 浙江省特種設(shè)備檢驗(yàn)研究院,杭州 310020;2.寧波市特種設(shè)備檢驗(yàn)研究院,寧波 315048)

采用凸面超聲相控陣探頭對(duì)管狀設(shè)備進(jìn)行內(nèi)檢測(cè),針對(duì)水-鋼二層介質(zhì)的結(jié)構(gòu)特征,理論分析了相控陣聲場(chǎng)的聚焦性能。根據(jù)近場(chǎng)球面波束形成模型,按聲束路徑求得的延時(shí)時(shí)間作為各陣元的提前激發(fā)時(shí)間,設(shè)計(jì)了水-鋼二層介質(zhì)相控陣聚焦法則。利用CIVA軟件對(duì)相控陣聲場(chǎng)進(jìn)行建模仿真研究,確定了陣元數(shù)目、陣元尺寸、超聲波頻率等相控陣探頭參數(shù)。最后,根據(jù)優(yōu)化后的相控陣參數(shù)研制了512大陣元數(shù)雙層結(jié)構(gòu)的圓環(huán)形探頭,再通過(guò)帶有人工缺陷的管狀模擬試樣進(jìn)行檢測(cè)試驗(yàn)。結(jié)果表明:提出的方法對(duì)腐蝕、氣孔缺陷具有更好的檢測(cè)靈敏度和效率,能夠?qū)崿F(xiàn)管狀設(shè)備的腐蝕和埋藏缺陷的內(nèi)檢測(cè)。

管狀設(shè)備;內(nèi)檢測(cè);超聲檢測(cè);相控陣;超聲仿真

儲(chǔ)氣井和壓力管道等管狀設(shè)備通常是埋地和架空敷設(shè),在外面對(duì)其進(jìn)行檢測(cè)十分困難。如何實(shí)現(xiàn)不搭腳手架、不拆保溫層、不開挖地面的內(nèi)檢測(cè),是對(duì)此類設(shè)備實(shí)行定期檢驗(yàn)的關(guān)鍵。管狀設(shè)備內(nèi)檢測(cè)常用的無(wú)損檢測(cè)技術(shù)包括:針對(duì)鐵磁性材料的漏磁通檢測(cè)法、超聲檢測(cè)法、基于激光與紅外等光源的成像和視頻技術(shù)[1-3]。超聲檢測(cè)具有檢測(cè)精度高、檢測(cè)結(jié)果可定量、內(nèi)外腐蝕可直接分辨、檢測(cè)結(jié)果不受管材影響等優(yōu)點(diǎn)[4-5]。筆者針對(duì)管狀設(shè)備的結(jié)構(gòu)特點(diǎn),提出了圓柱形凸面相控陣檢測(cè)技術(shù),通過(guò)水浸耦合,使探頭產(chǎn)生的聲束能夠覆蓋整個(gè)管道的橫截面積。針對(duì)凸面相控陣探頭在水-鋼二層介質(zhì)中的結(jié)構(gòu)特征,理論分析了相控陣聲場(chǎng)的聚焦性能,并根據(jù)近場(chǎng)球面波束形成模型,以聲束路徑求得的延時(shí)時(shí)間作為各陣元的提前激發(fā)時(shí)間,設(shè)計(jì)了水-鋼二層介質(zhì)相控陣聚焦法則。在此基礎(chǔ)上,利用CIVA軟件對(duì)相控陣聲場(chǎng)進(jìn)行建模仿真研究,以波束的主瓣尖銳、較低的旁瓣、消除柵瓣等為相控陣探頭的設(shè)計(jì)準(zhǔn)則,對(duì)相控陣參數(shù)進(jìn)行優(yōu)化分析,根據(jù)優(yōu)化后的相控陣參數(shù)研制了512大陣元圓柱形凸面相控陣探頭,通過(guò)帶有人工缺陷的管道模擬試樣進(jìn)行檢測(cè)。結(jié)果表明:該方法對(duì)腐蝕、氣孔具有更好的檢測(cè)靈敏度和效率,能夠?qū)崿F(xiàn)管狀設(shè)備的自動(dòng)內(nèi)檢測(cè),具有很大的應(yīng)用潛力。

1 凸面超聲相控陣內(nèi)檢測(cè)技術(shù)

采用圓環(huán)形凸面相控陣探頭伸入管狀設(shè)備內(nèi)部進(jìn)行檢測(cè),通過(guò)充水耦合,使圓環(huán)形探頭產(chǎn)生的聲束能夠覆蓋整個(gè)管道的橫截面積。凸面相控陣檢測(cè)原理示意如圖1所示,以同時(shí)激發(fā)16個(gè)陣元為例,采用相控陣電子線掃描方式,先激發(fā)第1~16陣元組,形成在鋼中聚焦點(diǎn)1的聚焦聲束;完成掃描后以同樣的聚焦法則激發(fā)第2~17陣元組, 形成在鋼中聚焦點(diǎn)2的聚焦聲束;依此類推,直到完成整個(gè)圓周向的電子線掃描[6-7]。

圖1 凸面相控陣檢測(cè)原理示意

(1)

式中:c1為超聲在鋼中的聲速;c2為超聲在水中的聲速。

圖2 相控陣延時(shí)時(shí)間簡(jiǎn)易計(jì)算模型

在圖2所示的平面坐標(biāo)中,相控陣探頭通過(guò)延時(shí)發(fā)射,使得原半徑為R的凸面陣元,形成為曲率半徑為F(AD段長(zhǎng)度)的凹面新波陣面,由此可以求得第m號(hào)陣元的延時(shí)時(shí)間τm,設(shè)AC段長(zhǎng)度為L(zhǎng)AC。

(2)

同樣可以求得第m-1號(hào)陣元的延時(shí)時(shí)間τm-1為:

(3)

用同樣方法可以求得第m-2、m-3、…、-m號(hào)陣元的延時(shí)時(shí)間τm-2、τm-3、…、τ-m,形成聚焦法則激發(fā)第m至-m陣元組,完成掃描后以同樣的聚焦法則激發(fā)第m-1至-(m+1)陣元組,依此類推,直到完成整個(gè)圓周向的電子線掃描。

圖3 檢測(cè)工作原理示意

檢測(cè)工作原理示意如圖3所示:探頭激發(fā)的超聲脈沖從水中傳播到達(dá)管壁,管道內(nèi)外壁的脈沖回波經(jīng)過(guò)水被探頭接收,由此便可以對(duì)壓力管道進(jìn)行檢測(cè)。實(shí)際上是采用超聲測(cè)厚原理,超聲波T1通過(guò)水介質(zhì)沿徑向入射到管道內(nèi)壁,在水-鋼界面會(huì)產(chǎn)生一個(gè)強(qiáng)的反射波即R1,該反射波為探頭接收,形成水-鋼界面回波B1;其余透射能量形成透射波T2繼續(xù)前行,到達(dá)管道外壁,同樣會(huì)產(chǎn)生一個(gè)反射波R2,此反射波又一次經(jīng)鋼/水界面透射至水中,被探頭接收,形成外壁界面回波B2;由于鋼中聲速固定,壁厚減少,B1和B2回波之間的距離也隨之減少,從而達(dá)到壁厚檢測(cè)的目的。

2 凸面相控陣CIVA聲場(chǎng)仿真

CIVA軟件是一款能夠?qū)崿F(xiàn)超聲檢測(cè)仿真的商用計(jì)算軟件,以往通過(guò)CIVA實(shí)例[8-9]對(duì)仿真可靠性進(jìn)行了驗(yàn)證。筆者針對(duì)國(guó)內(nèi)常用的規(guī)格(外徑×厚度)為φ177.8 mm×10.36 mm的埋地高壓儲(chǔ)氣井管狀設(shè)備,利用CIVA軟件建立相控陣模型并計(jì)算相控陣探頭在工件內(nèi)部形成的超聲場(chǎng),并對(duì)相控陣的不同激發(fā)陣元個(gè)數(shù)、陣元尺寸、頻率等參數(shù)進(jìn)行仿真??紤]到圓柱形凸面相控陣探頭通過(guò)水浸耦合,為保證探頭在管內(nèi)的通過(guò)性能和第二次避免水-鋼界面回波對(duì)檢測(cè)的干擾,初步設(shè)計(jì)圓柱形凸面相控陣探頭直徑為φ130.4 mm,以512陣元數(shù)計(jì)算,陣元間距為0.8 mm。

2.1 激發(fā)陣元個(gè)數(shù)對(duì)聚焦性能的影響 根據(jù)圖1的檢測(cè)原理,初步選取頻率為5 MHz、陣元寬度為0.7 mm、間距為0.8 mm、長(zhǎng)度為10 mm、偏轉(zhuǎn)角為0°的相控陣參數(shù),焦點(diǎn)落在規(guī)格為φ177.8 mm×10.36 mm的儲(chǔ)氣井外壁,分別對(duì)4,8,16,32個(gè)激發(fā)陣元的聲場(chǎng)聚焦性能進(jìn)行仿真,激發(fā)陣元個(gè)數(shù)對(duì)聚焦性能的影響如圖4所示。

圖4 激發(fā)陣元個(gè)數(shù)對(duì)聚焦性能的影響

由圖4可知:激發(fā)陣元個(gè)數(shù)越多,在焦點(diǎn)處的聚焦效果越好,但超聲波在水中越容易產(chǎn)生柵瓣和旁瓣。針對(duì)5 MHz的檢測(cè)頻率,超聲波在水中的波長(zhǎng)僅為0.29 mm,而陣元間距為0.8 mm,約為波長(zhǎng)的2.8倍,因此,激發(fā)陣元個(gè)數(shù)為32時(shí),即使聲束不偏轉(zhuǎn)也會(huì)出現(xiàn)柵瓣;激發(fā)陣元個(gè)數(shù)為16時(shí),在水中產(chǎn)生嚴(yán)重的旁瓣;激發(fā)陣元個(gè)數(shù)為8時(shí),在水中產(chǎn)生一定程度的旁瓣;激發(fā)陣元個(gè)數(shù)為4時(shí),只在水中產(chǎn)生輕微的旁瓣??紤]到凸面相控陣陣元偏轉(zhuǎn)角度更大,為保證水中聲場(chǎng)不出現(xiàn)嚴(yán)重的柵瓣和旁瓣,采用4~8個(gè)激發(fā)陣元較為合適,具體選擇時(shí)根據(jù)工藝試驗(yàn)確定。

2.2 陣元尺寸對(duì)聚焦性能的影響 陣元尺寸包括陣元間距d、陣元寬度a、陣元長(zhǎng)度L。其中對(duì)相控陣聚焦性能影響較大的是陣元間距d和陣元寬度a,通常情況下為了抑制旁瓣,應(yīng)使陣元寬度a盡量接近陣元間距d。筆者針對(duì)陣元數(shù)N為8,12,16的不同激發(fā)陣元,對(duì)a為0.3 mm、d為0.4 mm,a為0.7 mm、d為0.8 mm,a為1.1 mm、d為1.2 mm的三組不同陣元尺寸條件進(jìn)行CIVA聲場(chǎng)仿真。不同陣元尺寸的CIVA聲場(chǎng)仿真結(jié)果如圖5所示。

由圖5仿真結(jié)果可得: 在3組不同尺寸的相控陣列中,只有a為0.3 mm、d為0.4 mm陣元尺寸的相控陣列的超聲波能夠?qū)崿F(xiàn)控制旁瓣和消除柵瓣的要求,其他兩種陣元尺寸的相控陣列在激發(fā)陣元數(shù)不小于8時(shí),在水中均不可避免地會(huì)產(chǎn)生較嚴(yán)重的旁瓣,甚至出現(xiàn)柵瓣。通過(guò)上述比較,為保證水中聲場(chǎng)不出現(xiàn)嚴(yán)重的柵瓣和旁瓣,宜采用a為0.3 mm、d為0.4 mm的陣元尺寸。具體選擇激發(fā)陣元數(shù)N時(shí),應(yīng)綜合考慮超聲波在水中和工件中的聲場(chǎng)聚焦效果。雖然N為8時(shí)超聲波在水中的聚焦效果較好,但在工件中的聚焦效果較差;在N為12,16時(shí),超聲波在工件中設(shè)定焦點(diǎn)處的聚焦效果較好,且在水中只產(chǎn)生輕微的旁瓣,因此激發(fā)陣元數(shù)宜選用12~16較為合適,具體選擇時(shí)根據(jù)工藝試驗(yàn)確定。

2.3 檢測(cè)頻率對(duì)聚焦性能的影響

在超聲相控陣參數(shù)設(shè)計(jì)時(shí),超聲波頻率的選取至關(guān)重要。頻率對(duì)超聲檢測(cè)的影響有多個(gè)方面。頻率越高,相控陣探頭的焦點(diǎn)越小,聚焦能力和分辨力提高,但容易產(chǎn)生旁瓣和柵瓣。

針對(duì)陣元寬度為0.7 mm、間距為0.8 mm、長(zhǎng)度為10 mm、偏轉(zhuǎn)角為0°的相控陣參數(shù),陣元數(shù)為16、焦點(diǎn)落在φ177.8 mm×10.36 mm的儲(chǔ)氣井外壁,頻率分別選取2.5,5,10 MHz進(jìn)行仿真分析,不同頻率的CIVA聲場(chǎng)仿真結(jié)果如圖6所示。

圖5 不同陣元尺寸的CIVA聲場(chǎng)仿真結(jié)果

圖6 不同頻率的CIVA聲場(chǎng)仿真結(jié)果

從圖6的仿真結(jié)果可以得出:?jiǎn)渭儚乃芯劢剐Ч鸵种古园旰蜄虐昕紤],選擇頻率為2.5 MHz較為合適;選擇頻率為10 MHz時(shí),在水中會(huì)不可避免地出現(xiàn)柵瓣;選擇頻率為5 MHz時(shí),也會(huì)在水中產(chǎn)生不同程度的旁瓣,激發(fā)陣元數(shù)越多,旁瓣越嚴(yán)重。但頻率為2.5 MHz的相控陣探頭在工件中的聚焦能力較差,縱向和橫向分辨力也較低,因此在選取頻率時(shí)應(yīng)綜合考慮聚焦能力、分辨力和超聲波衰減等多種因素,如果頻率選取2.5 MHz,激發(fā)陣元數(shù)可適當(dāng)多一些,如果頻率選取5 MHz,為抑止水中聲場(chǎng)的旁瓣,激發(fā)陣元數(shù)不宜超過(guò)8個(gè)。

3 檢測(cè)過(guò)程

根據(jù)上述的CIVA仿真分析,為滿足φ177.8 mm×10.36 mm的埋地高壓儲(chǔ)氣井這類較大直徑的管狀設(shè)備全截面檢測(cè)的要求,選擇陣元尺寸a為0.3 mm、d為0.4 mm,能夠?qū)崿F(xiàn)控制旁瓣和消除柵瓣的要求。然而選擇a為0.3 mm、d為0.4 mm陣元尺寸時(shí),全截面檢測(cè)需要1 024個(gè)陣元相控陣儀器,研發(fā)成本很高。而選擇a為0.7 mm、d為0.8 mm陣元尺寸時(shí),全截面檢測(cè)只需要512個(gè)陣元相控陣儀器,研發(fā)成本大幅降低。筆者通過(guò)CIVA仿真分析得出,如果選取頻率為5 MHz,在激發(fā)陣元數(shù)N為4的條件下的聲場(chǎng)聚焦效果幾乎與激發(fā)陣元數(shù)N為8的條件下,a為0.3 mm、d為0.4 mm陣元尺寸的聲場(chǎng)相同,只是旁瓣稍嚴(yán)重一些。因此筆者采用512個(gè)陣元的圓周凸面超聲相控陣檢測(cè)技術(shù),并委托汕頭超聲儀器研究所有限公司定制了512個(gè)陣元的相控陣模塊和陣元寬度為0.7 mm、間距為0.8 mm、陣元長(zhǎng)度為10 mm、頻率為5 MHz的圓柱形相控陣探頭裝置,探頭裝置由4個(gè)128陣元的凸形探頭組合而成,圓柱形相控陣探頭裝置和相控陣模塊實(shí)物如圖7所示。

圖7 圓柱形相控陣探頭裝置和相控陣模塊實(shí)物

圖8 樣管尺寸示意和實(shí)物

為驗(yàn)證提出的檢測(cè)方法的可靠性,筆者加工了2個(gè)試驗(yàn)樣管,樣管尺寸示意和實(shí)物如圖8所示。檢測(cè)時(shí)將樣管內(nèi)充滿水,將探頭通過(guò)探頭架放入水中進(jìn)行檢測(cè),超聲波通過(guò)水耦合進(jìn)入工件中,再以反射波通過(guò)水回到探頭中,探頭組件和編碼器用于采集信號(hào)輸入,組合相控陣模塊將采集信號(hào)轉(zhuǎn)換成數(shù)字信號(hào),并將數(shù)字信號(hào)傳輸至超聲波檢測(cè)儀。通過(guò)成像軟件形成實(shí)時(shí)A掃描、B掃描、C掃描圖像。如只需檢測(cè)缺陷,則將信號(hào)門閥套在壓力管道內(nèi)部,檢測(cè)并記錄壓力管道內(nèi)部缺陷的位置信息和超聲信息;如只需檢測(cè)壁厚,則將門閥套住水-鋼界面回波和管子外壁回波,并以回波的時(shí)間差(或距離)作為分析信號(hào),并賦以一系列的顏色編碼,以色彩方式顯示管道壁厚值,可以通過(guò)在 C 掃描圖上顏色的變化,顯示整根管道的壁厚情況。

樣管1缺陷檢測(cè)的超聲圖像如圖9所示,樣管2檢測(cè)壁厚得到的C掃描超聲圖像如圖10所示。從圖9,10的試驗(yàn)結(jié)果可以得出:該方法能夠較為精確且有效地檢測(cè)出缺陷的形狀、尺寸和位置,并能顯示管道內(nèi)缺陷與腐蝕情況,例如圖9中2號(hào)與3號(hào)探頭檢測(cè)到的缺陷分別距離外壁大約2 mm和4 mm,且缺陷環(huán)向位置也與實(shí)際情況相符;圖10中C掃描厚度成像能清晰地顯示腐蝕的位置和形狀,對(duì)壁厚的測(cè)量誤差在0.15 mm范圍內(nèi)。

圖9 樣管1缺陷檢測(cè)的超聲圖像

圖10 樣管2檢測(cè)壁厚得到的C掃描超聲圖像

4 結(jié)論

(1) 針對(duì)水-鋼二層介質(zhì)的管狀設(shè)備凸面相控陣檢測(cè),采用CIVA軟件的仿真分析可知,如果陣元間距大于0.8 mm,則在水中容易產(chǎn)生柵瓣和旁瓣,激發(fā)陣元個(gè)數(shù)越多,柵瓣和旁瓣效應(yīng)越嚴(yán)重。因此,從消除柵瓣和抑止旁瓣的角度考慮,宜采用較小的陣元間距和較少的激發(fā)陣元個(gè)數(shù)。

(2) 通過(guò)CIVA仿真分析,針對(duì)規(guī)格為φ177.8 mm×10.36 mm的管狀設(shè)備,采用512個(gè)陣元的圓周凸面超聲相控陣檢測(cè)技術(shù),選擇合適的相控陣參數(shù),能夠在水中消除柵瓣和抑止旁瓣,并在鋼中獲得較好的聚焦效果。

(3) 對(duì)帶有人工缺陷的樣管進(jìn)行檢測(cè)試驗(yàn)結(jié)果表明,提出的圓周凸面相控陣檢測(cè)技術(shù),能夠較為精確且有效地檢測(cè)出管狀設(shè)備缺陷的形狀、尺寸和位置,并能顯示管道內(nèi)缺陷。

[1] 楊理踐,張國(guó)光,劉剛. 管道漏磁內(nèi)檢測(cè)技術(shù)[M].北京: 化學(xué)工業(yè)出版社, 2014.

[2] 劉慧芳,張鵬,周俊杰,等. 油氣管道內(nèi)腐蝕檢測(cè)技術(shù)的現(xiàn)狀與發(fā)展趨勢(shì)[J]. 管道技術(shù)與設(shè)備, 2008 (5): 46-48.

[3] 李響. 管道視覺檢測(cè)技術(shù)的研究[D].沈陽(yáng):沈陽(yáng)工業(yè)大學(xué), 2009.

[4] 肖雄. 換熱管內(nèi)旋轉(zhuǎn)超聲檢測(cè)技術(shù)[J]. 壓力容器, 2013, 30(9): 60-64.

[5] 趙晶. 管道腐蝕裂紋超聲內(nèi)檢測(cè)信號(hào)處理研究[D].北京:北京化工大學(xué), 2008.

[6] 鄭暉,林樹青. 超聲檢測(cè)[M].北京:新華出版社,2008.

[7] 施克仁,郭寓岷. 相控陣超聲成像檢測(cè)[M]. 北京:高等教育出版社,2010.

[8] 錢盛杰,郭偉燦.換熱器管板角焊縫相控陣超聲聲場(chǎng)CIVA仿真和檢測(cè)[J].無(wú)損檢測(cè),2015,37(1):11-15.

[9] 郭偉燦,錢盛杰.厚壁筒形件內(nèi)表面徑向缺陷檢測(cè)的超聲聚焦技術(shù)[J]. 無(wú)損檢測(cè), 2016,38(12):55-58,61.

歡迎網(wǎng)上投稿 歡迎訂閱 歡迎刊登廣告

Convex Phased Array Simulation of Water-Steel Interface Acoustic Field

GUO Weican1, QIAN Shengjie2, MIAO Cunjian1

(1.Zhejiang Provincial Special Equipment Inspection and Research Institute, Hangzhou 310020, China;2.Ningbo Special Equipment Inspection and Research Institute, Ningbo 315048, China)

The convex phased array probe was adopted in this paper for performing internal inspection of the tubular equipment. According to the structure characteristics of two layer medium, the paper theoretically analyzed the focusing performance of phased array acoustic field. Based on the model of near field spherical beam forming, the delay time can be obtained by the route of acoustic beam which can be used as stimulation time for each array elements. Phased array focusing law of the water-steel two layer medium is designed in this paper. The simulation on the phased array acoustic field with CIVA software is completed to determine the phased array probe parameters, including array elements, array element size, ultrasonic frequency and so on. At last,the double layer circular ring probe of 512 array elements is manufactured according to the phased array parameters which have been optimized,the test on tubular samples which contain artificial flaws is performed. Experiments indicate that this technique can detect the corrosion and the porosity flaw,with high sensitivity and efficiency. The corrosion and embedded flaws of the tubular equipment can be detected by internal inspection.

tubular equipment; internal inspection; ultrasonic testing; phased array; ultrasonic simulation

2016-08-23

浙江省質(zhì)量技術(shù)監(jiān)督系統(tǒng)科研計(jì)劃資助項(xiàng)目(20160123)

郭偉燦(1968-),男,博士,教授級(jí)高工,主要從事承壓設(shè)備檢驗(yàn)與無(wú)損檢測(cè)工作

郭偉燦,gwcndt@126.com

10.11973/wsjc201705015

TG115.28

A

1000-6656(2017)05-0070-06

猜你喜歡
柵瓣凸面管狀
利用雙載頻方向圖相乘的稀疏均勻陣柵瓣抑制方法
活動(dòng)期油源斷裂輸導(dǎo)油氣有利凸面脊分布預(yù)測(cè)方法
淺談管狀骨架噴涂方法
大間距非周期陣列柵瓣抑制研究*
無(wú)限管狀區(qū)域上次調(diào)和函數(shù)的邊界性質(zhì)
全息離子束擺動(dòng)刻蝕凸面閃耀光柵制備技術(shù)
用于門、百葉窗或等同物的鉸鏈
科技資訊(2016年19期)2016-11-15 10:46:23
彩寶藝術(shù)大師 Bulgari
鳳凰生活(2016年8期)2016-08-10 12:38:53
一種超寬帶相控陣天線陣方向圖柵瓣抑制方法
基于柵瓣等波動(dòng)陣列設(shè)計(jì)方法
電子科技(2015年7期)2015-03-06 01:30:54
恭城| 平安县| 徐州市| 鱼台县| 许昌县| 根河市| 安龙县| 炎陵县| 汽车| 上高县| 清涧县| 永清县| 安阳县| 昭觉县| 赤城县| 东海县| 沅江市| 吉木乃县| 扬州市| 长垣县| 安陆市| 德化县| 辽中县| 临夏县| 广东省| 新沂市| 台东市| 隆安县| 胶南市| 旅游| 正阳县| 贺兰县| 大冶市| 台湾省| 驻马店市| 华坪县| 北川| 外汇| 新平| 绥阳县| 台山市|