• <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ā)的多智能體輸入飽和一致性控制
    18禁动态无遮挡网站| 久久久久久久大尺度免费视频| 一级片免费观看大全| 一边摸一边做爽爽视频免费| 国产成人午夜福利电影在线观看| 国产精品久久久久久久久免| 亚洲第一青青草原| 美女高潮到喷水免费观看| av网站免费在线观看视频| 久久精品久久精品一区二区三区| 国产精品一二三区在线看| 国产精品久久久久久人妻精品电影 | 亚洲成av片中文字幕在线观看| 大片免费播放器 马上看| 国产精品国产三级专区第一集| 日本欧美国产在线视频| 国产亚洲一区二区精品| 99re6热这里在线精品视频| av片东京热男人的天堂| a级片在线免费高清观看视频| 最近手机中文字幕大全| 日韩av免费高清视频| 久久久精品区二区三区| 青春草国产在线视频| av国产久精品久网站免费入址| 一区二区日韩欧美中文字幕| 99久久99久久久精品蜜桃| 欧美精品人与动牲交sv欧美| 欧美av亚洲av综合av国产av | 妹子高潮喷水视频| 男的添女的下面高潮视频| 十八禁高潮呻吟视频| 国产精品无大码| 精品一区二区三卡| 9色porny在线观看| 在线看a的网站| 免费黄色在线免费观看| 18禁裸乳无遮挡动漫免费视频| 国产男女超爽视频在线观看| 嫩草影视91久久| 欧美精品亚洲一区二区| 亚洲国产av新网站| 宅男免费午夜| 国产国语露脸激情在线看| 丰满饥渴人妻一区二区三| 香蕉丝袜av| 欧美av亚洲av综合av国产av | 五月开心婷婷网| 赤兔流量卡办理| 亚洲国产最新在线播放| 亚洲成色77777| 超碰97精品在线观看| 亚洲欧美一区二区三区久久| 国产精品三级大全| 制服丝袜香蕉在线| 久久婷婷青草| 乱人伦中国视频| 亚洲五月色婷婷综合| 国产片特级美女逼逼视频| 综合色丁香网| 啦啦啦在线观看免费高清www| 国产黄频视频在线观看| 啦啦啦啦在线视频资源| 亚洲,一卡二卡三卡| 日韩精品有码人妻一区| 国产淫语在线视频| 午夜福利影视在线免费观看| 亚洲精品,欧美精品| 91aial.com中文字幕在线观看| 夫妻性生交免费视频一级片| 欧美 日韩 精品 国产| 91精品伊人久久大香线蕉| 婷婷色综合大香蕉| 在线观看免费日韩欧美大片| 国产成人精品无人区| 久久人妻熟女aⅴ| a级片在线免费高清观看视频| 亚洲欧美成人精品一区二区| 精品国产超薄肉色丝袜足j| 伦理电影免费视频| 国产一区二区在线观看av| 亚洲国产日韩一区二区| 男的添女的下面高潮视频| 激情视频va一区二区三区| 日本av手机在线免费观看| 又粗又硬又长又爽又黄的视频| 日本wwww免费看| 日本av免费视频播放| 日韩人妻精品一区2区三区| 黑丝袜美女国产一区| 亚洲国产日韩一区二区| 人妻 亚洲 视频| 韩国高清视频一区二区三区| 久久影院123| 日韩制服骚丝袜av| 成人毛片60女人毛片免费| 久久精品国产亚洲av涩爱| 少妇被粗大的猛进出69影院| 国产亚洲av片在线观看秒播厂| www日本在线高清视频| 久久青草综合色| 中文字幕色久视频| 一区在线观看完整版| 国产精品国产三级专区第一集| 精品福利永久在线观看| 国产免费现黄频在线看| 亚洲av日韩精品久久久久久密 | 精品一区二区三卡| 嫩草影院入口| 成人影院久久| 日日撸夜夜添| 成年女人毛片免费观看观看9 | 精品久久久久久电影网| 老司机影院毛片| 久久精品亚洲熟妇少妇任你| 国产激情久久老熟女| 99精国产麻豆久久婷婷| 久久精品aⅴ一区二区三区四区| 国产一级毛片在线| 国产爽快片一区二区三区| 久久精品亚洲熟妇少妇任你| 午夜日韩欧美国产| 夜夜骑夜夜射夜夜干| 天堂8中文在线网| 亚洲免费av在线视频| 亚洲免费av在线视频| 精品少妇一区二区三区视频日本电影 | 夫妻性生交免费视频一级片| 赤兔流量卡办理| 亚洲精品在线美女| 毛片一级片免费看久久久久| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产一区二区精华液| 国产精品无大码| 免费在线观看黄色视频的| 国产精品国产三级专区第一集| 国产 一区精品| 亚洲av国产av综合av卡| 欧美人与善性xxx| 人体艺术视频欧美日本| 亚洲第一区二区三区不卡| 国产欧美日韩一区二区三区在线| 午夜老司机福利片| av线在线观看网站| 久久久久国产精品人妻一区二区| 无限看片的www在线观看| xxxhd国产人妻xxx| 亚洲av电影在线观看一区二区三区| 精品一区在线观看国产| 久久久久精品国产欧美久久久 | 亚洲精品久久久久久婷婷小说| 天堂中文最新版在线下载| 在线看a的网站| 一级毛片黄色毛片免费观看视频| 国产在线视频一区二区| 少妇被粗大的猛进出69影院| 久久国产精品男人的天堂亚洲| 搡老岳熟女国产| 三上悠亚av全集在线观看| 精品国产超薄肉色丝袜足j| 激情视频va一区二区三区| 另类精品久久| 国产亚洲av高清不卡| 最近的中文字幕免费完整| 亚洲av日韩在线播放| 国产成人精品福利久久| 亚洲国产日韩一区二区| 国产毛片在线视频| 人妻一区二区av| 十八禁网站网址无遮挡| 99re6热这里在线精品视频| 纯流量卡能插随身wifi吗| 免费久久久久久久精品成人欧美视频| 中国三级夫妇交换| 国产成人91sexporn| 七月丁香在线播放| 免费av中文字幕在线| 大片电影免费在线观看免费| 丁香六月欧美| 另类精品久久| 国产色婷婷99| 另类亚洲欧美激情| av视频免费观看在线观看| 午夜福利乱码中文字幕| 亚洲欧美激情在线| 久久久久久免费高清国产稀缺| 国产亚洲av片在线观看秒播厂| 深夜精品福利| 欧美最新免费一区二区三区| 成年美女黄网站色视频大全免费| 亚洲国产欧美一区二区综合| 亚洲国产精品一区二区三区在线| 一级片免费观看大全| 久久精品亚洲av国产电影网| 少妇的丰满在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 成人国产麻豆网| 精品久久蜜臀av无| svipshipincom国产片| 我的亚洲天堂| 我的亚洲天堂| 大片免费播放器 马上看| 高清欧美精品videossex| 黄色视频在线播放观看不卡| 亚洲熟女毛片儿| 老司机在亚洲福利影院| 性色av一级| 嫩草影视91久久| 国产xxxxx性猛交| 久久天堂一区二区三区四区| 精品少妇黑人巨大在线播放| 国产在线视频一区二区| 国产精品二区激情视频| 夫妻性生交免费视频一级片| 亚洲欧美精品自产自拍| 成人毛片60女人毛片免费| 欧美亚洲 丝袜 人妻 在线| 欧美日韩福利视频一区二区| 亚洲国产av新网站| 交换朋友夫妻互换小说| 欧美精品一区二区大全| 久久精品国产亚洲av高清一级| 久久99精品国语久久久| 久久国产精品大桥未久av| 久久毛片免费看一区二区三区| av网站在线播放免费| 欧美另类一区| 狂野欧美激情性xxxx| 国产成人91sexporn| 欧美精品一区二区免费开放| 永久免费av网站大全| 搡老乐熟女国产| 欧美日韩一级在线毛片| 十八禁网站网址无遮挡| 十八禁高潮呻吟视频| 韩国av在线不卡| 大码成人一级视频| 老司机深夜福利视频在线观看 | 亚洲国产精品999| 亚洲av综合色区一区| 久久精品国产综合久久久| 美女扒开内裤让男人捅视频| 一级毛片黄色毛片免费观看视频| 国产在线免费精品| 日韩中文字幕欧美一区二区 | 国产成人a∨麻豆精品| 80岁老熟妇乱子伦牲交| 2021少妇久久久久久久久久久| 在线观看免费视频网站a站| av一本久久久久| 18禁观看日本| 两个人免费观看高清视频| 色吧在线观看| 青春草国产在线视频| 久久久国产精品麻豆| 久久久久精品性色| 飞空精品影院首页| 久久精品国产a三级三级三级| 亚洲第一av免费看| 国产伦理片在线播放av一区| 国产成人91sexporn| 亚洲精品国产色婷婷电影| www.自偷自拍.com| 国产成人精品久久久久久| 亚洲美女视频黄频| 一区二区三区精品91| 国精品久久久久久国模美| 肉色欧美久久久久久久蜜桃| 国产精品一国产av| 亚洲伊人久久精品综合| 国产片内射在线| 国产精品一区二区在线不卡| 国产精品人妻久久久影院| 国产一区亚洲一区在线观看| 黄色怎么调成土黄色| 国产一区有黄有色的免费视频| 亚洲国产精品999| 国产极品天堂在线| 久久ye,这里只有精品| 99久久综合免费| 一二三四在线观看免费中文在| 国产又爽黄色视频| 国产成人精品福利久久| 欧美精品人与动牲交sv欧美| 国产精品亚洲av一区麻豆 | 啦啦啦在线观看免费高清www| 亚洲情色 制服丝袜| 在线观看国产h片| 久久韩国三级中文字幕| 飞空精品影院首页| 男女无遮挡免费网站观看| 啦啦啦中文免费视频观看日本| 成年动漫av网址| 不卡视频在线观看欧美| 免费人妻精品一区二区三区视频| 日韩中文字幕欧美一区二区 | 女人爽到高潮嗷嗷叫在线视频| 丝袜喷水一区| 尾随美女入室| 美女国产高潮福利片在线看| 亚洲精品美女久久av网站| 中文精品一卡2卡3卡4更新| 免费在线观看完整版高清| 黄色怎么调成土黄色| 汤姆久久久久久久影院中文字幕| 黄色 视频免费看| 最近的中文字幕免费完整| 各种免费的搞黄视频| 欧美日韩亚洲综合一区二区三区_| 国产精品二区激情视频| 亚洲精品久久成人aⅴ小说| 波多野结衣av一区二区av| 亚洲图色成人| 青春草国产在线视频| 久久综合国产亚洲精品| 悠悠久久av| 亚洲免费av在线视频| 免费人妻精品一区二区三区视频| 精品酒店卫生间| 亚洲精品一区蜜桃| 欧美日韩视频高清一区二区三区二| 精品卡一卡二卡四卡免费| 18在线观看网站| 精品少妇一区二区三区视频日本电影 | 国产精品一二三区在线看| 日本欧美视频一区| 国产不卡av网站在线观看| 51午夜福利影视在线观看| 国产在视频线精品| 777米奇影视久久| 国产成人精品久久二区二区91 | 亚洲成人手机| 一级片'在线观看视频| 久久青草综合色| 少妇人妻久久综合中文| 一级毛片 在线播放| 国产一区二区在线观看av| 色综合欧美亚洲国产小说| 丰满乱子伦码专区| 天美传媒精品一区二区| 国产1区2区3区精品| av国产久精品久网站免费入址| 国产免费一区二区三区四区乱码| 欧美 亚洲 国产 日韩一| 亚洲欧美日韩另类电影网站| 久热爱精品视频在线9| 韩国av在线不卡| 婷婷色综合www| www日本在线高清视频| 国产精品嫩草影院av在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲三区欧美一区| 曰老女人黄片| 中文字幕亚洲精品专区| 日韩电影二区| h视频一区二区三区| 丝袜喷水一区| 蜜桃在线观看..| 精品亚洲成国产av| 热re99久久国产66热| 色婷婷久久久亚洲欧美| 五月天丁香电影| 在线 av 中文字幕| 一区二区三区乱码不卡18| 国产精品嫩草影院av在线观看| 国产在视频线精品| 高清不卡的av网站| 两个人看的免费小视频| 卡戴珊不雅视频在线播放| 91老司机精品| xxxhd国产人妻xxx| 精品一区在线观看国产| 黑人猛操日本美女一级片| 大片免费播放器 马上看| 国产成人a∨麻豆精品| av电影中文网址| 丝袜脚勾引网站| 美女视频免费永久观看网站| 国产精品久久久久成人av| 9191精品国产免费久久| 黄色怎么调成土黄色| 亚洲一级一片aⅴ在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 夫妻午夜视频| 91精品伊人久久大香线蕉| 99九九在线精品视频| 国产精品久久久久久精品电影小说| 老司机亚洲免费影院| 国产成人av激情在线播放| 国产99久久九九免费精品| 中文字幕高清在线视频| 欧美国产精品va在线观看不卡| 久久影院123| 亚洲 欧美一区二区三区| 人人妻人人添人人爽欧美一区卜| 青青草视频在线视频观看| 丝袜美足系列| 99久久综合免费| 下体分泌物呈黄色| 少妇 在线观看| 国产精品蜜桃在线观看| 丝袜脚勾引网站| 成年女人毛片免费观看观看9 | 久久综合国产亚洲精品| 黄片播放在线免费| svipshipincom国产片| 一个人免费看片子| 色网站视频免费| 一区在线观看完整版| 男女边吃奶边做爰视频| 久久 成人 亚洲| 中文字幕最新亚洲高清| 又黄又粗又硬又大视频| 免费看av在线观看网站| 久热这里只有精品99| 日本黄色日本黄色录像| 一级毛片黄色毛片免费观看视频| 欧美日韩亚洲综合一区二区三区_| 一级a爱视频在线免费观看| 亚洲自偷自拍图片 自拍| 国产在线一区二区三区精| 日本vs欧美在线观看视频| 91老司机精品| 精品第一国产精品| 亚洲综合精品二区| 大片免费播放器 马上看| 91精品三级在线观看| 高清视频免费观看一区二区| 秋霞伦理黄片| 久久久久精品人妻al黑| 国产日韩欧美视频二区| 两个人看的免费小视频| 丰满乱子伦码专区| 天美传媒精品一区二区| 久久久久精品性色| 亚洲精品一区蜜桃| 新久久久久国产一级毛片| 欧美日韩一区二区视频在线观看视频在线| 极品人妻少妇av视频| 老司机深夜福利视频在线观看 | 亚洲人成77777在线视频| 亚洲av日韩在线播放| 免费观看a级毛片全部| 欧美精品一区二区免费开放| 免费观看av网站的网址| 久久99精品国语久久久| 亚洲精品日韩在线中文字幕| 日韩制服骚丝袜av| 人人澡人人妻人| 久久久久网色| 97人妻天天添夜夜摸| 飞空精品影院首页| 国产午夜精品一二区理论片| 久久人妻熟女aⅴ| 日本av免费视频播放| svipshipincom国产片| 天天躁日日躁夜夜躁夜夜| 侵犯人妻中文字幕一二三四区| 午夜激情av网站| 久久狼人影院| 波野结衣二区三区在线| 欧美日韩成人在线一区二区| 亚洲四区av| 日韩中文字幕视频在线看片| 99热全是精品| 亚洲人成77777在线视频| 免费久久久久久久精品成人欧美视频| 晚上一个人看的免费电影| 成人亚洲欧美一区二区av| 母亲3免费完整高清在线观看| 99久久精品国产亚洲精品| 我的亚洲天堂| xxxhd国产人妻xxx| 校园人妻丝袜中文字幕| 精品少妇一区二区三区视频日本电影 | 水蜜桃什么品种好| 久久人人爽人人片av| 嫩草影视91久久| 在线观看免费视频网站a站| 国产一区二区激情短视频 | 看免费av毛片| 色吧在线观看| 国产日韩一区二区三区精品不卡| 中文欧美无线码| 色综合欧美亚洲国产小说| 久久综合国产亚洲精品| 久热爱精品视频在线9| 国产亚洲精品第一综合不卡| 黄片小视频在线播放| 午夜免费男女啪啪视频观看| 日韩欧美一区视频在线观看| 国产一区二区三区综合在线观看| 久久av网站| 亚洲,欧美精品.| 日韩,欧美,国产一区二区三区| 亚洲欧美激情在线| 天天影视国产精品| 啦啦啦 在线观看视频| 纵有疾风起免费观看全集完整版| 国产97色在线日韩免费| 在线观看免费高清a一片| 中国国产av一级| 国产精品一区二区精品视频观看| 自线自在国产av| 精品人妻熟女毛片av久久网站| 大香蕉久久成人网| 国产一区二区在线观看av| 最新的欧美精品一区二区| 欧美国产精品一级二级三级| 久久久久国产一级毛片高清牌| av在线播放精品| 欧美日韩av久久| 如何舔出高潮| 国产成人精品在线电影| 成年人免费黄色播放视频| 久久精品aⅴ一区二区三区四区| 精品久久蜜臀av无| 亚洲国产成人一精品久久久| 你懂的网址亚洲精品在线观看| 国产精品偷伦视频观看了| 99久久综合免费| 麻豆乱淫一区二区| 波野结衣二区三区在线| 午夜老司机福利片| 人妻人人澡人人爽人人| 久久久久视频综合| 中文字幕色久视频| 丰满乱子伦码专区| 欧美成人午夜精品| 考比视频在线观看| 成年美女黄网站色视频大全免费| 大码成人一级视频| 免费观看人在逋| 天天影视国产精品| 国产一区二区激情短视频 | 午夜福利免费观看在线| 人人妻人人澡人人看| 啦啦啦在线免费观看视频4| 美女福利国产在线| 悠悠久久av| 欧美激情极品国产一区二区三区| 午夜老司机福利片| 啦啦啦在线免费观看视频4| 韩国精品一区二区三区| 久久国产精品大桥未久av| 亚洲 欧美一区二区三区| 久久天堂一区二区三区四区| 熟女少妇亚洲综合色aaa.| 国产淫语在线视频| 亚洲欧洲日产国产| 80岁老熟妇乱子伦牲交| 制服人妻中文乱码| 欧美97在线视频| 精品一区二区三区四区五区乱码 | 国产熟女午夜一区二区三区| 狂野欧美激情性bbbbbb| 天天操日日干夜夜撸| 久久精品亚洲熟妇少妇任你| 看免费成人av毛片| 成人影院久久| 另类亚洲欧美激情| 成人18禁高潮啪啪吃奶动态图| 一区二区日韩欧美中文字幕| 中文精品一卡2卡3卡4更新| av在线老鸭窝| 国产精品亚洲av一区麻豆 | 别揉我奶头~嗯~啊~动态视频 | 男女之事视频高清在线观看 | 美女大奶头黄色视频| 在现免费观看毛片| 国产免费福利视频在线观看| 成人18禁高潮啪啪吃奶动态图| 午夜福利乱码中文字幕| 亚洲成国产人片在线观看| 亚洲伊人色综图| 国产精品成人在线| 青春草国产在线视频| 在线看a的网站| 亚洲自偷自拍图片 自拍| 一本一本久久a久久精品综合妖精| 最近最新中文字幕免费大全7| 亚洲精品成人av观看孕妇| 叶爱在线成人免费视频播放| 最近中文字幕高清免费大全6| 日韩电影二区| 色婷婷av一区二区三区视频| 亚洲欧美日韩另类电影网站| 99精国产麻豆久久婷婷| 一区二区三区四区激情视频| 久久狼人影院| 午夜老司机福利片| 婷婷成人精品国产| 大香蕉久久成人网| 国产99久久九九免费精品| 久久精品亚洲熟妇少妇任你| 日韩熟女老妇一区二区性免费视频| 久久精品国产综合久久久| 亚洲精品一二三| 亚洲国产av新网站| 日韩一区二区三区影片| 亚洲欧美一区二区三区黑人| 中文字幕色久视频| 国产免费一区二区三区四区乱码| 久久久欧美国产精品| 亚洲精品,欧美精品| 日韩成人av中文字幕在线观看| 欧美97在线视频| 美女主播在线视频| 国产在视频线精品| 亚洲精品一二三| 成年女人毛片免费观看观看9 | 婷婷色麻豆天堂久久| 热re99久久国产66热| 男女床上黄色一级片免费看| 国产一区二区 视频在线|