两个人的电影免费视频_国产精品久久久久久久久成人_97视频在线观看播放_久久这里只有精品777_亚洲熟女少妇二三区_4438x8成人网亚洲av_内谢国产内射夫妻免费视频_人妻精品久久久久中国字幕

用于地震數(shù)據(jù)采集的自適應濾波裝置和方法

文檔序號:6017456閱讀:233來源:國知局
專利名稱:用于地震數(shù)據(jù)采集的自適應濾波裝置和方法
技術領域
本發(fā)明大體上涉及地質勘測,特別地,涉及一種用于采集和處理地震數(shù)據(jù)的裝置和方法。
背景技術
用于地震勘探的傳統(tǒng)地震檢波器、加速計和水下地震檢波器系統(tǒng)通常有多個傳感器,用來產(chǎn)生表示通常為地震波形式的聲能的信號。地震波通常由能源,例如振動車、爆炸物或者在水下地震檢波器系統(tǒng)情況下的氣槍產(chǎn)生。然后這些地震信號傳輸?shù)讲杉?轉換電路。對來自一個或多個遠程地震傳感器(水下地震檢波器、地震檢波器或其它地震傳感器)的模擬信號進行采樣,并通過采集/轉換電路轉換成一系列數(shù)值。通??梢耘渲迷摬杉?轉換電路以實現(xiàn),例如,調(diào)節(jié)采樣率、更改任何數(shù)字濾波或其它數(shù)字信號處理參數(shù)或者進行診斷。
將一個或多個這樣的采集/轉換電路連接到一個數(shù)據(jù)收集單元。每個數(shù)據(jù)收集單元收集來自所有地震傳感器的一系列數(shù)值,這些地震傳感器連接到所有連在該數(shù)據(jù)收集單元上的采集/轉換單元。該數(shù)據(jù)收集單元將這些數(shù)據(jù)傳輸給一個可以包括地震記錄設備的系統(tǒng)控制器。
聲地震波包括多種能量信號。它們包括通過區(qū)分到達檢測器的時間(區(qū)分速度)來識別的P & Shear信號。對于任意給定的地質事件,該信號包含不同的頻率。因此,它們在時限的垂直和水平時間窗口之間不互相關聯(lián)。
其它地震信號是不希望得到的信號,如噪聲。噪聲的一種形式是地滾(ground roll)。地滾噪聲,或簡稱地滾,在相同時間以相同頻率到達垂直和水平分量。這些信號在90度互相極化。因此,噪聲在時限的垂直和水平記錄之間互相關聯(lián)。
當聲源開始活動時產(chǎn)生地滾。從淺阻抗以及目標深度的阻抗反射聲能。與從目標深度的反射相比,淺反射有較大的振幅,并且有時同時到達表面檢測器。噪聲的存在嚴重影響了地震數(shù)據(jù)采集,甚至許多傳統(tǒng)系統(tǒng)關于噪聲特性而不是地震信號的需要來設計。例如,通常采用陣列配置地震檢波器的應用以盡量消除地滾的影響,其中是根據(jù)相比于信號事件(至少頻帶限制)這些噪聲信號頻率不同、有較大振幅并且是分散的。單獨使用陣列以減少地滾在消除地滾方面只能夠獲得部分成功。此外,如果使用傳統(tǒng)濾波方法例如頻率波數(shù)(FK)濾波來消除地滾,它緩慢的速度要求嚴格的空間采樣。
極化濾波器是一大類基于地滾橢圓微粒運動的檢測和分離的濾波器。由于地滾的分散特性以及許多假定因當?shù)氐刭|條件而失效,一些極化濾波器只獲得了有限的成功。
自適應系統(tǒng)通常特征在于能夠以時變的方式使用動態(tài)噪聲預計輸入來調(diào)節(jié)性能。已知的自適應濾波技術不包括從單個多軸傳感器使用垂直和水平分量的能力。而且,普通的濾波方法無法從其它形式的能量中區(qū)分出地滾,并且可能在過濾地滾的同時過濾所需的信號。最后,由于對空間采樣和接收機陣列的依賴,普通系統(tǒng)不能提供逐站減少地滾的能力。

