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

    柱狀成層地層中儀器偏心對陣列感應(yīng)測井響應(yīng)的影響

    2016-09-29 08:11:31魏寶君王成園常欣莉曹景強(qiáng)
    地球物理學(xué)報(bào) 2016年8期
    關(guān)鍵詞:絕緣層井眼偏心

    魏寶君,王成園,常欣莉,曹景強(qiáng)

    1 中國石油大學(xué)(華東)理學(xué)院,青島 266580 2 山東省高校新能源物理與材料科學(xué)重點(diǎn)實(shí)驗(yàn)室(中國石油大學(xué)),青島 266580

    ?

    柱狀成層地層中儀器偏心對陣列感應(yīng)測井響應(yīng)的影響

    魏寶君1,2,王成園1,常欣莉1,曹景強(qiáng)1

    1 中國石油大學(xué)(華東)理學(xué)院,青島266580 2 山東省高校新能源物理與材料科學(xué)重點(diǎn)實(shí)驗(yàn)室(中國石油大學(xué)),青島266580

    采用柱狀成層介質(zhì)的并矢Green函數(shù)高效計(jì)算井眼中偏心發(fā)射源產(chǎn)生的電磁波,并利用該算法模擬了陣列感應(yīng)測井儀器在井眼中偏心時(shí)的響應(yīng).在模擬時(shí)既考慮到了金屬心軸的存在也考慮到了線圈系所處絕緣層的存在,并重點(diǎn)分析了絕緣層的影響.數(shù)值模擬結(jié)果表明,對相對低泥漿電導(dǎo)率情況而言,無論線圈距長短,線圈系所處絕緣層對儀器響應(yīng)的影響均較小,可忽略不計(jì).而對于相對高泥漿電導(dǎo)率情況,線圈系所處絕緣層對儀器響應(yīng)的影響明顯,不可忽略,且線圈距越短影響越大,儀器偏心距越大影響越大.另由數(shù)值模擬結(jié)果可以看出,在相同泥漿電導(dǎo)率情況下,儀器線圈距越短偏心對線圈系響應(yīng)的影響越大.

    陣列感應(yīng)測井;偏心;并矢Green函數(shù);柱狀成層地層;金屬心軸;絕緣層

    The algorithm was based on the dyadic Green′s functions of cylindrical multilayered medium,in which the coils are supposed to be positioned in the 0-th layer,i.e.the insulating protection layer.The borehole is taken as the 1-th layer in the algorithm,which was used to be taken as the 0-th layer where the coils are supposed to be positioned in conventional analytical methods.Thus the new model adopted is more suitable for describing real electromagnetic well-logging tools,whose responses can not be directly simulated by conventional analytical methods.In the new developed analytical algorithm,the continuity condition of the tangential electromagnetic fields across the boundary between the insulating protection layer and the borehole fluid was used,thus the fields in the borehole can be expressed by the fields in the insulating layer.The relationship between the electromagnetic fields under the cylindrical coordinates relative to the borehole and to the tool was then used and finally the new algorithm was developed.

    The responses of centric and eccentric array induction logging tools in a borehole were simulated by the algorithm,and the difference between the responses of the tool with or without an insulating protection layer was compared.The observations include:(1)Given the borehole diameter and the conductivity of the drilling fluid and formation,the relationship between the receiver coil′s response and the eccentric distance was obtained for different arrays.The results showed that the effect of the insulating layer on the tool′s response is small and can be neglected when the fluid′s conductivity is relatively low,whether the transmitter-receiver spacing is long or short.The effect is big for the condition of relatively high fluid′s conductivity,and the shorter the transmitter-receiver spacing the bigger the effect,the bigger the eccentric distance the bigger the effect.(2)Given the conductivity of the drilling fluid and formation,the relationship among the coil′s geometric factor and the borehole′s diameter was obtained for different arrays with the tool being in the borehole center.Since the tool was positioned in the borehole center,the effect of the insulating layer was relatively small.(3)Given the conductivity of the drilling fluid and formation,the relationship among the coil′s geometric factor,borehole′s diameter and eccentricity ratio for different arrays was obtained.The results can be used to analyze the radial detection depth with different borehole diameter and eccentricity ratio.

    From all the results we can conclude that it is better to consider all the specific constructions of the tool and coils when simulating the responses of electromagnetic well-logging tools in a borehole,and it is better to adopt a model as much as closer to real tools,thus the simulated results will be much more precise.

    1 引言

    偏心對電磁測井儀器在井眼中響應(yīng)的影響一直是測井分析家極為關(guān)注的問題之一,對于測量數(shù)據(jù)的精確解釋具有重要價(jià)值(Gianzero,1978).Gianzero(1978)分析了偏心對傳統(tǒng)軸向感應(yīng)線圈響應(yīng)特性的影響,其給出的解析表達(dá)式包含兩層介質(zhì),可適用于僅存在井眼和均質(zhì)地層的模型,沒有考慮金屬心軸和絕緣保護(hù)層等因素的影響.另外,他在分析時(shí)將線圈系近似為磁偶極子,忽略了線圈具體尺寸,故模擬結(jié)果與真實(shí)結(jié)果存在差別.Hagiwara等(2005)基于Lovell和Chew(1987)的方法模擬了含心軸、井眼和侵入帶的地層模型中傾斜線圈系在井眼內(nèi)居中時(shí)的響應(yīng).Wang等(2008)利用給出的計(jì)算柱狀成層各向異性介質(zhì)(橫向各向同性)中任意位置點(diǎn)源響應(yīng)的解析表達(dá)式模擬了多分量感應(yīng)測井儀器在井眼中居中和偏心時(shí)的響應(yīng)特性,模擬時(shí)亦將線圈系近似為磁偶極子,且忽略了線圈具體尺寸以及金屬心軸和絕緣保護(hù)層等因素.Lovell和Chew(1990)、Hue和Teixeira(2006)基于Lovell和Chew(1987)給出的計(jì)算柱狀成層介質(zhì)中任意位置點(diǎn)源響應(yīng)的解析表達(dá)式,推導(dǎo)出了含心軸電磁測井儀器在井眼中偏心時(shí)的響應(yīng),并分別分析了偏心對傳統(tǒng)軸向感應(yīng)線圈系和傾斜線圈系響應(yīng)特性的影響.上述推導(dǎo)假設(shè)線圈系位于儀器外側(cè)的井眼中,忽略了金屬心軸外面線圈系所處絕緣保護(hù)層的影響.對于實(shí)際電磁測井儀器如陣列感應(yīng)測井儀器而言,其金屬心軸外面一般含有一層絕緣介質(zhì)以保護(hù)線圈系.這層絕緣介質(zhì)占據(jù)了井眼中部分泥漿所占據(jù)的空間,其對儀器在井眼中居中和偏心時(shí)的響應(yīng)均將產(chǎn)生影響.盡管通過在線圈系內(nèi)側(cè)增加一層介質(zhì)的方式亦可直接利用上述推導(dǎo)模擬含絕緣保護(hù)層的儀器的響應(yīng),但仍須假設(shè)線圈系位于絕緣層外的井眼泥漿中而不是絕緣層內(nèi),與儀器實(shí)際結(jié)構(gòu)仍存在差距.本文基于均勻介質(zhì)和柱狀成層介質(zhì)中并矢Green函數(shù)的矢量本征函數(shù)展開式(戴振鐸和魯述,1996;魏寶君,2007;魏寶君等,2010,2012,2015;俞燕明等,2014),經(jīng)推導(dǎo)得到了計(jì)算含心軸和絕緣保護(hù)層的電磁測井儀器在井眼中偏心時(shí)響應(yīng)的解析表達(dá)式,并利用該解析方法模擬了陣列感應(yīng)測井儀器在井眼中偏心時(shí)的響應(yīng),重點(diǎn)分析了絕緣層對儀器響應(yīng)的影響.與Lovell和Chew(1990)、Hue和Teixeira(2006)假設(shè)線圈系位于儀器外側(cè)的井眼泥漿中相比,本文推導(dǎo)假設(shè)線圈系位于儀器內(nèi)部的絕緣層中,更接近于儀器實(shí)際結(jié)構(gòu).

    2 基本理論

    2.1均勻介質(zhì)中譜域并矢Green函數(shù)

    采用圓柱坐標(biāo)系并假設(shè)發(fā)射源隨時(shí)間的變化關(guān)系為exp(-iωt),其中ω為角頻率.戴振鐸和魯述(1996)、魏寶君(2007)、魏寶君等(2010,2012)給出了圓柱坐標(biāo)系下均勻介質(zhì)中頻率域并矢Green函數(shù)的矢量本征函數(shù)展開式,它們均可以表示為譜域并矢Green函數(shù)積分的形式.例如對電流源的電型和磁型并矢Green函數(shù)GEJ、GHJ,有:

    (1a)

    (1b)

    設(shè)電流源的體電流密度為JT,則介質(zhì)各處的電磁場可表示為:

    (2a)

    (2b)

    定義如下矩陣(Lovell and Chew,1990):

    (3a)

    (3b)

    (4a)

    (4b)

    (4c)

    (5a)

    (5b)

    (5c)

    結(jié)合式(1),頻率域內(nèi)不同類型和尺寸的源產(chǎn)生的電磁場可表示為

    (6a)

    (6b)

    (6c)

    2.2柱狀成層介質(zhì)中譜域并矢Green函數(shù)

    對于沿ρ方向柱狀成層介質(zhì),設(shè)源在第0層,介質(zhì)從里向外編號為-m,-m+1,…,0,1,…,n,半徑分別為ρ-m,ρ-m+1,…,ρ0,ρ1,…,ρn-1,類比式(5)可寫出任意第i層介質(zhì)中譜域并矢Green函數(shù)的z分量、φ分量和ρ分量為

    (7a)

    (7b)

    (7c)

    柱狀成層介質(zhì)各處頻率域內(nèi)不同類型和尺寸的源產(chǎn)生的電磁場可表示為

    (8a)

    (8b)

    (8c)

    2.3儀器在井眼中偏心時(shí)的譜域并矢Green函數(shù)

    圖1給出了電磁測井儀器軸心在井眼內(nèi)偏離井軸距離為d、角度為φE時(shí)的示意圖.不帶撇的坐標(biāo)系以井軸為軸心,帶撇的坐標(biāo)系以儀器軸為軸心.與Lovell和Chew(1990)、Hue和Teixeira (2006)假設(shè)線圈系位于儀器外側(cè)的井眼泥漿中、忽略線圈系所處絕緣保護(hù)層的地層模型相比,圖1中的發(fā)射線圈和接收線圈均在第0層介質(zhì)即絕緣層中(圖中虛線所示),而將井眼作為第1層介質(zhì).

    在帶撇的儀器坐標(biāo)系中,設(shè)發(fā)射源徑向坐標(biāo)為

    圖1 井眼內(nèi)儀器偏心示意圖Fig.1 Illustration of eccentric borehole geometry

    (9a)

    (9b)

    (9c)

    頻率域內(nèi)不同類型和尺寸的源產(chǎn)生的電磁場可表示為

    (10a)

    (10b)

    (10c)

    (11a)

    (11c)

    頻率域內(nèi)不同類型和尺寸的源產(chǎn)生的電磁場可表示為

    (12a)

    (12b)

    (12c)

    在帶撇的儀器坐標(biāo)系中,第1層介質(zhì)即井眼中譜域并矢Green函數(shù)的z分量、φ′分量和ρ′分量可表示為如下通式:

    (13a)

    (13b)

    (13c)

    式中a1ν′和b1ν′為待定系數(shù).

    頻率域內(nèi)不同類型和尺寸的源產(chǎn)生的電磁場可表示為

    (14a)

    (14b)

    (14c)

    在交界面ρ′=ρ0處,電磁場z分量和φ′分量連續(xù),由式(12a)和式(14a)相等、式(12b)和式(14b)相等得到

    (15a)

    (15b)

    (16a)

    (16b)

    在不帶撇的地層坐標(biāo)系中,根據(jù)式(7),在第1層介質(zhì)即井眼中半徑大于發(fā)射源的位置處,譜域并矢Green函數(shù)的z分量、φ分量和ρ分量可表示為

    (17a)

    (17b)

    (17c)

    頻率域內(nèi)不同類型和尺寸的源產(chǎn)生的電磁場可表示為

    (18a)

    (18b)

    (18c)

    根據(jù)針對Bessel函數(shù)的Graf加法定理(Abramowitz and Stegun,1970),

    將上述兩式與式(18a)對比,有:

    (19)

    當(dāng)發(fā)射源和場點(diǎn)均在第0層介質(zhì)中時(shí),由于電磁場包含背景項(xiàng),式(19)的收斂速度較慢,為加快式(19)的收斂速度,令:

    (20)

    (ν=-νmax,-νmax+1,…,νmax),

    (21)

    考慮到式(20),式(10)或式(12)扣除背景項(xiàng)后的散射部分表達(dá)式相同,為

    (22a)

    (22b)

    (22c)

    2.4陣列感應(yīng)測井儀器響應(yīng)的計(jì)算

    利用式(22)并結(jié)合接收線圈的形狀和尺寸,可以計(jì)算含心軸和絕緣保護(hù)層的電磁測井儀器在井眼中偏心時(shí)的感應(yīng)電動(dòng)勢,儀器的線圈系可以是任意結(jié)構(gòu),也可以是任意方向.本文以陣列感應(yīng)測井儀器為例計(jì)算其偏心響應(yīng).設(shè)陣列感應(yīng)測井儀器發(fā)射電流強(qiáng)度為IT,軸向發(fā)射線圈匝數(shù)為NT、中心點(diǎn)軸向坐標(biāo)為zT、軸向?qū)挾葹棣T、線圈半徑為ρT,軸向接收線圈匝數(shù)為NR、中心點(diǎn)軸向坐標(biāo)為zR、軸向?qū)挾葹棣R、線圈半徑為ρR.在儀器坐標(biāo)系中發(fā)射線圈電流源沿eφ′方向,表達(dá)式為:

    (23)

    結(jié)合式(1a)、式(2a),有

    (24)

    接收線圈處感應(yīng)電動(dòng)勢的散射項(xiàng)為:

    (26)

    3 數(shù)值算例及討論

    下面利用上述解析方法模擬陣列感應(yīng)測井儀器不同子陣列的線圈系在井眼內(nèi)居中和偏心時(shí)的響應(yīng),重點(diǎn)對比忽略絕緣層和考慮絕緣層時(shí)線圈系的響應(yīng).采用4層介質(zhì)模型,-1層為金屬心軸,0層為絕緣層,1層為井眼,2層為均質(zhì)地層.

    取井眼直徑dh=0.20 m、泥漿電導(dǎo)率σ1=1.0 S·m-1、均質(zhì)地層電導(dǎo)率σ2=0.001 S·m-1,考慮絕緣層的存在,表1和表2分別給出了偏心率Φ為0.45和0.9時(shí)本文解析方法得到的視電導(dǎo)率與商用有限元軟件計(jì)算結(jié)果的對比.由表1和表2可以看出,兩種方法的計(jì)算結(jié)果是一致的,說明了本文解析方法的正確性.需指出的是,在同一計(jì)算機(jī)上運(yùn)算,本文解析方法比有限元軟件的計(jì)算速度快一個(gè)數(shù)量級以上,體現(xiàn)了該算法的高效性.

    表1 偏心率為0.45時(shí)解析方法與有限元計(jì)算結(jié)果的對比(單位:S·m-1)

    表2 偏心率為0.9時(shí)解析方法與有限元計(jì)算結(jié)果的對比(單位:S·m-1)

    取井眼直徑dh=0.20 m,圖2為泥漿電導(dǎo)率σ1=1.0 S·m-1、均質(zhì)地層電導(dǎo)率σ2=0.01 S·m-1時(shí)接收線圈的感應(yīng)電動(dòng)勢隨偏心距離d的變化關(guān)系.圖3為泥漿電導(dǎo)率σ1=0.001 S·m-1、均質(zhì)地層電導(dǎo)率σ2=1.0 S·m-1時(shí)接收線圈的感應(yīng)電動(dòng)勢隨偏心距離d的變化關(guān)系.兩幅圖均對無絕緣層和有絕緣層兩種情況下的響應(yīng)進(jìn)行了對比.由圖2可以看出,對于相對高泥漿電導(dǎo)率情況而言,絕緣層對線圈系響應(yīng)的影響明顯,且儀器偏心距越大影響越大,線圈距越短影響越大.由圖3可以看出,對于相對低泥漿電導(dǎo)率情況而言,無論線圈距長短,線圈系所處絕緣層對儀器響應(yīng)的影響均較小,可忽略不計(jì).另通過圖2和圖3可以看出,相同泥漿電導(dǎo)率情況下線圈距越短偏心對線圈系響應(yīng)的影響越大,同一組線圈系在相對高泥漿電導(dǎo)率井眼中偏心對其響應(yīng)的影響更大.

    取泥漿電導(dǎo)率σ1=1.0 S·m-1、地層電導(dǎo)率σ2=0.01 S·m-1,圖4為儀器在井眼中居中時(shí)線圈系幾何因子g隨井眼直徑dh的變化關(guān)系.正如前面所述,由于儀器居中,無論井眼直徑是大是小,線圈系所處絕緣層對儀器響應(yīng)的影響均較小.除線圈距最短的第一組線圈系響應(yīng)區(qū)別較明顯外,絕緣層對另外兩組線圈距較長的線圈系的響應(yīng)可忽略不計(jì).利用圖4可以分析儀器在井眼中居中時(shí)不同線圈系的徑向探測深度.

    圖2 泥漿電導(dǎo)率σ1=1.0 S·m-1、地層電導(dǎo)率σ2=0.01 S·m-1時(shí)接收線圈響應(yīng)隨偏心距離的變化關(guān)系(a)第一組;(b)第二組;(c)第三組.Fig.2 Relationship between receiver coil′s response and eccentric distance when the drilling fluid′s conductivity σ1=1.0 S·m-1 and the formation′s conductivity σ2=0.01 S·m-1

    圖3 泥漿電導(dǎo)率σ1=0.001 S·m-1、地層電導(dǎo)率σ2=1.0 S·m-1時(shí)接收線圈響應(yīng)隨偏心距離的變化關(guān)系(a)第一組;(b)第二組;(c)第三組.Fig.3 Relationship between receiver coil′s response and eccentric distance when the drilling fluid′s conductivity σ1=0.001 S·m-1 and the formation′s conductivity σ2=1.0 S·m-1

    圖4 儀器在井眼中居中時(shí)線圈系幾何因子隨井徑變化關(guān)系A(chǔ)1n—第一組線圈無絕緣層,A1y—第一組線圈有絕緣層,其余類推.Fig.4 Relationship between coil′s geometric factor and borehole diameter with tool being in the borehole center

    仍取泥漿電導(dǎo)率σ1=1.0 S·m-1、地層電導(dǎo)率σ2=0.01 S·m-1,圖5為第一組線圈系幾何因子g隨井徑dh及偏心率Φ的變化關(guān)系,圖6為第二組線圈系幾何因子隨井徑及偏心率的變化關(guān)系.由于在正常的井徑范圍內(nèi)線圈距最長的第三組線圈系的幾何因子接近于零,受井徑及偏心率的影響可忽略不計(jì),故沒有給出.對比圖5a、5b可以發(fā)現(xiàn),對于第一組短線圈距線圈系,在偏心率較小時(shí),無論井徑是大是小,無絕緣層時(shí)的幾何因子比有絕緣層時(shí)的幾何因子數(shù)值大.隨著偏心率的增加,兩種情況下的幾何因子數(shù)值均逐漸減小.在井徑較小時(shí),無絕緣層幾何因子減小速度比有絕緣層幾何因子減小速度要慢,但隨著井徑的增加,無絕緣層幾何因子減小速度逐漸加快,最終高于有絕緣層幾何因子的減小速度.總體上,第一組短線圈距線圈系有絕緣層幾何因子比無絕緣層幾何因子隨井徑及偏心率的變化要平緩一些.對比圖6a、6b可以發(fā)現(xiàn),絕緣層對第二組線圈系幾何因子的影響要小于對第一組線圈系的影響.在井徑較小時(shí),幾何因子隨偏心率的增加而增加.但隨著井徑的增加,偏心率較小時(shí)幾何因子增加速度更快,偏心率越大幾何因子隨井徑增加而增加的速度越慢甚至出現(xiàn)減小現(xiàn)象.因而在井徑較大時(shí),幾何因子隨偏心率的增加而減小.利用圖5、圖6可以分析儀器在不同直徑井眼中偏心程度不同時(shí)不同線圈系的徑向探測深度.

    4 結(jié)論

    (1)本文基于柱狀成層介質(zhì)的并矢Green函數(shù)

    圖5 第一組線圈系幾何因子隨井徑及偏心率的變化關(guān)系(a)無絕緣層;(b)有絕緣層.Fig.5 Relationship between coil′s geometric factor,borehole diameter and eccentricity ratio for the first coils

    圖6 第二組線圈系幾何因子隨井徑及偏心率變化關(guān)系(a)無絕緣層;(b)有絕緣層.Fig.6 Relationship between coil′s geometric factor,borehole diameter and eccentricity ratio for the second coils

    得到解析表達(dá)式,可高效模擬電磁測井儀器在井眼中偏心時(shí)的電磁場分布及響應(yīng),并且在模擬時(shí)既可考慮金屬心軸的存在也可考慮線圈系所處絕緣保護(hù)層的存在.由于在推導(dǎo)過程中假設(shè)線圈系位于儀器內(nèi)部的絕緣層內(nèi)而非儀器外側(cè)的井眼中,所采用模型更接近于儀器實(shí)際結(jié)構(gòu).

    (2)對于相對高泥漿電導(dǎo)率情況而言,絕緣層對線圈系響應(yīng)的影響明顯,且儀器偏心距越大影響越大,線圈距越短影響越大.對于相對低泥漿電導(dǎo)率情況而言,無論線圈距長短,線圈系所處絕緣層對儀器響應(yīng)的影響均較小,可忽略不計(jì).相同泥漿電導(dǎo)率情況下線圈距越短偏心對線圈系響應(yīng)的影響越大,同一組線圈系在相對高泥漿電導(dǎo)率井眼中偏心對其響應(yīng)的影響更大.

    References

    Abramowitz M,Stegun I A.1970.Handbook of Mathematical Functions.New York:Dover Publications.

    Dai Z D,Lu S.1996.Dyadic Green′s Functions in Electromagnetic Theory (in Chinese).Wuhan:Wuhan University Press.

    Gianzero S C.1978.Effect of sonde eccentricity on responses of conventional induction-logging tools.IEEE Transactions on Geoscience Electronics,16(4):332-339.

    Hagiwara T,Banning E J,Ostermeier R M,et al.2005.Effects of mandrel,borehole,and invasion for tilt-coil antennas.//SPE Reservoir Evaluation &Engineering.SPE.

    Hue Y K,Teixeira F L.2006.Analysis of tilted-coil eccentric borehole antennas in cylindrical multilayered formations for well-logging applications.IEEE Trans.Antennas Propag.,54(4):1058-1064.

    Lovell J R,Chew W C.1987.Response of a point source in a multicylindrcally layered medium.IEEE Trans.Geosci.Remote Sens.,GE-25(6):850-858.

    Lovell J R,Chew W C.1990.Effect of tool eccentricity on some electrical well-logging tools.IEEE Trans.Geosci.Remote Sens.,28(1):127-136.

    Wang H N,So P,Yang S W,et al.2008.Numerical modeling of multicomponent induction well-logging tools in the cylindrically stratified anisotropic media.IEEE Trans.Geosci.Remote Sens.,46(4):1134-1147.

    Wei B J.2007.Response and calibration of a new logging-while-drilling resistivity tool.Chinese Journal of Geophysics (in Chinese),50(2):632-641.

    Wei B J,Tian K,Zhang X,et al.2010.Evaluating influence of eccentricity on response of electromagnetic wave resistivity logging-while-drilling by vector eigenfunction expansion formulae for dyadic Green′s functions.Journal of China University of Petroleum (Edition of Natural Science)(in Chinese),34(5):57-62.

    Wei B J,Ou Y F,Wu Y,et al.2012.Simulation of tilted coil′s response in cylindrically stratified media and its application for electromagnetic propagation measurement-while-drilling.Journal of China University of Petroleum (Edition of Natural Science)(in Chinese),36(5):72-79.

    Wei B J,Wang C Y,Yu Y M,et al.2015.Simulating multi-component induction logging tool′s calibration by vector eigenfunction expansion formulae for dyadic Green′s functions.Journal of China University of Petroleum (Edition of Natural Science)(in Chinese),39(2):39-47.

    Yu Y M,Xiao J Q,Wei B J,et al.2014.Evaluating influence of metal mandrel on response of multi-component induction logging by vector eigenfunction expansion formulae for dyadic Green′s functions.Journal of China University of Petroleum (Edition of Natural Science)(in Chinese),38(4):57-64.

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

    戴振鐸,魯述.1996.電磁理論中的并矢格林函數(shù).武漢:武漢大學(xué)出版社.

    魏寶君.2007.一種新型隨鉆電阻率測井儀器的響應(yīng)和刻度.地球物理學(xué)報(bào),50(2):632-641.

    魏寶君,田坤,張旭等.2010.用并矢Green函數(shù)的矢量本征函數(shù)展開式評價(jià)偏心對隨鉆電磁波電阻率測井響應(yīng)的影響.中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,34(5):57-62.

    魏寶君,歐永峰,武楊等.2012.柱狀成層介質(zhì)中傾斜線圈響應(yīng)的模擬及其在電磁波傳播隨鉆測量中的應(yīng)用.中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,36(5):72-79.

    魏寶君,王成園,俞燕明等.2015.用并矢Green函數(shù)的矢量本征函數(shù)展開式模擬多分量感應(yīng)測井儀器的刻度.中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,39(2):39-47.

    俞燕明,肖加奇,魏寶君等.2014.用并矢Green函數(shù)的矢量本征函數(shù)展開式評價(jià)金屬心軸對多分量感應(yīng)測井響應(yīng)的影響.中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,38(4):57-64.

    (本文編輯何燕)

    2016年(首屆)“中國區(qū)域地震學(xué)參考模型”學(xué)術(shù)交流會(huì)紀(jì)要

    2016年6月27日至29日,國家自然科學(xué)基金委員會(huì)“中國區(qū)域地震學(xué)參考模型工作組”(以下簡稱“工作組”)在中國科學(xué)院地質(zhì)與地球物理研究所成功舉辦首屆“中國區(qū)域地震學(xué)參考模型”學(xué)術(shù)交流會(huì).會(huì)議匯聚了來自35個(gè)國內(nèi)外高校、科研院所和地震局系統(tǒng)的專家學(xué)者313人.《地球物理學(xué)報(bào)》編輯部何燕參加了會(huì)議.

    國家自然科學(xué)基金委員會(huì)地球物理和空間物理學(xué)科主任于晟教授致開幕辭,提出建立中國區(qū)域地震學(xué)參考模型的任務(wù)和培養(yǎng)地震學(xué)專業(yè)科研人才的期望.國內(nèi)外共24位科研一線的地震學(xué)專家學(xué)者作了大會(huì)報(bào)告.大會(huì)報(bào)告分別是美國地質(zhì)調(diào)查局(USGS)地震學(xué)家Walter Mooney的“Global and regional crustal and lithospheric models:a synthesis”,法國巴黎地球物理研究所(IPGP)地震學(xué)家Nikolai Shapiro的“Tomography based on cross-correlations of the ambient seismic noise:accounting for inhomogeneous source distribution”,美國密蘇里科技大學(xué)(MUST)Steve Gao的“Crustal and mantle seismic anisotropy beneath the Tibetan Plateau and adjacent areas from shear-wave splitting analyses”,中國地質(zhì)大學(xué)(武漢)Lupei Zhu的“Receiver function studies of crustal structure of China—a review and prospect”,中國地震局地震預(yù)測研究所高原的“利用近場記錄剪切波分裂研究中國大陸中上地殼地震各向異性”,北京大學(xué)陳永順的“地震學(xué)參考模型探討:以東北地區(qū)為例—NECESSArray”,中國科學(xué)院地質(zhì)與地球物理研究所趙亮的“中國東部上地幔速度結(jié)構(gòu)與變形特征”,中國地震局地球物理研究所李永華的“中國東北地區(qū)地殼上地幔模型研究”,中國科學(xué)院測量與地球物理研究所儲(chǔ)日升的“青藏高原巖石圈結(jié)構(gòu)研究”,中國地震局地球物理勘探中心王夫運(yùn)的“基于人工地震測深構(gòu)建中國大陸地殼三維模型”,中國石油大學(xué)(北京)唐有彩的“鄂爾多斯地區(qū)地殼結(jié)構(gòu)與演化的地震學(xué)證據(jù)”,中國科學(xué)院地質(zhì)與地球物理研究所陳凌的“華北巖石圈結(jié)構(gòu)及其全球?qū)Ρ取?,中國科學(xué)院地質(zhì)與地球物理研究所李娟的“東北亞地區(qū)地幔過渡帶結(jié)構(gòu)研究、進(jìn)展及思考”,中國科學(xué)院地質(zhì)與地球物理研究所趙連鋒的“中國區(qū)域地震波衰減模型研究進(jìn)展”,中國地震局地殼應(yīng)力研究所雷建設(shè)的“青藏高原東緣殼幔波速結(jié)構(gòu)與深層動(dòng)力學(xué)過程:騰沖火山與強(qiáng)震”,中國地質(zhì)大學(xué)(北京)李紅誼的“青藏高原及鄰區(qū)地殼上地幔速度結(jié)構(gòu)研究進(jìn)展”,成都理工大學(xué)梁春濤的“汶川和蘆山地震之間地震空區(qū)研究”,中國科學(xué)技術(shù)大學(xué)姚華建的“基于臺(tái)陣數(shù)據(jù)聯(lián)合反演的地震學(xué)速度模型構(gòu)建”,云南大學(xué)劉魁的“Global Rayleigh wave phase velocity maps from finite-frequency tomography”,中國地震局地球物理研究所鄭秀芬的“國家測震臺(tái)網(wǎng)地震波形數(shù)據(jù)共享”,中國地震局地球物理研究所丁志峰的“南北地震帶的地震臺(tái)陣探測及其深部結(jié)構(gòu)模型”,中國科學(xué)院地質(zhì)與地球物理研究所鄭天愉的“華北地區(qū)地殼-上地幔地震波速度結(jié)構(gòu)模型”,中國地震局地質(zhì)所劉啟元的“China Crust 0.5:關(guān)于中國大陸地殼模型研究項(xiàng)目的建議”,中國科學(xué)技術(shù)大學(xué)溫聯(lián)星的“中國區(qū)域地震學(xué)參考模型:現(xiàn)狀”.所有報(bào)告材料和影像資料均可在中國區(qū)域地震學(xué)參考模型的網(wǎng)站(http://chinageorefmodel.org)下載.報(bào)告內(nèi)容從數(shù)據(jù)采集、整合、處理、分析到結(jié)果,從地殼到上地幔,再到核幔邊界和地核,從地震波速度到各向異性和衰減研究,從不同區(qū)域的地震學(xué)參考模型到全球模型,報(bào)告精彩紛呈,會(huì)場討論熱烈.除大會(huì)報(bào)告外,還有110個(gè)展板報(bào)告.展板內(nèi)容豐富、敘述和討論嚴(yán)謹(jǐn),展示了地震學(xué)研究的新進(jìn)展.最后,中國科學(xué)技術(shù)大學(xué)溫聯(lián)星教授主持了“工作組”總結(jié)會(huì)議.

    通過這次會(huì)議報(bào)告和交流,“工作組”總結(jié)了中國區(qū)域地震學(xué)參考模型的研究現(xiàn)狀,確定了建立中國區(qū)域地震學(xué)參考模型的近期目標(biāo)是現(xiàn)有模型匯總、公開和檢驗(yàn),工作重點(diǎn)是數(shù)據(jù)整合、方法優(yōu)化和參考模型的逐步建立.這次會(huì)議對于凝聚國內(nèi)地震學(xué)研究力量、形成利用集群優(yōu)勢共同攻關(guān)的學(xué)科戰(zhàn)略起到了重要的推動(dòng)作用.同時(shí),搭建了旨在建立中國區(qū)域地震學(xué)參考模型的學(xué)術(shù)交流平臺(tái),有益于促進(jìn)科研、教育和相關(guān)文化事業(yè)的發(fā)展,帶動(dòng)地震學(xué)研究和科研文化的創(chuàng)新.

    (中國科學(xué)院地球與行星物理重點(diǎn)實(shí)驗(yàn)室趙連鋒供稿)

    Effect of tool eccentricity on responses of array induction logging tools in cylindrical multilayered formations

    WEI Bao-Jun1,2,WANG Cheng-Yuan1,CHANG Xin-Li1,CAO Jing-Qiang1

    1 College of Science,China University of Petroleum,Qingdao 266580,China 2 Key Laboratory of New Energy Physics &Material Science in Universities of Shandong (China University of Petroleum), Qingdao 266580,China

    The effect of eccentricity on the responses of electromagnetic well-logging tools in a borehole is of great value for precisely interpreting the measured data.For real electromagnetic well-logging tools such as array induction-logging tools,usually both a metal mandrel and an insulating protection layer are included in the tools.Only the existence of the metal mandrel was considered for simulating and analyzing the effect of tool eccentricity in conventional analytical methods,and the influence of the insulating protection layer has been neglected.Since the insulating protection layer has occupied part of the space of the drilling fluid,its influence on the responses should also be considered whether the tool is centric or eccentric in a borehole.An analytical algorithm was developed that can simulate the responses of centric or eccentric electromagnetic well-logging tools with both a metal mandrel and an insulating protection layer in a borehole.

    Array induction logging;Eccentricity;Dyadic Green′s functions;Cylindrical multilayered formations;Metal mandrel;Insulating layer

    魏寶君,王成園,常欣莉等.2016.柱狀成層地層中儀器偏心對陣列感應(yīng)測井響應(yīng)的影響.地球物理學(xué)報(bào),59(8):3121-3132,

    10.6038/cjg20160834.

    Wei B J,Wang C Y,Chang X L,et al.2016.Effect of tool eccentricity on responses of array induction logging tools in cylindrical multilayered formations.Chinese J.Geophys.(in Chinese),59(8):3121-3132,doi:10.6038/cjg20160834.

    中國石油天然氣集團(tuán)公司科學(xué)研究與技術(shù)開發(fā)項(xiàng)目(2014D-4018);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(15CX05047A)資助.

    魏寶君,男,1969年生,教授,博士,教育部“新世紀(jì)優(yōu)秀人才支持計(jì)劃”入選者,研究方向?yàn)殡姶艤y井理論及應(yīng)用.

    E-mail:weibj@upc.edu.cn

    10.6038/cjg20160834

    P631

    2015-03-31,2015-06-23收修定稿

    猜你喜歡
    絕緣層井眼偏心
    剪切滑移裂縫對井眼聲波傳播的影響
    云南化工(2021年10期)2021-12-21 07:33:46
    伊拉克H 油田Sadi 油藏魚骨井井眼布置方案研究
    媽媽不偏心
    快樂語文(2018年7期)2018-05-25 02:32:15
    復(fù)合絕緣層漆包線熱老化壽命數(shù)據(jù)的統(tǒng)計(jì)分析
    電線電纜(2018年1期)2018-03-16 02:19:08
    一種水電容式電纜絕緣層缺陷的檢測方法研究
    電線電纜(2017年4期)2017-07-25 07:49:49
    長慶油田儲(chǔ)氣庫水平井大井眼鉆井技術(shù)
    偏心的母親
    小說月刊(2015年9期)2015-04-23 08:48:23
    受井眼約束帶接頭管柱的縱橫彎曲分析
    麥加輕軌站臺(tái)2.5 m絕緣層施工技術(shù)
    巧妙應(yīng)對老師的“偏心”
    人生十六七(2015年5期)2015-02-28 13:08:24
    欧美国产日韩亚洲一区| 又黄又粗又硬又大视频| 亚洲专区国产一区二区| 久久天堂一区二区三区四区| 久久亚洲真实| 欧美另类亚洲清纯唯美| 亚洲av片天天在线观看| 久久久国产欧美日韩av| 嫩草影院精品99| 欧美日韩福利视频一区二区| 亚洲精品av麻豆狂野| 少妇熟女aⅴ在线视频| 乱人伦中国视频| 19禁男女啪啪无遮挡网站| 一区二区三区精品91| 色尼玛亚洲综合影院| 色婷婷久久久亚洲欧美| av在线播放免费不卡| 免费观看人在逋| 久久香蕉激情| 国产激情欧美一区二区| 神马国产精品三级电影在线观看 | 中亚洲国语对白在线视频| 亚洲第一欧美日韩一区二区三区| 亚洲一区中文字幕在线| 男女做爰动态图高潮gif福利片 | 免费在线观看视频国产中文字幕亚洲| 美女大奶头视频| 国产精品自产拍在线观看55亚洲| 十分钟在线观看高清视频www| 中文字幕高清在线视频| 久久欧美精品欧美久久欧美| 亚洲av日韩精品久久久久久密| 国产欧美日韩一区二区精品| 丝袜美足系列| 色在线成人网| 国产成人精品久久二区二区91| 久久久久国产一级毛片高清牌| 美国免费a级毛片| 成人av一区二区三区在线看| 露出奶头的视频| av福利片在线| 在线观看舔阴道视频| 日韩精品青青久久久久久| 亚洲中文av在线| 欧美日韩福利视频一区二区| 国产伦一二天堂av在线观看| 日韩精品免费视频一区二区三区| www.www免费av| 黄色丝袜av网址大全| 一二三四在线观看免费中文在| 在线观看免费视频网站a站| 高清黄色对白视频在线免费看| www国产在线视频色| 国产成+人综合+亚洲专区| 亚洲精品在线观看二区| 国产高清有码在线观看视频 | 国产xxxxx性猛交| 国产成人系列免费观看| 日韩欧美一区视频在线观看| 他把我摸到了高潮在线观看| 99热只有精品国产| 国产成人精品在线电影| a级毛片在线看网站| 一级作爱视频免费观看| 精品日产1卡2卡| 人人澡人人妻人| netflix在线观看网站| 老熟妇乱子伦视频在线观看| 亚洲一区二区三区色噜噜| 亚洲成人精品中文字幕电影| 国产一区二区三区综合在线观看| 狠狠狠狠99中文字幕| 欧美成人免费av一区二区三区| 欧美性长视频在线观看| 一区二区日韩欧美中文字幕| 国产亚洲av嫩草精品影院| 精品第一国产精品| ponron亚洲| 国产一区二区激情短视频| 人人妻人人澡人人看| 久久久久久久久中文| 亚洲av成人不卡在线观看播放网| 欧美激情高清一区二区三区| svipshipincom国产片| 又紧又爽又黄一区二区| 亚洲精品美女久久久久99蜜臀| 真人一进一出gif抽搐免费| 啦啦啦观看免费观看视频高清 | 99久久精品国产亚洲精品| 欧美精品亚洲一区二区| 一区福利在线观看| 日韩高清综合在线| 中文字幕高清在线视频| 日韩精品免费视频一区二区三区| 真人一进一出gif抽搐免费| 99精品在免费线老司机午夜| 村上凉子中文字幕在线| 一边摸一边抽搐一进一小说| 国产精品野战在线观看| 亚洲欧美日韩无卡精品| 久久亚洲真实| 精品一区二区三区av网在线观看| 女人被躁到高潮嗷嗷叫费观| 国产精品精品国产色婷婷| 亚洲欧美激情综合另类| 黑人操中国人逼视频| 国产野战对白在线观看| 51午夜福利影视在线观看| 久久草成人影院| 岛国在线观看网站| 亚洲第一青青草原| 国产成人一区二区三区免费视频网站| www日本在线高清视频| 亚洲国产看品久久| 欧美成狂野欧美在线观看| 精品国产乱子伦一区二区三区| 亚洲一区二区三区色噜噜| 啦啦啦 在线观看视频| 99香蕉大伊视频| 亚洲中文日韩欧美视频| 亚洲国产精品sss在线观看| 亚洲精品中文字幕一二三四区| 成人三级做爰电影| 一区二区三区激情视频| 国产亚洲欧美在线一区二区| 亚洲第一欧美日韩一区二区三区| 午夜两性在线视频| 女生性感内裤真人,穿戴方法视频| 国产麻豆成人av免费视频| 韩国av一区二区三区四区| 丁香欧美五月| 桃红色精品国产亚洲av| 色播在线永久视频| 久久狼人影院| 精品欧美国产一区二区三| 国产激情欧美一区二区| 国产精品久久久av美女十八| 国产成人影院久久av| 淫妇啪啪啪对白视频| 纯流量卡能插随身wifi吗| 久久久久久久午夜电影| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩黄片免| 视频区欧美日本亚洲| 欧美色欧美亚洲另类二区 | cao死你这个sao货| 国内精品久久久久精免费| 日本免费一区二区三区高清不卡 | 日本在线视频免费播放| 老司机靠b影院| 久久久久久大精品| 国产片内射在线| 日本黄色视频三级网站网址| 叶爱在线成人免费视频播放| 国产精品一区二区在线不卡| 午夜影院日韩av| 中亚洲国语对白在线视频| 一级毛片女人18水好多| av在线播放免费不卡| 国产亚洲av嫩草精品影院| 亚洲av成人一区二区三| 久久婷婷人人爽人人干人人爱 | 女人高潮潮喷娇喘18禁视频| 嫁个100分男人电影在线观看| 黄色毛片三级朝国网站| 久久婷婷人人爽人人干人人爱 | 免费少妇av软件| 一边摸一边抽搐一进一小说| 亚洲一区二区三区不卡视频| 男女床上黄色一级片免费看| 成人亚洲精品av一区二区| 精品一区二区三区视频在线观看免费| 一夜夜www| 国产日韩一区二区三区精品不卡| 男女午夜视频在线观看| 欧美成狂野欧美在线观看| 咕卡用的链子| 91精品三级在线观看| 欧美丝袜亚洲另类 | 日韩欧美免费精品| www.熟女人妻精品国产| 99精品欧美一区二区三区四区| 午夜老司机福利片| 成人特级黄色片久久久久久久| 久久久久久大精品| 午夜成年电影在线免费观看| 老熟妇仑乱视频hdxx| 搡老妇女老女人老熟妇| xxx96com| 国产真人三级小视频在线观看| 国产色视频综合| 国产精品亚洲美女久久久| 亚洲片人在线观看| 香蕉久久夜色| 国产亚洲av高清不卡| 一本综合久久免费| 国产精品98久久久久久宅男小说| 麻豆久久精品国产亚洲av| 涩涩av久久男人的天堂| 一二三四社区在线视频社区8| 亚洲av熟女| 亚洲av熟女| 精品久久久久久久毛片微露脸| 叶爱在线成人免费视频播放| 久久亚洲真实| 国产国语露脸激情在线看| 国产激情久久老熟女| 91精品国产国语对白视频| 激情在线观看视频在线高清| 国产成人欧美| 午夜福利一区二区在线看| 黑丝袜美女国产一区| 欧美乱妇无乱码| 色综合欧美亚洲国产小说| 日日摸夜夜添夜夜添小说| 女同久久另类99精品国产91| 日本欧美视频一区| 日本免费a在线| 中文字幕人妻熟女乱码| 欧美人与性动交α欧美精品济南到| 一区福利在线观看| 精品国产一区二区三区四区第35| 波多野结衣av一区二区av| 精品久久久久久,| 亚洲精品美女久久av网站| 亚洲精品中文字幕一二三四区| 后天国语完整版免费观看| 熟女少妇亚洲综合色aaa.| 人妻丰满熟妇av一区二区三区| 午夜福利影视在线免费观看| 欧美日韩黄片免| 精品人妻在线不人妻| 叶爱在线成人免费视频播放| 国产欧美日韩综合在线一区二区| 亚洲在线自拍视频| 亚洲专区国产一区二区| 亚洲欧美日韩无卡精品| 操出白浆在线播放| 成人特级黄色片久久久久久久| 香蕉丝袜av| 两个人看的免费小视频| avwww免费| 国产极品粉嫩免费观看在线| 亚洲电影在线观看av| 国内毛片毛片毛片毛片毛片| 啦啦啦韩国在线观看视频| 日韩精品免费视频一区二区三区| 国产97色在线日韩免费| 天天躁夜夜躁狠狠躁躁| 亚洲色图 男人天堂 中文字幕| 老司机深夜福利视频在线观看| 人人妻人人爽人人添夜夜欢视频| 欧美中文综合在线视频| 久久久久亚洲av毛片大全| 大香蕉久久成人网| 国产精品一区二区精品视频观看| 少妇裸体淫交视频免费看高清 | 狠狠狠狠99中文字幕| 欧美成人性av电影在线观看| 在线观看免费日韩欧美大片| 一区二区三区精品91| 精品一区二区三区av网在线观看| 亚洲av五月六月丁香网| 国产高清视频在线播放一区| 午夜福利欧美成人| 深夜精品福利| 亚洲五月天丁香| 久久草成人影院| 一本综合久久免费| 日本五十路高清| 亚洲无线在线观看| 国产欧美日韩综合在线一区二区| 91老司机精品| 黄色 视频免费看| 国产一区二区三区综合在线观看| 国产欧美日韩一区二区精品| 怎么达到女性高潮| 国产一区二区三区视频了| 欧美日本中文国产一区发布| 一区二区三区精品91| 高潮久久久久久久久久久不卡| 亚洲欧美日韩另类电影网站| 国产精品,欧美在线| e午夜精品久久久久久久| 999精品在线视频| 最近最新免费中文字幕在线| 亚洲va日本ⅴa欧美va伊人久久| 女生性感内裤真人,穿戴方法视频| 午夜久久久久精精品| 国产在线观看jvid| 亚洲成a人片在线一区二区| 一边摸一边做爽爽视频免费| 母亲3免费完整高清在线观看| 亚洲国产欧美网| 亚洲久久久国产精品| 多毛熟女@视频| 精品久久久久久,| 国产精品美女特级片免费视频播放器 | 好看av亚洲va欧美ⅴa在| 成人国语在线视频| 成人三级做爰电影| 成人18禁在线播放| 国产99白浆流出| 亚洲第一青青草原| 一区二区三区精品91| 91精品国产国语对白视频| 免费看a级黄色片| 免费人成视频x8x8入口观看| 国内精品久久久久精免费| av中文乱码字幕在线| av视频在线观看入口| 人人澡人人妻人| 此物有八面人人有两片| 91av网站免费观看| 天天躁夜夜躁狠狠躁躁| 成人亚洲精品一区在线观看| 久久香蕉激情| 中文字幕久久专区| 一级黄色大片毛片| 99在线人妻在线中文字幕| 午夜福利视频1000在线观看 | 999精品在线视频| 1024香蕉在线观看| 日韩欧美免费精品| 99久久精品国产亚洲精品| 两个人看的免费小视频| 久久久久国产精品人妻aⅴ院| 久久人人精品亚洲av| 成在线人永久免费视频| 亚洲一区高清亚洲精品| 天天一区二区日本电影三级 | 免费在线观看亚洲国产| 韩国av一区二区三区四区| 无遮挡黄片免费观看| 欧美成人性av电影在线观看| 国产成人啪精品午夜网站| 黑人巨大精品欧美一区二区mp4| 欧美黄色片欧美黄色片| 免费在线观看黄色视频的| 国产成+人综合+亚洲专区| 日本 av在线| 又大又爽又粗| 久久午夜综合久久蜜桃| 日本免费a在线| 日本一区二区免费在线视频| 国产精品一区二区免费欧美| 黄网站色视频无遮挡免费观看| 亚洲性夜色夜夜综合| 亚洲精品av麻豆狂野| 久久精品影院6| 精品一区二区三区av网在线观看| 十八禁网站免费在线| 黑人欧美特级aaaaaa片| 亚洲,欧美精品.| 国产片内射在线| 两个人免费观看高清视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产野战对白在线观看| 狂野欧美激情性xxxx| 成人特级黄色片久久久久久久| 久久天堂一区二区三区四区| 母亲3免费完整高清在线观看| 日日爽夜夜爽网站| 亚洲精品国产精品久久久不卡| 人妻丰满熟妇av一区二区三区| 免费在线观看完整版高清| 午夜久久久久精精品| 国产精品久久久久久精品电影 | 日韩大尺度精品在线看网址 | 黑人操中国人逼视频| 一个人免费在线观看的高清视频| 啦啦啦观看免费观看视频高清 | 精品国产乱子伦一区二区三区| 免费在线观看完整版高清| 成人18禁高潮啪啪吃奶动态图| 在线永久观看黄色视频| 久久午夜综合久久蜜桃| 亚洲久久久国产精品| 夜夜躁狠狠躁天天躁| 中文字幕人妻丝袜一区二区| 亚洲色图综合在线观看| 国产成人精品在线电影| 久久精品亚洲精品国产色婷小说| 久久人人精品亚洲av| 操美女的视频在线观看| 97人妻精品一区二区三区麻豆 | 级片在线观看| 99热只有精品国产| 99精品在免费线老司机午夜| 90打野战视频偷拍视频| 国产免费男女视频| 男女下面插进去视频免费观看| 又黄又爽又免费观看的视频| 色播亚洲综合网| 一二三四在线观看免费中文在| 欧美成人一区二区免费高清观看 | av福利片在线| 精品人妻1区二区| 欧美丝袜亚洲另类 | 丰满人妻熟妇乱又伦精品不卡| 色综合亚洲欧美另类图片| 高清在线国产一区| 人人妻人人爽人人添夜夜欢视频| 日日夜夜操网爽| 女性生殖器流出的白浆| 国产成+人综合+亚洲专区| 亚洲av第一区精品v没综合| 日韩精品中文字幕看吧| 视频在线观看一区二区三区| 久久久久久免费高清国产稀缺| 一区在线观看完整版| 999久久久国产精品视频| 国产成人一区二区三区免费视频网站| 午夜福利18| www国产在线视频色| 亚洲精品在线观看二区| 久久精品成人免费网站| 日韩欧美国产在线观看| 欧美色欧美亚洲另类二区 | 亚洲伊人色综图| 国产一区二区三区在线臀色熟女| 日韩欧美在线二视频| 99久久久亚洲精品蜜臀av| 搡老妇女老女人老熟妇| 国产成人av激情在线播放| 国产亚洲精品第一综合不卡| 亚洲人成电影观看| 亚洲国产日韩欧美精品在线观看 | 757午夜福利合集在线观看| 亚洲欧美精品综合久久99| 天天躁夜夜躁狠狠躁躁| 性欧美人与动物交配| 亚洲国产看品久久| 精品第一国产精品| videosex国产| 欧美性长视频在线观看| 这个男人来自地球电影免费观看| 99久久久亚洲精品蜜臀av| 午夜老司机福利片| 日本欧美视频一区| 欧美大码av| 啦啦啦 在线观看视频| 9191精品国产免费久久| 欧美在线一区亚洲| 精品一区二区三区视频在线观看免费| 黄色片一级片一级黄色片| 男女午夜视频在线观看| 亚洲男人的天堂狠狠| 一进一出好大好爽视频| 丁香六月欧美| 国产亚洲精品久久久久5区| 曰老女人黄片| 香蕉国产在线看| 国产欧美日韩一区二区三区在线| 亚洲成国产人片在线观看| 一区二区日韩欧美中文字幕| 怎么达到女性高潮| 日韩视频一区二区在线观看| 9热在线视频观看99| 大香蕉久久成人网| 欧美亚洲日本最大视频资源| 人人妻人人澡人人看| 色在线成人网| 精品久久蜜臀av无| 男女午夜视频在线观看| 亚洲成a人片在线一区二区| 国产xxxxx性猛交| 亚洲欧美精品综合久久99| 亚洲少妇的诱惑av| 两人在一起打扑克的视频| 一区福利在线观看| 我的亚洲天堂| 国产精品九九99| 国产成+人综合+亚洲专区| av网站免费在线观看视频| 一级黄色大片毛片| bbb黄色大片| 亚洲美女黄片视频| 多毛熟女@视频| 亚洲精品美女久久久久99蜜臀| 国产亚洲精品av在线| 日韩av在线大香蕉| 色哟哟哟哟哟哟| 国产麻豆69| 亚洲国产精品sss在线观看| 亚洲第一电影网av| 国产精品九九99| 国产精品乱码一区二三区的特点 | 国产成人欧美| 少妇的丰满在线观看| 精品国产美女av久久久久小说| 丁香六月欧美| 人成视频在线观看免费观看| 亚洲国产欧美日韩在线播放| www国产在线视频色| 亚洲成人国产一区在线观看| 一级毛片高清免费大全| 两个人视频免费观看高清| 精品一区二区三区视频在线观看免费| 亚洲成a人片在线一区二区| 午夜福利一区二区在线看| 成人欧美大片| 日韩av在线大香蕉| 女人被躁到高潮嗷嗷叫费观| 欧美成人免费av一区二区三区| 人人妻人人爽人人添夜夜欢视频| 午夜福利影视在线免费观看| 嫁个100分男人电影在线观看| 波多野结衣高清无吗| 日韩欧美在线二视频| 久久香蕉国产精品| 亚洲av五月六月丁香网| 午夜精品久久久久久毛片777| 久久国产乱子伦精品免费另类| 制服人妻中文乱码| 亚洲av电影在线进入| 777久久人妻少妇嫩草av网站| 精品熟女少妇八av免费久了| 国产精品香港三级国产av潘金莲| 亚洲精品av麻豆狂野| 高清黄色对白视频在线免费看| 老汉色∧v一级毛片| 成年女人毛片免费观看观看9| 午夜福利,免费看| 国产亚洲精品一区二区www| 国产成人影院久久av| 国产av在哪里看| 高清在线国产一区| 淫秽高清视频在线观看| 久热爱精品视频在线9| 无人区码免费观看不卡| 在线永久观看黄色视频| 91在线观看av| 亚洲五月婷婷丁香| 午夜福利18| 亚洲在线自拍视频| 久久精品国产亚洲av高清一级| 麻豆久久精品国产亚洲av| 精品午夜福利视频在线观看一区| 久久久水蜜桃国产精品网| 国产精品国产高清国产av| 精品熟女少妇八av免费久了| 一区二区日韩欧美中文字幕| 色播在线永久视频| 一本久久中文字幕| 亚洲国产毛片av蜜桃av| 精品国产乱子伦一区二区三区| 久久精品国产亚洲av高清一级| 免费在线观看影片大全网站| 熟女少妇亚洲综合色aaa.| 久久久久久国产a免费观看| 欧美日韩精品网址| x7x7x7水蜜桃| 黄片小视频在线播放| 国产伦一二天堂av在线观看| tocl精华| 美女午夜性视频免费| 日日爽夜夜爽网站| 日本精品一区二区三区蜜桃| 热re99久久国产66热| 这个男人来自地球电影免费观看| 男女午夜视频在线观看| 精品一区二区三区av网在线观看| 亚洲七黄色美女视频| 一区二区三区高清视频在线| 亚洲精品国产一区二区精华液| 国产一区二区在线av高清观看| 精品福利观看| 99久久精品国产亚洲精品| 亚洲av美国av| 色婷婷久久久亚洲欧美| 一进一出抽搐动态| 午夜免费观看网址| 搞女人的毛片| 在线观看免费午夜福利视频| av视频免费观看在线观看| 精品卡一卡二卡四卡免费| 日本 欧美在线| 涩涩av久久男人的天堂| √禁漫天堂资源中文www| 欧美黄色片欧美黄色片| 成人精品一区二区免费| 啪啪无遮挡十八禁网站| 国产精品久久久久久亚洲av鲁大| 人妻久久中文字幕网| 久久精品国产清高在天天线| 看片在线看免费视频| 18禁裸乳无遮挡免费网站照片 | 欧美午夜高清在线| 中文字幕精品免费在线观看视频| 99香蕉大伊视频| 他把我摸到了高潮在线观看| 老司机深夜福利视频在线观看| 欧美日韩福利视频一区二区| 免费av毛片视频| 久久精品亚洲熟妇少妇任你| 欧美日韩福利视频一区二区| 中文字幕人成人乱码亚洲影| 色综合婷婷激情| 国产色视频综合| 久久亚洲精品不卡| 午夜福利影视在线免费观看| 欧美日韩一级在线毛片| 免费av毛片视频| 亚洲最大成人中文| 国产麻豆69| 亚洲精品美女久久久久99蜜臀| 黄色成人免费大全| 中文字幕色久视频| 久久久久久大精品| 身体一侧抽搐| 亚洲精品国产色婷婷电影| 亚洲七黄色美女视频| 乱人伦中国视频|