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

    面向滲漏探測的堤防磁場分布特性分析

    2024-01-01 00:00:00宋俊磊周丹肖國強周華敏董凱鋒晉芳莫文琴惠亞娟

    摘要:快速確定堤防滲漏隱患的位置信息,對堤防的治理具有重要意義。然而在滲漏通道弱磁數(shù)據(jù)處理過程中,堤防正常場計算往往需要花費較長時間,嚴重制約著滲漏隱患的定位效率。為解決上述問題,本文基于有限單元法開展了磁電阻率法探測堤防滲漏的建模仿真工作,探究了不同尺寸堤防內(nèi)穩(wěn)定電流場及其外部磁場的分布特性,并構(gòu)建了一定數(shù)量的模型數(shù)據(jù)庫。研究結(jié)果表明,磁感應(yīng)強度y分量在測域中部的整體波動相較于x分量和z分量更大,其值隨堤防上寬的增加呈下降趨勢,而隨斜坡角度的增大或堤防主體高度的增加呈上升趨勢。依據(jù)本文總結(jié)的分布特性,可以預(yù)測數(shù)據(jù)庫中缺失的堤防正常場,堤防模型的預(yù)測平均相對誤差為9.63%。

    關(guān)鍵詞:堤防滲漏;磁電阻率法;穩(wěn)定電流場分布;磁場分布

    doi:10.13278/j.cnki.jjuese.20230096

    中圖分類號:P631

    文獻標志碼:A

    Supported by the National Key Research and Development Program of China (2017YFC1502600)

    Analysis of Magnetic Field Distribution Characteristics of

    Embankment for Leakage Detection

    Song Junlei1, 2, 3, Zhou Dan1, 2, 3, Xiao Guoqiang4, Zhou Huamin4,

    Dong Kaifeng1, 2, 3, Jin Fang1, 2, 3, Mo Wenqin1, 2, 3, Hui Yajuan1, 2, 3

    1. School of Automation, China University of Geosciences, Wuhan 430074, China

    2. Hubei Key Laboratory of Advanced Control and Intelligent Automation for Complex Systems, Wuhan 430074, China

    3. Engineering Research Center of Intelligent Technology for Geo-Exploration, Ministry of Education, Wuhan 430074, China

    4. Changjiang River Scientific Research Institute of Changjiang Water Resources Commission, Wuhan 430019, China

    Abstract: It is of great significance to quickly determine the location information of seepage hazard for the treatment of embankments. However, in the process of magnetic weakening data processing of leakage channel, it often takes a long time to calculate the normal field of an embankment, which seriously restricts the locating efficiency of potential leakage hazards. In order to solve the above problems, based on the finite element method," the modeling and simulation of the detection of embankment leakage by magnetometric resistivity method

    was carried out, and" the distribution characteristics of the internal stable current field and the external magnetic field of the embankment with different sizes were explored, and" a certain number of model databases were built. The results show that the overall fluctuation value of" y component of magnetic induction intensity in the middle part of the measurement domain is larger than that of

    x component and z

    component, and its value decreases with the increase of the upper width of the embankment, and increases with the increase of the slope angle or the height of the embankment. According to the distribution characteristics summarized in this paper," the missing normal field of embankment in the database can be predicted, and the mean relative error of the embankment model is 9.63%.

    Key words: embankment leakage; magnetometric resistivity method;

    stable current field distribution; magnetic field distribution

    0 引言

    堤防是人民生活和社會生產(chǎn)的重要保障[1],在我國防洪體系中起著舉足輕重的作用。但大量堤防受建設(shè)時期社會經(jīng)濟和技術(shù)水平的限制,以及地質(zhì)災(zāi)害、環(huán)境因素和后期運管維護情況的綜合影響,存在裂縫、松動、滲漏和管涌等安全隱患[24]。

    磁電阻率法是一種測量由人工電流源在待測區(qū)所激發(fā)磁場的分布反演地質(zhì)結(jié)構(gòu)的物探技術(shù)[5]?;诖烹娮杪史?,本科研團隊與長江科學(xué)院開展合作,已研制了一套堤防滲漏通道探測裝備,包括軟、硬件系統(tǒng)和數(shù)據(jù)處理方法,并成功應(yīng)用于江西省余干縣鄱陽湖流域堤防和湖北省陽新縣長江干堤等堤防滲漏通道探測,實現(xiàn)了堤防滲漏通道和滲漏入水口位置和走向的探測[6]。

    使用磁測儀獲取的待測區(qū)域總磁場包含該區(qū)域的正常場與滲漏通道產(chǎn)生的磁場,因此需消除目標區(qū)域的正常場才能獲得滲漏通道產(chǎn)生的磁場分布特征。為得到待測區(qū)域的正常場,需對所測堤防進行建模和數(shù)據(jù)分析。堤防建模的主要目的是獲得堤防內(nèi)激勵電流的分布特性,依據(jù)電流分布特性得到待測區(qū)域正常場。該過程花費時間較長,在這期間堤防滲漏狀況可能會惡化,顯然無法滿足汛期堤防應(yīng)急搶險的需求。

    在防汛期間,排查效率是重中之重。若能掌握不同尺寸堤防內(nèi)穩(wěn)定電流場的分布特性,就能根據(jù)堤防的具體尺寸快速獲得堤防的正常場,大大減少每次依據(jù)具體堤防尺寸建模分析的時間,提升堤防滲漏的探測效率。因此,探究并獲得堤防內(nèi)部穩(wěn)定電流場及其外部磁場的分布特性是實現(xiàn)堤防快速建模的基礎(chǔ),為快速確定滲漏通道深度和位置信息提供數(shù)據(jù)支撐,對滿足堤防應(yīng)急搶險需要具有重要意義。

    目前關(guān)于堤防內(nèi)穩(wěn)定電流場分布規(guī)律的研究較少,且現(xiàn)有研究往往僅停留在對電流場的定性解釋上,沒有對其分布規(guī)律進行定量的分析與解釋。如:王朋[7]運用二維電阻率層析成像技術(shù)和不同電極裝置對土石壩滲漏進行探測,研究涉及了堤防二維電場的分布應(yīng)用,但未對二維電場分布的規(guī)律進行研究和總結(jié);張欣[8]通過構(gòu)建有限差分模型,采用MATLAB編程探究了含隱患土石堤壩的電場分布規(guī)律,但并未對規(guī)律進行理論分析與解釋;歐元超[9]通過試驗研究獲得了角度及埋深影響下滲漏通道的電場分布規(guī)律,但對影響因素的探討不足。綜上,目前采用電法手段進行堤防滲漏診斷后,診斷結(jié)論大多是通過人為判斷獲得,缺乏針對堤防三維電流場分布規(guī)律的理論支撐。因此,本文針對不同尺寸的堤防進行了有限元電流場仿真,根據(jù)仿真現(xiàn)象與仿真數(shù)據(jù)對堤防內(nèi)電流場的分布特性進行分析,并在電流分布特性的基礎(chǔ)上進一步開展磁場分布特性研究。

    1 穩(wěn)定電流場模擬方法對比

    電法勘探模擬方法可分為三類:物理模擬、解析計算和數(shù)值模擬。其中數(shù)值模擬方法的主流方法包括:有限差分法、積分方程法、邊界單元法和有限單元法等。每種模擬方法都有其優(yōu)缺點[10](表1)。數(shù)值模擬法能夠很好地解決復(fù)雜場源以及具有不規(guī)則輪廓地下構(gòu)造的模擬問題,并且徹底擺脫實物模型和解析公式,其中有限差分法和邊界單元法的適應(yīng)性不如有限單元法。計算機技術(shù)發(fā)展至今,有限單元法對計算機的要求已較易實現(xiàn),故該方法已經(jīng)成為現(xiàn)代地球物理研究的重要手段,應(yīng)用范圍很廣,本文亦采用此方法。

    2 堤防建模與內(nèi)部電流場仿真分析

    2.1 穩(wěn)定電流場基本理論

    假設(shè)在電阻率為ρ的均勻各向同性無限介質(zhì)中,有一位于地表的點電源A(圖1),其電流大小為I,則與A點距離為r的M點的電流密度值為

    J=I2πr2。(1)

    根據(jù)歐姆定律微分形式可知,M點的電場強度為

    E=Jρ=Iρ2πr2。(2)

    由式(2)可求出M點的電壓為

    U=∫SymboleB@rEdr=Iρ2πr。(3)

    由式(3)可知,在以r為半徑的半球面上電壓相等。

    當有兩個相距一定距離的電極置于地表時,兩個電極所形成的電流場都會影響近地表任何一點的電位。圖2為雙電極供電時電流線與等位線的分布圖。

    本文使用的磁電阻率法采用的是兩個異性點電源的供電方式,其電流場分布應(yīng)該服從上述理論,在此基礎(chǔ)上觀察地形對其電流場分布的影響。

    2.2 堤防內(nèi)電流場仿真計算

    本文作者所在科研團隊已在鄱陽湖流域堤防進行了堤防滲漏探測[6],本文基于此堤防的實際尺寸建立基礎(chǔ)仿真模型,并在此基礎(chǔ)上改變堤防的不同參數(shù),探究不同尺寸堤防內(nèi)穩(wěn)定電流場和磁場的分布特性。

    圖3為磁電阻率法測量示意圖,本文即是對此過程進行仿真分析。

    本模型采用的物理場類型為AC/DC(交流/直流),為了將電場和磁場結(jié)合起來選用了多物理場耦合。堤防基礎(chǔ)模型如圖4所示,其中長方體為大地,棱臺為堤防主體部分。大地的整體寬度(x方向)設(shè)置為180 m,長度(y方向)為100 m,高度(z方向)為20 m。同時考慮到實際情況中大地向無限遠處延伸的屬性,在模型周圍建立了無限元域表征大地?zé)o限遠的屬性。本文的研究內(nèi)容是探究棱臺參數(shù)的變化對堤防內(nèi)穩(wěn)定電流場與外部磁場分布的影響,即探究不同尺寸堤防內(nèi)穩(wěn)定電流場和磁場的分布特性。

    使用磁電阻率法進行測量時需要施加人工激勵電流,電極的位置與埋深會影響電流的分布[12],因此考慮實際測量需求將兩電極位置分別固定在點A(-60, 50, 0)與點B(60, 50, 0)處,其中A處為負電極,B處為正電極,電極施加的電流大小為0.5 A。堤防上寬為5 m,堤防斜坡的坡度為19°,坡長為24 m。令堤防寬度變化值為M,基礎(chǔ)高度為H,變化值為N;考慮堤防的材質(zhì)以黃土為主,故將堤防整體電阻率設(shè)置為200 Ω·m。

    完成幾何建模與物理場設(shè)置之后,需對模型進行網(wǎng)格剖分。網(wǎng)格剖分分為物理場控制網(wǎng)格和用戶控制網(wǎng)格兩種:物理場控制網(wǎng)格能自動完成網(wǎng)格剖分,比較方便,但是不一定能滿足用戶需求;而用戶控制網(wǎng)格可以根據(jù)需求自定義網(wǎng)格。自由四面體網(wǎng)格對幾何的適應(yīng)性最好,因此本文采用的網(wǎng)格類型是自由四面體網(wǎng)格。

    在劃分網(wǎng)格大小時,網(wǎng)格劃分越細致,計算結(jié)果越精確,但也意味著更長的求解時間和更多的內(nèi)存需求。本文將對重點研究對象堤防主體網(wǎng)格進行細化剖分,對堤防邊界無限元域網(wǎng)格進行粗化剖分。這樣可以在所關(guān)注的區(qū)域獲得更加精細化的結(jié)果,保證數(shù)據(jù)的可靠性。網(wǎng)格剖分結(jié)果如圖5所示。

    構(gòu)建好堤防模型后,對模型的相關(guān)電場參數(shù)進行求解,圖6為堤防模型電壓等值線結(jié)果圖。由圖6可知,在正電極處電壓最高,負電極處電壓最低,電壓等值線為以電極為中心的圓,符合式(3)。

    圖7是堤防模型電場模結(jié)果。由圖7a、b可知,兩個電極所在處電場模較大,越遠離電極,電場模越小。因為是均勻介質(zhì),故可以認為堤防非主體部分電場線分布符合圖2所示的雙電極供電時的電流分布。由圖7c可知,堤防主體部分對大地電流的分布產(chǎn)生了影響,導(dǎo)致部分電流向上分散,電場線向上彎曲。下文即針對堤防不同參數(shù)的變化,分析堤防主體分走的電流占比,初步探究其內(nèi)部電流分布特性。

    在模型中取一組對稱的截面(y=30, 70 m,圖8a)與截線(x=40, -40 m,z=-5 m,圖8b),探究電流分布的對稱性。圖8a中截點在x與z方向間距均為5 m,x方向布置范圍為[-100, 100] m,z方向覆蓋范圍為[-20, 10] m。由圖8a可知,兩個關(guān)于y=50 m截面對稱的截面上對應(yīng)測點的電流密度相差較小,在電極x坐標位置(x=-60, 60 m)處電流密度最大,越遠離電極電流密度越小,電流密度分布關(guān)于x=0對稱。圖8b可以進一步驗證該結(jié)論,兩根關(guān)于x=0截線對稱的截線上電流密度相等,且均在y=50 m截面上電流密度最大。

    2.3 堤防尺寸變化時電流分布特性

    以圖4模型為基礎(chǔ),通過改變堤防的上寬、斜坡角度和高度來建立一定數(shù)量的模型數(shù)據(jù)庫,研究堤防尺寸變化對電流分布的影響,并在電流場分布的基礎(chǔ)上研究堤防背景磁場的分布特性。

    基于基礎(chǔ)模型參數(shù),以2 m為跨度,建立上寬為3~15 m的堤防模型,其他主要參數(shù)保持不變。對各個模型中的電流密度進行仿真計算,求得堤防主體部分在某個具體截面中的占比關(guān)系,結(jié)果如圖9a所示。由圖9a可知:隨著堤防上寬的增加,堤防主體電流密度占比整體呈上升趨勢;截面越遠離中心面(y=50 m),相同條件下占比越大。

    基于基礎(chǔ)模型參數(shù),保持其他參數(shù)不變,建立斜坡角度為12°~25°的堤防模型,斜坡角度變化時堤防主體部分電流密度在截面中的占比關(guān)系結(jié)果如圖9b所示。由圖9b可知,隨著斜坡角度的增加,堤防主體電流密度占比整體呈下降趨勢,且截面越遠離中心面,相同條件下占比越大。

    基于基礎(chǔ)模型參數(shù),以0.5 m為跨度,建立高度為5.35~10.35 m的堤防模型,其他主要參數(shù)保持不變。圖9c為在不同堤防高度下,堤防主體電流密度在對應(yīng)截面中的占比變化。由圖9c可知,隨著堤防高度的增加,占比整體呈上升趨勢,且截面越遠離中心面,相同條件下占比越大。

    為了進一步了解堤防尺寸對電流分布的影響,本文探究了電流密度x、y、z三分量的分布特性。圖10以y=30 m、z=-5 m截線為例,繪制了電流密度三分量在該截線上的分布特性。由圖10可知:電流密度x分量是軸對稱的,其值在兩個電極向外達到了最大值,在靠近電極處下降,堤防中心處有極大值;電流密度y分量在正電極處(x=60 m)附近下降且為最小值,在負電極處(x=-60 m)附近上升且為最大值;電流密度z分量的分布特性與y分量一致,但其絕對值比y分量小。

    基于基礎(chǔ)模型改變堤防參數(shù),探究堤防尺寸變化對電流密度三分量的變化。圖11分別為堤防上寬、斜坡角度和高度變化時,在y=30 m、z=-5 m截線處電流密度x、y、z三分量的值。由圖11可知,在堤防尺寸變化的情況下,電流密度y分量曲線重合度最高,即受影響最??;電流密度z分量曲線波動稍大,受影響較大;電流密度x分量受影響程度處于兩者之間。不難發(fā)現(xiàn),在模型中段,即x∈[-20, 20] m區(qū)間內(nèi),三分量受堤防尺寸影響較小,曲線重合度較高且相對平滑;而在電極處曲線波動較大。

    3 堤防背景磁場仿真分析

    3.1 磁電阻率法基本理論

    在本實驗方案中,外部電流通過一對電極施加到地面,磁場為人工施加的電流通過地下導(dǎo)電介質(zhì)而激發(fā)。因為所施加的激勵電流為直流電,故頻域中一般麥克斯韋方程組的解可以簡化為穩(wěn)態(tài)問題,即:

    式中:E為電場強度;B為磁感應(yīng)強度,B=μH,H為磁場強度;μ為真空磁導(dǎo)率,μ=4π×10-7 T·m/A;J為電流密度。

    電流密度由與信號源相關(guān)的一次電流與待研究區(qū)域內(nèi)的傳導(dǎo)電流組成[13],即

    J=σE+Js。(7)

    式中:Js為信號源的電流密度;σ為電導(dǎo)率;σE為

    待測區(qū)域電流密度。E可以寫為[14]

    E=-

    式中,φ為電勢??梢缘贸鯷13]

    式中:r為某一點的位置矢量;r+A與r-B為激勵電極A和B的位置。

    Jx、Jy、Jz分別為電流密度x、y、z分量。

    為了求解磁場強度,需引入磁矢勢(A)[14],使得

    μH=SymbolQC@×A。(10)

    對磁矢勢施加經(jīng)典規(guī)范條件SymbolQC@·A=0,避免由磁矢勢定義引起的不確定性,即可將方程(5)寫成SymbolQC@A=μJ。又因為A(|r|→∞)=0,則可以得出測量點的通解[15]:

    A(r)=μ4π∫ΩJ(r′)r-r′dτ′(11)

    式中:r′為源鄰域中積分點的位置;dτ′為圍繞源點周圍的體積元。

    對于磁電阻率法來說,可以先由電流密度計算磁矢勢,然后通過計算磁矢勢的旋度得到磁感應(yīng)強度。從而畢奧-薩伐爾定律可以表示為[15]

    B(r)=

    μ4π∫[Js(r′)-σ(r′)SymbolQC@φ(r′)]r-r′r-r′3dτ′。(12)

    3.2 堤防尺寸變化時磁場分布特性

    在研究堤防內(nèi)電流分布特性的基礎(chǔ)上,開展了堤防正常場分布特性的研究工作。圖12為z=10 m截面堤防正常場分布圖。由圖12可知,在兩個電極處磁感應(yīng)強度最大,且電極附近磁場梯度較大,而在堤防中部磁場變化較小。

    圖13以y=30 m、z=10 m截線為例,繪制了磁感應(yīng)強度三分量在該截線上的分布特性。由圖13可知,磁感應(yīng)強度x分量在正電極處最小,在負電極處最大;磁感應(yīng)強度y分量整體上在兩個電極中部即x∈[-60, 60] m區(qū)間內(nèi)最大;磁感應(yīng)強度z分量較小,在x∈[-20, 20] m區(qū)間內(nèi)其曲線較為平滑,在其他區(qū)間內(nèi)波動較大。

    基于基礎(chǔ)模型改變堤防參數(shù),探究堤防尺寸變化對磁感應(yīng)強度三分量的變化。圖14分別為堤防上寬、斜坡角度和高度變化時,y=30 m、z=10 m截線處磁感應(yīng)強度x、y、z三分量的值。由圖14可知,在堤防尺寸變化的情況下,磁感應(yīng)強度三分量在電極附近變化較大,在曲線中段x∈[-20, 20] m區(qū)間內(nèi)受堤防尺寸影響較小,與電流規(guī)律吻合。即使曲線中段磁感應(yīng)強度受堤防尺寸影響較小,但是磁感應(yīng)強度y分量曲線中段整體波動相較于另外兩個分量更大。在曲線中段:隨著堤防上寬的增加,磁感應(yīng)強度y分量整體呈下降趨勢(圖14b);隨著斜坡角度的增大,整體上呈上升趨勢(圖14e);隨著堤防主體高度的增加,整體上亦呈上升趨勢(圖14h)。在曲線中段外的其他區(qū)域,磁感應(yīng)強度x分量相較于y分量變化波動更大,其變化規(guī)律在堤防上寬和斜坡角度的影響下與y分量一致,堤防高度增加時在兩個電極處與之相反。磁感應(yīng)強度z分量的變化規(guī)律不明晰,且其值相較于x和y分量更小,對堤防總磁感應(yīng)強度影響最小。

    在對數(shù)據(jù)庫中缺失的堤防模型進行預(yù)測時,可以依據(jù)此分布特性進行。

    3.3 堤防測點數(shù)據(jù)預(yù)測

    根據(jù)實地測量需求,在仿真模型中布置測點如圖15所示。測點布置范圍為100 m×60 m,測點之間在x與y方向上均間距10 m。電流密度測點布置在z=-5 m平面上(圖15a),磁感應(yīng)強度測點布置在z=10 m平面上(圖15b)。

    完成測點的布置后,分別對電流密度和磁感應(yīng)強度進行仿真研究,部分仿真結(jié)果如圖16所示。由圖16a可知:在曲線中段x∈[-20, 20] m區(qū)間內(nèi),電流密度隨x的增大變化梯度較小,在[-40, 20] m與[20, 40] m區(qū)間內(nèi)變化較大;同時,測點所在截面越靠近中心面(y=50 m),電流密度越大,越靠近兩端的測點,電流密度也越大。所得磁場分布趨勢(圖16b、c)與電場分布趨一致。

    通過建立一定數(shù)量的模型數(shù)據(jù)庫,可以獲得較多的測點數(shù)據(jù),以達到對數(shù)據(jù)庫中缺失的模型尺寸進行預(yù)測的目的。例如現(xiàn)需要對上寬為6 m、斜坡角度為17°、高度為8.35 m的堤防進行滲漏通道探測工作,使用磁電阻率法對該堤防開展作業(yè)前,首先需要獲取其正常場,即獲得如圖17所示測點處的磁感應(yīng)強度,進而可以依據(jù)現(xiàn)有數(shù)據(jù)庫以及堤防磁場分布特性對測點處的磁感應(yīng)強度進行預(yù)測,預(yù)測結(jié)果如圖18a所示。為了驗證預(yù)測結(jié)果的準確性,在仿真軟件中搭建了該尺寸的模型,仿真獲得該尺寸下堤防的實際磁感應(yīng)強度數(shù)據(jù),結(jié)果如圖18b所示。

    分別計算每條預(yù)測曲線與實際曲線對應(yīng)的均方誤差以及模型整體的均方誤差,其結(jié)果如表2所示。

    由圖18及表2可知,模型的均方誤差較小,從整體趨勢來看,越靠近電極所在y平面,預(yù)測數(shù)據(jù)的均方誤差越大。

    另外計算得知,此假設(shè)模型中,預(yù)測數(shù)據(jù)的平均相對誤差為9.63%,整體預(yù)測效果較好。

    4 結(jié)論

    本文針對滲漏通道弱磁數(shù)據(jù)處理過程中堤防正常場計算時間長的問題,開展了堤防磁場分布特性分析的研究工作,得到以下結(jié)論。

    1)在堤防中段,磁感應(yīng)強度y分量最大,在實際測量過程中,中段的測量區(qū)域覆蓋了堤防自身所在的地域,故為研究中更關(guān)心的區(qū)域。

    2)磁感應(yīng)強度y分量在測域中部的整體波動相較于另外兩個分量更大,其值隨著堤防上寬的增加呈下降趨勢,而隨著斜坡角度的增大或堤防主體高度的增加呈上升趨勢。在對數(shù)據(jù)庫中缺失的堤防模型進行預(yù)測時,可以依據(jù)此分布特性進行。

    3)依據(jù)本文總結(jié)的分布特性和現(xiàn)有的模型數(shù)據(jù)庫能夠較好地預(yù)測與之尺寸相近堤防的正常場,對遠離電極的平面預(yù)測準確度較高,堤防模型預(yù)測的平均相對誤差為9.63%。

    本研究可以在探測現(xiàn)場快速得到待測堤防正常場,大大縮減了堤防正常場精細化建模所需花費的時間。綜上,本研究工作對于提升堤防滲漏探測的效率具有一定的實際工程意義。

    參考文獻(References):

    [1] 蘇懷智,周仁練. 土石堤壩滲漏病險探測模式和方法研究進展[J]. 水利水電科技進展, 2022, 42(1): 110.

    Su Huaizhi, Zhou Renlian.Research Progress and Prospect of Earth-Rockfill Dam Leakage Detection Modes and Method[J]. Advances in Science and Technology of Water Resources, 2022, 42(1): 110.

    [2] 王少偉,蘇懷智,付啟民. 病險水利工程除險加固效果評價研究進展[J]. 水利水電科技進展, 2018, 38(6): 7785.

    Wang Shaowei, Su Huaizhi, Fu Qimin.Advances in Reinforcement Effect Evaluation for Dangerous Hydraulic Projects[J]. Advances in Science and Technology of Water Resources, 2018, 38(6): 7785.

    [3] 杜新強,何立瀅,任思睿,等. 中國北方地區(qū)水資源演變和供水構(gòu)成變化特征[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版), 2023, 53(2): 566577.

    Du Xinqiang, He Liying, Ren Sirei, et al. Water Resources Evolution and Water Supply Composition Changes in Northern China[J]. Journal of Jilin University (Earth Science Edition), 2023, 53(2): 566577.

    [4] 王福剛,盤惠林,李勝偉,等. 含膏巖地層地下水環(huán)境對混凝土建筑侵蝕實驗[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版), 2023, 53(4):11851194.

    Wang Fugang, Pan Huilin, Li Shengwei, et al. Experiment on Concrete Building Erosion by Underground Water in Gypsum Rock Formation[J]. Journal of Jilin University (Earth Science Edition), 2023, 53(4):11851194.

    [5] Abderrezak B, Bernard G, Michel A. Down-Hole Magnetometric Resistivity Inversion for Zinc and Lead Lenses Localization at Tobermalug, County Limerick, Ireland[J]. Journal of Applied Geophysics, 2017, 137: 2533.

    [6] 周華敏,鄔愛清,肖國強,等. 堤防管涌通道微弱磁場信號檢測技術(shù)[J]. 長江科學(xué)院院報, 2023, 40(8): 112119.

    Zhou Huamin, Wu Aiqing, Xiao Guoqiang, et al.Weak Magnetic Field Signal Detection Technology for Dike Piping Channel[J]. Journal of Changjiang River Scientific Research Institute, 2023, 40(8): 112119.

    [7] 王朋. 二維電阻率層析成像技術(shù)在土石壩滲漏診斷中的應(yīng)用[D]. 重慶:重慶交通大學(xué), 2009.

    Wang Peng. The Applied Research of Two-Dimensional Rsistivity Tomography Technology in the Earth Dam Leakage[D]. Chongqing: Chongqing Jiaotong University, 2009.

    [8] 張欣. 含隱患土石堤壩三維電場特征分析及成像識別技術(shù)研究[D]. 重慶:重慶交通大學(xué), 2017.

    Zhang Xin. Analysis of 3D Electric Field Characteristics and Imaging Recognition Technology of Earth-Rock Dam with Hidden Trouble[D]. Chongqing: Chongqing Jiaotong University, 2017.

    [9] 歐元超. 土石壩水庫滲漏并行電法探查工程應(yīng)用研究[D]. 淮南:安徽理工大學(xué), 2017.

    Ou Yuanchao. Study on the Application of Parallel Electrical Prospecting Method for Seepage of Earth Rock Dam[D]. Huainan:Anhui University of Science and Technology, 2017.

    [10] 宋昊翔. 起伏地形下高密度電法地形響應(yīng)研究:以松潘地區(qū)為例[D]. 北京:中國地質(zhì)大學(xué)(北京),2017.

    Song Haoxiang.A Research on Topographic Response of High-Density Resistivity Method in Rugged Area-Illustrated by the Case of Songpan Area[D]. Beijing: China University of Geosciences (Beijing), 2017.

    [11] 李鵬慧. 露天礦臺階內(nèi)電流場分布對高密度電法勘探的影響研究[D].包頭:內(nèi)蒙古科技大學(xué),2021.

    Li Penghui.Research on the Influence of Current Field Distribution in Strip Mine Steps on High-Density Electrical Exploration[D]. Baotou: Inner Mongolia University of Science and Technology, 2021.

    [12] 葉青,王曉,杜學(xué)彬,等. 中國地震井下地電阻率研究進展[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版), 2022, 52(3): 669683.

    Ye Qing, Wang Xiao, Du Xuebin, et al. Research Progress of Seismic Underground Geo-Resistivity in China[J]. Journal of Jilin University (Earth Science Edition), 2022, 52(3): 669683.

    [13] 萬玲,葉睿,汪滿滿,等. 基于磁電阻率法的堤壩滲漏路徑磁異常分布特征分析[J]. 中南大學(xué)學(xué)報(自然科學(xué)版),2022,53(8):30313039.

    Wan Ling, Ye Rui, Wang Manman, et al. Analysis of Magnetic Anomaly Distribution Characteristics Caused by Seepage Path in Dam Based on Magnetometric Resistivity[J]. Journal of Central South University (Science and Technology), 2022, 53(8): 30313039.

    [14] Chen J, Haber E, Oldenburg D W. Three-Dimensional Numerical Modelling and Inversion of Magnetometric Resistivity Data[J]. Geophysical Journal International, 2002, 149(3): 679697.

    [15] Jessop M, Jardani A, Revil A, et al. Magnetometric Resistivity: A New Approach and Its Application to the Detection of Preferential Flow Paths in Mine Waste Rock Dumps[J]. Geophysical Journal International, 2018, 215(1): 222239.

    国产亚洲精品综合一区在线观看 | 久久中文字幕一级| 久久精品91无色码中文字幕| 老司机福利观看| 啦啦啦观看免费观看视频高清 | 亚洲五月婷婷丁香| 国产精品亚洲av一区麻豆| 九色国产91popny在线| 黑人操中国人逼视频| 久久香蕉精品热| 中文字幕人成人乱码亚洲影| 久久中文字幕一级| 国产单亲对白刺激| 真人做人爱边吃奶动态| 少妇粗大呻吟视频| 久久久水蜜桃国产精品网| 色精品久久人妻99蜜桃| 日本三级黄在线观看| 又黄又粗又硬又大视频| 人人妻人人澡欧美一区二区 | 搡老熟女国产l中国老女人| 999精品在线视频| 国产免费av片在线观看野外av| 精品国产乱码久久久久久男人| 国产欧美日韩一区二区精品| 国内精品久久久久精免费| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利18| 琪琪午夜伦伦电影理论片6080| 嫁个100分男人电影在线观看| 精品午夜福利视频在线观看一区| 国产伦一二天堂av在线观看| 神马国产精品三级电影在线观看 | 波多野结衣高清无吗| 国产精品一区二区三区四区久久 | 99精品久久久久人妻精品| 午夜久久久久精精品| 在线十欧美十亚洲十日本专区| 久久人人97超碰香蕉20202| 国产在线精品亚洲第一网站| 亚洲欧美激情在线| 国产私拍福利视频在线观看| 男人的好看免费观看在线视频 | 黄色片一级片一级黄色片| 黑人欧美特级aaaaaa片| 妹子高潮喷水视频| 一二三四社区在线视频社区8| 一级,二级,三级黄色视频| 精品午夜福利视频在线观看一区| 97人妻天天添夜夜摸| 国产精品久久久av美女十八| 一夜夜www| av片东京热男人的天堂| 咕卡用的链子| 成人三级黄色视频| 99国产精品一区二区蜜桃av| 日本三级黄在线观看| 色老头精品视频在线观看| 男女午夜视频在线观看| 99国产综合亚洲精品| 一级作爱视频免费观看| 中文字幕高清在线视频| 国产男靠女视频免费网站| 琪琪午夜伦伦电影理论片6080| 69精品国产乱码久久久| 亚洲精品av麻豆狂野| 国产一区二区三区视频了| 91麻豆精品激情在线观看国产| 两个人视频免费观看高清| 成人国语在线视频| 色婷婷久久久亚洲欧美| 午夜福利一区二区在线看| 狠狠狠狠99中文字幕| 操美女的视频在线观看| 精品久久久精品久久久| 搡老岳熟女国产| 欧美日韩亚洲综合一区二区三区_| 99久久国产精品久久久| 国产亚洲精品综合一区在线观看 | 欧美精品亚洲一区二区| 亚洲第一青青草原| 午夜精品久久久久久毛片777| 久久午夜综合久久蜜桃| 亚洲中文字幕日韩| 十八禁网站免费在线| 给我免费播放毛片高清在线观看| 波多野结衣巨乳人妻| 在线观看66精品国产| 午夜视频精品福利| 亚洲成人久久性| 麻豆一二三区av精品| 波多野结衣巨乳人妻| 真人做人爱边吃奶动态| 亚洲av成人一区二区三| 女人爽到高潮嗷嗷叫在线视频| 99国产精品一区二区三区| 在线十欧美十亚洲十日本专区| 免费无遮挡裸体视频| 国产精品亚洲一级av第二区| 国产一级毛片七仙女欲春2 | av中文乱码字幕在线| 午夜福利18| 男女做爰动态图高潮gif福利片 | 国产单亲对白刺激| 欧美性长视频在线观看| 丁香欧美五月| 国产成人精品久久二区二区91| 成熟少妇高潮喷水视频| 黄网站色视频无遮挡免费观看| 欧美激情 高清一区二区三区| 这个男人来自地球电影免费观看| 久久人人爽av亚洲精品天堂| 亚洲专区中文字幕在线| 国产精品爽爽va在线观看网站 | 色av中文字幕| 少妇 在线观看| 18禁国产床啪视频网站| 免费在线观看黄色视频的| 国产精品野战在线观看| 亚洲欧美日韩无卡精品| 这个男人来自地球电影免费观看| 久久久久久久久免费视频了| 丝袜在线中文字幕| 天堂√8在线中文| 欧美日韩一级在线毛片| 午夜日韩欧美国产| 久久国产乱子伦精品免费另类| 国产成人精品久久二区二区91| 久久伊人香网站| 女同久久另类99精品国产91| 久久久久久久午夜电影| 久久久国产成人精品二区| 一级毛片女人18水好多| 曰老女人黄片| 动漫黄色视频在线观看| 日韩免费av在线播放| 亚洲国产看品久久| 中文字幕人妻丝袜一区二区| 免费观看人在逋| 久久精品91蜜桃| 亚洲第一电影网av| 亚洲电影在线观看av| 国产欧美日韩精品亚洲av| 午夜老司机福利片| 国产精品乱码一区二三区的特点 | 国产高清videossex| 韩国av一区二区三区四区| 国产伦一二天堂av在线观看| 亚洲精品国产区一区二| 久久人人97超碰香蕉20202| 欧美绝顶高潮抽搐喷水| 国产成人一区二区三区免费视频网站| 一区在线观看完整版| 黄色片一级片一级黄色片| 狂野欧美激情性xxxx| 欧洲精品卡2卡3卡4卡5卡区| 国产成人免费无遮挡视频| 精品福利观看| 成人av一区二区三区在线看| 久久国产乱子伦精品免费另类| 1024香蕉在线观看| 亚洲精品国产精品久久久不卡| 亚洲aⅴ乱码一区二区在线播放 | 精品国产亚洲在线| 国产又色又爽无遮挡免费看| 母亲3免费完整高清在线观看| 欧美中文日本在线观看视频| 国产精品亚洲一级av第二区| 欧美午夜高清在线| 真人一进一出gif抽搐免费| 制服丝袜大香蕉在线| 制服丝袜大香蕉在线| 国产成人影院久久av| 成人三级做爰电影| 日韩国内少妇激情av| 神马国产精品三级电影在线观看 | 精品久久久久久久人妻蜜臀av | 午夜福利,免费看| 男人的好看免费观看在线视频 | 国产一区二区激情短视频| 别揉我奶头~嗯~啊~动态视频| 首页视频小说图片口味搜索| 一本综合久久免费| 婷婷精品国产亚洲av在线| av有码第一页| 亚洲人成伊人成综合网2020| 757午夜福利合集在线观看| 国产av一区二区精品久久| 男女下面插进去视频免费观看| 大码成人一级视频| 啦啦啦观看免费观看视频高清 | av福利片在线| 老司机深夜福利视频在线观看| 国产成人系列免费观看| 欧美日本中文国产一区发布| 国产精品久久久久久人妻精品电影| 国产成人精品无人区| 色播亚洲综合网| 免费在线观看亚洲国产| 欧美日韩一级在线毛片| 女人被狂操c到高潮| 免费久久久久久久精品成人欧美视频| 亚洲精品av麻豆狂野| 免费在线观看影片大全网站| 淫秽高清视频在线观看| 久久精品国产99精品国产亚洲性色 | 美女高潮喷水抽搐中文字幕| 给我免费播放毛片高清在线观看| 亚洲中文日韩欧美视频| 一进一出抽搐gif免费好疼| 深夜精品福利| 午夜福利成人在线免费观看| 在线观看66精品国产| 丰满的人妻完整版| 精品熟女少妇八av免费久了| 欧美绝顶高潮抽搐喷水| 亚洲精品粉嫩美女一区| 精品国产一区二区久久| 国产精品亚洲美女久久久| 日本 av在线| 亚洲一码二码三码区别大吗| 免费一级毛片在线播放高清视频 | 一卡2卡三卡四卡精品乱码亚洲| 久久香蕉精品热| 国产一区二区在线av高清观看| 精品国产一区二区久久| 国产精品一区二区精品视频观看| 国产精品一区二区三区四区久久 | 老司机午夜福利在线观看视频| 久久国产精品影院| 欧美日韩亚洲综合一区二区三区_| 国产精品久久电影中文字幕| 悠悠久久av| 91av网站免费观看| 中文字幕高清在线视频| 99热只有精品国产| 久久久国产成人精品二区| 久久国产精品影院| 日本a在线网址| 日韩国内少妇激情av| 久久精品国产亚洲av香蕉五月| 两性夫妻黄色片| 久久精品亚洲精品国产色婷小说| 亚洲,欧美精品.| 黄色片一级片一级黄色片| 九色亚洲精品在线播放| 国产一区二区三区综合在线观看| 欧美成人免费av一区二区三区| 亚洲精品在线美女| 久久久久久久久久久久大奶| 国产亚洲精品av在线| 欧美av亚洲av综合av国产av| 精品欧美国产一区二区三| av电影中文网址| 亚洲片人在线观看| aaaaa片日本免费| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av日韩精品久久久久久密| 中出人妻视频一区二区| 午夜福利高清视频| 免费搜索国产男女视频| 中文字幕精品免费在线观看视频| 精品电影一区二区在线| АⅤ资源中文在线天堂| 麻豆一二三区av精品| 黄色 视频免费看| 色综合站精品国产| 亚洲一卡2卡3卡4卡5卡精品中文| 国内久久婷婷六月综合欲色啪| 日韩精品中文字幕看吧| 精品一品国产午夜福利视频| 天堂√8在线中文| 99国产精品一区二区蜜桃av| 国产成人av教育| 精品熟女少妇八av免费久了| 国产精品秋霞免费鲁丝片| 亚洲av美国av| 老鸭窝网址在线观看| 国产亚洲欧美精品永久| 激情在线观看视频在线高清| 亚洲第一av免费看| 久久人妻福利社区极品人妻图片| x7x7x7水蜜桃| 亚洲av成人一区二区三| 丁香欧美五月| 18禁美女被吸乳视频| 国产97色在线日韩免费| 12—13女人毛片做爰片一| 最近最新中文字幕大全免费视频| 亚洲 欧美一区二区三区| 亚洲国产看品久久| 亚洲免费av在线视频| 亚洲人成网站在线播放欧美日韩| 国内毛片毛片毛片毛片毛片| 免费观看精品视频网站| av在线播放免费不卡| av片东京热男人的天堂| 午夜视频精品福利| 一级,二级,三级黄色视频| 丝袜美足系列| av有码第一页| 免费观看精品视频网站| 欧美黄色片欧美黄色片| 婷婷丁香在线五月| 国产男靠女视频免费网站| 国产av又大| 一级a爱视频在线免费观看| 天天添夜夜摸| 亚洲精品av麻豆狂野| 男女下面进入的视频免费午夜 | 一级毛片高清免费大全| 村上凉子中文字幕在线| 久久人妻熟女aⅴ| 久99久视频精品免费| 国产av一区二区精品久久| 制服丝袜大香蕉在线| 国产国语露脸激情在线看| 波多野结衣巨乳人妻| 身体一侧抽搐| 国产午夜精品久久久久久| 99在线视频只有这里精品首页| 精品国内亚洲2022精品成人| 久久国产精品影院| 亚洲午夜精品一区,二区,三区| 好看av亚洲va欧美ⅴa在| 天天一区二区日本电影三级 | 久久香蕉精品热| 中文字幕人成人乱码亚洲影| 长腿黑丝高跟| 久久人妻熟女aⅴ| 两个人视频免费观看高清| 午夜久久久在线观看| 美女扒开内裤让男人捅视频| 天堂影院成人在线观看| 91老司机精品| 韩国av一区二区三区四区| 色综合婷婷激情| 精品高清国产在线一区| 又黄又粗又硬又大视频| 亚洲欧美激情在线| 最好的美女福利视频网| 色在线成人网| 免费在线观看完整版高清| 国产一区二区三区在线臀色熟女| 色精品久久人妻99蜜桃| 亚洲精品av麻豆狂野| 国产精品久久久久久亚洲av鲁大| 欧美日韩亚洲综合一区二区三区_| 精品午夜福利视频在线观看一区| 曰老女人黄片| 久久亚洲精品不卡| 69av精品久久久久久| 日本在线视频免费播放| 亚洲欧美激情在线| 变态另类丝袜制服| 亚洲欧美激情综合另类| 日韩精品中文字幕看吧| 久久中文看片网| 久久香蕉激情| 精品高清国产在线一区| 村上凉子中文字幕在线| 午夜福利在线观看吧| 久久久久久久久免费视频了| 亚洲专区中文字幕在线| 国产一区二区三区在线臀色熟女| 欧美丝袜亚洲另类 | 欧美日本中文国产一区发布| 欧美久久黑人一区二区| 国产私拍福利视频在线观看| 久久人妻av系列| 丁香欧美五月| 老司机午夜福利在线观看视频| 女人高潮潮喷娇喘18禁视频| 国产熟女午夜一区二区三区| 亚洲成人免费电影在线观看| 日韩高清综合在线| 国产精品一区二区精品视频观看| 亚洲国产精品成人综合色| 国产成+人综合+亚洲专区| 90打野战视频偷拍视频| 欧美日本视频| 国产精品日韩av在线免费观看 | 丝袜美腿诱惑在线| 亚洲午夜理论影院| 长腿黑丝高跟| 国产激情欧美一区二区| 亚洲av成人av| 亚洲人成电影观看| 久久国产亚洲av麻豆专区| 日本欧美视频一区| 日韩高清综合在线| 精品久久久久久久人妻蜜臀av | 男女之事视频高清在线观看| 法律面前人人平等表现在哪些方面| 日本黄色视频三级网站网址| 日韩三级视频一区二区三区| 久久精品国产99精品国产亚洲性色 | 久久午夜综合久久蜜桃| 亚洲国产高清在线一区二区三 | 男人舔女人的私密视频| 国产熟女午夜一区二区三区| 国产成人欧美| 亚洲五月婷婷丁香| 少妇粗大呻吟视频| 禁无遮挡网站| 亚洲人成网站在线播放欧美日韩| 中亚洲国语对白在线视频| 少妇 在线观看| av网站免费在线观看视频| 无限看片的www在线观看| 亚洲狠狠婷婷综合久久图片| 国产免费男女视频| 亚洲欧美一区二区三区黑人| 韩国av一区二区三区四区| 久久久精品欧美日韩精品| 亚洲av成人一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 日日摸夜夜添夜夜添小说| 国产av一区二区精品久久| 亚洲国产精品久久男人天堂| 1024香蕉在线观看| 国产高清videossex| 无遮挡黄片免费观看| 夜夜夜夜夜久久久久| 91成人精品电影| 久久这里只有精品19| 亚洲精品中文字幕在线视频| 操出白浆在线播放| 成人国产一区最新在线观看| 在线观看免费午夜福利视频| 午夜视频精品福利| 桃色一区二区三区在线观看| 久久久国产欧美日韩av| 狂野欧美激情性xxxx| 国产欧美日韩一区二区三| 国产欧美日韩一区二区精品| 中文亚洲av片在线观看爽| 久久天躁狠狠躁夜夜2o2o| 又紧又爽又黄一区二区| av网站免费在线观看视频| 欧美日韩瑟瑟在线播放| 国产成人精品久久二区二区91| 天天一区二区日本电影三级 | 免费观看精品视频网站| 国产一区二区在线av高清观看| 国内精品久久久久精免费| 村上凉子中文字幕在线| 给我免费播放毛片高清在线观看| 国产欧美日韩一区二区三区在线| 国产国语露脸激情在线看| 嫩草影视91久久| 99精品久久久久人妻精品| 国产一区二区三区视频了| 欧美日韩福利视频一区二区| 亚洲av成人不卡在线观看播放网| 国产精品美女特级片免费视频播放器 | 1024香蕉在线观看| 国产aⅴ精品一区二区三区波| 国产高清激情床上av| 亚洲男人的天堂狠狠| 国产高清激情床上av| 久久人人爽av亚洲精品天堂| 国产亚洲欧美在线一区二区| 久久人妻av系列| 欧美黑人欧美精品刺激| 9191精品国产免费久久| 男女之事视频高清在线观看| 亚洲av五月六月丁香网| 午夜福利视频1000在线观看 | 国产区一区二久久| 国产真人三级小视频在线观看| 欧美另类亚洲清纯唯美| 午夜免费观看网址| 国产在线精品亚洲第一网站| 国产午夜精品久久久久久| 精品国产国语对白av| 久久人妻熟女aⅴ| 国产精品国产高清国产av| 久久精品成人免费网站| 国产高清视频在线播放一区| 夜夜夜夜夜久久久久| 丁香欧美五月| 日韩欧美在线二视频| 亚洲av美国av| 国产1区2区3区精品| 欧美激情极品国产一区二区三区| 午夜福利影视在线免费观看| 免费女性裸体啪啪无遮挡网站| 欧美 亚洲 国产 日韩一| 欧美在线一区亚洲| 一级黄色大片毛片| 制服诱惑二区| 成人国产一区最新在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 他把我摸到了高潮在线观看| 精品久久久久久久久久免费视频| 在线天堂中文资源库| 无限看片的www在线观看| 成人三级黄色视频| 欧美人与性动交α欧美精品济南到| 久久欧美精品欧美久久欧美| 久久草成人影院| 日韩欧美一区二区三区在线观看| 免费在线观看日本一区| 91av网站免费观看| 九色亚洲精品在线播放| 亚洲成av人片免费观看| 黑人巨大精品欧美一区二区蜜桃| 人妻久久中文字幕网| 国产欧美日韩一区二区三| 国产高清videossex| 黄色视频不卡| 国产亚洲精品av在线| 高清黄色对白视频在线免费看| 亚洲专区国产一区二区| 国产成人一区二区三区免费视频网站| 国产又爽黄色视频| 国产精品一区二区免费欧美| 久久久国产精品麻豆| 亚洲伊人色综图| 国产区一区二久久| 丰满的人妻完整版| 亚洲在线自拍视频| 国产一区二区三区在线臀色熟女| 看黄色毛片网站| 国产伦人伦偷精品视频| 又紧又爽又黄一区二区| 视频在线观看一区二区三区| 亚洲精品中文字幕在线视频| 亚洲人成伊人成综合网2020| 一a级毛片在线观看| 精品国内亚洲2022精品成人| 国内毛片毛片毛片毛片毛片| 国语自产精品视频在线第100页| 久久久久九九精品影院| 女性生殖器流出的白浆| 日韩精品中文字幕看吧| 国产av又大| 亚洲国产毛片av蜜桃av| 精品久久久久久久久久免费视频| 欧美日韩乱码在线| 亚洲一区二区三区色噜噜| 少妇的丰满在线观看| 99精品欧美一区二区三区四区| 国产亚洲精品综合一区在线观看 | 久久中文字幕人妻熟女| 夜夜看夜夜爽夜夜摸| 性少妇av在线| 欧美激情 高清一区二区三区| 韩国av一区二区三区四区| 国内精品久久久久久久电影| 国产主播在线观看一区二区| 窝窝影院91人妻| 亚洲黑人精品在线| 嫩草影视91久久| 丰满人妻熟妇乱又伦精品不卡| 国产成人精品在线电影| 在线观看免费视频日本深夜| 两人在一起打扑克的视频| 午夜福利欧美成人| 国产99久久九九免费精品| 亚洲aⅴ乱码一区二区在线播放 | 精品久久久精品久久久| 日本 欧美在线| 亚洲伊人色综图| 精品久久久久久成人av| 免费高清视频大片| 一边摸一边抽搐一进一出视频| 两性夫妻黄色片| 久久国产精品男人的天堂亚洲| 亚洲 欧美一区二区三区| 高清黄色对白视频在线免费看| 午夜视频精品福利| 久久伊人香网站| 国产一区二区在线av高清观看| 777久久人妻少妇嫩草av网站| 久久久精品国产亚洲av高清涩受| 国产av一区二区精品久久| 一级毛片高清免费大全| www.www免费av| 久久欧美精品欧美久久欧美| 午夜久久久久精精品| 大陆偷拍与自拍| 国产精品99久久99久久久不卡| 亚洲人成伊人成综合网2020| 精品乱码久久久久久99久播| 咕卡用的链子| 后天国语完整版免费观看| 巨乳人妻的诱惑在线观看| 脱女人内裤的视频| 人成视频在线观看免费观看| 波多野结衣高清无吗| 欧美国产日韩亚洲一区| 俄罗斯特黄特色一大片| 久久婷婷成人综合色麻豆| 色av中文字幕| 亚洲av电影在线进入| 岛国在线观看网站| 一区在线观看完整版| 日韩精品免费视频一区二区三区| 丁香六月欧美| 中国美女看黄片| 在线视频色国产色| 精品国内亚洲2022精品成人| 久久人妻福利社区极品人妻图片| 亚洲精品在线美女| 美女扒开内裤让男人捅视频| 成人亚洲精品一区在线观看| 亚洲五月色婷婷综合| 成年女人毛片免费观看观看9| 国产一区二区三区综合在线观看| 手机成人av网站| 非洲黑人性xxxx精品又粗又长| 国产一卡二卡三卡精品| 熟妇人妻久久中文字幕3abv|