發(fā)明內(nèi)容
本發(fā)明通過提供一種有多個區(qū)域單元的系統(tǒng),致力于解決上述在傳統(tǒng)地震數(shù)據(jù)采集系統(tǒng)中遇到的難題,其中所述區(qū)域單元包括一個或多個地震源、一個或多個用于接收由震源產(chǎn)生并且在地層中反射的聲能的多軸地震傳感器、一個數(shù)據(jù)采集單元、一個接收機/發(fā)射機和一個遠離一個或多個區(qū)域單元的中心控制處理器。處理器根據(jù)初始濾波器系數(shù),使用自適應濾波算法來處理垂直和水平傳感器輸出,其中所述濾波器系數(shù)根據(jù)逐個采樣而改變作為噪聲變化的估計值。在首選實施例中,初始濾波器系數(shù)保存在控制器的存儲設備中。
依照本發(fā)明的方法包括兩個步驟。第一步是濾波器設計,其中根據(jù)垂直和水平分量之間的瞬間相位差異來定義偏移依賴應用窗口。如果需要,可以使用帶通濾波器來選擇頻帶。第二步是僅對偏移窗口之內(nèi)的數(shù)據(jù)和進出窗口的時間錐形應用最小均方(LMS)自適應濾波算法。
一方面,本發(fā)明用于減少或取消地震檢波器陣列以減弱地滾。本發(fā)明的系統(tǒng)不依賴與許多接收機陣列設計相關的方位設計參數(shù)。本發(fā)明提供了使用所需信號作為主要設計因素而不是噪聲信號來設計系統(tǒng)的能力。
本發(fā)明提供了一種有效的地滾濾波器而無需考慮陣列、空間采樣,同時減少對噪聲極化的依賴。提供過濾的地震信號的裝置包括一個接地以接收聲波的檢測器,該檢測器提供一個輸出,包含表明接收聲波的方向分量的第一信號和第二信號。將一個處理器連接到該檢測器以接收所述第一和第二信號,所述處理器處理對第一信號應用自適應濾波算法的命令,所述自適應濾波包括一個表明結合了第二信號的第一信號的選擇部分的反饋信號。
檢測器最好是MEMS多軸加速計,可選的帶通濾波器可用于過濾自適應濾波輸出。作為替代,帶通濾波器可以是其它任意的有用的通帶濾波器,例如帶阻濾波器。


