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

    基于區(qū)間變權(quán)和未知測度理論的煤層底板突水危險(xiǎn)性評價(jià)

    2023-10-13 06:14:48崔俊超郝建卿任君豪孟建勇王心義
    地下水 2023年5期
    關(guān)鍵詞:變權(quán)隔水層突水

    崔俊超,平 宇,郝建卿,任君豪,孟建勇,王心義

    (1.河南平煤神馬梁北二井煤業(yè)有限公司,河南 禹州 461670,2.河南理工大學(xué)資源環(huán)境學(xué)院,河南 焦作 454000;3.煤炭安全生產(chǎn)與清潔高效利用省部共建協(xié)同創(chuàng)新中心,河南 焦作 454000;4.晉能控股集團(tuán)麻家梁煤業(yè)有限責(zé)任公司,山西 朔州 038500)

    0 引言

    煤炭是我國能源結(jié)構(gòu)的主體,2021年我國能源生產(chǎn)結(jié)構(gòu)中煤炭占67.0%,煤炭消費(fèi)占一次能源消費(fèi)總量的比重為56.0%[1]。華北型煤田主采石炭-二疊紀(jì)煤層,奧陶系或寒武系厚層灰?guī)r含水層普遍發(fā)育,富水性較強(qiáng)且不均一,礦井開采過程中普遍受到底板灰?guī)r含水層水害威脅[2]。因此,客觀準(zhǔn)確地評價(jià)煤層底板突水危險(xiǎn)性,有利于礦井指定針對性的防治水措施,是保證礦井安全生產(chǎn)的基礎(chǔ)。

    近年來,國內(nèi)外相關(guān)學(xué)者針對煤層底板突水危險(xiǎn)性評價(jià)開展了大量研究,取得了一系列研究成果。Liu[3]根據(jù)呂溝煤礦的地質(zhì)條件、水文地質(zhì)特征、開采情況三個方面,確定12個指標(biāo)因子,基于層次分析法(AHP)及模糊綜合評價(jià)法,構(gòu)建了充填開采條件下的突水潰砂危險(xiǎn)性評價(jià)模型。李博[4]從巖溶發(fā)育程度、含水層巖性組合、地質(zhì)構(gòu)造、含水層水文地質(zhì)參數(shù)和地形地貌5個方面,確定了12項(xiàng)評價(jià)指標(biāo),構(gòu)建了較全面的西南地區(qū)煤層頂板巖溶含水層富水性評價(jià)指標(biāo)體系。并結(jié)合導(dǎo)水裂隙帶發(fā)育高度,確定研究區(qū)頂板突水危險(xiǎn)性等級。Li[5]利用主成分邏輯回歸分析(PCLRA)和改進(jìn)的層次分析法(IAHP)建立了綜合評價(jià)模型,以確定每個突水評價(jià)指標(biāo)的綜合權(quán)重,并利用GIS疊加法確定了突水風(fēng)險(xiǎn)區(qū)劃。結(jié)果表明此方法相較于突水系數(shù)法具有更高的擬合精度和分辨率。Ruan[6]在綜合分析礦山突水風(fēng)險(xiǎn)因素的基礎(chǔ)上,采用層次分析法對Dempster-Shafer(D-S)綜合規(guī)則進(jìn)行改進(jìn),建立了突水風(fēng)險(xiǎn)評價(jià)體系,并對王家?guī)X煤礦突水危險(xiǎn)性進(jìn)行了預(yù)測,預(yù)測結(jié)果符合實(shí)際情況。姚輝[7]以河北省華北型煤田東歡坨礦為研究對象,選取含水層性能、隔水層性能、地質(zhì)條件、煤層條件的評價(jià)因素集,綜合考慮10個評價(jià)因素,建立主客觀組合賦權(quán)-加權(quán)秩和比評價(jià)模型,利用GIS系統(tǒng)對研究區(qū)底板突水危險(xiǎn)性進(jìn)行了分區(qū)。施龍青[8]以現(xiàn)場實(shí)際數(shù)據(jù)為輸入樣本,通過灰狼優(yōu)化算法(GWO)得到Elman神經(jīng)網(wǎng)絡(luò)優(yōu)化的最佳權(quán)重和閾值,建立GWO-Elman神經(jīng)網(wǎng)絡(luò)底板突水預(yù)測模型,并在此基礎(chǔ)上通過測試樣本輸入模型驗(yàn)證,結(jié)果準(zhǔn)確率達(dá)到100%。

    上述研究將多種理論和技術(shù)引入煤層底板突水危險(xiǎn)性評價(jià)中,為水害預(yù)測技術(shù)的發(fā)展提供了支持。但這些評價(jià)預(yù)測方法也存在一些問題,如實(shí)際決策系統(tǒng)中指標(biāo)值存在“短板效應(yīng)”,即指標(biāo)值變化時(shí)其權(quán)重亦會發(fā)生變化;在評價(jià)體系建立時(shí),因勘探資料的不足,常忽略了煤層底板復(fù)合巖體的巖性組合特征和巖體質(zhì)量指標(biāo);構(gòu)造復(fù)雜程度界定時(shí)僅考慮構(gòu)造軌跡而未涉及構(gòu)造交叉、斷層落差和構(gòu)造性質(zhì)等要素。

    針對上述問題,本文基于收集和實(shí)測的地質(zhì)、水文地質(zhì)、工程地質(zhì)資料及鉆孔數(shù)據(jù),結(jié)合組合賦權(quán)法和變權(quán)理論確定指標(biāo)權(quán)重,基于未知測度理論建立煤層底板突水危險(xiǎn)性評價(jià)模型。評價(jià)結(jié)果相較常權(quán)-未知測度理論、區(qū)間變權(quán)-物元可拓兩種模型更加符合礦井水文地質(zhì)條件,并得到實(shí)際開采的良好驗(yàn)證。研究成果可以幫助礦井快速、準(zhǔn)確地識別底板突水危險(xiǎn)性,并采取針對性的防治水措施。

    1 理論與方法

    1.1 常權(quán)確定

    1.1.1 主觀權(quán)重

    層次分析法(AHP)是一種解決多目標(biāo)復(fù)雜問題的定性和定量相結(jié)合進(jìn)行計(jì)算主觀權(quán)重的研究方法[9]。采用5級標(biāo)度法構(gòu)建判斷矩陣,計(jì)算確定指標(biāo)主觀權(quán)重W1,并對判斷矩陣進(jìn)行一致性檢驗(yàn)[10]。定義一致性指標(biāo)CI:

    CI=(λmax-n)/(n-1)

    (1)

    式中:λmax為判斷矩陣的最大特征根;n為判斷矩陣的階數(shù)。

    CR=CI/RI

    (2)

    式中:RI為判斷矩陣的平均隨機(jī)一致性指標(biāo),RI取值根據(jù)查表確定。

    1.1.2 客觀權(quán)重

    灰色關(guān)聯(lián)分析法[11]通過計(jì)算比較序列與參考序列之間的關(guān)聯(lián)度來確定指標(biāo)因子的客觀權(quán)重值。首先按式(3)、式(4)對指標(biāo)進(jìn)行歸一化處理:

    越小越優(yōu)型指標(biāo):

    (3)

    越大越優(yōu)型指標(biāo):

    (4)

    設(shè)各指標(biāo)最大值為參考序列Z=(zj),鉆孔數(shù)據(jù)為比較序列X=(xij)。計(jì)算比較序列與參考序列之間的關(guān)聯(lián)度矩陣D=[dij],其中dij:

    (5)

    式中:Δij=xij-zj;ρ為分辨系數(shù),ρ∈[0,1],一般情況下,ρ=0.5。

    關(guān)聯(lián)度向量γ=(γj)由關(guān)聯(lián)度矩陣D每一列的均值組成,將γ歸一化處理后可得客觀權(quán)重向量W2。

    1.1.3 組合賦權(quán)

    層次分析法主要取決于決策者的主觀判斷,而灰色關(guān)聯(lián)分析法主要受指標(biāo)值的統(tǒng)計(jì)特征影響。根據(jù)最小相對信息熵原理,建立優(yōu)化組合賦權(quán)模型,利用拉格朗日乘子法[12]確定主、客觀權(quán)重的組合權(quán)重W0,可以一定程度上削弱兩種方法的弊端,計(jì)算公式為:

    (6)

    1.2 區(qū)間變權(quán)

    常權(quán)權(quán)重是現(xiàn)在多數(shù)評價(jià)模型中使用的權(quán)重。但在實(shí)際情況中,指標(biāo)因子的突變會改變其對于決策目標(biāo)的影響程度,即“短板效應(yīng)”。為客觀反映指標(biāo)因子異常對評價(jià)結(jié)果的影響程度,變權(quán)理論應(yīng)運(yùn)而生[13]。

    根據(jù)變權(quán)理論,按照指標(biāo)值大小將指標(biāo)區(qū)間劃分為“懲罰”、“不激勵不懲罰”、“初激勵”和“強(qiáng)激勵”區(qū)間,當(dāng)指標(biāo)量化值位于其中某一區(qū)間時(shí),其權(quán)重值會得到相應(yīng)的增大或減小。區(qū)間變權(quán)向量S=[sij]如下式[14]:

    (7)

    則變權(quán)權(quán)重wij:

    (8)

    1.3 未知測度理論

    未確知測度理論是一種不同于灰色信息、隨機(jī)信息和模糊信息的不確定性信息。測度函數(shù)的表達(dá)形式有很多種,主要包括直線型、指數(shù)函數(shù)型和二次函數(shù)型等,上述函數(shù)表達(dá)式都建立在“非負(fù)、歸一、可加”的基礎(chǔ)上[15]。本文采用應(yīng)用廣泛、計(jì)算簡單的直線型未確知測度函數(shù),計(jì)算樣本的單指標(biāo)測度評價(jià)矩陣和多指標(biāo)測度向量,最后依照置信度準(zhǔn)則進(jìn)行綜合評價(jià)。

    設(shè)某評價(jià)對象U有m個評價(jià)指標(biāo),分別用x1,x2,…,xm表示,則指標(biāo)空間記為X={x1,x2,…,xm},xj表示第j個評價(jià)指標(biāo)的測量值。對xj有k(1,2,…,p)個評價(jià)等級C1,C2,…,Cp,評價(jià)等級空間記為V={C1,C2,…,Cp}。設(shè)第k級比k+1級強(qiáng),即Ck>Ck+1,如果存在C1>C2>…>Cp,則稱{C1,C2,…,Cp}是評價(jià)空間V的一個有序分割類。

    1.3.1 單指標(biāo)測度函數(shù)

    單指標(biāo)未確知測度是整個評價(jià)指標(biāo)體系中不可再分的因素,若令μk=μ(x∈Ck)表示評價(jià)指標(biāo)x屬于第k個評價(jià)等級Ck的程度,且μ滿足非負(fù)性、歸一性、可加性,則稱μjk為單指標(biāo)的未確知測度,矩陣(μk)m×p為單指標(biāo)評價(jià)矩陣。

    各指標(biāo)的未確知測度值,需要通過具體的測度函數(shù)和該指標(biāo)的實(shí)際測量值計(jì)算得出,直線測度函數(shù)的表達(dá)式為:

    (9)

    (10)

    式中:x為某單指標(biāo)的實(shí)際測量值;μk(x)、μk+1(x)分別為該測量值x屬于Ck和Ck+1等級的測度;ak、ak+1分別為Ck、Ck+1等級上的中間值。特殊的,當(dāng)x≤a1或x≥ap時(shí),μk(x)=1。

    1.3.2 多指標(biāo)綜合測度向量

    令μk=μ(U∈Ck)表示評價(jià)對象U屬于第k個評價(jià)等級的程度,則有式(11)。

    μk=∑wj×μjk

    (11)

    式中:wj為評價(jià)指標(biāo)xj在評價(jià)指標(biāo)體系中的權(quán)重。

    μjk滿足0≤μjk≤1及Σμjk=1未確知測度條件,式(11)為評價(jià)對象U的多指標(biāo)綜合測度向量。

    1.3.3 置信度識別準(zhǔn)則

    為了給評價(jià)對象進(jìn)行最后的結(jié)果評價(jià),引入置信度識別準(zhǔn)則[16],可以克服最大隸屬度原則的缺點(diǎn):設(shè)λ為置信度(λ>0.5,常取λ=0.6或0.7),若C1>C2>…>Cp,且令:

    (12)

    則認(rèn)為評價(jià)對象U屬于第k0個評價(jià)等級Ck0。

    2 工程應(yīng)用

    2.1 研究區(qū)概況

    研究礦井為我國某典型華北型煤田礦井,目前主采二疊系山西組4號煤層,產(chǎn)量1 000萬t/a。井田內(nèi)地質(zhì)構(gòu)造以斷層為主,煤層下伏充水水源為奧陶系上馬家溝組巖溶水,水壓在3.32~6.03 MPa之間,平均水壓4.67 MPa。隔水層主要為砂-泥巖互層的復(fù)合巖體,厚度107.50~148.82 m,平均121.03 m。研究礦井的工作面布置、構(gòu)造展布、鉆孔分布如圖1,綜合地質(zhì)柱狀圖如圖2。

    圖1 礦井采掘平面示意圖

    圖2 礦井綜合地質(zhì)柱狀圖

    圖3 突水危險(xiǎn)性評價(jià)體系

    2.2 評價(jià)體系建立

    為綜合考慮含水層水文地質(zhì)特征、隔水層工程地質(zhì)特征及軟弱結(jié)構(gòu)面發(fā)育特征,本文選取水壓、單位涌水量、隔水層厚度、隔水層巖體質(zhì)量、層理面數(shù)量、泥巖占比、構(gòu)造復(fù)雜程度量化值7個因子作為評價(jià)預(yù)測底板突水危險(xiǎn)性的主控指標(biāo)。

    2.2.1 含水層水壓

    煤層底板含水層水壓是底板突水的動力來源。根據(jù)研究區(qū)水文地質(zhì)鉆孔資料及奧灰含水層長期觀測孔資料,奧陶系灰?guī)r含水層水壓呈現(xiàn)西南部較高,東北部較低的分布特征,如圖4a。

    圖4 指標(biāo)因子量化值分布趨勢

    2.2.2 單位涌水量

    單位涌水量反映含水層的富水性,單位涌水量越大,含水層富水性越強(qiáng),補(bǔ)給條件越好。根據(jù)礦井歷次抽水試驗(yàn)資料,繪制的奧灰含水層單位涌水量等值線圖如圖4b。

    2.2.3 有效隔水層厚度

    隔水層厚度指煤層底板與含水層頂板之間的距離,是反映隔水層隔水性能的關(guān)鍵指標(biāo)。開采活動會對底板巖體造成擾動破壞,有效隔水層厚度是煤層底板隔水層厚度減去底板擾動破壞深度的值??紤]麻家梁煤礦4號煤層傾角較緩,為5°~30°,工作面斜長一般小于200 m,且采高為5 m左右,故選取底板擾動破壞深度經(jīng)驗(yàn)公式[17]為:

    h=0.700 7+0.107 9L

    (13)

    式中:h為底板擾動破壞深度,m;L為工作面斜長,m。

    經(jīng)計(jì)算,研究區(qū)范圍有效隔水層厚度分布如圖4c。

    2.2.4 泥巖占比

    泥巖是具有良好隔水性能的塑性巖石,隔水層巖體中的泥巖占比一定程度上影響著隔水層的隔水性能[18]。根據(jù)井田地質(zhì)資料及鉆探成果,研究礦井4號煤隔水層主要由砂巖和泥巖組成,其泥巖占比分布趨勢如圖4d。

    2.2.5 巖體質(zhì)量

    按照《煤礦床水文地質(zhì)、工程地質(zhì)及環(huán)境地質(zhì)勘察評價(jià)標(biāo)準(zhǔn)》[19]規(guī)定,隔水層巖體質(zhì)量常用指標(biāo)M法來表征。M值越大,巖體越不容易發(fā)生塑性破壞,隔水性能越強(qiáng)。M值計(jì)算式為:

    (14)

    式中:Rc指巖塊飽和軸向抗壓強(qiáng)度,MPa/cm2;RQD指巖石質(zhì)量指標(biāo),%。

    RQD是指鉆探過程中,每次進(jìn)尺中等于或大于10 cm的柱狀巖芯累計(jì)長度與每個鉆進(jìn)回次進(jìn)尺之比,即:

    (15)

    式中:n指大于等于10 cm的巖芯數(shù)量;li指大于等于10 cm的巖芯長度,m;L指鉆孔或掃描線長度,m。

    研究礦井共20個鉆孔施工工程中進(jìn)行了取芯及巖樣物理力學(xué)參數(shù)測定,據(jù)此繪制的M值等值線如圖4e。

    2.2.6 層理面數(shù)量

    層理面是隔水層巖體中的原生軟弱結(jié)構(gòu)面,已有研究[20]顯示,在相同厚度相同傾角的砂泥巖復(fù)合巖體中,層理面越多,層厚越小,巖體越易發(fā)生形變與破壞。根據(jù)鉆孔資料,隔水層層理面數(shù)量等值線圖如圖4f。

    2.2.7 構(gòu)造復(fù)雜程度

    地質(zhì)構(gòu)造在煤層底板突水中起著關(guān)鍵作用,巖體中斷層或褶皺越復(fù)雜,其隔水能力越差,發(fā)生底板突水的可能性越大。以往構(gòu)造復(fù)雜程度常根據(jù)構(gòu)造軌跡,利用分形維數(shù)[21]來定量評價(jià)。但由于分形維數(shù)計(jì)算時(shí)未考慮構(gòu)造性質(zhì)、交叉復(fù)合特性、斷層落差等要素,致使所界定的構(gòu)造復(fù)雜程度難以精準(zhǔn)展現(xiàn)構(gòu)造的影響程度。

    這里將研究區(qū)劃分為87個1 000 m×1 000 m的單元(圖5),融合單元中的構(gòu)造交點(diǎn)數(shù)F1、斷層平均落差F2、分形維數(shù)F3、導(dǎo)水性能F4,來評估構(gòu)造復(fù)雜程度F。

    圖5 研究區(qū)單元劃分及編號

    不同性質(zhì)構(gòu)造的導(dǎo)水性能F4可利用模糊語氣算子理論[22]進(jìn)行定義。逆斷層因受水平方向最大主應(yīng)力的擠壓作用[23],其破碎帶充填程度及泥質(zhì)含量往往高于正斷層,導(dǎo)水性能一般較正斷層更弱,且區(qū)內(nèi)褶皺均為兩翼傾角2°~9°的寬緩褶皺,軸部巖體破壞程度更弱。因此確定正斷層、逆斷層、褶皺的相對導(dǎo)水性能語氣算子為{‘Same’,‘Slightly’,‘More obviously’},相對導(dǎo)水性能量化值為K=(1.0,0.667,0.429),F4即為塊段內(nèi)各構(gòu)造導(dǎo)水性能量化值之和。

    利用層次分析法[9]確定F1、F2、F3、F4的權(quán)重向量為Wf=(0.237,0.099,0.446,0.218),線性加權(quán)后得到塊段的構(gòu)造復(fù)雜程度F其分布趨勢如圖4g。

    2.3 權(quán)重確定

    結(jié)合專家意見及研究區(qū)地質(zhì)特征,確定比較矩陣R:

    經(jīng)過主觀、客觀權(quán)重耦合計(jì)算,得到指標(biāo)常權(quán)權(quán)重如表1。

    表1 常權(quán)計(jì)算結(jié)果 %

    利用SPSS軟件的K-均值聚類法[24]將各指標(biāo)值分為4類,以相鄰類別臨界值之間的均值作為變權(quán)區(qū)間閾值,即“懲罰”、“不激勵不懲罰”、“初激勵”、“強(qiáng)激勵”區(qū)間,如表2。

    表2 變權(quán)區(qū)間閾值

    鉆孔S1的指標(biāo)向量為Xs1=(0.664 2,0.040 4,0.554 9,0.442 2,0.363 0,0.128 2,0.086 1),確定其理想變權(quán)權(quán)重為W=(0.231 9,0.055 3,0.169 2,0.093 4,0.072 4,0.220 2,0.157 6)。將Xs1、W0及W代入式(7)和式(8),利用MATLAB軟件計(jì)算得調(diào)權(quán)參數(shù)c=0.502 5,a1=0.083 7,a2=0.841 5,a3=0.813 6。再將各鉆孔指標(biāo)向量進(jìn)行循環(huán)運(yùn)算,得到各鉆孔變權(quán)權(quán)重,如表3。

    表3 部分鉆孔變權(quán)權(quán)重

    2.4 模型構(gòu)建與評價(jià)

    針對研究礦井水文地質(zhì)特征,將其突水危險(xiǎn)性分為4個等級,級別越高,突水危險(xiǎn)性越強(qiáng),即Y={Ⅰ級,Ⅱ級,Ⅲ級,Ⅳ級}={弱,較弱,較強(qiáng),強(qiáng)},各指標(biāo)的分級區(qū)間閾值沿用表2分析結(jié)果,即懲罰區(qū)間對應(yīng)Ⅰ級,不激勵不懲罰區(qū)間對應(yīng)Ⅱ級,初激勵區(qū)間對應(yīng)Ⅲ級,強(qiáng)激勵區(qū)間對應(yīng)Ⅳ級。

    根據(jù)未知測度理論,首先確定等級中點(diǎn)值矩陣A=[ajk]:

    根據(jù)單指標(biāo)測度基本理論中的式(9)、式(10)和矩陣A構(gòu)建單指標(biāo)測度函數(shù),如圖6。

    圖6 單指標(biāo)測度函數(shù)

    以鉆孔S1為例,將指標(biāo)向量Xs1帶入到單指標(biāo)測度函數(shù)中,得到其單指標(biāo)測度矩陣為:

    則據(jù)式(11),鉆孔S1的多指標(biāo)綜合測度向量為μk=(0.464 9,0.215 1,0.290 7,0.029 4)。根據(jù)置信度識別原則,置信度λ取0.6,對鉆孔S1有C1+C2=0.679 9>0.6,k0=2,即鉆孔S1處的底板突水危險(xiǎn)性為Ⅱ級,危險(xiǎn)性較弱。同理計(jì)算可得其它鉆孔處的多指標(biāo)綜合測度向量和突水危險(xiǎn)性評價(jià)結(jié)果。

    3 結(jié)果對比與分析

    根據(jù)評價(jià)結(jié)果繪制研究區(qū)煤層底板突水危險(xiǎn)性分區(qū)圖,如圖7。同時(shí)繪制常權(quán)-未知測度理論、變權(quán)-物元可拓法[25]的評價(jià)分區(qū)結(jié)果,如圖8、圖9。

    圖7 區(qū)間變權(quán)-未知測度評價(jià)分區(qū)圖

    圖8 常權(quán)-未知測度評價(jià)分區(qū)圖

    圖9 區(qū)間變權(quán)-物元可拓評價(jià)分區(qū)圖

    結(jié)合鉆探資料分析:一采區(qū)14101-1工作面西側(cè)的63141鉆孔、63218鉆孔處有效隔水層厚度分別為84.52 m、86.16 m,遠(yuǎn)低于礦井均值97.35 m且接近最小值82.63 m;兩鉆孔隔水層層理面數(shù)量分別為79、80,遠(yuǎn)大于礦區(qū)均值50.54;加之此處斷層密集發(fā)育,構(gòu)造復(fù)雜程度較高,隔水層阻隔水能力遭到嚴(yán)重削弱,更易發(fā)生底板突水。礦井東南角水壓整體偏大,而有效隔水層厚度整體偏小,如B4005鉆孔處水壓達(dá)5.95 MPa(礦區(qū)最大值6.0 MPa),有效隔水層厚度為次小值83.67 m,致使此處較其它區(qū)域具有更強(qiáng)的突水危險(xiǎn)性。

    實(shí)際開采過程中,2020年4月開始回采的一采區(qū)14104-1工作面,出現(xiàn)過煤層底板奧陶系灰?guī)r含水層涌水量突然增大現(xiàn)象。其主要原因是工作面西側(cè)發(fā)育有DF9、DF8、F4等多個斷層,其中DF9斷層落差達(dá)40 m,導(dǎo)通了下伏奧灰含水層。斷層的集中發(fā)育削弱了底板巖體的穩(wěn)定性及隔水性能,隔水層在下伏含水層水壓的作用下發(fā)生破壞,形成了導(dǎo)水通道。

    由圖中可以看出,常權(quán)-未知測度理論模型在一采區(qū)14101-1工作面附近構(gòu)造密集發(fā)育、隔水層巖體完整性較差的情況下,不能夠調(diào)整隔水層厚度、層理面數(shù)量和構(gòu)造復(fù)雜程度的指標(biāo)權(quán)重,造成對14101-1工作面附近的突水危險(xiǎn)性造成了誤判;而變權(quán)-物元可拓模型中,指標(biāo)因子對各隸屬區(qū)間的隸屬關(guān)系為“非此即彼”的關(guān)系,即相較于未知測度理論不能體現(xiàn)指標(biāo)因子的模糊性與不確定性,因此造成評價(jià)結(jié)果整體偏高。區(qū)間變權(quán)-未知測度模型的評價(jià)結(jié)果與實(shí)際開采情況基本相符,且精準(zhǔn)判別了14104-1工作面的底板突水危險(xiǎn)性,評價(jià)體系的構(gòu)建、指標(biāo)權(quán)重的確定及評價(jià)結(jié)果更加合理、準(zhǔn)確。

    無論是開采實(shí)際還是模型評價(jià)結(jié)果,都指示出研究區(qū)內(nèi)煤層底板突水危險(xiǎn)性主要受到構(gòu)造的影響。斷層、褶皺等軟弱結(jié)構(gòu)面的發(fā)育會導(dǎo)致隔水層巖體中發(fā)育大量裂隙,誘發(fā)次生破壞,直接或間接地提供煤層底板突水導(dǎo)水通道,在礦井生產(chǎn)過程中應(yīng)特別注意斷層的探查工作,尤其在斷層集中發(fā)育的區(qū)域。

    區(qū)間變權(quán)-未知測度理論模型的評價(jià)結(jié)果中,突水危險(xiǎn)性為弱的區(qū)域面積最小,占比為2.26%,主要集中在礦井中部的3406、J1、J2孔附近,以及礦井東北角小部分區(qū)域。突水危險(xiǎn)性較弱的區(qū)域面積最大,占比54.81%,主要分布于一采區(qū)北部、三采區(qū)、四采區(qū)東部、五采區(qū)西部、六采區(qū)、七采區(qū)東部。突水危險(xiǎn)性較強(qiáng)的區(qū)域占比為38.05%,主要分布于一采區(qū)南部,二采區(qū)、四采區(qū)西部、五采區(qū)中部和七采區(qū)西部。突水危險(xiǎn)性為強(qiáng)的區(qū)域占比為4.88%,主要集中在一采區(qū)14101-1工作面西南側(cè)和礦井東南角。

    4 結(jié)語

    (1)結(jié)合層次分析法、灰色關(guān)聯(lián)分析及最組合賦權(quán)法,確定評價(jià)指標(biāo)組合常權(quán)W0=(0.194 2 0.057 5,0.179 6,0.099 2,0.076 5,0.228 4,0.164 7),基于區(qū)間變權(quán)理論確定各指標(biāo)因子的變權(quán)權(quán)重。

    (2)將4號煤底板突水危險(xiǎn)性分為弱,較弱,較強(qiáng),強(qiáng)四個等級,構(gòu)建區(qū)間變權(quán)-未知測度綜合評價(jià)模型,對研究區(qū)內(nèi)各鉆孔處的底板突水危險(xiǎn)性進(jìn)行定量分析與評價(jià),并繪制分區(qū)圖。

    (3)與常權(quán)-未知測度理論、區(qū)間變權(quán)-物元可拓兩種模型的評價(jià)結(jié)果對比顯示,區(qū)間變權(quán)-未知測度理論的權(quán)重確定、隸屬度計(jì)算更加合理,評價(jià)結(jié)果更加符合礦井地質(zhì)條件及開采實(shí)際。

    (4)評價(jià)結(jié)果顯示,研究區(qū)內(nèi)底板突水危險(xiǎn)性為弱、較弱、較強(qiáng)、強(qiáng)的區(qū)域面積占比分別為2.26%,54.81%,38.05%,4.88%,底板突水危險(xiǎn)性總體上較弱,危險(xiǎn)性強(qiáng)的區(qū)域集中一采區(qū)14101-1工作面西南側(cè)和礦井東南角,礦井生產(chǎn)過程中應(yīng)特別注意斷層的探查工作,尤其在斷層集中發(fā)育的區(qū)域。

    猜你喜歡
    變權(quán)隔水層突水
    礦井突水水源的判別方法
    礦井滯后突水機(jī)理探討
    變權(quán)空間權(quán)重構(gòu)造及空間效應(yīng)分析
    新疆生產(chǎn)建設(shè)兵團(tuán)城鎮(zhèn)化水平的變權(quán)組合預(yù)測
    李雅莊礦滯后突水主控因素分析
    西藏阿里結(jié)則茶卡湖西隔水層的賦存狀態(tài)及滲透性研究
    西部資源(2016年4期)2016-10-18 07:11:00
    基于變權(quán)的物流資源公平分配方法
    基于黃金分割法優(yōu)選的中長期負(fù)荷變權(quán)組合預(yù)測
    電測與儀表(2016年3期)2016-04-12 00:27:46
    小屯煤礦上組煤水文地質(zhì)特征及類型劃分
    科技資訊(2015年5期)2016-01-14 17:44:01
    大相嶺隧道高壓突水機(jī)理與預(yù)測分析
    午夜日本视频在线| 91在线精品国自产拍蜜月| 王馨瑶露胸无遮挡在线观看| 国产av国产精品国产| 一级二级三级毛片免费看| 欧美精品一区二区大全| 女人久久www免费人成看片| 日韩视频在线欧美| 国精品久久久久久国模美| 国产一区亚洲一区在线观看| 国内揄拍国产精品人妻在线| 最新中文字幕久久久久| 久久99精品国语久久久| 少妇丰满av| 精品国产乱码久久久久久小说| 亚洲天堂av无毛| 午夜av观看不卡| 哪个播放器可以免费观看大片| 91aial.com中文字幕在线观看| 视频区图区小说| 亚洲美女视频黄频| 人人妻人人澡人人看| 大片电影免费在线观看免费| 亚洲欧洲日产国产| 毛片一级片免费看久久久久| 极品少妇高潮喷水抽搐| 三级国产精品欧美在线观看| 久久久国产一区二区| 免费观看的影片在线观看| 又爽又黄a免费视频| 亚洲中文av在线| 亚洲熟女精品中文字幕| 国产成人免费无遮挡视频| 亚洲人与动物交配视频| 国产成人一区二区在线| 如日韩欧美国产精品一区二区三区 | 国产一区有黄有色的免费视频| 国产白丝娇喘喷水9色精品| 哪个播放器可以免费观看大片| 成人特级av手机在线观看| 曰老女人黄片| 国产精品国产av在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲内射少妇av| 热99国产精品久久久久久7| 中文在线观看免费www的网站| 久久久国产精品麻豆| 亚洲真实伦在线观看| 青春草亚洲视频在线观看| 一级黄片播放器| 一级毛片aaaaaa免费看小| 久久99一区二区三区| 97超视频在线观看视频| 99久久精品国产国产毛片| 久久久a久久爽久久v久久| av.在线天堂| 最近最新中文字幕免费大全7| videossex国产| av线在线观看网站| 国产精品久久久久久久电影| 在线观看三级黄色| 亚洲欧美一区二区三区国产| 人妻一区二区av| 国产精品成人在线| 亚洲色图综合在线观看| 十八禁高潮呻吟视频 | 乱系列少妇在线播放| a级一级毛片免费在线观看| 久久热精品热| 综合色丁香网| 高清毛片免费看| 国产精品.久久久| 亚洲欧美成人综合另类久久久| 啦啦啦视频在线资源免费观看| 欧美日韩亚洲高清精品| 精品国产一区二区三区久久久樱花| 亚洲怡红院男人天堂| 卡戴珊不雅视频在线播放| 自线自在国产av| 亚洲成色77777| www.色视频.com| 性色avwww在线观看| 91成人精品电影| 寂寞人妻少妇视频99o| 777米奇影视久久| 精品一区二区三区视频在线| 中国美白少妇内射xxxbb| 中文欧美无线码| 精品少妇久久久久久888优播| 一级av片app| 免费大片黄手机在线观看| 国产一区二区在线观看日韩| 欧美 亚洲 国产 日韩一| 久久国产精品男人的天堂亚洲 | 黄色一级大片看看| 日韩三级伦理在线观看| 麻豆精品久久久久久蜜桃| 色5月婷婷丁香| 国产一区二区三区av在线| 久久青草综合色| 三上悠亚av全集在线观看 | 熟女人妻精品中文字幕| 欧美激情国产日韩精品一区| 国内揄拍国产精品人妻在线| 国产片特级美女逼逼视频| 赤兔流量卡办理| 只有这里有精品99| 视频中文字幕在线观看| 亚洲国产精品国产精品| 91久久精品电影网| 久久久久久久精品精品| av在线app专区| 天堂中文最新版在线下载| 国产亚洲午夜精品一区二区久久| 午夜激情福利司机影院| 卡戴珊不雅视频在线播放| 99热这里只有精品一区| 国产免费福利视频在线观看| 18禁在线播放成人免费| 国产深夜福利视频在线观看| 久久热精品热| 卡戴珊不雅视频在线播放| 日本91视频免费播放| 国产免费一级a男人的天堂| 久久久国产一区二区| 99热全是精品| 色网站视频免费| 国产成人精品婷婷| 日本黄色片子视频| 大陆偷拍与自拍| 丰满人妻一区二区三区视频av| 亚洲国产色片| 大又大粗又爽又黄少妇毛片口| 欧美激情极品国产一区二区三区 | 国产探花极品一区二区| 亚洲一级一片aⅴ在线观看| 欧美日韩av久久| 欧美老熟妇乱子伦牲交| www.色视频.com| 麻豆成人午夜福利视频| 亚洲欧美日韩东京热| 亚洲国产精品专区欧美| 卡戴珊不雅视频在线播放| 美女视频免费永久观看网站| 久久久久久久大尺度免费视频| 国产精品久久久久久精品古装| 国产在线一区二区三区精| 午夜精品国产一区二区电影| 国产淫语在线视频| 色94色欧美一区二区| 国产真实伦视频高清在线观看| 多毛熟女@视频| 国产精品熟女久久久久浪| av在线app专区| 色94色欧美一区二区| 69精品国产乱码久久久| 日韩强制内射视频| 最近2019中文字幕mv第一页| 欧美97在线视频| 久久人妻熟女aⅴ| 99久久精品热视频| 亚洲精品久久午夜乱码| 啦啦啦啦在线视频资源| 国产又色又爽无遮挡免| 草草在线视频免费看| 日韩亚洲欧美综合| 亚洲国产日韩一区二区| 日本猛色少妇xxxxx猛交久久| 国产亚洲一区二区精品| 熟女电影av网| 视频中文字幕在线观看| 亚洲丝袜综合中文字幕| 精品久久久久久电影网| 国产69精品久久久久777片| 国产精品一区二区在线观看99| 国产国拍精品亚洲av在线观看| 热re99久久精品国产66热6| 久久国产精品大桥未久av | 午夜福利,免费看| 亚洲综合色惰| 国产精品伦人一区二区| 三级国产精品欧美在线观看| 热re99久久国产66热| h视频一区二区三区| 国产真实伦视频高清在线观看| 性色avwww在线观看| 国产成人午夜福利电影在线观看| 精品久久久久久电影网| 一级av片app| 成年美女黄网站色视频大全免费 | 久热久热在线精品观看| 中文字幕人妻丝袜制服| 久久青草综合色| 日本欧美视频一区| 免费播放大片免费观看视频在线观看| a级毛色黄片| 国产精品一区二区在线观看99| 国产亚洲午夜精品一区二区久久| 乱码一卡2卡4卡精品| 亚洲一区二区三区欧美精品| 亚洲综合精品二区| 少妇高潮的动态图| 我要看日韩黄色一级片| 亚洲美女搞黄在线观看| 国产欧美亚洲国产| 成年av动漫网址| 香蕉精品网在线| 免费少妇av软件| 精品熟女少妇av免费看| 国产午夜精品久久久久久一区二区三区| 久久久精品94久久精品| 国产免费一区二区三区四区乱码| 少妇人妻精品综合一区二区| a级毛色黄片| 色94色欧美一区二区| 国产成人精品婷婷| 美女视频免费永久观看网站| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品国产色婷婷电影| 亚洲欧美中文字幕日韩二区| 欧美xxxx性猛交bbbb| 国产毛片在线视频| 最黄视频免费看| 一区二区三区免费毛片| 国产视频内射| 国产av一区二区精品久久| 成人毛片60女人毛片免费| 日韩三级伦理在线观看| 亚洲人与动物交配视频| 国产探花极品一区二区| 亚洲精品日韩av片在线观看| 久久久久人妻精品一区果冻| 三级国产精品片| 久热这里只有精品99| 日日啪夜夜撸| 纵有疾风起免费观看全集完整版| 人人妻人人看人人澡| 久久国产精品大桥未久av | 一本色道久久久久久精品综合| 男人舔奶头视频| 国产精品99久久久久久久久| 一级,二级,三级黄色视频| 99九九线精品视频在线观看视频| 久久人人爽人人片av| 乱码一卡2卡4卡精品| 亚洲久久久国产精品| 日韩伦理黄色片| 卡戴珊不雅视频在线播放| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品999| 日韩伦理黄色片| 免费观看的影片在线观看| 国产精品麻豆人妻色哟哟久久| 在线观看av片永久免费下载| 美女xxoo啪啪120秒动态图| 婷婷色av中文字幕| 日韩,欧美,国产一区二区三区| 国产黄频视频在线观看| 综合色丁香网| 国产精品久久久久久久久免| 国产成人精品福利久久| 如日韩欧美国产精品一区二区三区 | 国内精品宾馆在线| 久久久久久久大尺度免费视频| 黄色日韩在线| 国产国拍精品亚洲av在线观看| 丰满饥渴人妻一区二区三| 精品酒店卫生间| 精品久久久久久久久av| 欧美 日韩 精品 国产| 91久久精品电影网| 91久久精品国产一区二区成人| 精品一品国产午夜福利视频| 久久久久精品久久久久真实原创| 人妻人人澡人人爽人人| 欧美人与善性xxx| 精品午夜福利在线看| 久久久久久久大尺度免费视频| 久久久久久久久久久久大奶| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩精品一区二区| 少妇丰满av| 少妇人妻精品综合一区二区| 精品国产一区二区三区久久久樱花| 久久综合国产亚洲精品| 日本欧美视频一区| 久久久久精品久久久久真实原创| 少妇 在线观看| 熟女电影av网| 搡老乐熟女国产| 亚洲图色成人| av网站免费在线观看视频| 在线观看美女被高潮喷水网站| 十八禁高潮呻吟视频 | 国产美女午夜福利| 岛国毛片在线播放| 亚洲自偷自拍三级| 久久精品久久久久久久性| 老司机影院毛片| 国产亚洲最大av| 亚洲av在线观看美女高潮| 边亲边吃奶的免费视频| 男女国产视频网站| 我要看黄色一级片免费的| 国产伦在线观看视频一区| 少妇人妻久久综合中文| 国产视频内射| 国产有黄有色有爽视频| 久久午夜综合久久蜜桃| 黑人巨大精品欧美一区二区蜜桃 | kizo精华| www.色视频.com| 在线观看www视频免费| 久久久久久久久久人人人人人人| 激情五月婷婷亚洲| 如何舔出高潮| 国产精品久久久久久久电影| 成人影院久久| 丁香六月天网| 街头女战士在线观看网站| 好男人视频免费观看在线| 国产在线男女| 精品国产一区二区久久| 少妇人妻精品综合一区二区| 精品酒店卫生间| 有码 亚洲区| 女的被弄到高潮叫床怎么办| 三级国产精品欧美在线观看| 9色porny在线观看| 国产乱人偷精品视频| 午夜福利网站1000一区二区三区| 亚洲,一卡二卡三卡| 欧美变态另类bdsm刘玥| 日本色播在线视频| 日本欧美国产在线视频| 免费看不卡的av| 看非洲黑人一级黄片| 成人无遮挡网站| 两个人免费观看高清视频 | 精品人妻一区二区三区麻豆| 久久99精品国语久久久| 高清欧美精品videossex| 99国产精品免费福利视频| 五月玫瑰六月丁香| 国产欧美日韩一区二区三区在线 | av在线观看视频网站免费| 我要看黄色一级片免费的| 人人妻人人看人人澡| av免费在线看不卡| 51国产日韩欧美| 黑丝袜美女国产一区| 国国产精品蜜臀av免费| 国产精品久久久久久久久免| 欧美精品人与动牲交sv欧美| 91久久精品国产一区二区成人| 日日摸夜夜添夜夜添av毛片| 久久免费观看电影| 国产亚洲精品久久久com| 人妻 亚洲 视频| 国产一区二区三区av在线| 国产精品99久久99久久久不卡 | 大香蕉久久网| 国产男女内射视频| 自线自在国产av| freevideosex欧美| 久久综合国产亚洲精品| 国产免费又黄又爽又色| 日韩中字成人| 韩国av在线不卡| 九九在线视频观看精品| 色视频在线一区二区三区| 国产色爽女视频免费观看| freevideosex欧美| 最后的刺客免费高清国语| 亚洲无线观看免费| 777米奇影视久久| 久久精品国产亚洲av天美| 国产 一区精品| 少妇人妻精品综合一区二区| 你懂的网址亚洲精品在线观看| 在线观看美女被高潮喷水网站| 黄色欧美视频在线观看| 亚洲天堂av无毛| 插阴视频在线观看视频| 十八禁网站网址无遮挡 | 纯流量卡能插随身wifi吗| 精品国产一区二区三区久久久樱花| 天堂8中文在线网| 国产日韩欧美在线精品| 免费观看a级毛片全部| 中文乱码字字幕精品一区二区三区| 成人免费观看视频高清| 午夜免费鲁丝| 国产av一区二区精品久久| 七月丁香在线播放| 国产亚洲一区二区精品| 国产免费又黄又爽又色| 在线播放无遮挡| 亚洲av成人精品一二三区| 麻豆精品久久久久久蜜桃| 寂寞人妻少妇视频99o| 成人影院久久| 国产乱人偷精品视频| 国产精品国产三级专区第一集| 一本久久精品| 精品少妇黑人巨大在线播放| 亚洲欧洲日产国产| 亚洲欧美清纯卡通| 观看免费一级毛片| 亚洲精品日韩在线中文字幕| 精品少妇内射三级| 尾随美女入室| 91aial.com中文字幕在线观看| 亚洲欧美日韩卡通动漫| 秋霞伦理黄片| 天天躁夜夜躁狠狠久久av| 日韩成人av中文字幕在线观看| 国产淫语在线视频| 国产成人免费观看mmmm| 亚洲国产成人一精品久久久| 亚洲精品乱码久久久久久按摩| 美女cb高潮喷水在线观看| av有码第一页| 在线观看国产h片| 久久婷婷青草| 18禁裸乳无遮挡动漫免费视频| 嘟嘟电影网在线观看| 汤姆久久久久久久影院中文字幕| 久久久久精品性色| 免费观看在线日韩| 一二三四中文在线观看免费高清| 免费大片黄手机在线观看| 美女内射精品一级片tv| 久久人妻熟女aⅴ| 黄色视频在线播放观看不卡| 久久99精品国语久久久| 成年av动漫网址| 日韩大片免费观看网站| 日韩欧美一区视频在线观看 | 久久久久久久久久久丰满| 久久久久久伊人网av| 欧美97在线视频| 亚洲精品国产色婷婷电影| 一个人看视频在线观看www免费| 91久久精品国产一区二区成人| 插逼视频在线观看| 精品人妻熟女av久视频| 街头女战士在线观看网站| 熟妇人妻不卡中文字幕| xxx大片免费视频| 久久毛片免费看一区二区三区| 黄色怎么调成土黄色| 国产欧美亚洲国产| 免费观看的影片在线观看| 天天躁夜夜躁狠狠久久av| 欧美3d第一页| 极品人妻少妇av视频| 日韩视频在线欧美| 青春草亚洲视频在线观看| www.av在线官网国产| 亚洲综合精品二区| 亚洲一区二区三区欧美精品| 三上悠亚av全集在线观看 | 亚洲欧美成人精品一区二区| 激情五月婷婷亚洲| 亚洲精华国产精华液的使用体验| 天美传媒精品一区二区| 国产有黄有色有爽视频| 国产中年淑女户外野战色| 五月伊人婷婷丁香| 99热这里只有是精品50| 91久久精品国产一区二区三区| 国产探花极品一区二区| 亚洲欧美清纯卡通| 日产精品乱码卡一卡2卡三| 午夜精品国产一区二区电影| 久久久久视频综合| 亚洲性久久影院| 王馨瑶露胸无遮挡在线观看| 高清黄色对白视频在线免费看 | 极品人妻少妇av视频| 日韩电影二区| 久久精品久久精品一区二区三区| 中文字幕人妻丝袜制服| 国产精品国产三级国产专区5o| 久久久久精品久久久久真实原创| 中文字幕久久专区| 欧美日韩视频高清一区二区三区二| 黄色日韩在线| 国产淫片久久久久久久久| 亚洲欧美精品专区久久| 在线观看人妻少妇| 我要看日韩黄色一级片| 欧美最新免费一区二区三区| 国产成人aa在线观看| 哪个播放器可以免费观看大片| 特大巨黑吊av在线直播| 婷婷色综合大香蕉| 18禁动态无遮挡网站| 亚洲av国产av综合av卡| 亚州av有码| 日韩 亚洲 欧美在线| 亚洲av国产av综合av卡| 日韩一区二区三区影片| 日本色播在线视频| 国产爽快片一区二区三区| 69精品国产乱码久久久| 久久精品国产a三级三级三级| 少妇人妻 视频| 中文字幕人妻熟人妻熟丝袜美| 极品人妻少妇av视频| 亚洲精品乱码久久久v下载方式| 亚洲欧美精品专区久久| 色视频www国产| 2022亚洲国产成人精品| 最近的中文字幕免费完整| 亚洲精品视频女| 亚洲人与动物交配视频| freevideosex欧美| 国产视频内射| 国产精品无大码| 女的被弄到高潮叫床怎么办| 少妇被粗大猛烈的视频| 内地一区二区视频在线| 少妇人妻 视频| 日日啪夜夜爽| 99久久中文字幕三级久久日本| 精品人妻一区二区三区麻豆| 色吧在线观看| tube8黄色片| 国产成人精品婷婷| 久久这里有精品视频免费| 欧美另类一区| 中国美白少妇内射xxxbb| 欧美日韩视频高清一区二区三区二| 久久久久久久久久成人| 人妻人人澡人人爽人人| av播播在线观看一区| 一个人免费看片子| 亚洲情色 制服丝袜| 成人亚洲欧美一区二区av| 少妇熟女欧美另类| 亚洲国产精品专区欧美| 国产伦理片在线播放av一区| 性色avwww在线观看| 久久99蜜桃精品久久| 久久久久久久亚洲中文字幕| 久久人人爽av亚洲精品天堂| 国产熟女欧美一区二区| 欧美3d第一页| 午夜激情久久久久久久| 久久久久国产网址| 欧美 亚洲 国产 日韩一| 欧美日韩视频精品一区| 国产淫语在线视频| 毛片一级片免费看久久久久| 日日啪夜夜撸| 国产色爽女视频免费观看| 天美传媒精品一区二区| 亚洲精品aⅴ在线观看| 精品99又大又爽又粗少妇毛片| 黑人猛操日本美女一级片| 大香蕉久久网| 亚洲精品第二区| 女人精品久久久久毛片| 精品少妇黑人巨大在线播放| 精品一区二区三区视频在线| 亚洲欧洲国产日韩| 蜜桃在线观看..| 婷婷色综合www| 亚洲美女黄色视频免费看| 又黄又爽又刺激的免费视频.| 嫩草影院新地址| 91久久精品电影网| 久久精品国产鲁丝片午夜精品| 丰满迷人的少妇在线观看| 欧美日韩一区二区视频在线观看视频在线| 国产精品伦人一区二区| 精品久久国产蜜桃| 大陆偷拍与自拍| 成人二区视频| 久久韩国三级中文字幕| a级毛片在线看网站| 极品少妇高潮喷水抽搐| 色视频在线一区二区三区| 天天操日日干夜夜撸| 丝袜在线中文字幕| 欧美另类一区| 国产精品久久久久久久电影| 欧美成人精品欧美一级黄| 欧美国产精品一级二级三级 | 日日摸夜夜添夜夜爱| 亚洲国产精品一区二区三区在线| 国产爽快片一区二区三区| 国产高清国产精品国产三级| 午夜福利在线观看免费完整高清在| 在线亚洲精品国产二区图片欧美 | 久久人人爽人人片av| 久久久久视频综合| av有码第一页| 三级国产精品片| 午夜激情久久久久久久| 夫妻性生交免费视频一级片| 亚洲精品中文字幕在线视频 | 国产色爽女视频免费观看| 亚洲欧美日韩卡通动漫| 婷婷色麻豆天堂久久| 久久久久久久亚洲中文字幕| 少妇猛男粗大的猛烈进出视频| 久久99精品国语久久久| 青春草视频在线免费观看| 91精品伊人久久大香线蕉| 2022亚洲国产成人精品| 成人综合一区亚洲| 亚洲高清免费不卡视频|