• <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
    亚洲人成77777在线视频| 永久免费av网站大全| 久久精品国产自在天天线| 97在线视频观看| 黄色毛片三级朝国网站| 国产深夜福利视频在线观看| av有码第一页| 哪个播放器可以免费观看大片| 男女高潮啪啪啪动态图| 韩国av在线不卡| 久久狼人影院| 日韩人妻高清精品专区| 欧美另类一区| 日本av免费视频播放| 国产高清国产精品国产三级| 建设人人有责人人尽责人人享有的| 久久婷婷青草| 两个人的视频大全免费| 大码成人一级视频| 久久99热这里只频精品6学生| 午夜久久久在线观看| 日韩在线高清观看一区二区三区| 久久这里有精品视频免费| 精品一区二区三卡| 国产午夜精品一二区理论片| 国产精品人妻久久久久久| 人妻一区二区av| 熟女电影av网| 22中文网久久字幕| 欧美日韩视频精品一区| 99热网站在线观看| 亚洲内射少妇av| 中国国产av一级| 综合色丁香网| 午夜老司机福利剧场| 青春草国产在线视频| 80岁老熟妇乱子伦牲交| av播播在线观看一区| 国产视频首页在线观看| 亚洲av综合色区一区| 91精品一卡2卡3卡4卡| av网站免费在线观看视频| 国产成人精品婷婷| 国产黄色视频一区二区在线观看| 久久毛片免费看一区二区三区| 22中文网久久字幕| 亚洲欧洲日产国产| 久久久久国产精品人妻一区二区| 欧美+日韩+精品| 国产精品 国内视频| 有码 亚洲区| 曰老女人黄片| kizo精华| 大香蕉久久成人网| 国产成人a∨麻豆精品| 亚洲人成网站在线播| 日韩熟女老妇一区二区性免费视频| 国产欧美日韩综合在线一区二区| 一级毛片我不卡| 精品亚洲乱码少妇综合久久| 久久人妻熟女aⅴ| 一本色道久久久久久精品综合| 丰满少妇做爰视频| 美女中出高潮动态图| 寂寞人妻少妇视频99o| 日韩三级伦理在线观看| 高清欧美精品videossex| tube8黄色片| 久久人妻熟女aⅴ| 自拍欧美九色日韩亚洲蝌蚪91| 日本欧美国产在线视频| 日韩视频在线欧美| 18在线观看网站| 亚洲精品国产色婷婷电影| 最近的中文字幕免费完整| av视频免费观看在线观看| 国产永久视频网站| 国产亚洲最大av| 亚洲精品成人av观看孕妇| 国产精品人妻久久久久久| 免费人妻精品一区二区三区视频| 一二三四中文在线观看免费高清| 国产高清三级在线| 国产精品久久久久久久电影| 亚洲精品亚洲一区二区| 亚洲av.av天堂| 人妻人人澡人人爽人人| av在线老鸭窝| 2022亚洲国产成人精品| 黑人欧美特级aaaaaa片| 亚洲图色成人| 你懂的网址亚洲精品在线观看| 国产精品国产三级国产av玫瑰| 国产精品久久久久成人av| 中文字幕久久专区| 免费高清在线观看视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 久久久精品免费免费高清| 亚洲一区二区三区欧美精品| 日本爱情动作片www.在线观看| 久久久久久久国产电影| 午夜福利在线观看免费完整高清在| 在线观看三级黄色| 天堂俺去俺来也www色官网| 久久精品国产亚洲av涩爱| 免费观看a级毛片全部| 久久韩国三级中文字幕| 精品一区二区三卡| 亚洲精品日韩av片在线观看| 一级毛片我不卡| 日韩视频在线欧美| 午夜av观看不卡| 国产免费视频播放在线视频| 大又大粗又爽又黄少妇毛片口| 国产永久视频网站| 啦啦啦视频在线资源免费观看| 高清毛片免费看| 狠狠婷婷综合久久久久久88av| 国产成人精品婷婷| 美女视频免费永久观看网站| 中文字幕人妻丝袜制服| 9色porny在线观看| 乱人伦中国视频| 99国产综合亚洲精品| 日韩伦理黄色片| 大香蕉久久成人网| a级毛片在线看网站| 搡女人真爽免费视频火全软件| 亚洲内射少妇av| 美女内射精品一级片tv| 免费黄色在线免费观看| 老女人水多毛片| 成人免费观看视频高清| 18在线观看网站| 九色成人免费人妻av| av国产久精品久网站免费入址| 久久99热6这里只有精品| 日韩视频在线欧美| 亚洲成人av在线免费| 日本色播在线视频| 激情五月婷婷亚洲| 99久国产av精品国产电影| 精品国产乱码久久久久久小说| 日韩强制内射视频| 久久久久网色| 亚洲国产毛片av蜜桃av| 欧美日本中文国产一区发布| 免费观看av网站的网址| 国产视频内射| 99热6这里只有精品| av黄色大香蕉| 久久鲁丝午夜福利片| 成人黄色视频免费在线看| 搡老乐熟女国产| 亚洲欧美一区二区三区国产| 日本欧美视频一区| 精品国产露脸久久av麻豆| 最近最新中文字幕免费大全7| 黄色怎么调成土黄色| 建设人人有责人人尽责人人享有的| av专区在线播放| 亚洲精品亚洲一区二区| 久久影院123| 久久99热6这里只有精品| 少妇的逼水好多| 成人亚洲精品一区在线观看| 国产精品一国产av| 国产日韩一区二区三区精品不卡 | 男女国产视频网站| 夜夜骑夜夜射夜夜干| 亚洲欧美色中文字幕在线| 亚洲精品久久久久久婷婷小说| 黄色毛片三级朝国网站| 国产av一区二区精品久久| 免费观看a级毛片全部| 婷婷色综合www| 亚洲精华国产精华液的使用体验| 成人国语在线视频| 午夜激情福利司机影院| 美女脱内裤让男人舔精品视频| 精品一区在线观看国产| 涩涩av久久男人的天堂| 欧美人与性动交α欧美精品济南到 | 国产又色又爽无遮挡免| 亚洲av男天堂| 一边摸一边做爽爽视频免费| 亚洲第一区二区三区不卡| 久久ye,这里只有精品| 嘟嘟电影网在线观看| 国产有黄有色有爽视频| 亚洲精品久久午夜乱码| 亚洲人成网站在线观看播放| 各种免费的搞黄视频| 久久人人爽人人片av| 亚洲经典国产精华液单| 一级二级三级毛片免费看| 亚洲国产精品国产精品| 成人黄色视频免费在线看| 久久国内精品自在自线图片| 国产在线免费精品| 国产国拍精品亚洲av在线观看| 亚洲精品中文字幕在线视频| 成人毛片a级毛片在线播放| 日本爱情动作片www.在线观看| 又黄又爽又刺激的免费视频.| av播播在线观看一区| 夜夜骑夜夜射夜夜干| 人人妻人人爽人人添夜夜欢视频| 久久久国产一区二区| 如何舔出高潮| 这个男人来自地球电影免费观看 | 满18在线观看网站| 国产精品久久久久久av不卡| 高清视频免费观看一区二区| 高清av免费在线| 18禁观看日本| 国产黄频视频在线观看| 午夜91福利影院| 狠狠精品人妻久久久久久综合| 人人妻人人澡人人爽人人夜夜| 亚洲欧美成人综合另类久久久| 51国产日韩欧美| 国产成人精品在线电影| 黄片播放在线免费| 欧美97在线视频| 69精品国产乱码久久久| 亚洲精品第二区| 九九在线视频观看精品| 天美传媒精品一区二区| av视频免费观看在线观看| 菩萨蛮人人尽说江南好唐韦庄| 秋霞在线观看毛片| 日本-黄色视频高清免费观看| 天堂8中文在线网| 国产高清三级在线| 极品人妻少妇av视频| 建设人人有责人人尽责人人享有的| av.在线天堂| 天美传媒精品一区二区| 日韩精品免费视频一区二区三区 | 国产成人精品婷婷| 亚洲不卡免费看| 男女啪啪激烈高潮av片| 成年人午夜在线观看视频| 五月开心婷婷网| 欧美国产精品一级二级三级| 亚洲国产精品国产精品| 日韩伦理黄色片| 最近的中文字幕免费完整| 99国产综合亚洲精品| 观看美女的网站| 在线观看免费视频网站a站| 日本wwww免费看| 韩国高清视频一区二区三区| a级毛片在线看网站| 亚洲国产最新在线播放| 亚洲精品美女久久av网站| 国国产精品蜜臀av免费| 在线观看一区二区三区激情| 亚洲av男天堂| 久久精品国产鲁丝片午夜精品| 在线观看免费日韩欧美大片 | 国产精品一国产av| 国产精品国产三级国产专区5o| 午夜福利在线观看免费完整高清在| 亚洲欧美日韩卡通动漫| 亚洲在久久综合| 精品久久国产蜜桃| 91成人精品电影| 国产熟女午夜一区二区三区 | 26uuu在线亚洲综合色| 一级黄片播放器| 成年美女黄网站色视频大全免费 | 亚洲国产精品一区二区三区在线| av国产久精品久网站免费入址| 亚洲,欧美,日韩| 97超视频在线观看视频| 久久久久久久久久久丰满| 日本与韩国留学比较| 在线 av 中文字幕| 欧美日韩一区二区视频在线观看视频在线| 黄色毛片三级朝国网站| 看免费成人av毛片| 久久青草综合色| 少妇的逼水好多| 五月玫瑰六月丁香| 又黄又爽又刺激的免费视频.| 视频区图区小说| 一级a做视频免费观看| 久久精品国产亚洲av天美| 亚洲激情五月婷婷啪啪| 亚洲国产色片| 亚洲欧美一区二区三区黑人 | 日本爱情动作片www.在线观看| 人妻 亚洲 视频| 亚洲一级一片aⅴ在线观看| 熟妇人妻不卡中文字幕| 精品久久久久久久久亚洲| 99热6这里只有精品| 国产成人精品无人区| 丰满少妇做爰视频| 18+在线观看网站| 欧美日韩精品成人综合77777| 日韩成人av中文字幕在线观看| 不卡视频在线观看欧美| 精品午夜福利在线看| 热99国产精品久久久久久7| 国产精品人妻久久久久久| 草草在线视频免费看| 99热全是精品| 黄色一级大片看看| 97在线视频观看| 国产爽快片一区二区三区| 交换朋友夫妻互换小说| 人人澡人人妻人| 午夜福利视频在线观看免费| 国产 一区精品| 免费日韩欧美在线观看| 18+在线观看网站| 亚洲av成人精品一二三区| 丝袜脚勾引网站| 日韩,欧美,国产一区二区三区| 日韩三级伦理在线观看| 另类精品久久| 午夜福利网站1000一区二区三区| 视频中文字幕在线观看| 男人添女人高潮全过程视频| 国产精品一国产av| 美女cb高潮喷水在线观看| 人人妻人人爽人人添夜夜欢视频| 七月丁香在线播放| 人妻系列 视频| 最近中文字幕2019免费版| 日韩强制内射视频| 精品久久久噜噜| 久久精品久久久久久噜噜老黄| 国产 精品1| 亚洲精品美女久久av网站| 国产男女内射视频| 中文字幕免费在线视频6| 寂寞人妻少妇视频99o| 国产黄片视频在线免费观看| 欧美精品高潮呻吟av久久| 日韩精品有码人妻一区| 黄色视频在线播放观看不卡| 色94色欧美一区二区| 国国产精品蜜臀av免费| 在线亚洲精品国产二区图片欧美 | 免费av中文字幕在线| 亚洲精品国产av蜜桃| 天天躁夜夜躁狠狠久久av| 国产免费一区二区三区四区乱码| 国产午夜精品一二区理论片| 国产一区二区在线观看av| 51国产日韩欧美| 黑丝袜美女国产一区| 夫妻午夜视频| 99热这里只有精品一区| 欧美日韩一区二区视频在线观看视频在线| 搡老乐熟女国产| 在线观看人妻少妇| 成人手机av| 亚洲欧美清纯卡通| 人妻一区二区av| 又粗又硬又长又爽又黄的视频| 晚上一个人看的免费电影| 欧美xxⅹ黑人| 一区二区三区四区激情视频| 亚洲三级黄色毛片| 69精品国产乱码久久久| 国产精品熟女久久久久浪| 99re6热这里在线精品视频| 各种免费的搞黄视频| 久久精品国产亚洲网站| 国产av码专区亚洲av| 国产日韩一区二区三区精品不卡 | 一级毛片黄色毛片免费观看视频| 丰满迷人的少妇在线观看| 久久精品国产a三级三级三级| 亚洲精品乱码久久久久久按摩| 久久午夜福利片| 日韩成人伦理影院| 女人精品久久久久毛片| 亚洲精品成人av观看孕妇| 久久久国产精品麻豆| 国产高清国产精品国产三级| 日本91视频免费播放| 日韩伦理黄色片| 99热6这里只有精品| 免费观看av网站的网址| 久久青草综合色| 男人添女人高潮全过程视频| 欧美人与性动交α欧美精品济南到 | 中文字幕制服av| 久久久亚洲精品成人影院| 美女福利国产在线| 蜜桃在线观看..| 18禁裸乳无遮挡动漫免费视频| 热re99久久精品国产66热6| 免费高清在线观看日韩| 丰满迷人的少妇在线观看| 高清在线视频一区二区三区| 国产精品成人在线| 国产亚洲精品久久久com| 免费看光身美女| 中文字幕最新亚洲高清| 午夜福利视频在线观看免费| 午夜av观看不卡| 九色亚洲精品在线播放| 女性生殖器流出的白浆| 狂野欧美激情性xxxx在线观看| 久久久国产欧美日韩av| 哪个播放器可以免费观看大片| 简卡轻食公司| 不卡视频在线观看欧美| √禁漫天堂资源中文www| 亚洲图色成人| 国产精品无大码| 精品一区二区三区视频在线| 国产伦理片在线播放av一区| 黑人高潮一二区| av线在线观看网站| 黄片无遮挡物在线观看| 丝袜脚勾引网站| 插逼视频在线观看| 亚洲色图综合在线观看| 精品人妻在线不人妻| 亚洲精品日韩av片在线观看| av女优亚洲男人天堂| 国产免费现黄频在线看| 青青草视频在线视频观看| 欧美97在线视频| 国产成人91sexporn| 建设人人有责人人尽责人人享有的| 天堂中文最新版在线下载| 久久国产亚洲av麻豆专区| 色5月婷婷丁香| 国产精品免费大片| 精品亚洲乱码少妇综合久久| 日本色播在线视频| 狂野欧美激情性bbbbbb| 寂寞人妻少妇视频99o| 黄色毛片三级朝国网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 蜜臀久久99精品久久宅男| 日本欧美国产在线视频| 五月开心婷婷网| 国产精品女同一区二区软件| 一级爰片在线观看| 国产国语露脸激情在线看| 午夜激情福利司机影院| 天天操日日干夜夜撸| 日韩免费高清中文字幕av| 欧美激情国产日韩精品一区| 亚洲av电影在线观看一区二区三区| 99久久精品一区二区三区| 亚州av有码| 涩涩av久久男人的天堂| 精品久久久久久电影网| 五月玫瑰六月丁香| 18禁在线播放成人免费| 日韩欧美一区视频在线观看| 黑人高潮一二区| 欧美成人午夜免费资源| 哪个播放器可以免费观看大片| 国产精品一二三区在线看| 18在线观看网站| 97精品久久久久久久久久精品| .国产精品久久| 久久青草综合色| 99国产精品免费福利视频| 亚洲情色 制服丝袜| 热99久久久久精品小说推荐| 亚洲国产欧美在线一区| 国产色爽女视频免费观看| 亚洲国产av影院在线观看| 丝袜在线中文字幕| av播播在线观看一区| 国产精品一二三区在线看| 超色免费av| 欧美精品亚洲一区二区| 午夜福利在线观看免费完整高清在| 啦啦啦中文免费视频观看日本| 赤兔流量卡办理| 少妇的逼好多水| videos熟女内射| 久久久久久伊人网av| 熟妇人妻不卡中文字幕| 国产国语露脸激情在线看| 天堂俺去俺来也www色官网| 男女免费视频国产| 人妻人人澡人人爽人人| 久久国产精品大桥未久av| 国精品久久久久久国模美| 日日摸夜夜添夜夜添av毛片| 狂野欧美激情性xxxx在线观看| 国产又色又爽无遮挡免| av一本久久久久| 一级毛片aaaaaa免费看小| 高清av免费在线| 在线观看免费视频网站a站| 一级片'在线观看视频| 久久久久久久亚洲中文字幕| 美女国产视频在线观看| av线在线观看网站| 看非洲黑人一级黄片| 免费黄网站久久成人精品| 少妇丰满av| 中文字幕免费在线视频6| 美女主播在线视频| 久久人人爽人人片av| 亚洲激情五月婷婷啪啪| 午夜激情av网站| 国产永久视频网站| 97精品久久久久久久久久精品| 国产黄色免费在线视频| 水蜜桃什么品种好| 男女啪啪激烈高潮av片| 欧美精品亚洲一区二区| 天天躁夜夜躁狠狠久久av| av女优亚洲男人天堂| 熟女电影av网| 热99国产精品久久久久久7| 久久99一区二区三区| 国产老妇伦熟女老妇高清| 国产永久视频网站| 九色亚洲精品在线播放| 18在线观看网站| 青春草国产在线视频| 在线观看免费日韩欧美大片 | 狠狠婷婷综合久久久久久88av| 春色校园在线视频观看| 一级毛片黄色毛片免费观看视频| 久久国内精品自在自线图片| 大话2 男鬼变身卡| a 毛片基地| 视频在线观看一区二区三区| 午夜激情福利司机影院| 久久久午夜欧美精品| a级毛片黄视频| 日韩一本色道免费dvd| 秋霞在线观看毛片| 国产精品国产三级专区第一集| 97精品久久久久久久久久精品| 亚洲激情五月婷婷啪啪| 建设人人有责人人尽责人人享有的| av在线观看视频网站免费| 日本黄色片子视频| 欧美人与性动交α欧美精品济南到 | 精品人妻在线不人妻| 制服人妻中文乱码| 亚洲性久久影院| 欧美性感艳星| 人妻制服诱惑在线中文字幕| 免费观看性生交大片5| 色94色欧美一区二区| 99九九在线精品视频| 国产av精品麻豆| 婷婷色综合www| 亚洲精品aⅴ在线观看| 午夜激情av网站| 欧美日韩成人在线一区二区| 欧美亚洲日本最大视频资源| 少妇被粗大猛烈的视频| 久久久久久久国产电影| av免费观看日本| 丝袜在线中文字幕| 国产极品粉嫩免费观看在线 | 五月开心婷婷网| 99九九在线精品视频| 欧美最新免费一区二区三区| 日韩免费高清中文字幕av| 亚洲欧美精品自产自拍| 亚洲欧美一区二区三区国产| 国产精品偷伦视频观看了| 美女视频免费永久观看网站| 少妇人妻 视频| 在线观看人妻少妇| 一级a做视频免费观看| 亚洲精品一二三| 人人妻人人添人人爽欧美一区卜| www.av在线官网国产| 午夜91福利影院| 国产免费一级a男人的天堂| 下体分泌物呈黄色| av播播在线观看一区| 美女视频免费永久观看网站| 草草在线视频免费看| 欧美精品高潮呻吟av久久| 丰满乱子伦码专区| 国产探花极品一区二区| 午夜福利视频精品| 精品熟女少妇av免费看| 99热6这里只有精品| 婷婷色综合大香蕉| 国产在线一区二区三区精| 五月玫瑰六月丁香| 热99久久久久精品小说推荐| 国产综合精华液| 国产免费现黄频在线看| 国产精品嫩草影院av在线观看| 这个男人来自地球电影免费观看 | 成人综合一区亚洲| 成人影院久久| 国产探花极品一区二区| 久久 成人 亚洲| 免费观看av网站的网址| 少妇人妻精品综合一区二区| 日产精品乱码卡一卡2卡三| 女的被弄到高潮叫床怎么办| 少妇人妻久久综合中文| 日韩电影二区| 18禁在线播放成人免费| 色94色欧美一区二区| 中文乱码字字幕精品一区二区三区|