通過附圖,以及考慮下述的描述,將能夠最好地理解本發(fā)明的新穎特性以及發(fā)明本身,其中相似的附圖標記表示相似的部分,附圖包括圖1是一種依照本發(fā)明的地震數(shù)據(jù)采集系統(tǒng)的實施例的系統(tǒng)示意圖;圖2是用于如圖1所示系統(tǒng)中的多軸傳感器;圖3是圖2的傳感器模塊的實施例;圖4是圖3的傳感器模塊的替代實施例;
圖5A表示一種傳統(tǒng)的自適應濾波器;圖5B是一種依照本發(fā)明的自適應濾波器;圖6是一種依照本發(fā)明的方法的流程圖;以及圖7是一種依照本發(fā)明的替代方法的流程圖。
具體實施例方式
圖1是一種依照本發(fā)明的地震勘測系統(tǒng)的實施例的系統(tǒng)示意圖。系統(tǒng)100包括一個用于向系統(tǒng)100中的一個或多個遠程單元發(fā)送命令并從其接收數(shù)據(jù)信號的主控制和記錄單元102。如圖所示,主控制和記錄單元102是個如對技術熟悉的人所知的“車”或“儀器車”的移動單元。作為替代,主控制和記錄單元也可以放置在任何車輛或載體上,或者半永久性地安裝在預定位置。
主控制和記錄單元102(在下文中簡稱為中心控制器102)連接到一個或多個接收機單元104。接收機單元可以與中心控制器102一起放置,或者放置在例如擴展的區(qū)域中。通過雙向通信鏈路124將接收機104連接到多個數(shù)據(jù)采集單元(“DAUs”)106a-106b等。數(shù)據(jù)采集單元106包括用于接收來自地震傳感器108a的陣列108的信號的電路。每個地震傳感器最好能夠提供一個表示感應到的地震活動的數(shù)字輸出信號。在替代實施例中,傳感器輸出信號是模擬的,在DAU中執(zhí)行的A/D轉換。
系統(tǒng)100還包括一個地震源110,例如一個或多個振動車的全體,爆炸物(例如炸藥),氣槍或類似物。震源110和中心控制器102經(jīng)由雙向通信鏈路122互相連接。通信鏈路122可以是模擬或者數(shù)字的,通信的數(shù)據(jù)最好包括控制設置、發(fā)射命令、位置(相關的和/或全面的)狀態(tài)信息、震源可用性、震源等級和任何其它在管理勘測方面有用的信息。
在操作上,震源110將聲源傳進地球表面114,從而設置了大體向下輻射進地球的聲波116。聲波116將在阻抗變化例如地層邊界120處部分反射,從而產(chǎn)生大體上傳播回表面114的反射波118。
通過獨立傳感器108a在檢測器陣列108處檢測反射波118。傳感器可以是適合于檢測聲能的任意型號傳感器,例如地震檢波器、擴音器、基于加速計的檢測器、在濕地條件下的水下地震檢波器以及類似物。作為替代,陣列可以是單個元件或者包括許多元件。檢測器108a可以提供模擬或數(shù)字輸出信號,可以組合陣列108中的檢測器,例如剛才描述的檢測器型號和/或輸出信號的組合。
在首選實施例中,檢測器108a是使用微機電系統(tǒng)(MEMS)加速計的多軸檢測器。這些MEMS檢測器在提交于2001年9月12日的U.S.申請09/936,634中有詳細討論,該申請是基于提交于2000年3月17日的PCT申請PCT/US00/07310。因此,在此引用了這兩個申請的全部內(nèi)容。
來自陣列108的輸出信號可以包括來自獨立檢測器108a的信號或者包括每個檢測器一個信號。在實施例中,使用在每個數(shù)據(jù)采集單元106a-106b中的處理器(未單獨畫出)來預處理信號。然后將數(shù)據(jù)傳輸?shù)浇邮諜C104。接收機104簡化地表示成單個接收機。熟悉技術的人將能理解實踐中可以配置任意數(shù)量的接收機以接收來自傳感器陣列108的數(shù)據(jù)。可以依照本實施例來處理信號,或者在系統(tǒng)100的任何適當位置,例如在接收機104或中心控制器102之內(nèi)使用依照本發(fā)明的自適應濾波器,來處理信號。
在一個實施例中,中心控制器102和接收機104之間的連接112,接收機104和DAUs106之間以及每個DAU和傳感器陣列108之間的連接124,和中心控制器102和震源110之間的連接122都是導線112a和112b、122a、122b、124a、124b。這些連接統(tǒng)稱為連接112/124,或者當需要時,將使用各自的附圖標記來引用單條連接。選擇包括連接112/124在內(nèi)的這些導線以充分傳輸特定應用所要求的電源和數(shù)據(jù)。作為替代,連接112/124可以是能夠傳輸電源和數(shù)據(jù)信號的任何適當連接。在完整的信號路徑上,一條連接,例如連接112a,不必是與任何其它連接例如112b相同的連接類型。典型連接可以是導線、射頻、光纖或電磁的任意組合。作為替代的,當向任意特殊單元例如接收機中的電池或發(fā)電機提供電源作為震源110時,連接112/124也可以包括一個無線電鏈路。在這些替代實施例中,收發(fā)機(接收機和發(fā)射機的組合)安裝在連接的單元中。中心控制器102和接收機104之間的連接112/124用虛線表示以標明該點。
參見圖2,傳感器裝置108a的實施例包括一個連接到第一封端帽210、第二封端帽215和傳感器模塊305的外殼205。外殼205通過連接元件例如機械夾具310、315、320和325連接到封端帽210和215。
密封元件330a-d密封外殼205和第一封端帽210之間的連接,密封元件335a-d密封外殼205和第二封端帽215之間的連接。密封元件330a-d和335a-d可以是能夠被壓縮在預定位置以進行密封的彈性環(huán)。
外殼205最好包括一個空腔340和平坦表面345。在實施例中,外殼205是由高強度材料制造的金屬管以提供堅固的壓力容器。
傳感器模塊305在外殼205的空腔340中由平坦表面345支撐,并通過PC板連接355連接到第一封端帽210。
參見圖3,傳感器模塊305最好包括傳感器組件405a-c和襯底410。每個傳感器組件405a-c包括對應的感應軸415、420和425。每個感應軸415、420和425分別與x軸、y軸和z軸大致平行。
參見圖4,傳感器模塊305的一個替代實施例包括多個傳感器組件405、襯底410和單片組件1902。傳感器組件405連接到單片組件1902。在實施例中,傳感器模塊305包括第一傳感器組件405a、第二傳感器組件405b和第三傳感器組件405c。第一傳感器組件405a包括一個感應軸415。該感應軸415與x軸大致平行。第一傳感器組件405a連接到單片電路組件1902以保持感應軸415平行于x軸。第二傳感器組件405b包括一個感應軸420。該感應軸420與y軸大致平行。第二傳感器組件405b連接到單片電路組件1902以保持感應軸420平行于y軸。第三傳感器組件405c包括一個感應軸425。該感應軸425與z軸大致平行。第三傳感器組件405c連接到單片電路組件1902以保持感應軸425平行于z軸。
本發(fā)明的裝置提供了一種自適應濾波器以在信號處理之前消除不需要的信號,例如地滾。圖5A表示一種傳統(tǒng)的自適應濾波器。類似的自適應濾波器用于消除來自地震信號的不一致的分散噪聲,例如機械噪聲。在Crews等人的U.S.專利4,890,264中描述了這種濾波器應用的例子。傳統(tǒng)濾波器的另一應用是消除地滾噪聲。在Widrow的U.S.專利4,556,962中可以看到這種地滾濾波器應用的例子。Widrow設備利用一個濾波輸入信號n’(k),其來自用于檢測地滾的獨立地震檢波器或者直接來自聲源的信號。術語n’(k)代表噪聲輸入,x(k)是自適應濾波器輸入的普通表示,并且其輸出信號用y(k)表示。第二系統(tǒng)輸入是所需的信號d(k),所需信號d(k)是噪聲n(k)和已知聲傳感器檢測到的聲信號s(k)的組合。然后加總輸入d(k)和自適應濾波器轉換的輸出y(k)以產(chǎn)生表示信號誤差的輸出e(k)。然后將誤差信號e(k)反饋回自適應濾波器以更新濾波輸出y(k)來進行連續(xù)的數(shù)據(jù)處理。在這一典型系統(tǒng)中,通過將噪聲輸入與來自地震檢波器陣列的平均信號相加,消除了地滾。在所描述的所有例子中,地震檢波器是單向的檢測器。
圖5B是一種依照本發(fā)明的自適應濾波器500。該自適應濾波器包括向濾波器提供垂直和水平(或放射狀radial)分量的多軸傳感器輸入信號。地滾表現(xiàn)為橢圓逆向微粒運動,這導致能量到達如上面所述以及圖2-4所示的那些三分量地震檢測器的垂直和水平分量。而且,地滾的極化特性允許地滾和可能與水平和垂直分量之間相關的其它信號和/或噪聲之間有區(qū)別。
地滾的問題在于它速度緩慢并且分散。這些因素共同使消除地滾的處理變得麻煩。依照本發(fā)明的自適應濾波器提供了利用分量之間的關聯(lián)的單個接收臺濾波器,以使用非固定濾波器來識別和消除噪聲。
當?shù)貪L的振幅大到足以成為地震數(shù)據(jù)的考慮因素時,在發(fā)射記錄中可以清楚地看到作為高振幅、低頻率、分散的噪聲序列。然而,還有其它的波類型,例如能夠幾乎同時到達并且被誤認為地滾的純剪切折射到達波和樂甫波(Love waves)。本發(fā)明的濾波器設計在濾波之前及時隔離了地滾以及頻率,以便只過濾被地滾明確污染的地震數(shù)據(jù)。
圖5B表示檢測的垂直數(shù)據(jù)使用信號加噪聲s(k)+n(k)作為系統(tǒng)輸入d(k)。檢測放射狀波,輸入是對系統(tǒng)自適應濾波器的噪聲信號n’(k)。在本發(fā)明中,通過產(chǎn)生一個使用時間參數(shù)來隔離特定放射狀數(shù)據(jù)的“窗口”,來進一步凈化濾波器輸出??蛇x的帶通(BP)濾波器可用于使數(shù)據(jù)頻帶變小。因此,在使用d(k)加總之前,轉換的輸出y(k)更加精確。這導致了反饋回自適應濾波器的更小的系統(tǒng)誤差輸出e(k)。
通過顯示來自一個單個三分量傳感器的垂直和水平數(shù)據(jù)的瞬間相位之間的不同區(qū)域來完成地滾的區(qū)分。產(chǎn)生的差異部分提供了簡易識別極化噪聲的能力。
上述的和圖1-5B所示的系統(tǒng)100可用于本發(fā)明的首選方法中。聲源可用于向地球施加聲能信號。聲能也會產(chǎn)生與從所需或目標深度反射回的聲信號基本同時到達多軸傳感器的地滾噪聲信號。多軸檢測器將接收到的聲信號轉換成垂直和水平輸出信號。將來自多軸檢測器108a的垂直和水平信號發(fā)送到放置在例如數(shù)據(jù)采集單元中的處理器中。數(shù)據(jù)采集單元包括用于處理應用于垂直和水平輸出信號的濾波算法的存儲指令。濾波算法是一個自適應濾波器,通過該自適應濾波器來根據(jù)改變輸入條件來改變輸入噪聲系數(shù)。初始濾波系數(shù)可以保存在存儲設備中以供處理器應用于自適應濾波器。存儲設備和處理器無需放在一起。在實施例中,處理器和存儲設備位于中心控制器中。作為替代,系數(shù)可以放在DAU的存儲設備中或者任何其它所需的位置。同樣,處理器無需放在中心控制器中。處理器可以放在任何所需的位置,例如DAU中。下面參照圖6和7進一步描述該方法。
現(xiàn)在參見圖6,顯示的是依照本發(fā)明采用自適應濾波器以從選擇信號中減弱地滾的首選方法600的流程。在首選實施例中,本方法的指令包含在上述的和圖1所示的系統(tǒng)中的存儲器中。由系統(tǒng)中的處理器來執(zhí)行這些指令以從包括所需地震反射波以及地滾噪聲的垂直信號中消除地滾。
如上所述,使用從多軸地震傳感器的垂直和水平輸出接收的信號確定窗口參數(shù)602。該參數(shù)表示有特殊地滾特征的那些信號。例如,當垂直和水平信號相位差異大約為90度時。
向自適應濾波器發(fā)送傳感器水平輸出信號604,應用可選帶通濾波器606以選擇特殊的頻帶。然后,比較水平信號與窗口參數(shù)608,將那些符合窗口參數(shù)標準的頻率與包括地滾噪聲的垂直傳感器信號相加610。產(chǎn)生的輸出重復地發(fā)送回自適應濾波器以進行進一步濾波,以便最終輸出的垂直信號基本上沒有地滾噪聲。
現(xiàn)在參見圖7,顯示了依照本發(fā)明的替代的自適應濾波方法700。該方法基本上與圖6的方法相似,除了反轉了垂直和水平傳感器輸出信號。當需要時,該方法提供了過濾的水平信號。使用從多軸地震傳感器的垂直和水平輸出接收的信號來確定窗口參數(shù)702。該參數(shù)表示有特殊地滾特征的那些信號。例如,當垂直和水平信號相位差異大約為90度時。
向自適應濾波器發(fā)送傳感器垂直輸出信號704,應用可選帶通濾波器706以選擇特殊的頻帶。然后,比較垂直信號與窗口參數(shù)708,將那些符合窗口參數(shù)標準的頻率與包括地滾噪聲的水平傳感器信號相加710。產(chǎn)生的輸出重復地發(fā)送回自適應濾波器以進行進一步濾波,以便最終輸出的水平信號基本上沒有地滾噪聲。
前述的描述是出于說明和解釋的目的而針對于本發(fā)明的特定實施例。然而,對于熟悉技術的人來說,在不背離本發(fā)明范圍的條件下,有可能對上述提出的實施例進行修改和變化是顯而易見的。附上的權利要求書意在解釋包括所有這種修改和變化。
權利要求
1.一種在地震勘測系統(tǒng)中過濾地震信號的方法,該方法包括a)使用一個檢測器來接收一個聲波;b)從所述檢測器發(fā)送第一信號和第二信號,所述第一和第二信號表示接收到的聲波的方向分量;以及c)對所述第一信號應用一個自適應濾波算法,所述自適應濾波包括一個反饋信號,該反饋信號表明結合了第二信號的第一信號的一個選中部分。
2.權利要求1的方法,進一步包括,至少部分地通過對所述自適應濾波算法的輸出應用一套預定窗口參數(shù)來確定所述第一信號的選中部分。
3.權利要求1的方法,進一步包括,至少部分地通過對所述自適應濾波算法的輸出應用頻帶濾波器來確定所述第一信號的選中部分。
4.權利要求3的方法,其中所述頻帶濾波器是一個帶通濾波器。
5.權利要求1的方法,其中所述第一和第二信號包括有大約90度相位差異的信號。
6.權利要求1的方法,其中所述第一信號是水平信號,所述第二信號是垂直信號。
7.權利要求1的方法,其中所述算法用于從所述地震信號中消除地滾。
8.權利要求2的方法,其中所述預定窗口參數(shù)包括相位和定時參數(shù)。
9.權利要求1的方法,其中應用所述自適應濾波器進一步包括,至少部分地根據(jù)初始濾波器系數(shù)來應用自適應濾波器。
10.權利要求9的方法,進一步包括,當向連續(xù)接收的聲波應用自適應濾波時,更新所述初始濾波器系數(shù)。
11.權利要求2的方法,其中確定窗口參數(shù)包括至少部分地根據(jù)聲波的垂直和水平分量之間的相位差異來選擇偏移依賴應用窗口。
12.權利要求1的方法,其中所述檢測器包括一個多軸檢測器。
13.權利要求12的方法進一步包括使用一個MEMS多軸檢測器來接收聲波。
14.一種用于提供過濾后的地震信號的裝置,包括a)一個連接到地球以接收聲波的檢測器,該檢測器提供一個包括第一和第二信號的輸出,所述第一和第二信號表示接收到的聲波的方向分量;b)一個連接到所述檢測器以接收第一和第二信號的處理器,所述處理器處理指令以向第一信號應用自適應濾波算法,所述自適應濾波包括一個反饋信號,該反饋信號表明結合了第二信號的第一信號的選中部分。
15.根據(jù)權利要求14的裝置,其中所述檢測器包括一個多軸檢測器。
16.根據(jù)權利要求15的裝置,其中所述多軸檢測器包括一個多軸MEMS加速計。
17.根據(jù)權利要求14的裝置進一步包括一個頻帶濾波器以過濾自適應濾波器的輸出。
18.根據(jù)權利要求17的裝置,其中所述頻帶濾波器包括一個帶通濾波器。
19.根據(jù)權利要求14的裝置,其中至少部分地通過對所述自適應濾波算法的輸出應用一套預定窗口參數(shù)來確定所述第一信號的選中部分。
20.根據(jù)權利要求14的裝置,其中所述第一和第二信號包括有大約90度的相位差異的信號。
21.根據(jù)權利要求14的裝置,其中所述第一信號是水平信號,所述第二信號是垂直信號。
22.根據(jù)權利要求14的裝置進一步包括一個裝置輸出信號,其中應用所述算法以從所述地震信號中消除地滾。
23.根據(jù)權利要求19的裝置,其中所述預定窗口參數(shù)包括相位和定時參數(shù)。
24.根據(jù)權利要求14的裝置進一步包括,一個其中保存了至少部分用于自適應濾波器應用的初始濾波器系數(shù)的存儲設備。
25.根據(jù)權利要求24的裝置,其中當向連續(xù)接收的聲波應用自適應濾波器時,所述處理器更新所述初始濾波器系數(shù)。
26.根據(jù)權利要求19的裝置,其中所述窗口參數(shù)包括至少部分地根據(jù)聲波的垂直和水平分量之間相位差異的偏移依賴應用窗口。
全文摘要
本發(fā)明提供了一種基本上消除了在地震數(shù)據(jù)采集時遇到的地滾噪聲的自適應濾波方法(圖5B)。(圖1)還提供了一種包括一個連接到一個多軸傳感器設備的自適應濾波電路(圖5B)的用于采集地震數(shù)據(jù)的裝置 (305) 。
文檔編號G01V1/18GK1646938SQ03808905
公開日2005年7月27日 申請日期2003年3月19日 優(yōu)先權日2002年3月20日
發(fā)明者拉塞爾·A·卡皮烏斯, 加里·A·格雷斯 申請人:輸入輸出公司
網(wǎng)友詢問留言 已有0條留言
  • 還沒有人留言評論。精彩留言會獲得點贊!
1
正镶白旗| 德昌县| 晋城| 抚州市| 兰西县| 达孜县| 永嘉县| 会东县| 疏附县| 紫阳县| 高平市| 尼勒克县| 曲阳县| 淮安市| 厦门市| 遂溪县| 壤塘县| 武定县| 申扎县| 宜阳县| 象州县| 静宁县| 舟曲县| 华阴市| 广安市| 淅川县| 新蔡县| 淄博市| 行唐县| 于田县| 贞丰县| 西充县| 商河县| 合川市| 丹凤县| 北海市| 厦门市| 淮阳县| 额尔古纳市| 高清| 临清市|