• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于加權(quán)一致性的蟻群算法在斷層檢測中的應(yīng)用

    2016-11-08 03:04:22劉財(cái)劉海燕彭沖張營劉瓊勾福巖
    地球物理學(xué)報(bào) 2016年10期
    關(guān)鍵詞:蟻群梯度一致性

    劉財(cái), 劉海燕 *, 彭沖, 張營, 劉瓊, 勾福巖

    1 吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院, 長春 130026 2 山西師范大學(xué)地理科學(xué)學(xué)院, 山西臨汾 041004 3 中國石化東北油氣分公司, 長春 130062 4 中國石油吐哈油田公司勘探開發(fā)研究院, 新疆哈密 839009

    ?

    基于加權(quán)一致性的蟻群算法在斷層檢測中的應(yīng)用

    劉財(cái)1, 劉海燕1 *, 彭沖2, 張營3, 劉瓊4, 勾福巖1

    1 吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院, 長春130026 2 山西師范大學(xué)地理科學(xué)學(xué)院, 山西臨汾041004 3 中國石化東北油氣分公司, 長春130062 4 中國石油吐哈油田公司勘探開發(fā)研究院, 新疆哈密839009

    突出地震數(shù)據(jù)振幅空間不連續(xù)性的斷層增強(qiáng)屬性體是斷層自動解釋的基礎(chǔ),而如何壓制噪聲、地層殘余響應(yīng)、角度不整合接觸等偽斷層的影響是斷層自動檢測的關(guān)鍵問題之一.在指紋圖像邊緣檢測處理中,一致性是對指示局部方向場信息各向異性強(qiáng)度的一個度量,本文將其作為一種新的斷層增強(qiáng)屬性引入到地震數(shù)據(jù)處理中.為了更有效地壓制噪聲,本文通過對一致性進(jìn)行加權(quán)處理,再結(jié)合蟻群算法,提出了新的斷層自動檢測方法——基于加權(quán)一致性的蟻群算法.理論地震模型和實(shí)際地震數(shù)據(jù)的處理結(jié)果表明,對于相同的斷層檢測過程,對比常規(guī)的基于方差屬性的蟻群算法和基于C3相干屬性的蟻群算法,本文所提方法不僅更有效地壓制了噪聲和地層殘余響應(yīng)產(chǎn)生的偽斷層信息,而且更完整地展現(xiàn)了斷層的延伸長度,同時具有更高的穩(wěn)定性和計(jì)算效率.

    一致性; 加權(quán)一致性; 蟻群算法; 斷層檢測

    1 引言

    在地震數(shù)據(jù)中,斷層主要表現(xiàn)為反射層位的位移或中斷、橫向上振幅不連續(xù),斷層兩側(cè)的反射振幅存在明顯差異,這些特征在地震剖面上可以被肉眼識別,這也是人工識別斷層的基礎(chǔ)(嚴(yán)哲,2010).但是人工識別斷層效率低,主觀性強(qiáng),解釋能力局限.

    所以,學(xué)者們針對斷裂系統(tǒng)的精確解釋,提出并應(yīng)用了很多斷層的自動或半自動檢測技術(shù).Bahorich和Farmer(1995)提出了在三維斷層解釋技術(shù)中占據(jù)重要地位的基于互相關(guān)的第一代相干算法;Marfurt等(1998)提出了基于多道相似性的第二代相干算法;Gersztenkorn和Marfurt(1999)提出了基于特征結(jié)構(gòu)的第三代相干算法.相干體技術(shù)在很大程度上提高了斷層的解釋效率,大部分?jǐn)鄬幼詣訖z測技術(shù)也主要是以這三代相干體技術(shù)為核心所進(jìn)行的改進(jìn).然而,傳統(tǒng)的基于滑動窗口的相干體屬性對處理數(shù)據(jù)信息中的局部非平穩(wěn)特征有一定的難度(劉洋等,2014).為了解決以上問題,學(xué)者們陸續(xù)開發(fā)了一些新的斷層自動檢測技術(shù),其中Randen等(2001)提出了“人工螞蟻”斷層自動檢測技術(shù),并由Pedersen等(2002,2003)進(jìn)行了更深層次的分析闡述,該技術(shù)利用人工螞蟻對三維地震數(shù)據(jù)體生成的突出振幅空間不連續(xù)性的斷層增強(qiáng)屬性體(傾角方位屬性、混沌屬性、方差屬性等)進(jìn)行處理,從而提取斷層信息,該技術(shù)在2002年以Ant tracking模塊被引入到Petrel軟件中,得到廣泛商用.Dorn等(2005)提出基于相干體數(shù)據(jù)的斷層自動提取技術(shù),該技術(shù)融合了地質(zhì)先驗(yàn)知識的信號處理技術(shù),后被應(yīng)用到商業(yè)軟件Epos中,但該技術(shù)只是對相干體斷層響應(yīng)進(jìn)行了加強(qiáng),對斷層響應(yīng)的連續(xù)性改善作用不大.趙偉(2009)利用蟻群算法對相干切片圖像進(jìn)行聚類處理,實(shí)現(xiàn)了“簡單”斷層線的提取和重構(gòu),但結(jié)果依賴相干體的斷層檢測能力,且計(jì)算效率低.李香臣(2010)利用蟻群優(yōu)化聚類算法和基于圖搜索式蟻群優(yōu)化算法對構(gòu)造特征進(jìn)行了增強(qiáng).嚴(yán)哲(2010)提出方向約束的蟻群追蹤算法,進(jìn)一步改進(jìn)之后,提出了基于蟻群算法的自動斷層追蹤的方法(Yan et al.,2013),該方法是將相干體數(shù)據(jù)作為斷層增強(qiáng)屬性體,再利用蟻群算法對相干體數(shù)據(jù)進(jìn)行處理以達(dá)到檢測斷層的目的,與Petrel軟件中Ant tracking模塊的斷層追蹤結(jié)果對比,該方法檢測到的斷層連續(xù)性高,但結(jié)果好壞也直接依賴于相干體算法對斷層的檢測能力.

    斷層增強(qiáng)屬性體會沿著斷層面產(chǎn)生響應(yīng),這是斷層自動檢測的基礎(chǔ),但也夾雜了噪聲、地層殘余響應(yīng)、角度不整合接觸等假象.因此,為了提高斷層檢測方法的信噪比,本文引入了一致性的概念.一致性由Kass和Witkin(1987)提出,是指在對指紋圖像邊緣檢測處理中,對指示局部方向場信息的各向異性強(qiáng)度的一個度量,其計(jì)算過程不涉及滑動窗口,計(jì)算結(jié)果受人為選取參數(shù)約束小,且不受類似地震同相軸的指紋脊谷線形狀的影響,即若將其引用到地震數(shù)據(jù)處理,其計(jì)算結(jié)果將不會受地層傾角的影響.因此,本文將其作為一種新的斷層增強(qiáng)屬性應(yīng)用于地震數(shù)據(jù)斷層檢測,并提出加權(quán)一致性的概念,來提高一致性屬性的抗噪性,再結(jié)合蟻群算法構(gòu)成了新的斷層檢測方法.在理論模型和實(shí)際數(shù)據(jù)試算階段,對于相同的斷層檢測過程,本文所提方法在抗噪性、收斂性及壓制地層殘余響應(yīng)等方面均明顯優(yōu)于常規(guī)的基于方差屬性的蟻群算法和基于C3相干屬性的蟻群算法.

    2 方法原理

    2.1一致性

    在指紋圖像中,灰度沿指紋脊谷線垂向方向變化最快,沿走向方向變化最慢,這種特征被稱為各向異性(Kass and Witkin,1987).該特征可用于求取指紋脊谷線的局部方向場信息,Kass和Witkin(1987)將衡量指示局部方向場的各向異性強(qiáng)度的度量稱為“一致性”.

    地震數(shù)據(jù)實(shí)際是對地下沉積地層的反射波成像,地震記錄中的同相軸類似于指紋圖像中的脊谷線,因此借鑒一致性算法,對地震數(shù)據(jù)中指示局部地層走向的各向異性強(qiáng)度作以衡量.一致性的求取過程如下:

    (1)

    在極坐標(biāo)系下,公式(1)可表示為:

    (2)

    其中,Gρ為極徑,表示梯度向量的模長,Gφ為極角,表示梯度向量的方位角.則直角坐標(biāo)系下的梯度向量又可表示為:

    (3)

    (4)

    然后,提取一個以點(diǎn)(x,y)為中心的包含n×n個數(shù)據(jù)點(diǎn)的窗口B,將該數(shù)據(jù)塊的平方梯度向量表示為:

    (5)

    其中,

    類似指紋圖像,地震記錄中每個數(shù)據(jù)點(diǎn)的梯度向量總是指向灰度變化最快的方向,并與地震同相軸走向垂直.同一地震同相軸具有兩側(cè)邊界線,分布在兩側(cè)邊界線上數(shù)據(jù)點(diǎn)的梯度向量呈對稱分布,即指向相反.平方梯度向量方向角是梯度向量方向角的兩倍,這樣指向相反的向量變成了指向相同的向量.所以對于同一條同相軸,兩側(cè)邊界上點(diǎn)的平方梯度向量指向一致.

    為了衡量數(shù)據(jù)塊B的各向異性強(qiáng)度,即數(shù)據(jù)塊內(nèi)各點(diǎn)的平方梯度向量方向的集中程度,Kass和Witkin(1987)提出一致性概念:

    (6)

    定義CohB為數(shù)據(jù)塊B中心點(diǎn)(x,y)的一致性,其中,||表示向量的模.

    如果數(shù)據(jù)塊B內(nèi)各點(diǎn)的平方梯度向量指向都相同,則該局部范圍內(nèi)各向異性指示性很強(qiáng),與平方梯度向量方向垂直的地震同相軸走向也一致,表明地層連續(xù),此時,所有向量模的和等于所有向量和的模,CohB等于1;如果數(shù)據(jù)塊B內(nèi)各點(diǎn)的平方梯度向量指向在(0-2π)范圍內(nèi)均勻分布,則該局部范圍內(nèi)各向異性特征極不明顯,很難獲取地震同相軸走向信息,表明地層不連續(xù)或存在強(qiáng)隨機(jī)噪聲,此時,所有向量模的和遠(yuǎn)大于所有向量和的模,CohB趨向于0;其他的情況,CohB介于0和1之間.本文給出的圖件中一致性均用Coh表示.

    對圖1a無噪聲理論合成地震記錄計(jì)算一致性,結(jié)果如圖1c所示,可以看出,在地層連續(xù)處,一致性數(shù)值較高,而在斷層處及地層不整合接觸位置,一致性數(shù)值較低,這樣就突出了地層的不連續(xù)性,該結(jié)果表明一致性數(shù)據(jù)可以作為新的斷層增強(qiáng)屬性體應(yīng)用于斷層檢測.對加入強(qiáng)隨機(jī)噪聲的合成地震記錄(圖1b)計(jì)算一致性,結(jié)果如圖1d所示,可以看出,一致性數(shù)據(jù)在強(qiáng)噪聲干擾下,對地層的橫向不連續(xù)性識別效果非常差.所以,為了提高抗噪性,本文提出了加權(quán)一致性的概念.

    2.2加權(quán)一致性

    接下來計(jì)算目標(biāo)數(shù)據(jù)塊V的一致性.從圖2可以看出,目標(biāo)數(shù)據(jù)塊V在相鄰的四個鄰域D中都有被包含,只是在不同的數(shù)據(jù)塊D中所處的位置不同,分別處在四個不同方向,每個基礎(chǔ)數(shù)據(jù)塊Γ的平方梯度向量由公式(5)決定,一致性由公式(6)決定.則數(shù)據(jù)塊Dj的平方梯度向量為:

    圖1 理論地震數(shù)據(jù)的一致性(a) 二維合成地震記錄; (b) 加入隨機(jī)噪聲之后的合成地震記錄;(c) 圖1a的一致性結(jié)果; (d) 圖1b的一致性結(jié)果.Fig.1 Coherence for synthetic seismic data(a) 2D synthetic data; (b) Synthetic data with random noise added; (c) Coherence for Fig.1a; (d) Coherence for Fig.1b.

    圖2 3×3個基礎(chǔ)數(shù)據(jù)塊Γ構(gòu)成的數(shù)據(jù)塊M:中心數(shù)據(jù)塊V為目標(biāo)數(shù)據(jù)塊,D1,D2,D3,D4分別是以數(shù)據(jù)塊V為中心的四個相互重疊的鄰域(Wang et al.,2005);;;Fig.2 A block of data M consist of 3×3 base data block Γ:the centered data block V is the target block,D1,D2,D3,D4 are the four overlapping neighborhoods(Wang et al.,2005)

    (7)

    (8)

    (9)

    理想情況下,在數(shù)據(jù)塊M內(nèi),如果反映局部地層信息的同相軸連續(xù)且相互平行,那么相互重疊的四個數(shù)據(jù)塊D1,D2,D3,D4就具有非常高的相關(guān)性.此時,在這四個數(shù)據(jù)塊內(nèi),所有數(shù)據(jù)點(diǎn)的平方梯度向量指向一致,那么由公式(7)可知,這四個數(shù)據(jù)塊的平方梯度向量指向也一致,衡量數(shù)據(jù)塊各向異性強(qiáng)度的度量值CohDj均相等,并且與目標(biāo)數(shù)據(jù)塊V的一致性CohV相等.在數(shù)字圖像處理中,利用鄰域信息進(jìn)行冗余計(jì)算可壓制噪聲干擾,并能有效恢復(fù)噪聲圖像(Wangetal.,2005).類似的,如果目標(biāo)數(shù)據(jù)塊V內(nèi)噪聲干擾大,那么數(shù)據(jù)塊V內(nèi)各數(shù)據(jù)點(diǎn)的平方梯度向量方向就比較散亂,CohV值相應(yīng)會變小,不能指示出真正的地層走向信息.此時,就可利用數(shù)據(jù)塊D1,D2,D3,D4的平方梯度向量信息來估算目標(biāo)數(shù)據(jù)塊V的一致性,從而壓制噪聲干擾,提高信噪比.

    將數(shù)據(jù)塊D1,D2,D3,D4組成新的數(shù)據(jù)塊D*,那么數(shù)據(jù)塊D*的平方梯度向量為:

    (10)

    再結(jié)合公式(8)、(9)可知,公式(10)可用基礎(chǔ)數(shù)據(jù)塊的平方梯度向量表示:

    (11)

    其中,M={Ⅰ,Ⅱ,Ⅲ,Ⅳ,Ⅴ,Ⅵ,Ⅶ,Ⅷ,Ⅸ},rΓ表示基礎(chǔ)數(shù)據(jù)塊Γ參與計(jì)算的次數(shù),顯而易見,其值等于該數(shù)據(jù)塊Γ被數(shù)據(jù)塊Dj包含的次數(shù).將系數(shù)rΓ進(jìn)行歸一化:

    (12)

    其中k=∑rΓ.定義ωΓ為覆蓋率(Wangetal.,2005),令ωΓ替換公式(11)中的rΓ,則有:

    (13)

    可見,目標(biāo)數(shù)據(jù)塊V的覆蓋率是最高的.覆蓋率不同表明每個基礎(chǔ)數(shù)據(jù)塊的重要程度不同.將公式(13)定義為數(shù)據(jù)塊M的加權(quán)平方梯度向量.如果任意數(shù)據(jù)塊?!蔒的平方梯度向量的方向和大小都相同,即

    (14)

    則有:

    (15)

    由公式(15)可知,當(dāng)任意數(shù)據(jù)塊?!蔒的平方梯度向量的方向和大小都相同時,目標(biāo)數(shù)據(jù)塊V的平方梯度向量等于數(shù)據(jù)塊M的加權(quán)平方梯度向量.因此,目標(biāo)數(shù)據(jù)塊V的一致性CohV就可根據(jù)基礎(chǔ)數(shù)據(jù)塊Γ的平方梯度向量加權(quán)而得:

    (16)

    當(dāng)公式(14)成立時,有:

    =1.

    (17)

    定義CohwV為目標(biāo)數(shù)據(jù)塊V中心點(diǎn)的加權(quán)一致性.

    公式(17)表明:如果M內(nèi)各基礎(chǔ)數(shù)據(jù)塊Γ的平方梯度向量的方向和大小都相同,CohwV等于1,表明該局部范圍內(nèi)地震同相軸走向一致,說明地層連續(xù).如果M內(nèi)各基礎(chǔ)數(shù)據(jù)塊Γ的平方梯度向量的方向和大小不盡相同,CohwV介于0和1之間.CohwV接近于0時,表明該局部范圍內(nèi)地震同相軸信息很弱,可能是強(qiáng)噪聲的影響,也可能存在不連續(xù)地層.本文給出的圖件中加權(quán)一致性均用Cohw表示.

    圖3a和圖3b分別是對圖1a和圖1b計(jì)算得到的加權(quán)一致性數(shù)據(jù).可見,加權(quán)一致性數(shù)據(jù)相較于一致性數(shù)據(jù),不僅能在一定程度上壓制噪聲干擾,還能弱化地層不整合接觸和地層殘余響應(yīng)所產(chǎn)生的偽斷層信息,使地層的橫向不連續(xù)性更為突出、準(zhǔn)確.

    為了驗(yàn)證加權(quán)一致性數(shù)據(jù)的優(yōu)越性,對實(shí)際地震數(shù)據(jù)(圖4a)計(jì)算一致性(圖4b)和加權(quán)一致性(圖4c).由圖可見,對于實(shí)際地震數(shù)據(jù),加權(quán)一致性數(shù)據(jù)相較于一致性數(shù)據(jù)能更好地壓制隨機(jī)噪聲和地層殘余響應(yīng),突出地層的橫向不連續(xù)性,提高信噪比.但是作為新的斷層增強(qiáng)屬性體,無論是一致性數(shù)據(jù)還是加權(quán)一致性數(shù)據(jù),都是基于一定窗口范圍內(nèi)數(shù)據(jù)塊的計(jì)算,對不連續(xù)性存在著水平方向或垂直方向的彌散作用(嚴(yán)哲,2010),使識別到的地層不連續(xù)性表現(xiàn)為一定寬度且不連續(xù)的條帶,而且加權(quán)一致性并不能徹底地消除噪聲干擾及地層殘余響應(yīng).所以,為了進(jìn)一步提高斷層檢測的質(zhì)量,本文將加權(quán)一致性屬性與蟻群追蹤算法相結(jié)合,提出了新的斷層檢測方法——基于加權(quán)一致性的蟻群算法.

    圖3 理論地震數(shù)據(jù)的加權(quán)一致性(a) 圖1a的加權(quán)一致性結(jié)果; (b) 圖1b的加權(quán)一致性結(jié)果.Fig.3 Weighted coherence of synthetic seismic data(a) Weighted coherence for Fig.1a; (b) Weighted coherence for Fig.1b.

    圖4 實(shí)際地震數(shù)據(jù)一致性與加權(quán)一致性的對比(a) 實(shí)際地震數(shù)據(jù); (b) 圖4a的一致性結(jié)果; (c) 圖4a的加權(quán)一致性結(jié)果.Fig.4 Comparison between coherence and weighted coherence for real seismic data(a) Real seismic data; (b) Coherence for Fig.4a; (c) Weighted coherence for Fig.4a.

    2.3基于加權(quán)一致性的蟻群算法在斷層檢測中的應(yīng)用

    基于加權(quán)一致性的蟻群追蹤算法進(jìn)行斷層檢測的過程主要分為以下幾步:

    (1) 初始化螞蟻位置

    將二維加權(quán)一致性數(shù)據(jù)分為若干個互不相交的m×m的相鄰數(shù)據(jù)塊,在數(shù)據(jù)塊中,螞蟻的初始位置由第i個數(shù)據(jù)點(diǎn)的被選擇概率Pi決定:

    (18)

    其中,Cohwi表示數(shù)據(jù)塊內(nèi)第i個數(shù)據(jù)點(diǎn)的加權(quán)一致性,加權(quán)一致性值越低的點(diǎn)被選中的概率越大,可降低螞蟻在其他區(qū)進(jìn)行無意義追蹤的概率.

    (2) 斷層追蹤

    (19)

    (20)

    可見,在初始信息素相同的情況下,加權(quán)一致性值越低的點(diǎn)被螞蟻選中的概率越大,可使螞蟻的搜索路徑盡快收斂到斷層上.

    (3)判斷終止條件

    螞蟻?zhàn)粉櫟慕K止條件是:

    (21)

    其中,終止條件S(0

    (4)更新信息素

    若當(dāng)前迭代中所有螞蟻的追蹤被終止,即按(22)式對數(shù)據(jù)點(diǎn)(i,j)處的信息素濃度進(jìn)行實(shí)時更新(段海濱,2005):

    (22)

    (23)

    (24)

    c為信息素更新系數(shù).

    當(dāng)前迭代結(jié)束時,記錄所有螞蟻的當(dāng)前所在位置,使其作為下一次迭代時的初始位置.

    (5) 輸出結(jié)果

    不斷重復(fù)(2)—(4)步驟,直到所有迭代次數(shù)完成.常規(guī)蟻群算法是直接輸出信息素剖面圖,但因螞蟻?zhàn)粉櫟碾S機(jī)性,致使不同參數(shù)設(shè)置下所輸出的信息素濃度范圍不同,易導(dǎo)致色標(biāo)不一致.因此,本文對追蹤結(jié)果做歸一化改進(jìn)后再輸出:

    (25)

    其中,τn表示第n次迭代之后的信息素濃度矩陣,C表示利用蟻群算法進(jìn)行斷層追蹤的最終輸出結(jié)果.

    基于加權(quán)一致性的蟻群算法進(jìn)行斷層檢測的流程如圖5所示.其中kmax表示進(jìn)行斷層追蹤的螞蟻總數(shù),nmax表示最大迭代次數(shù).

    圖5 基于加權(quán)一致性的蟻群算法進(jìn)行斷層檢測的流程圖Fig.5 Flow diagram of fault detection using ant colony algorithm based on weighted coherence

    3 理論模型分析

    為了驗(yàn)證基于加權(quán)一致性的蟻群算法進(jìn)行斷層檢測的可靠性,選取三維“qdome”模型(Claerbout and Fomel, 2012)中的主測線1.5 km位置的聯(lián)絡(luò)測線垂直剖面進(jìn)行測試(圖1a),該模型包含頂部的水平地層、中部的高斯型地層、底部的傾斜地層、不同類型地層區(qū)域邊界的水平不整合面及一個正斷層.圖1b是加入正態(tài)分布隨機(jī)噪聲以后的合成地震記錄.圖1a和圖1b的地震數(shù)據(jù)的時間深度為500個采樣點(diǎn),采樣間隔為0.01 s,測線長度為150道地震道,道間距為0.01 km.

    用基于加權(quán)一致性的蟻群算法分別對圖1a和圖1b進(jìn)行斷層檢測.參數(shù)設(shè)置為:分割地震數(shù)據(jù)的小數(shù)據(jù)塊大小為2×2;終止條件Smin=0.2、Smax=0.5;信息素濃度初始值為1,信息素更新系數(shù)c=1;信息素和啟發(fā)函數(shù)的重要程度因子分別是α=1、β=3;信息素?fù)]發(fā)系數(shù)ρ=0.01;迭代次數(shù)為5次.檢測結(jié)果如圖6a和圖6b所示,可見,無論是理論無噪聲地震數(shù)據(jù)(圖1a),還是加入強(qiáng)隨機(jī)噪聲的地震數(shù)據(jù)(圖1b),本文所提方法都能準(zhǔn)確、清晰且完整地檢測出斷層信息,抗噪性強(qiáng).

    選取基于方差屬性的蟻群算法(圖6c、圖6d)和基于C3相干屬性的蟻群算法(圖6e、圖6f)對圖1a和圖1b進(jìn)行斷層檢測.對比可見,對于理論無噪聲地震數(shù)據(jù)(圖1a),三種方法都能準(zhǔn)確、清晰、完整地檢測出斷層位置,但基于C3相干屬性的蟻群算法對大傾角地層很敏感,檢測結(jié)果伴有明顯的地層殘余響應(yīng).而對于加入強(qiáng)隨機(jī)噪聲的地震數(shù)據(jù)(圖1b),后兩種方法對噪聲壓制效果差,嚴(yán)重干擾對斷層的識別.圖6斷層檢測結(jié)果表明本文所提方法在理論模型中進(jìn)行斷層檢測是可行且有效的,而且本文所提方法的抗噪性最優(yōu).

    圖6 理論地震數(shù)據(jù)的斷層檢測結(jié)果(a) 基于加權(quán)一致性的蟻群算法對圖1a斷層檢測結(jié)果; (b) 基于加權(quán)一致性的蟻群算法對圖1b斷層檢測結(jié)果; (c) 基于方差屬性的蟻群算法對圖1a斷層檢測結(jié)果; (d) 基于方差屬性的蟻群算法對圖1b斷層檢測結(jié)果; (e) 基于C3相干屬性的蟻群算法對圖1a斷層檢測結(jié)果; (f) 基于C3相干屬性的蟻群算法對圖1b斷層檢測結(jié)果.Fig.6 Results of fault detection on synthetic seismic data(a) Fault detection for Fig.1a by ant colony algorithm based on weighted coherence; (b) Fault detection for Fig.1b by ant colony algorithm based on weighted coherence; (c) Fault detection for Fig.1a by ant colony algorithm based on variance attribute; (d) Fault detection for Fig.1b by ant colony algorithm based on variance attribute; (e) Fault detection for Fig.1a by ant colony algorithm based on C3 coherence attribute; (f) Fault detection for Fig.1b by ant colony algorithm based on C3 coherence attribute.

    4 實(shí)際數(shù)據(jù)分析

    選取墨西哥灣某地區(qū)的二維時間偏移剖面對本文所提方法進(jìn)行測試(Claerbout, 2010)(圖4a).該地震數(shù)據(jù)的時間深度范圍為0.4~1.8 s,有350個采樣點(diǎn),采樣間隔為0.004 s,測線長度為200道地震道.

    用基于加權(quán)一致性的蟻群算法對圖4a進(jìn)行斷層檢測.參數(shù)設(shè)置為:分割地震數(shù)據(jù)的小數(shù)據(jù)塊大小為2×2;終止條件Smin=0.2、Smax=0.5;各數(shù)據(jù)點(diǎn)的信息素濃度初始值為1,信息素更新系數(shù)c=0.5;信息素和啟發(fā)函數(shù)的重要程度因子分別是α=1、β=3;信息素?fù)]發(fā)系數(shù)ρ=0.1;迭代次數(shù)為20次.檢測結(jié)果如圖7a所示,結(jié)果能清晰且較完整地指示出斷層位置及延伸長度,并較好地壓制了噪聲和地層殘余響應(yīng)這兩類偽斷層信息的影響.為了進(jìn)一步比較,給出由基于方差屬性的蟻群算法和基于C3相干屬性的蟻群算法對圖4a計(jì)算而得到的斷層信息(圖7b、圖7c).因?yàn)榉讲钏惴ê虲3相干算法對噪聲和地層殘余響應(yīng)的壓制效果差,致使蟻群追蹤到的斷層信息較模糊且不完整,對于地層信息匱乏的斷層區(qū)域,這兩種方法并不能檢測出應(yīng)有的斷層信息.圖7的對比結(jié)果證明了本文所提方法在實(shí)際地震數(shù)據(jù)中進(jìn)行斷層檢測的可行性和優(yōu)越性.

    圖7 實(shí)際地震數(shù)據(jù)(圖4a)的斷層檢測結(jié)果(a) 基于加權(quán)一致性的蟻群算法對圖4a的斷層檢測結(jié)果; (b) 基于方差屬性的蟻群算法對圖4a的斷層檢測結(jié)果; (c) 基于C3相干屬性的蟻群算法對圖4a的斷層檢測結(jié)果.Fig.7 Results of fault detection on real seismic data (Fig.4a)(a) Fault detection for Fig.4a by ant colony algorithm based on weighted coherence; (b) Fault detection for Fig.4a by ant colony algorithm based on variance attribute; (c) Fault detection for Fig.4a by ant colony algorithm based on C3 coherence attribute.

    5 計(jì)算效率分析

    為了分析本文所提方法進(jìn)行斷層檢測的計(jì)算效率,本文列出基于加權(quán)一致性的蟻群算法、基于方差屬性的蟻群算法和基于C3相干屬性的蟻群算法對相同地震數(shù)據(jù)(圖1a,圖1b,圖4a)進(jìn)行斷層檢測的運(yùn)行時間表(表1).本文運(yùn)行程序所用電腦配置為:操作系統(tǒng):Windows 10 專業(yè)版64位;處理器:Inter(R) Core(TM) i3-3220 CPU@ 3.30GHz(4CPUs),~3.3 GHz;內(nèi)存:4096 MB RAM.三種方法在蟻群追蹤環(huán)節(jié)的參數(shù)設(shè)置相同,由表1可見,本文所提方法的計(jì)算效率最高.

    表1 不同方法對相同地震數(shù)據(jù)進(jìn)行斷層檢測的運(yùn)行時間(單位:s)Table 1 Operation times consumed by different methods for fault detection on same seismic data (unit: s)

    6 結(jié)論

    本文提出了新的斷層自動檢測方法——基于加權(quán)一致性的蟻群算法,得出了以下結(jié)論:

    (1) 一致性是對指紋圖像中指示局部方向場信息的各向異性強(qiáng)度的一個度量,將其應(yīng)用于地震數(shù)據(jù)處理,可指示地層的不連續(xù)性.

    (2) 提出加權(quán)一致性的概念,與一致性相比,前者具有更高的抗噪性能,更大程度地消除了不整合面和地層殘余響應(yīng)這兩類偽斷層信息,更能突顯地層的橫向不連續(xù)性.

    (3) 對于相同的斷層檢測過程,對比常規(guī)的基于方差屬性的蟻群算法及基于C3相干屬性的蟻群算法,本文所提方法具有以下優(yōu)點(diǎn):抗噪性更高,對偽斷層信息壓制效果更好,斷層邊界收斂能力更強(qiáng);在增加斷層連續(xù)性的基礎(chǔ)之上,對斷層的延伸長度展現(xiàn)的更完整;穩(wěn)定性和計(jì)算效率更高.

    Bahorich M, Farmer S. 1995. 3-D seismic discontinuity for faults and stratigraphic features: The coherence cube.TheLeadingEdge, 14(10): 1053-1058. Claerbout J F. 2010.Basic Earth imaging.Stanford Exploration Project. http:∥sepwww.stanford.edu/sep/prof/bei11.2010.pdf.Claerbout J F,Fomel S.2012.Image estimation by example:Geophysical soundings image construction.Stanford Exploration Project. http://sepwww.stanford.edu/sep/prof/gee1-2012.pdf.

    Dorn G A,James H E,Evins L. Automatic fault extraction (AFE) in 3D seismic data.2005 CSEG National Convention,247-250.

    Duan H B. 2005. Ant Colony Algorithms: Theory and Applications (in Chinese). Beijing: Science Press. Gersztenkorn A, Marfurt K J. 1999. Eigenstructure-based coherence computations as an aid to 3-D structural and stratigraphic mapping.Geophysics, 64(5): 1468-1479. Kass M, Witkin A. 1987. Analyzing oriented patterns.Computer

    Vision,Graphics,andImageProcessing, 37(3): 362-385.

    Li X C. 2010. Ant colony optimization and application of 3-D seismic information for refined explanation of structure in mining section [Ph. D. thesis] (in Chinese). Xuzhou: China University of Ming and Technology. Liu W. 2013. Research on the ant colony optimization algorithm and its applications [Master′s thesis] (in Chinese). Beijing: Beijing University of Posts and Telecommunications.

    Liu Y, Wang D, Liu C, et al. 2014. Structure-oriented filtering and fault detection based on nonstationary similarity.ChineseJ.Geophys. (in Chinese), 57(4): 1177-1187, doi: 10.6038/cjg20140415.

    Marfurt K J, Kirlin R L, Farmer S L, et al. 1998. 3-D seismic attributes using a semblance-based coherency algorithm.Geophysics, 63(4): 1150-1165.

    Pedersen S I, Randen T, S?nneland L. 2002. Automatic fault extraction using artificial ants.SEGTechnicalProgramExpandedAbstracts, 21(1): 512-515.Pedersen S I, Skov T, Hetlelid A, et al. 2003. New paradigm of fault interpretation.SEGTechnicalProgramExpandedAbstracts, 22(1): 350-353.

    Randen T, Pedersen S I, S?nneland L. 2001. Automatic detection and extraction of faults from three-dimensional seismic data.SEGTechnicalProgramExpandedAbstracts, 20(1): 551-555. Wang Y, Hu J K, Schroder H. 2005. A gradient-based weighted averaging method for estimation of fingerprint orientation fields. ∥ Proceedings of the Digital Image Computing: Techniques and Applications. Piscataway, USA: IEEE.Yan Z. 2010. Fault automatic recognition and intelligent interpretation in 3-D seismic data [Ph. D. thesis] (in Chinese). Wuhan: China University of Geosciences. Yan Z, Gu H M, Cai C G. 2013. Automatic fault tracking based on ant colony algorithms.Computers&Geosciences, 51: 269-281.

    Zhao W. 2009. Research on 3D seismic fault identification method based on ant colony algorithm [Master′s thesis] (in Chinese). Nanjing: Nanjing University of Science & Technology.

    附中文參考文獻(xiàn)

    段海濱. 2005. 蟻群算法原理及其應(yīng)用. 北京: 科學(xué)出版社.

    李香臣. 2010. 采區(qū)構(gòu)造三維地震精細(xì)解釋的蟻群算法及其應(yīng)用[博士論文]. 徐州: 中國礦業(yè)大學(xué).

    劉聞. 2013. 蟻群算法及其應(yīng)用研究[碩士論文]. 北京: 北京郵電大學(xué).

    劉洋, 王典, 劉財(cái)?shù)? 2014. 基于非平穩(wěn)相似性系數(shù)的構(gòu)造導(dǎo)向?yàn)V波及斷層檢測方法. 地球物理學(xué)報(bào), 57(4): 1177-1187, doi: 10.6038/cjg20140415.

    嚴(yán)哲. 2010. 三維地震斷層自動識別與智能解釋[博士論文]. 武漢: 中國地質(zhì)大學(xué).

    趙偉. 2009. 基于蟻群算法的三維地震斷層識別方法研究[碩士論文]. 南京: 南京理工大學(xué).

    (本文編輯何燕)

    Application of the ant colony algorithm based on weighted coherence to fault detection

    LIU Cai1, LIU Hai-Yan1*, PENG Chong2, ZHANG Ying3, LIU Qiong4, GOU Fu-Yan1

    1CollegeofGeo-explorationScienceandTechnology,JilinUniversity,Changchun130026,China2CollegeofGeographicalSciences,ShanxiNormalUniversity,ShanxiLinfen041004,China>3NortheastOil&GasBranch,SINOPEC,Changchun130062,China4ResearchInstituteofExplorationandDevelopment,TuhaOilfieldCompany,PetroChina,XinjiangHami839009,China

    It is the basis of automatic interpretation of faults to highlight the amplitude spatial discontinuity of seismic data,while how to suppress the noise, the formation residual response, the angular unconformable contact, and other pseudo faults is one of the key issues in this approach.Coherence is a measure of the strength of anisotropy that indicates local direction field information in the detection process of fingerprint image edge,which is applied to seismic data processing as a new fault enhancement attribute in this paper.In order to suppress noise more effectively,the coherence is weighted firstly. Then we propose a new fault automatic detection method called the colony algorithm based on weighted coherence,which is formed by the weighted coherence and ant colony algorithm.Numerical tests on synthetic and real field seismic data show that the proposed method in this paper is superior to the conventional ant colony algorithm based on variance attribute and the other based on C3 coherence attribute for the same fault detection process. It can not only more effectively suppress the influence of pseudo fault information generated by the noise interference and formation residual response,but also indicate the extension length of the fault more completely. Meanwhile, this new algorithm has higher stability and computational efficiency than the existing ones aforementioned.

    Coherence; Weighted coherence; Ant colony algorithm; Fault detection

    10.6038/cjg20161029.

    國家自然科學(xué)基金重點(diǎn)基金(41430322),國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)課題(2013CB429805)聯(lián)合資助.

    劉財(cái),男,1963年生,教授,博士生導(dǎo)師,主要從事地震波場正反演理論、綜合地球物理等研究.E-mail:liucai@jlu.edu.cn

    劉海燕,女,1986年生,博士,主要從事地震數(shù)據(jù)處理與解釋等研究.E-mail:lhy509516@163.com

    10.6038/cjg20161029

    P631

    2016-02-01,2016-05-26收修定稿

    劉財(cái), 劉海燕, 彭沖等. 2016. 基于加權(quán)一致性的蟻群算法在斷層檢測中的應(yīng)用. 地球物理學(xué)報(bào),59(10):3859-3868,

    Liu C, Liu H Y, Peng C, et al. 2016. Application of the ant colony algorithm based on weighted coherence to fault detection.ChineseJ.Geophys. (in Chinese),59(10):3859-3868,doi:10.6038/cjg20161029.

    猜你喜歡
    蟻群梯度一致性
    關(guān)注減污降碳協(xié)同的一致性和整體性
    公民與法治(2022年5期)2022-07-29 00:47:28
    一個改進(jìn)的WYL型三項(xiàng)共軛梯度法
    注重教、學(xué)、評一致性 提高一輪復(fù)習(xí)效率
    IOl-master 700和Pentacam測量Kappa角一致性分析
    游戲社會:狼、猞猁和蟻群
    一種自適應(yīng)Dai-Liao共軛梯度法
    一類扭積形式的梯度近Ricci孤立子
    基于自適應(yīng)蟻群的FCM聚類優(yōu)化算法研究
    基于奇異值差分譜分析和蟻群算法的小波閾值降噪
    基于事件觸發(fā)的多智能體輸入飽和一致性控制
    淫妇啪啪啪对白视频| 久久亚洲真实| 男女视频在线观看网站免费| 久久久色成人| 老司机午夜福利在线观看视频| 狂野欧美激情性xxxx| 日本黄色片子视频| 99热这里只有精品一区 | 久久久国产成人免费| 午夜精品久久久久久毛片777| 久久久久性生活片| 成人18禁在线播放| 国产精品一区二区免费欧美| 国产精品乱码一区二三区的特点| 亚洲无线观看免费| 在线a可以看的网站| 国产视频一区二区在线看| 色哟哟哟哟哟哟| 婷婷精品国产亚洲av在线| 丁香欧美五月| 亚洲国产欧洲综合997久久,| 国产三级黄色录像| 亚洲精品中文字幕一二三四区| 亚洲精品久久国产高清桃花| 桃色一区二区三区在线观看| 亚洲在线观看片| 在线国产一区二区在线| 日韩精品中文字幕看吧| 日本黄色视频三级网站网址| 亚洲国产精品999在线| 久久中文看片网| 噜噜噜噜噜久久久久久91| 啪啪无遮挡十八禁网站| 日韩有码中文字幕| 欧美极品一区二区三区四区| 中文字幕高清在线视频| 少妇熟女aⅴ在线视频| 日本成人三级电影网站| 在线视频色国产色| 天堂动漫精品| 欧美在线一区亚洲| 国产成人影院久久av| 色视频www国产| 人妻久久中文字幕网| 91老司机精品| 久久久久精品国产欧美久久久| 不卡av一区二区三区| 999精品在线视频| 国产一区二区在线观看日韩 | 国产99白浆流出| 手机成人av网站| 久9热在线精品视频| 99国产精品一区二区蜜桃av| a级毛片在线看网站| 国产精品久久久久久精品电影| 男女午夜视频在线观看| 黄色女人牲交| 国产精品香港三级国产av潘金莲| 美女免费视频网站| 午夜免费成人在线视频| 色在线成人网| av国产免费在线观看| 91麻豆精品激情在线观看国产| 老汉色av国产亚洲站长工具| 国产又色又爽无遮挡免费看| 香蕉av资源在线| 日本精品一区二区三区蜜桃| 国产黄a三级三级三级人| 亚洲av电影在线进入| 丰满的人妻完整版| 久久精品国产亚洲av香蕉五月| 草草在线视频免费看| 色综合亚洲欧美另类图片| 一二三四社区在线视频社区8| 美女黄网站色视频| 国产午夜精品久久久久久| 国产精品一区二区三区四区久久| 日本 av在线| www.自偷自拍.com| 日韩免费av在线播放| 男女视频在线观看网站免费| 久久香蕉国产精品| 国产成年人精品一区二区| 亚洲激情在线av| 丰满人妻一区二区三区视频av | 久久久国产精品麻豆| 亚洲av熟女| 久久久久久久精品吃奶| 免费在线观看成人毛片| 黑人操中国人逼视频| 人人妻人人澡欧美一区二区| 免费搜索国产男女视频| 亚洲第一电影网av| 一区二区三区激情视频| 亚洲精品国产精品久久久不卡| 亚洲人成电影免费在线| 99久久99久久久精品蜜桃| 久久久国产欧美日韩av| 国产成人精品无人区| 精品久久久久久成人av| 欧美日韩综合久久久久久 | 夜夜夜夜夜久久久久| 中文字幕最新亚洲高清| 又黄又粗又硬又大视频| 久久久久久人人人人人| 国产高清激情床上av| 久久中文字幕人妻熟女| 成人av在线播放网站| 欧美国产日韩亚洲一区| 久久久久国产精品人妻aⅴ院| 97碰自拍视频| 亚洲中文av在线| 性色avwww在线观看| 老司机福利观看| aaaaa片日本免费| 亚洲在线自拍视频| 久久香蕉精品热| 一区福利在线观看| 观看美女的网站| 可以在线观看的亚洲视频| av黄色大香蕉| 色老头精品视频在线观看| 亚洲av五月六月丁香网| 村上凉子中文字幕在线| 天堂√8在线中文| 欧美黑人巨大hd| 亚洲一区高清亚洲精品| 非洲黑人性xxxx精品又粗又长| 欧美中文综合在线视频| 两个人看的免费小视频| 成人国产综合亚洲| 亚洲成人精品中文字幕电影| 黄色视频,在线免费观看| 成人av在线播放网站| 久久热在线av| 白带黄色成豆腐渣| 国产精品久久久av美女十八| 欧美日韩乱码在线| 日本撒尿小便嘘嘘汇集6| 国产又黄又爽又无遮挡在线| 亚洲电影在线观看av| 午夜福利18| 亚洲中文av在线| 亚洲五月婷婷丁香| 国产美女午夜福利| 成人性生交大片免费视频hd| 久久草成人影院| 美女被艹到高潮喷水动态| 999久久久精品免费观看国产| 亚洲美女视频黄频| 国产成人精品无人区| 少妇的丰满在线观看| 亚洲欧美精品综合久久99| 亚洲av成人一区二区三| 日韩欧美国产一区二区入口| 国产成人福利小说| 麻豆国产av国片精品| 中文字幕最新亚洲高清| 黄色成人免费大全| 国产蜜桃级精品一区二区三区| 日韩欧美免费精品| 久久亚洲真实| 日本在线视频免费播放| 成年版毛片免费区| 国产不卡一卡二| 免费大片18禁| 亚洲成a人片在线一区二区| 可以在线观看毛片的网站| 国产精品久久久久久久电影 | 一本综合久久免费| 亚洲真实伦在线观看| 女生性感内裤真人,穿戴方法视频| 给我免费播放毛片高清在线观看| 脱女人内裤的视频| 身体一侧抽搐| 久久久久久国产a免费观看| 精品欧美国产一区二区三| 久久午夜综合久久蜜桃| 国产精品一区二区精品视频观看| 亚洲国产精品999在线| 精品国产乱码久久久久久男人| 精品一区二区三区四区五区乱码| 亚洲 欧美一区二区三区| 999精品在线视频| 免费电影在线观看免费观看| 欧美成狂野欧美在线观看| 国产激情久久老熟女| 亚洲欧美日韩东京热| 国产精品影院久久| 午夜福利欧美成人| 最近最新免费中文字幕在线| 99久久成人亚洲精品观看| 老司机在亚洲福利影院| 亚洲美女视频黄频| 亚洲国产欧美人成| 一个人观看的视频www高清免费观看 | 中文字幕av在线有码专区| 亚洲成人免费电影在线观看| 美女高潮喷水抽搐中文字幕| 欧美乱妇无乱码| 高潮久久久久久久久久久不卡| 精品国产亚洲在线| 久久精品国产清高在天天线| 一级作爱视频免费观看| 精品久久蜜臀av无| 国产伦一二天堂av在线观看| 最好的美女福利视频网| 日本a在线网址| 成在线人永久免费视频| 亚洲一区高清亚洲精品| 免费在线观看视频国产中文字幕亚洲| 国产一级毛片七仙女欲春2| 日韩欧美精品v在线| 国产精品女同一区二区软件 | 中文字幕人妻丝袜一区二区| 午夜精品久久久久久毛片777| 国产精品亚洲av一区麻豆| 熟妇人妻久久中文字幕3abv| 无遮挡黄片免费观看| 亚洲九九香蕉| e午夜精品久久久久久久| 成人永久免费在线观看视频| 好男人电影高清在线观看| 五月玫瑰六月丁香| 三级毛片av免费| 欧美日本视频| 人人妻人人澡欧美一区二区| 男女之事视频高清在线观看| 亚洲第一欧美日韩一区二区三区| 国产精品乱码一区二三区的特点| 在线观看一区二区三区| 亚洲成av人片免费观看| 女人被狂操c到高潮| 无限看片的www在线观看| 成人国产综合亚洲| 88av欧美| 在线观看免费午夜福利视频| 欧美高清成人免费视频www| 午夜影院日韩av| 日日干狠狠操夜夜爽| 法律面前人人平等表现在哪些方面| 99国产精品一区二区三区| 国产精品香港三级国产av潘金莲| 性色avwww在线观看| 久久久精品大字幕| 国产高潮美女av| 国产成人福利小说| 久久久久性生活片| 日韩精品中文字幕看吧| 此物有八面人人有两片| 18禁美女被吸乳视频| 久久伊人香网站| 中文字幕久久专区| 亚洲精品一区av在线观看| 国产精品久久久人人做人人爽| 这个男人来自地球电影免费观看| 天堂网av新在线| 日韩 欧美 亚洲 中文字幕| 一级作爱视频免费观看| 日韩中文字幕欧美一区二区| 制服丝袜大香蕉在线| 五月玫瑰六月丁香| 国产高清videossex| 日韩人妻高清精品专区| 桃红色精品国产亚洲av| 岛国在线观看网站| 天天添夜夜摸| 日日摸夜夜添夜夜添小说| 亚洲美女视频黄频| 精品久久蜜臀av无| 国产男靠女视频免费网站| 久久九九热精品免费| 成人高潮视频无遮挡免费网站| 国产午夜精品久久久久久| 中文字幕最新亚洲高清| 男女下面进入的视频免费午夜| 午夜激情欧美在线| 亚洲电影在线观看av| 亚洲黑人精品在线| 特级一级黄色大片| 淫秽高清视频在线观看| 精品国产乱子伦一区二区三区| 一a级毛片在线观看| 热99re8久久精品国产| 精品国产超薄肉色丝袜足j| 岛国在线免费视频观看| 欧美高清成人免费视频www| 一二三四社区在线视频社区8| 欧美国产日韩亚洲一区| 国产精品香港三级国产av潘金莲| 中出人妻视频一区二区| 亚洲真实伦在线观看| 午夜精品久久久久久毛片777| 女同久久另类99精品国产91| 母亲3免费完整高清在线观看| 在线观看舔阴道视频| 日韩三级视频一区二区三区| 亚洲在线自拍视频| 美女黄网站色视频| 天堂影院成人在线观看| 在线观看一区二区三区| 特大巨黑吊av在线直播| 亚洲美女黄片视频| 久久这里只有精品19| 国产极品精品免费视频能看的| 精品久久久久久久末码| 国语自产精品视频在线第100页| 听说在线观看完整版免费高清| 国产av麻豆久久久久久久| 亚洲激情在线av| 国产成人啪精品午夜网站| 在线视频色国产色| 成人鲁丝片一二三区免费| 97碰自拍视频| x7x7x7水蜜桃| 久久精品人妻少妇| 国产淫片久久久久久久久 | 成熟少妇高潮喷水视频| 无遮挡黄片免费观看| 99热6这里只有精品| 国产高潮美女av| 禁无遮挡网站| 欧美日韩黄片免| 色精品久久人妻99蜜桃| av福利片在线观看| 久久天躁狠狠躁夜夜2o2o| 又黄又爽又免费观看的视频| 91久久精品国产一区二区成人 | 亚洲熟女毛片儿| 99久久无色码亚洲精品果冻| 日韩精品中文字幕看吧| 亚洲国产精品成人综合色| 日本熟妇午夜| 亚洲精品中文字幕一二三四区| 亚洲无线观看免费| 听说在线观看完整版免费高清| 少妇的逼水好多| 国产精品九九99| 国产成人av激情在线播放| 18禁国产床啪视频网站| 亚洲欧洲精品一区二区精品久久久| 一a级毛片在线观看| 狂野欧美激情性xxxx| 精品人妻1区二区| 国产成人aa在线观看| 黄色 视频免费看| 香蕉国产在线看| 精品国产超薄肉色丝袜足j| 午夜免费激情av| 亚洲精品一卡2卡三卡4卡5卡| 性色av乱码一区二区三区2| 国产亚洲精品一区二区www| 九九热线精品视视频播放| 99久久精品一区二区三区| 十八禁网站免费在线| 亚洲在线自拍视频| 九色成人免费人妻av| 91老司机精品| 久久香蕉国产精品| 欧美高清成人免费视频www| 一个人免费在线观看的高清视频| 好看av亚洲va欧美ⅴa在| 女生性感内裤真人,穿戴方法视频| 日韩欧美在线二视频| 日韩精品青青久久久久久| 一级作爱视频免费观看| 老熟妇仑乱视频hdxx| 欧美日韩一级在线毛片| 日韩免费av在线播放| 操出白浆在线播放| 久久久久亚洲av毛片大全| 女警被强在线播放| 国产三级中文精品| 精品熟女少妇八av免费久了| 久久欧美精品欧美久久欧美| 美女高潮的动态| 男人舔奶头视频| 狂野欧美白嫩少妇大欣赏| 日韩欧美 国产精品| 在线观看免费午夜福利视频| 午夜免费成人在线视频| 老司机午夜福利在线观看视频| 巨乳人妻的诱惑在线观看| 性欧美人与动物交配| 人人妻人人澡欧美一区二区| 婷婷丁香在线五月| av欧美777| 美女扒开内裤让男人捅视频| 色综合亚洲欧美另类图片| 级片在线观看| 天堂√8在线中文| 国产精品野战在线观看| av在线天堂中文字幕| 亚洲国产日韩欧美精品在线观看 | 在线观看免费午夜福利视频| 精品人妻1区二区| 又爽又黄无遮挡网站| 精品久久久久久,| 久久草成人影院| 国产又色又爽无遮挡免费看| 欧美大码av| 精品乱码久久久久久99久播| 日日摸夜夜添夜夜添小说| 国产精品爽爽va在线观看网站| 久久久久久久久中文| 色吧在线观看| 国产探花在线观看一区二区| 久久久精品欧美日韩精品| 欧美乱色亚洲激情| 嫩草影院入口| 国内精品美女久久久久久| 欧美日韩一级在线毛片| 男女午夜视频在线观看| 国产精品久久电影中文字幕| 1024香蕉在线观看| 一级毛片精品| 国产伦人伦偷精品视频| 亚洲专区国产一区二区| 后天国语完整版免费观看| 日韩人妻高清精品专区| 国产单亲对白刺激| 亚洲片人在线观看| 91老司机精品| 亚洲熟妇中文字幕五十中出| 成人欧美大片| 国产精品自产拍在线观看55亚洲| 不卡av一区二区三区| 欧美zozozo另类| 国语自产精品视频在线第100页| 十八禁人妻一区二区| 三级国产精品欧美在线观看 | 中文字幕人成人乱码亚洲影| 在线播放国产精品三级| 变态另类丝袜制服| 欧美乱码精品一区二区三区| 最近最新中文字幕大全免费视频| 午夜两性在线视频| 国产真实乱freesex| 国产精品久久久久久亚洲av鲁大| www国产在线视频色| 久久这里只有精品中国| 亚洲人成网站高清观看| 国产精品日韩av在线免费观看| 一本综合久久免费| 久久精品91无色码中文字幕| 国产高潮美女av| av片东京热男人的天堂| 国产不卡一卡二| 国产又黄又爽又无遮挡在线| 成人特级av手机在线观看| 国产熟女xx| 国产精品久久久久久人妻精品电影| 亚洲欧美精品综合一区二区三区| 精品一区二区三区视频在线 | 精品国产乱码久久久久久男人| 老熟妇乱子伦视频在线观看| 一本一本综合久久| 男女午夜视频在线观看| 成人18禁在线播放| 麻豆国产97在线/欧美| 一本综合久久免费| 一边摸一边抽搐一进一小说| 男人舔女人的私密视频| 成人三级做爰电影| 深夜精品福利| 日韩欧美在线乱码| 亚洲国产欧美一区二区综合| 成人av一区二区三区在线看| 亚洲精品456在线播放app | 国产亚洲欧美在线一区二区| av黄色大香蕉| 色综合站精品国产| 色老头精品视频在线观看| 嫁个100分男人电影在线观看| 男人舔女人的私密视频| 人妻久久中文字幕网| 99热只有精品国产| 亚洲中文日韩欧美视频| 色播亚洲综合网| 国产一级毛片七仙女欲春2| 很黄的视频免费| 91av网一区二区| 99久久成人亚洲精品观看| 国产成人福利小说| 757午夜福利合集在线观看| 桃红色精品国产亚洲av| av天堂在线播放| 人妻久久中文字幕网| 岛国视频午夜一区免费看| 99久久99久久久精品蜜桃| 99久久无色码亚洲精品果冻| 国产高清视频在线播放一区| 少妇的丰满在线观看| 1024香蕉在线观看| 黑人操中国人逼视频| 嫩草影院入口| 男人舔女人的私密视频| www.精华液| 国产成人福利小说| 久久久久国产精品人妻aⅴ院| 毛片女人毛片| 免费看a级黄色片| 法律面前人人平等表现在哪些方面| 亚洲精品美女久久av网站| 国产免费男女视频| 亚洲成av人片免费观看| 亚洲欧美日韩无卡精品| 亚洲人成网站高清观看| 我要搜黄色片| 亚洲精品在线美女| 又爽又黄无遮挡网站| 免费电影在线观看免费观看| 韩国av一区二区三区四区| 国产成人精品无人区| 我的老师免费观看完整版| 狂野欧美激情性xxxx| 97超视频在线观看视频| 久久精品亚洲精品国产色婷小说| 禁无遮挡网站| www日本黄色视频网| 亚洲在线观看片| 日本免费a在线| 精品国产超薄肉色丝袜足j| 美女高潮的动态| 黄片小视频在线播放| 日韩欧美在线乱码| 在线观看午夜福利视频| 亚洲精品在线观看二区| 亚洲天堂国产精品一区在线| 噜噜噜噜噜久久久久久91| av福利片在线观看| 免费在线观看成人毛片| 18禁国产床啪视频网站| 国产成+人综合+亚洲专区| 这个男人来自地球电影免费观看| 国产男靠女视频免费网站| 最新美女视频免费是黄的| 精品久久久久久久毛片微露脸| 精品福利观看| 国产美女午夜福利| 久久99热这里只有精品18| 亚洲国产欧美网| 日韩高清综合在线| 一区二区三区高清视频在线| 国产精品久久电影中文字幕| 亚洲无线在线观看| 黑人巨大精品欧美一区二区mp4| 日本撒尿小便嘘嘘汇集6| 国产私拍福利视频在线观看| 成人永久免费在线观看视频| tocl精华| 又粗又爽又猛毛片免费看| www.熟女人妻精品国产| 亚洲av免费在线观看| 欧美高清成人免费视频www| www.www免费av| 久久久久亚洲av毛片大全| 久久午夜综合久久蜜桃| 国产黄色小视频在线观看| 国产亚洲精品久久久久久毛片| 国产日本99.免费观看| 午夜a级毛片| 精品一区二区三区av网在线观看| 国产伦人伦偷精品视频| 一进一出抽搐gif免费好疼| 法律面前人人平等表现在哪些方面| 国产一区二区在线av高清观看| 黄色日韩在线| 亚洲成a人片在线一区二区| 亚洲精品色激情综合| www日本在线高清视频| 亚洲成人久久爱视频| 色综合亚洲欧美另类图片| 国产亚洲精品综合一区在线观看| 禁无遮挡网站| 91av网一区二区| 看免费av毛片| 久久久久久久精品吃奶| 999久久久国产精品视频| 高清在线国产一区| 国产精品自产拍在线观看55亚洲| 亚洲国产精品合色在线| 可以在线观看的亚洲视频| 男女下面进入的视频免费午夜| 国产成人精品久久二区二区免费| 桃红色精品国产亚洲av| 美女午夜性视频免费| 国产欧美日韩精品一区二区| 成人无遮挡网站| 国产三级在线视频| 91老司机精品| 免费搜索国产男女视频| 久久欧美精品欧美久久欧美| 日本一本二区三区精品| 亚洲国产精品成人综合色| 午夜免费成人在线视频| 午夜精品一区二区三区免费看| av视频在线观看入口| 日韩免费av在线播放| 欧美大码av| 成人午夜高清在线视频| 男人舔女人的私密视频| 九色国产91popny在线| 亚洲aⅴ乱码一区二区在线播放| 男人舔女人的私密视频| 九色国产91popny在线| 69av精品久久久久久| 国产av在哪里看| 一级毛片女人18水好多| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品999在线| 麻豆一二三区av精品| 丝袜人妻中文字幕| 美女午夜性视频免费| 91在线观看av|