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

    依據(jù)斷層摩擦準(zhǔn)則估算的地殼主應(yīng)力量值關(guān)系式及其參數(shù)分析

    2016-08-22 11:16:18曹慧靜崔效鋒樊文杰
    地震地質(zhì) 2016年2期

    曹慧靜 崔效鋒 樊文杰

    (中國地震局地殼應(yīng)力研究所, 地殼動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室, 北京 100085)

    ?

    依據(jù)斷層摩擦準(zhǔn)則估算的地殼主應(yīng)力量值關(guān)系式及其參數(shù)分析

    曹慧靜崔效鋒*樊文杰

    (中國地震局地殼應(yīng)力研究所, 地殼動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室, 北京100085)

    基于Zoback(1992b)提出的估算地殼構(gòu)造應(yīng)力量值的思路和2個(gè)假設(shè), 考慮了3個(gè)主應(yīng)力分別為垂直主應(yīng)力(對應(yīng)3種構(gòu)造應(yīng)力類型)的情況, 比較系統(tǒng)地建立了依據(jù)應(yīng)力形因子和斷層摩擦強(qiáng)度估算構(gòu)造主應(yīng)力量值的關(guān)系式。討論了摩擦系數(shù)、 孔隙壓力系數(shù)和應(yīng)力形因子對構(gòu)造主應(yīng)力量值的影響: 在最大主應(yīng)力為近水平的情況下(應(yīng)力類型為走滑型或逆斷型), 最大主應(yīng)力量值(或隨深度線性增加的斜率)隨摩擦系數(shù)取值的增大而增大, 隨孔隙壓力系數(shù)取值的增大而減小。在最小主應(yīng)力為近水平的情況下(應(yīng)力類型為走滑型或正斷型), 最小主應(yīng)力量值(或隨深度線性增加的斜率)隨摩擦系數(shù)取值的增大而減小, 隨孔隙壓力系數(shù)的增大而增大。并且, 3個(gè)參數(shù)對主應(yīng)力量值估算的影響很大, 如果摩擦系數(shù)取值較大, 同時(shí)孔隙壓力系數(shù)取值較小, 最大主應(yīng)力隨深度增加的斜率與最小主應(yīng)力隨深度增加的斜率之間的差異就大, 有可能導(dǎo)致不合理的結(jié)果出現(xiàn)。該方法只是對地殼巖石發(fā)生脆性破裂或摩擦滑動(dòng)時(shí)的構(gòu)造應(yīng)力量值進(jìn)行近似估算, 估算結(jié)果并不是地殼穩(wěn)定狀態(tài)下的構(gòu)造應(yīng)力量值。

    應(yīng)力量值摩擦系數(shù)孔隙壓力應(yīng)力形因子構(gòu)造應(yīng)力

    0 引言

    目前對地殼應(yīng)力量值的研究較少, 尤其是地殼深部應(yīng)力量值難以獲取。 雖然原地應(yīng)力測量可以直接獲取淺層地殼應(yīng)力狀態(tài), 得到地殼應(yīng)力的大小和方向, 但由于技術(shù)難度大、 資金耗費(fèi)高, 深部應(yīng)力量值數(shù)據(jù)較少。有些學(xué)者嘗試通過對淺層應(yīng)力量值的回歸分析推測地殼深部應(yīng)力大小, 并且取得了一定的進(jìn)展(王艷華等, 2012)?;瑒?dòng)方向擬合法的提出(Angelier, 1979; Gephartetal., 1984; 許忠淮, 1985), 使得人們不僅可以利用震源機(jī)制解、 斷層擦痕、 跨斷層形變觀測等數(shù)據(jù)資料對地殼構(gòu)造應(yīng)力的方向進(jìn)行反演計(jì)算, 而且可以得到反映主應(yīng)力相對大小的應(yīng)力形因子(謝富仁等, 1993, 2001; 崔效鋒等, 2006; Hardebeck等, 2006; 張浪平等, 2014)。應(yīng)力形因子為進(jìn)一步研究分析應(yīng)力量值提供了幫助。

    Zoback(1992b)利用滑動(dòng)擬合法給出的應(yīng)力形因子, 提出了根據(jù)優(yōu)勢方位斷層摩擦準(zhǔn)則對地殼應(yīng)力量值進(jìn)行估算的思路, 該思路是估算分析地殼巖石發(fā)生脆性破裂或摩擦滑動(dòng)時(shí)的構(gòu)造應(yīng)力量值的有效途徑之一, 之后有學(xué)者采用該思路對多個(gè)地區(qū)的構(gòu)造應(yīng)力量值進(jìn)行了估算分析(Plenefischetal., 1997; 趙建濤等, 2002)。該思路需要借助巖石力學(xué)參數(shù)(摩擦系數(shù)、 孔隙壓力)對構(gòu)造應(yīng)力量值進(jìn)行估算; 因此, 科學(xué)合理地選取巖石力學(xué)參數(shù)就成為該方法中的1個(gè)重要問題。本文利用Zoback(1992b)提出的構(gòu)造應(yīng)力量值估算思路, 比較系統(tǒng)地建立了依據(jù)應(yīng)力形因子和斷層摩擦強(qiáng)度估算主應(yīng)力量值的關(guān)系式; 在此基礎(chǔ)上, 對影響主應(yīng)力量值估算的巖石力學(xué)參數(shù)以及主應(yīng)力量值估算關(guān)系式的適用性進(jìn)行了初步討論。

    1 方法回顧

    滑動(dòng)方向擬合法是利用斷面上的剪應(yīng)力方向去擬合斷面上的滑動(dòng)方向, 從而獲取平均偏應(yīng)力張量的1種方法(Angelier, 1979; Gephartetal., 1984; 許忠淮, 1985)。滑動(dòng)方向擬合法在得到應(yīng)力張量的3個(gè)主應(yīng)力方向的同時(shí), 還可以得到1個(gè)主應(yīng)力量值之間比值的參數(shù)——應(yīng)力形因子Φ:

    (1)

    式(1)中, σ1、 σ2、 σ3分別為最大主應(yīng)力、 中間主應(yīng)力和最小主應(yīng)力。

    增加2個(gè)約束條件, 就可以借助應(yīng)力形因子Φ對3個(gè)主應(yīng)力量值進(jìn)行估算; 為此, Zoback提出了2個(gè)假設(shè), 建立了2個(gè)約束條件(Zoback, 1992b)。2個(gè)假設(shè)為:

    (1)若構(gòu)造應(yīng)力張量中有1個(gè)主應(yīng)力方向近垂直, 則近垂直的主應(yīng)力量值(用σv表示)近似于巖石靜壓力:

    (2)

    (2)地殼中最大應(yīng)力差受優(yōu)勢方位斷層摩擦強(qiáng)度所限(Byerlee, 1967; Jaegeretal., 1979):

    (3)

    式(2)、(3)中, ρ為巖石密度, g為重力加速度, μc為摩擦系數(shù), h為地殼深度, P0為孔隙壓力, σ1和σ3分別為最大主應(yīng)力和最小主應(yīng)力。

    全球大陸多數(shù)地區(qū)的應(yīng)力狀態(tài)為走滑應(yīng)力類型(Zoback, 1992a; 謝富仁等, 2003)。但也有一些地區(qū)存在不同的情況。以中國為例, 大量研究表明, 中國構(gòu)造應(yīng)力場在不同地區(qū)的應(yīng)力作用方向和構(gòu)造應(yīng)力類型不同, 既有最大主應(yīng)力為近水平的逆斷應(yīng)力類型和走滑應(yīng)力類型, 也有最大主應(yīng)力傾角很大的正斷應(yīng)力類型; 如山西地震帶和鄂爾多斯周邊地區(qū)多為正斷應(yīng)力類型, 而天山地區(qū)和龍門山地區(qū)則為逆斷應(yīng)力類型。因此根據(jù)不同的應(yīng)力類型, 以不同的主應(yīng)力作為垂直應(yīng)力, 運(yùn)用上述式(1)—(3)才可以對3個(gè)主應(yīng)力量值進(jìn)行合理的估算。

    從式(1)—(3)還可以看出, 除垂直應(yīng)力外, 摩擦系數(shù)、 孔隙壓力等巖石力學(xué)參數(shù)的大小也對主應(yīng)力量值的估算結(jié)果有很大的影響。 前人對這些巖石力學(xué)參數(shù)的影響沒有做過細(xì)致分析, 只是根據(jù)相應(yīng)的巖石力學(xué)參數(shù)研究結(jié)果, 選取一些特定的經(jīng)驗(yàn)值或理論值對主應(yīng)力量值進(jìn)行估算; 如趙建濤等(2002)選取摩擦系數(shù)為0.85和0.6, 對唐山震源區(qū)的主應(yīng)力量值進(jìn)行了估算分析。本文對摩擦系數(shù)、 孔隙壓力和應(yīng)力形因子的取值大小對主應(yīng)力量值的影響進(jìn)行了初步分析。

    2 主應(yīng)力量值的估算關(guān)系式

    上述估算方法的1個(gè)假設(shè)是構(gòu)造應(yīng)力張量中有1個(gè)主應(yīng)力近垂直, 因此, 我們分別考慮3種構(gòu)造應(yīng)力類型: 1)走滑應(yīng)力類型, 此時(shí)中間主應(yīng)力σ2近垂直; 2)正斷應(yīng)力類型, 此時(shí)最大主應(yīng)力σ1近垂直; 3)逆斷應(yīng)力類型, 此時(shí)最小主應(yīng)力σ3近垂直。通過(1)—(3)3個(gè)等量關(guān)系式, 分別以σ1、 σ2或σ3代替σv, 即可得到估算3個(gè)主應(yīng)力量值的關(guān)系式。

    為了便于推導(dǎo)和分析, 令

    (4)

    則式(3)可表示為

    (5)

    式(4)、 (5)中A為與摩擦系數(shù)μc有關(guān)的參數(shù), 并且A隨著摩擦系數(shù)的增加而單調(diào)增加。

    (1)應(yīng)力類型為走滑型。此時(shí)σ2為垂直應(yīng)力, 即σ2=σv; 其他2個(gè)主應(yīng)力水平, 利用式(1)—(3)和式(5)可以推出計(jì)算式

    (6)

    (7)

    (8)

    對于孔隙壓力, 可表示為1個(gè)比例系數(shù)乘以巖石靜壓力(陳颙等, 2009):

    (9)

    本文稱參數(shù)λ為孔隙壓力系數(shù)。 一般來說,λ的取值范圍為0~1, 其中幾種特殊的情況為: 當(dāng)λ=0時(shí), 孔隙壓力為零, 即相當(dāng)于地殼巖石孔隙中無水干燥的情況; 當(dāng)λ取0.4左右時(shí), 孔隙壓力近似為靜水壓力, 即地殼巖石中所有孔隙皆連通并且一直貫通至地面的情況。在不知道具體的孔隙壓力時(shí), 前人一般多用靜水壓力或相近值作為孔隙壓力的近似值(Plenefischetal., 1997; 趙建濤等, 2002)。當(dāng)λ=1時(shí), 孔隙壓力等于巖石靜壓力, 因此在實(shí)際物理狀態(tài)下,λ不可能為1。

    將式(9)帶入式(6)—(8)中, 則可以得到3個(gè)主應(yīng)力量值隨地殼深度變化的線性表達(dá)式:

    (10)

    (11)

    (12)

    從式(10)—(12)可以看出, 3個(gè)主應(yīng)力量值隨深度線性增加, 并受應(yīng)力形因子、 摩擦系數(shù)和孔隙壓力系數(shù)等3個(gè)參數(shù)的影響。另外, 當(dāng)摩擦系數(shù)取值增大時(shí),σ1值隨之變大,σ3值隨之變小; 當(dāng)孔隙壓力系數(shù)取值增大時(shí),σ1值隨之減小,σ3值隨之變大; 當(dāng)應(yīng)力形因子取值增大時(shí),σ1和σ3值隨之變小。

    (2)應(yīng)力類型為正斷型。此時(shí)σ1為垂直應(yīng)力, 即σ1=σv; 其他2個(gè)主應(yīng)力水平, 同樣可以根據(jù)式(1)—(3)、 (5)、 (9)推出3個(gè)主應(yīng)力量值隨地殼深度變化的線性表達(dá)式:

    (13)

    (14)

    (15)

    從式(13)—(15)可以看出, σ2和σ3與摩擦系數(shù)呈反比關(guān)系, 當(dāng)摩擦系數(shù)取值增大時(shí), σ2和σ3值隨之變小; σ2和σ3與孔隙壓力系數(shù)呈正比關(guān)系, 當(dāng)孔隙壓力系數(shù)取值增大時(shí), σ2和σ3值隨之變大; 應(yīng)力形因子只對σ2的變化有影響, 當(dāng)應(yīng)力形因子取值增大時(shí), σ2值隨之變大。

    (3)應(yīng)力類型為逆斷型。此時(shí)σ3為垂直應(yīng)力, 即σ3=σv; 其他2個(gè)主應(yīng)力水平, 同樣地根據(jù)式(1)—(3)、 (5)、 (9)不難推出3個(gè)主應(yīng)力量值隨地殼深度變化的線性表達(dá)式:

    (16)

    (17)

    (18)

    從式(16)—(18)可以看出, σ1和σ2與摩擦系數(shù)呈正比關(guān)系, 當(dāng)摩擦系數(shù)取值增大時(shí), σ1和σ2值隨之變大; σ1和σ2與孔隙壓力系數(shù)呈反比關(guān)系, 當(dāng)孔隙壓力系數(shù)取值增大時(shí), σ1和σ2值隨之減小; 應(yīng)力形因子只對σ2的變化有影響, 當(dāng)應(yīng)力形因子取值增大時(shí), σ2值隨之變大。

    雖然在式(1)—(3)中出現(xiàn)了3個(gè)參數(shù), 即摩擦系數(shù)、 孔隙壓力和應(yīng)力形因子; 但當(dāng)垂直應(yīng)力為最大主應(yīng)力或最小主應(yīng)力時(shí)(即上文(2)和(3) 2種情況), 最大主應(yīng)力值和最小主應(yīng)力值與應(yīng)力形因子取值的大小無關(guān), 中間主應(yīng)力值隨應(yīng)力形因子取值的增大而增大。

    3 討論

    3.1適用性討論

    本文基于Zoback(1992b)提出的估算地殼構(gòu)造應(yīng)力量值的思路, 推導(dǎo)了3種構(gòu)造應(yīng)力類型情況下的3個(gè)主應(yīng)力量值隨地殼深度變化的線性關(guān)系式。 該思路是建立在2個(gè)假設(shè)基礎(chǔ)上的, 因此有必要對上述方法給出的應(yīng)力量值做出說明。第1個(gè)假設(shè)是構(gòu)造應(yīng)力張量中的1個(gè)主應(yīng)力接近垂直, 并用巖石靜壓力估算垂直主應(yīng)力; 但實(shí)際情況中, 構(gòu)造應(yīng)力張量中的主應(yīng)力方向并不完全是水平和垂直的(石耀霖, 2004), 因此該方法只是對地殼構(gòu)造應(yīng)力量值進(jìn)行估算, 并不能得到準(zhǔn)確的地殼應(yīng)力量值。第2個(gè)假設(shè)是地殼中的最大差應(yīng)力受優(yōu)勢方位斷層摩擦強(qiáng)度所限, 此時(shí)巖石為脆性, 因此該方法只是對地殼所能承受的最大應(yīng)力量值進(jìn)行估算??偟膩碚f, 該方法只是對地殼巖石發(fā)生脆性破裂或摩擦滑動(dòng)時(shí)的構(gòu)造應(yīng)力量值進(jìn)行了近似估算, 估算結(jié)果并不是地殼穩(wěn)定狀態(tài)下的構(gòu)造應(yīng)力量值。

    圖1 巖石強(qiáng)度隨應(yīng)力變化示意圖(據(jù)Kohlstedt et al., 1995修改)Fig. 1 The change of rock strength with stress(adapted after Kohlstedt et al., 1995).BDT為脆韌性轉(zhuǎn)化處, BPT為脆塑性轉(zhuǎn)化處

    另外需要說明的是, 式(10)—(18)給出的3個(gè)主應(yīng)力量值隨地殼深度變化的線性表達(dá)式只對脆性的上地殼適用。地殼中影響巖石變形方式的因素有很多, 如圍壓、 溫度、 孔隙壓、 應(yīng)變速率、 巖性、 結(jié)構(gòu)等, 所以, 不同的學(xué)者對巖石變形的研究得到的結(jié)果不盡相同。在地殼中巖石干燥的情況下, 通常認(rèn)為巖石脆韌性的轉(zhuǎn)化發(fā)生在庫倫準(zhǔn)則和摩擦滑動(dòng)曲線的相交處(Kohlstedtetal., 1995; 圖1)。由于隨著地殼深度的增加, 溫度和壓力是上升的, 巖石狀態(tài)會(huì)發(fā)生轉(zhuǎn)化; 當(dāng)?shù)貧ど畈康膸r石發(fā)生流動(dòng), 巖石發(fā)生明顯的蠕變時(shí), 地殼應(yīng)力量值與蠕變速率及溫度分布有關(guān)(Kohlstedtetal., 1995; 圖2)。結(jié)合圖2 我們知道, 地殼中巖石脆性部分可能從地表延伸到十幾千米深處, 在這部分脆性破裂區(qū)域, 才可以運(yùn)用本文給出的關(guān)系式, 對地殼巖石脆性破裂或摩擦滑動(dòng)時(shí)的構(gòu)造應(yīng)力量值進(jìn)行近似的估算。

    圖2 大洋巖石圈和大陸巖石圈差應(yīng)力隨深度的變化圖(據(jù)Kohlstedt et al., 1995修改)Fig. 2 Change of differential stress with depth in oceanic lithosphere and continental lithosphere(adapted after Kohlstedt et al., 1995).

    3.2關(guān)于關(guān)系式中參數(shù)的討論

    由上節(jié)給出的關(guān)系式可以看出, 主應(yīng)力量值隨深度線性增加, 線性增加的斜率除了與巖石密度和重力加速度有關(guān)外, 還受摩擦系數(shù)、 孔隙壓力系數(shù)和應(yīng)力形因子取值的影響; 并且, 隨著摩擦系數(shù)、 孔隙壓力系數(shù)和應(yīng)力形因子取值的不同, 主應(yīng)力量值隨深度線性增加的斜率存在著較大的變化。

    3.2.1關(guān)于摩擦系數(shù)

    在最大主應(yīng)力σ1近水平的情況下(對應(yīng)的應(yīng)力類型為走滑型或逆斷型), 最大主應(yīng)力量值(或隨深度線性增加的斜率)隨摩擦系數(shù)取值的增大而增大(式(10)、 (16))。在最小主應(yīng)力σ3近水平的情況下(應(yīng)力類型為走滑型或正斷型), 最小主應(yīng)力量值(或隨深度線性增加的斜率)隨摩擦系數(shù)取值的增大而減小(式(12)、 (15))。當(dāng)中間主應(yīng)力σ2近水平時(shí)(應(yīng)力類型為正斷型或逆斷型), 摩擦系數(shù)對中間主應(yīng)力量值的影響趨勢與另一水平主應(yīng)力相同(式(14)—(17))。

    Byerlee(1978)在對巖石摩擦實(shí)驗(yàn)資料進(jìn)行分析后提出, 當(dāng)圍壓<200MPa時(shí), 巖石的摩擦系數(shù)為0.85; 當(dāng)圍壓≥200MPa時(shí), 巖石的摩擦系數(shù)為0.6。近年來很多研究也發(fā)現(xiàn), 某些地區(qū)斷層或巖石的摩擦系數(shù)可能較低, 甚至可以低至0.1(Plenefischetal., 1997; 單斌等, 2009; 何昌榮等, 2011)。如果摩擦系數(shù)可以在0.1~0.85范圍內(nèi)選取的話, 則最大主應(yīng)力量值(或隨深度線性增加的斜率)的變動(dòng)范圍就很大; 也就是說, 如果摩擦系數(shù)取值較大, 最大主應(yīng)力量值(或隨深度線性增加的斜率)就會(huì)很高。

    為了清晰地展現(xiàn)摩擦系數(shù)對主應(yīng)力量值估算的影響, 這里我們以走滑應(yīng)力類型的關(guān)系式為例, 固定應(yīng)力形因子和孔隙壓力系數(shù), 選取不同的摩擦系數(shù)估算主應(yīng)力量值。圖3 顯示了當(dāng)應(yīng)力形因子和孔隙壓力系數(shù)都取0.5, 摩擦系數(shù)分別取0.1、 0.4、 0.6和0.85時(shí), 最大主應(yīng)力和最小主應(yīng)力量值的估算結(jié)果。 其中, 圖3a清晰地展示了最大主應(yīng)力量值與摩擦系數(shù)呈正比, 即當(dāng)摩擦系數(shù)取值增大時(shí), 最大主應(yīng)力量值(隨深度線性增加的斜率)隨之變大; 圖3b顯示了最小主應(yīng)力量值與摩擦系數(shù)呈反比, 即當(dāng)摩擦系數(shù)取值增大時(shí), 最小主應(yīng)力量值(隨深度線性增加的斜率)隨之變小。

    圖3 摩擦系數(shù)對主應(yīng)力量值估算結(jié)果的影響Fig. 3 The magnitude of principal stresses influenced by friction coefficient.a 最大主應(yīng)力(σ1)量值估算結(jié)果; b 最小主應(yīng)力(σ3)量值估算結(jié)果

    3.2.2關(guān)于孔隙壓力系數(shù)

    孔隙壓力系數(shù)對主應(yīng)力量值的影響趨勢與摩擦系數(shù)相反; 即在最大主應(yīng)力σ1不是垂直應(yīng)力的情況下(應(yīng)力類型為走滑型或逆斷型), 最大主應(yīng)力量值(或隨深度線性增加的斜率)隨孔隙壓力系數(shù)取值的增大而減小(式(10)、 (16))。在最小主應(yīng)力σ3不是垂直應(yīng)力的情況下(應(yīng)力類型為走滑型或正斷型), 最小主應(yīng)力量值隨孔隙壓力系數(shù)的增大而增大(式(12)、 (15))。當(dāng)中間主應(yīng)力σ2為水平主應(yīng)力時(shí)(應(yīng)力類型為正斷型或逆斷型), 孔隙壓力系數(shù)對中間主應(yīng)力量值的影響趨勢與另一水平主應(yīng)力相同(式(14)—(17))。

    圖4 孔隙壓力系數(shù)對主應(yīng)力量值估算結(jié)果的影響Fig. 4 The magnitude of principal stresses influenced by pore pressure.a 最大主應(yīng)力(σ1)量值估算結(jié)果; b 最小主應(yīng)力(σ3)量值估算結(jié)果

    如前文所述, 如果孔隙壓力系數(shù)可以在0≤λ<1范圍內(nèi)選取的話, 則最大主應(yīng)力量值(或隨深度線性增加的斜率)的變動(dòng)范圍就很大。 也就是說, 如果孔隙壓力系數(shù)取值較大, 最大主應(yīng)力量值(或隨深度線性增加的斜率)就會(huì)小; 如果孔隙壓力系數(shù)取值較小, 最大主應(yīng)力量值(或隨深度線性增加的斜率)就會(huì)很大。圖4 顯示了當(dāng)應(yīng)力形因子和摩擦系數(shù)都取0.5, 孔隙壓力系數(shù)分別取0.1、 0.4、 0.6和0.9時(shí), 運(yùn)用走滑型關(guān)系式給出的最大主應(yīng)力和最小主應(yīng)力量值的估算結(jié)果, 直觀地展現(xiàn)了最大主應(yīng)力量值與孔隙壓力系數(shù)呈反比(圖4a), 最小主應(yīng)力量值與孔隙壓力系數(shù)呈正比(圖4b)。

    3.2.3關(guān)于應(yīng)力形因子

    根據(jù)應(yīng)力形因子的定義, 其取值范圍也在0~1之間。 當(dāng)應(yīng)力形因子為0時(shí), 中間主應(yīng)力等于最小主應(yīng)力; 當(dāng)應(yīng)力形因子為1時(shí), 中間主應(yīng)力等于最大主應(yīng)力; 當(dāng)應(yīng)力形因子等于0.5時(shí), 3個(gè)主應(yīng)力量值呈等差排列。

    圖5 應(yīng)力形因子對主應(yīng)力量值估算結(jié)果的影響Fig. 5 The magnitude of principal stresses influenced by stress form factor.

    除應(yīng)力形因子等于0和1這2種特殊情形外, 對于正斷和逆斷應(yīng)力類型來說, 中間主應(yīng)力σ2近水平, 應(yīng)力形因子只對中間主應(yīng)力隨深度線性增加的斜率有影響, 最大和最小主應(yīng)力隨深度線性增加的斜率與應(yīng)力形因子無關(guān)。對于走滑應(yīng)力類型, 中間主應(yīng)力σ2近垂直, 其量值用巖石靜壓力來估算; 這時(shí), 應(yīng)力形因子直接影響到最大和最小主應(yīng)力隨深度線性增加的斜率。

    同樣, 我們以走滑應(yīng)力類型的關(guān)系式為例, 取孔隙壓力系數(shù)為0.4(約為靜水壓力), 取摩擦系數(shù)為0.6, 應(yīng)力形因子分別取為0、 0.2、 0.4、 0.6、 0.8和1, 對主應(yīng)力量值進(jìn)行估算(圖5)。由圖5 可以看出, 最大主應(yīng)力量值和最小主應(yīng)力量值與應(yīng)力形因子呈反比, 隨著應(yīng)力形因子的增大, 最大主應(yīng)力量值和最小主應(yīng)力量值均減小。中間主應(yīng)力為垂直壓力, 用巖石靜壓力估算。

    3.2.4參數(shù)對主應(yīng)力量值的影響

    摩擦系數(shù)和孔隙壓力是地殼巖石力學(xué)參數(shù), 由上述分析不難看出, 它們的取值大小直接影響到主應(yīng)力量值的估算結(jié)果, 圖3 —5 只展現(xiàn)了部分參數(shù)取值情況下的估算結(jié)果。如果摩擦系數(shù)取值較大, 同時(shí)孔隙壓力系數(shù)取值較小, 最大主應(yīng)力隨深度增加斜率與最小主應(yīng)力隨深度增加斜率的差異就大, 在上地殼深部, 最大主應(yīng)力與最小主應(yīng)力的量值之差(即差應(yīng)力)就會(huì)很大; 反之, 最大主應(yīng)力與最小主應(yīng)力的斜率和量值之差就會(huì)小些。例如, 當(dāng)應(yīng)力形因子取0.1, 摩擦系數(shù)取0.8, 孔隙壓力系數(shù)取0.1時(shí), 由式(10)估算的最大主應(yīng)力量值在10km深處達(dá)800MPa, 此時(shí), 最大主應(yīng)力與最小主應(yīng)力的量值之差達(dá)595MPa, 這顯然偏高(參見3.1節(jié)); 巖石在這么高的應(yīng)力下已經(jīng)轉(zhuǎn)換為塑性(Kohlstedtetal., 1995; Xuetal., 2014)。

    這一方面說明, 只有對1個(gè)地區(qū)的地殼巖石力學(xué)參數(shù)和應(yīng)力形因子進(jìn)行深入細(xì)致的研究, 獲得比較可靠的摩擦系數(shù)、 孔隙壓力和應(yīng)力形因子等參數(shù), 才有可能給出研究區(qū)比較可信合理的應(yīng)力量值估算結(jié)果。另一方面也帶給我們1個(gè)啟示, 在不能很好地確定地殼巖石力學(xué)參數(shù)和應(yīng)力形因子的情況下, 如果能夠利用多學(xué)科的研究成果, 對摩擦系數(shù)、 孔隙壓力和應(yīng)力形因子等參數(shù)的取值范圍進(jìn)行限制, 就可以從宏觀上對地殼巖石發(fā)生脆性破裂或摩擦滑動(dòng)時(shí)的構(gòu)造應(yīng)力量值(即地殼能夠承受的最大應(yīng)力)范圍進(jìn)行估算, 進(jìn)而借助多學(xué)科方法開展地殼應(yīng)力積累變化研究, 以便獲得對地殼構(gòu)造應(yīng)力量值更科學(xué)、 更合理的認(rèn)識(shí)。

    4 結(jié)論

    (1)基于Zoback(1992b)提出的估算地殼構(gòu)造應(yīng)力量值的思路和2個(gè)假設(shè), 考慮了3個(gè)主應(yīng)力分別為垂直主應(yīng)力(對應(yīng)3種構(gòu)造應(yīng)力類型)的情況, 比較系統(tǒng)地建立了依據(jù)應(yīng)力形因子和斷層摩擦強(qiáng)度估算主應(yīng)力量值的關(guān)系式。

    (2)對影響主應(yīng)力量值估算的3個(gè)參數(shù)(摩擦系數(shù)、 孔隙壓力、 應(yīng)力形因子)進(jìn)行了初步分析。在最大主應(yīng)力為近水平的情況下(對應(yīng)的應(yīng)力類型為走滑型或逆斷型), 最大主應(yīng)力量值(或隨深度線性增加的斜率)隨摩擦系數(shù)取值的增大而增大, 隨孔隙壓力系數(shù)取值的增大而減小。在最小主應(yīng)力為近水平的情況下(對應(yīng)的應(yīng)力類型為走滑型或正斷型), 最小主應(yīng)力量值(或隨深度線性增加的斜率)隨摩擦系數(shù)取值的增大而減小, 隨孔隙壓力系數(shù)的增大而增大。當(dāng)中間主應(yīng)力為近水平時(shí), 摩擦系數(shù)和孔隙壓力系數(shù)對中間主應(yīng)力量值的影響趨勢與另1水平主應(yīng)力相同。3個(gè)參數(shù)對主應(yīng)力量值估算的影響很大, 如果摩擦系數(shù)取值較大, 同時(shí)孔隙壓力系數(shù)取值較小, 最大主應(yīng)力隨深度增加斜率與最小主應(yīng)力隨深度增加斜率的差異就大, 有可能導(dǎo)致不合理的結(jié)果出現(xiàn)。

    (3)本文給出的主應(yīng)力量值隨深度線性變化的估算關(guān)系式, 是建立在地殼中最大差應(yīng)力受優(yōu)勢方位斷層摩擦強(qiáng)度所限這一假設(shè)基礎(chǔ)上的, 此時(shí)巖石為脆性; 因此, 該方法只是對地殼巖石發(fā)生脆性破裂或摩擦滑動(dòng)時(shí)的構(gòu)造應(yīng)力量值進(jìn)行近似估算, 估算結(jié)果并不是地殼穩(wěn)定狀態(tài)下的構(gòu)造應(yīng)力量值。

    陳顒, 黃庭芳, 劉恩儒. 2001. 巖石物理學(xué) [M]. 北京: 北京大學(xué)出版社.

    CHEN Yong, HUANG Ting-fang, LIU En-ru. 2001. Rock Physics [M]. Peking University Press, Beijing(in Chinese).

    崔效鋒, 謝富仁, 張紅艷. 2006. 川滇地區(qū)現(xiàn)代構(gòu)造應(yīng)力場分區(qū)及動(dòng)力學(xué)意義 [J]. 地震學(xué)報(bào), 28(5): 451— 461.

    CUI Xiao-feng, XIE Fu-ren, ZHANG Hong-yan. 2006. Recent tectonic stress field zoning in Sichuan-Yunnan region and its dynamic interest [J]. Acta Seismologica Sinica, 28(5): 451— 461(in Chinese).

    何昌榮, Verberne B A, Spiers C J. 2011. 龍門山斷裂帶沉積巖和天然斷層泥的摩擦滑動(dòng)性質(zhì)與啟示 [J]. 巖石力學(xué)與工程學(xué)報(bào), 30(1): 113—131.

    HE Chang-rong, Verberne B A, Spiers C J. 2011. Frictional properties of sedimentary rocks and natural fault gouge from Longmenshan fault zone and their implications [J]. Chinese Journal of Rock Mechanics and Engineering, 30(1): 113—131(in Chinese).

    單斌, 熊熊, 鄭勇, 等. 2009. 2008年5月12日MW7.9汶川地震導(dǎo)致的周邊斷層應(yīng)力變化 [J]. 中國科學(xué)(D輯), 39(5): 537—545.

    SHAN Bin, XIONG Xiong, ZHENG Yong,etal. 2009. Stress changes on major faults caused by MW7.9 Whenchuan earthquakes, May 12, 2008 [J]. Sci China (Ser D), 39(5): 537—545(in Chinese).

    石耀霖. 2004. 地應(yīng)力主應(yīng)力的方位角求和與平均 [J]. 地震學(xué)報(bào), 26(1): 106—109.

    SHI Yao-lin. 2004. Summation and decomposition of principal stresses in the crust [J]. Acta Seismologica Sinica, 26(1): 106—109(in Chinese).

    王艷華, 崔效鋒, 胡幸平, 等. 2012. 基于原地應(yīng)力測量數(shù)據(jù)的中國大陸地殼上部應(yīng)力狀態(tài)研究 [J]. 地球物理學(xué)報(bào), 55(9): 3016—3027.

    WANG Yan-hua, CUI Xiao-feng, HU Xing-ping,etal. 2012. Study on the stress state in upper crust of China mainland based on in-situ stress measurements [J]. Chinese J Geophys, 55(9): 3016—3027(in Chinese).

    謝富仁, 崔效鋒, 趙建濤. 2003. 全球應(yīng)力場與構(gòu)造分析 [J]. 地學(xué)前緣, 10(增刊): 22—30.

    XIE Fu-ren, CUI Xiao-feng, ZHAO Jian-tao. 2003. Analysis of global tectonic stress field [J]. Earth Science Frontiers, 10(Suppl): 22—30(in Chinese).

    謝富仁, 蘇剛, 崔效鋒, 等. 2001. 滇西南地區(qū)現(xiàn)代構(gòu)造應(yīng)力場分析 [J]. 地震學(xué)報(bào), 23(1): 17—23.

    XIE Fu-ren, SU Gang, CUI Xiao-feng,etal. 2001. Modern tectonic stress field in southwestern Yunnan, China [J]. Acta Seismologica Sinica, 23(1): 17—23(in Chinese).

    謝富仁, 祝景忠, 粱海慶, 等. 1993. 中國西南地區(qū)現(xiàn)代構(gòu)造應(yīng)力場基本特征 [J]. 地震學(xué)報(bào), 15(4): 407— 417.

    XIE Fu-ren, ZHU Jing-zhong, LIANG Hai-qing,etal. 1993. Modern tectonic stress field in southwestern China [J]. Acta Seismologica Sinica, 15(4): 407— 417(in Chinese).

    許忠淮. 1985. 用滑動(dòng)方向擬合法反演唐山余震區(qū)的平均應(yīng)力場 [J]. 地震學(xué)報(bào), 7(4): 349—361.

    XU Zhong-huai. 1985. Mean stress filed in Tangshan aftershock area obtained from focal mechanism data by fitting slip directions [J]. Acta Seismologica Sinica, 7(4): 349—361(in Chinese).

    張浪平, 邵志剛, 李志海. 2014. 印度板塊與歐亞板塊在興都庫什—帕米爾地區(qū)相互俯沖的動(dòng)力作用分析 [J]. 地球物理學(xué)報(bào), 57(2): 459— 471.

    ZHANG Lang-ping, SHAO Zhi-gang, LI Zhi-hai. 2014. Dynamic action of mutual subduction between the Indian and the Eurasia plates in Hindu Kush-Pamir region [J]. Chinese J Geophys, 57(2): 459— 471(in Chinese).

    趙建濤, 崔效鋒, 謝富仁. 2002. 唐山地震震源區(qū)構(gòu)造應(yīng)力場強(qiáng)度的初步分析 [J]. 地震學(xué)報(bào), 24(3): 268—276.

    ZHAO Jian-tao, CUI Xiao-feng, XIE Fu-ren. 2002. Preliminary analysis of the tectonic stress intensity in the source region of Tangshan earthquake [J]. Acta Seismologica Sinica, 24(3): 268—276(in Chinese).

    Angelier J. 1979. Determination of the mean principal direction of stresses for a given fault population [J]. Tectonophysics, 56: 17—26.

    Byerlee J D. 1967. Frictional characteristics of granite under high confining pressure [J]. J Geophys Res, 72: 3639—3648.

    Byerlee J D. 1978. Friction of rocks [J]. Pure Appl Geophys, 116: 615— 626.

    Gephart J W, Forsyth D W. 1984. An improved method for determining the regional stress tensor using earthquake focal mechanism data: Application to the San Fernando earthquake sequence [J]. J Geophys Res, 89: 9305—9320.

    Hardebeck J L, Michael A J. 2006. Damped regional-scale stress inversions: Methodology and examples for southern California and the Coalinga aftershock sequence [J]. J Geophys Res, 111: B11310. doi: 10.1029/2005JB004144.

    Jaeger J C, Cook N G W. 1979. Fundamentals of Rock Mechanics [M]. Chapman and Hall, London.

    Xu Junshan, Wang Jianxin, Xu Congcong. 2014. The estimation of the maximum horizontal differential stress and its depth in the local lithosphere [C]∥Proceedings of the 48th US Rock Mechanics/Geomechanics Symposium, 2014, 2. 948—952.

    Kohlstedt D L, Brian Evans, Mackwell S J. 1995. Strength of the lithosphere: Constraints imposed by laboratory experiments [J]. J Geophysics Res, 100(B9): 17587—17602.

    Plenefisch T, Bonjer K P. 1997. The stress field in the Rhine Graben area inferred from earthquake focal mechanisms and estimation of frictional parameters [J]. Tectonophysics, 275: 71—97.

    Zoback M L. 1992a. First and second order patterns of stress in the lithosphere: The world stress map project [J]. J Geophysics Res, 97(B8): 11703—11728.

    Zoback M L. 1992b. Stress field constrains on intraplate seismicity in eastern North America [J]. J Geophysics Res, 97: 11761—11782.

    Abstract

    Based on Zoback’s method for estimating the tectonic stress magnitude and the two assumptions, we consider the conditions that three principal stresses are vertical principal stresses respectively(corresponding to three kinds of tectonic stress types). We deduced the formulae for estimating the tectonic stress magnitude by using the stress form factor and frictional strength of the fault and discussed the correlative influence of friction coefficient, pore pressure parameter and stress form factor on the stress value. When the maximum principal stress is approximately horizontal (when stress regime is strike-slip or reverse), the maximum principal stress (or the slope of stress increasing linearly with depth) is positively related with the friction coefficient and negatively related with the pore pressure coefficient. When the minimum principal stress is approximately horizontal (when stress regime is strike-slip or normal), the minimum principal stress (or the slope with depth) is negatively related to the friction coefficient, and positive to the pore pressure. Besides, these three parameters have great influence on the estimation of the tectonic stress magnitude. If the friction coefficient is too big and the pore pressure is too small, there could be a wide difference between the slope of the maximum principal stress increasing with depth and the slope of the minimum principal stress increasing with depth, which could lead to an unreasonable result. Our method is just an approximate estimation for the tectonic stress magnitude when crustal rocks have undergone brittle rupture or frictional sliding. The estimated results are not the tectonic stress magnitude when crust is in steady state.

    ESTIMATING THE MAGNITUDE OF TECTONIC STRESS BASED ON THE FRICTION CRITERIA OF FAULT AND ANALYSING THE PARAMETERS’ INFLUENCE

    CAO Hui-jingCUI Xiao-fengFAN Wen-jie

    (KeyLabofCrustalDynamics,InstituteofCrustalDynamics,ChinaEarthquakeAdministration,Beijing100085,China)

    stress magnitude, friction coefficient, pore pressure, stress form factor, tectonic stress

    2015-05-26收稿, 2016-02-23改回。

    國家科技支撐計(jì)劃課題(2012BAK19B03)、 深部探測技術(shù)與實(shí)驗(yàn)研究專項(xiàng)(SinoProbe0604)和中央級(jí)公益性科研院所基本科研業(yè)務(wù)專項(xiàng)(ZDJ201414)共同資助。
    *

    崔效鋒, 男, 研究員, E-mail: cuixfeng@sina.com。

    P315.2

    A

    0253-4967(2016)02-386-11

    曹慧靜, 女, 1989年生, 碩士研究生, 主要從事地殼應(yīng)力場研究, E-mail: doreen_jing@hotmail.com。

    doi:10.3969/j.issn.0253- 4967.2016.02.012

    免费无遮挡裸体视频| 日本 av在线| 国产蜜桃级精品一区二区三区| 亚洲av日韩精品久久久久久密| 午夜免费激情av| 亚洲国产中文字幕在线视频| 午夜日韩欧美国产| 国产一区二区三区视频了| 久久天堂一区二区三区四区| av超薄肉色丝袜交足视频| 老司机午夜福利在线观看视频| netflix在线观看网站| 亚洲中文字幕一区二区三区有码在线看 | 亚洲成人精品中文字幕电影| 国产精品爽爽va在线观看网站 | 国产单亲对白刺激| 波多野结衣巨乳人妻| 麻豆久久精品国产亚洲av| 1024香蕉在线观看| 一级片免费观看大全| 最新在线观看一区二区三区| 在线观看舔阴道视频| 亚洲全国av大片| 99热只有精品国产| 一边摸一边抽搐一进一小说| 色尼玛亚洲综合影院| 日韩精品中文字幕看吧| 亚洲国产毛片av蜜桃av| 欧美日韩亚洲综合一区二区三区_| bbb黄色大片| 国产精品电影一区二区三区| 12—13女人毛片做爰片一| 欧美黑人精品巨大| 色播亚洲综合网| 极品人妻少妇av视频| 日日爽夜夜爽网站| 午夜福利高清视频| 亚洲无线在线观看| 久久这里只有精品19| 中文字幕av电影在线播放| 国产99白浆流出| 精品乱码久久久久久99久播| 禁无遮挡网站| 最新美女视频免费是黄的| 黄色a级毛片大全视频| 一边摸一边抽搐一进一出视频| 免费少妇av软件| 欧美激情高清一区二区三区| 性欧美人与动物交配| 无限看片的www在线观看| 涩涩av久久男人的天堂| 久久中文字幕人妻熟女| 国产精品亚洲av一区麻豆| 一本综合久久免费| 美女大奶头视频| 精品免费久久久久久久清纯| 精品欧美一区二区三区在线| 在线十欧美十亚洲十日本专区| 91成人精品电影| 国产精品精品国产色婷婷| 操美女的视频在线观看| 不卡av一区二区三区| 欧美激情极品国产一区二区三区| 亚洲精品久久国产高清桃花| 美女扒开内裤让男人捅视频| 99国产精品99久久久久| 一夜夜www| 午夜免费观看网址| 在线观看免费视频网站a站| 国产私拍福利视频在线观看| 男女做爰动态图高潮gif福利片 | 露出奶头的视频| 久久天躁狠狠躁夜夜2o2o| 女警被强在线播放| 在线播放国产精品三级| 亚洲午夜理论影院| 日本黄色视频三级网站网址| 村上凉子中文字幕在线| 亚洲第一青青草原| 免费看美女性在线毛片视频| 欧美乱妇无乱码| www.精华液| 色婷婷久久久亚洲欧美| 一区二区日韩欧美中文字幕| 免费看十八禁软件| 9色porny在线观看| 精品电影一区二区在线| 亚洲av日韩精品久久久久久密| 老司机午夜十八禁免费视频| 中文字幕精品免费在线观看视频| 国产一区在线观看成人免费| 久久亚洲真实| xxx96com| 成人av一区二区三区在线看| 久久久国产精品麻豆| 亚洲欧美激情综合另类| 欧美性长视频在线观看| 国产av一区二区精品久久| 一本久久中文字幕| 国产激情欧美一区二区| 亚洲少妇的诱惑av| 两个人视频免费观看高清| 国产激情久久老熟女| 成在线人永久免费视频| 黄色a级毛片大全视频| 久久伊人香网站| 国产日韩一区二区三区精品不卡| 一区二区三区国产精品乱码| 女同久久另类99精品国产91| 国产成人系列免费观看| 美女高潮到喷水免费观看| 丝袜人妻中文字幕| 精品久久蜜臀av无| 91精品国产国语对白视频| 又紧又爽又黄一区二区| 人人妻人人澡欧美一区二区 | 国产精品亚洲av一区麻豆| 欧美最黄视频在线播放免费| 在线观看免费视频日本深夜| 欧美日本中文国产一区发布| 亚洲欧美激情综合另类| 手机成人av网站| 亚洲激情在线av| 香蕉国产在线看| 一级毛片高清免费大全| 久久国产精品男人的天堂亚洲| 久久天躁狠狠躁夜夜2o2o| 我的亚洲天堂| 亚洲少妇的诱惑av| 69av精品久久久久久| 伦理电影免费视频| 日韩免费av在线播放| 日韩欧美免费精品| 国产片内射在线| 国产成人欧美| 久久人妻福利社区极品人妻图片| √禁漫天堂资源中文www| 激情视频va一区二区三区| 天堂√8在线中文| 亚洲成a人片在线一区二区| 久久精品亚洲熟妇少妇任你| 国产精品免费一区二区三区在线| bbb黄色大片| 久久精品亚洲熟妇少妇任你| 美女高潮喷水抽搐中文字幕| 99久久精品国产亚洲精品| 又黄又爽又免费观看的视频| 他把我摸到了高潮在线观看| 国产极品粉嫩免费观看在线| 真人做人爱边吃奶动态| 亚洲第一av免费看| 香蕉国产在线看| 搡老熟女国产l中国老女人| 男女下面进入的视频免费午夜 | 97超级碰碰碰精品色视频在线观看| 天天一区二区日本电影三级 | 久久国产亚洲av麻豆专区| 涩涩av久久男人的天堂| 又大又爽又粗| 日日爽夜夜爽网站| 99re在线观看精品视频| 亚洲成av人片免费观看| 国产亚洲精品一区二区www| 久久亚洲精品不卡| 欧美成人性av电影在线观看| 久久午夜亚洲精品久久| 午夜成年电影在线免费观看| 一区二区三区激情视频| 国产一区二区三区在线臀色熟女| 多毛熟女@视频| 悠悠久久av| 免费看美女性在线毛片视频| 中出人妻视频一区二区| 日韩国内少妇激情av| 日本黄色视频三级网站网址| 国产又爽黄色视频| 亚洲九九香蕉| 97超级碰碰碰精品色视频在线观看| 黄色a级毛片大全视频| 亚洲一区中文字幕在线| 无遮挡黄片免费观看| 国产精品av久久久久免费| 黄色片一级片一级黄色片| 国产精品电影一区二区三区| 国产一区二区激情短视频| 久久久久久久精品吃奶| 日韩欧美国产在线观看| 亚洲精品一区av在线观看| 最近最新免费中文字幕在线| 9191精品国产免费久久| 国产主播在线观看一区二区| 久9热在线精品视频| 天天一区二区日本电影三级 | 一区二区三区国产精品乱码| avwww免费| 欧美日韩一级在线毛片| 亚洲一区中文字幕在线| 女性生殖器流出的白浆| 国产精品 欧美亚洲| 18禁观看日本| 国产精品 国内视频| 可以在线观看毛片的网站| 巨乳人妻的诱惑在线观看| 国产av又大| 中文字幕精品免费在线观看视频| 亚洲精品美女久久久久99蜜臀| 成人国语在线视频| 中文亚洲av片在线观看爽| 美女高潮到喷水免费观看| 国产成人啪精品午夜网站| 好男人电影高清在线观看| 亚洲aⅴ乱码一区二区在线播放 | 国产精品一区二区免费欧美| 欧美在线黄色| 日本一区二区免费在线视频| 在线观看免费视频日本深夜| 午夜亚洲福利在线播放| 欧美最黄视频在线播放免费| 国产高清激情床上av| 久久久精品国产亚洲av高清涩受| 亚洲欧美日韩另类电影网站| 国产av精品麻豆| 黑人欧美特级aaaaaa片| 国产成人系列免费观看| 国产极品粉嫩免费观看在线| 国产视频一区二区在线看| 嫩草影院精品99| 伦理电影免费视频| 97人妻精品一区二区三区麻豆 | 国产极品粉嫩免费观看在线| 精品不卡国产一区二区三区| 欧美日韩亚洲综合一区二区三区_| 亚洲情色 制服丝袜| 色在线成人网| 午夜精品久久久久久毛片777| 久久久久国产精品人妻aⅴ院| 亚洲熟妇熟女久久| 亚洲国产精品久久男人天堂| 欧美日本亚洲视频在线播放| 波多野结衣一区麻豆| 亚洲国产精品合色在线| 免费高清在线观看日韩| 91老司机精品| 两个人免费观看高清视频| 免费不卡黄色视频| 老汉色∧v一级毛片| 69av精品久久久久久| 在线国产一区二区在线| 99精品欧美一区二区三区四区| 成熟少妇高潮喷水视频| av中文乱码字幕在线| 欧美成狂野欧美在线观看| 色婷婷久久久亚洲欧美| 亚洲性夜色夜夜综合| 人妻丰满熟妇av一区二区三区| 男人舔女人下体高潮全视频| 女人精品久久久久毛片| 999久久久国产精品视频| 久久精品91无色码中文字幕| 在线十欧美十亚洲十日本专区| 美女高潮到喷水免费观看| 中文字幕最新亚洲高清| 丁香欧美五月| e午夜精品久久久久久久| 悠悠久久av| 午夜福利在线观看吧| 欧美乱色亚洲激情| 亚洲av电影在线进入| 人人澡人人妻人| 久久狼人影院| 十分钟在线观看高清视频www| 老司机靠b影院| 麻豆av在线久日| 两性午夜刺激爽爽歪歪视频在线观看 | 国产乱人伦免费视频| 天天添夜夜摸| 国产精品电影一区二区三区| 深夜精品福利| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美日韩无卡精品| 精品国产一区二区久久| 精品国产亚洲在线| 麻豆av在线久日| 少妇熟女aⅴ在线视频| 一个人免费在线观看的高清视频| 欧美最黄视频在线播放免费| 人人妻人人爽人人添夜夜欢视频| 国产人伦9x9x在线观看| 日韩欧美一区视频在线观看| 少妇的丰满在线观看| 亚洲专区中文字幕在线| 亚洲色图 男人天堂 中文字幕| 欧美日韩精品网址| 高清黄色对白视频在线免费看| netflix在线观看网站| 老汉色av国产亚洲站长工具| xxx96com| 国产精品免费一区二区三区在线| 高清在线国产一区| 亚洲片人在线观看| 精品日产1卡2卡| 人人妻,人人澡人人爽秒播| 欧美色视频一区免费| 久久香蕉精品热| 99国产精品一区二区蜜桃av| 亚洲最大成人中文| 91av网站免费观看| 亚洲一区二区三区不卡视频| 黑丝袜美女国产一区| 啦啦啦观看免费观看视频高清 | 精品卡一卡二卡四卡免费| 久久精品影院6| 最新美女视频免费是黄的| 91字幕亚洲| 露出奶头的视频| 国产人伦9x9x在线观看| 国产三级在线视频| 99香蕉大伊视频| 欧美精品亚洲一区二区| 变态另类丝袜制服| 国产精品 欧美亚洲| 欧美日韩瑟瑟在线播放| 日韩 欧美 亚洲 中文字幕| a级毛片在线看网站| 啦啦啦韩国在线观看视频| 村上凉子中文字幕在线| 亚洲五月天丁香| 窝窝影院91人妻| 亚洲国产看品久久| 真人一进一出gif抽搐免费| 国产亚洲精品一区二区www| 欧美午夜高清在线| 精品一区二区三区四区五区乱码| 777久久人妻少妇嫩草av网站| 国产亚洲欧美在线一区二区| 首页视频小说图片口味搜索| 一边摸一边做爽爽视频免费| 一个人免费在线观看的高清视频| 成人欧美大片| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品1区2区在线观看.| 啦啦啦 在线观看视频| 成人三级做爰电影| 久久婷婷成人综合色麻豆| 久久久久久久精品吃奶| 啦啦啦观看免费观看视频高清 | 日韩一卡2卡3卡4卡2021年| 亚洲五月色婷婷综合| 男女之事视频高清在线观看| 精品高清国产在线一区| 99精品欧美一区二区三区四区| 午夜久久久在线观看| 午夜免费成人在线视频| 久久九九热精品免费| 欧美日韩瑟瑟在线播放| 国产精品一区二区三区四区久久 | 在线观看日韩欧美| 欧美不卡视频在线免费观看 | 中文字幕最新亚洲高清| 97碰自拍视频| www.www免费av| 亚洲国产日韩欧美精品在线观看 | 午夜精品国产一区二区电影| 欧美成人性av电影在线观看| 一级毛片精品| 精品人妻在线不人妻| 国产精品秋霞免费鲁丝片| 女同久久另类99精品国产91| 亚洲狠狠婷婷综合久久图片| svipshipincom国产片| 日本免费一区二区三区高清不卡 | 老司机靠b影院| 99在线视频只有这里精品首页| 丝袜在线中文字幕| 国产成人av教育| 日韩欧美国产一区二区入口| 在线av久久热| 一a级毛片在线观看| 欧美另类亚洲清纯唯美| 一进一出好大好爽视频| 精品一区二区三区av网在线观看| 一二三四社区在线视频社区8| 一级毛片女人18水好多| 午夜福利免费观看在线| 国产成+人综合+亚洲专区| 精品国产乱码久久久久久男人| 久久精品亚洲精品国产色婷小说| 国产成人精品无人区| 成年版毛片免费区| 嫁个100分男人电影在线观看| 亚洲中文字幕日韩| 午夜精品久久久久久毛片777| 午夜精品国产一区二区电影| av片东京热男人的天堂| 啦啦啦免费观看视频1| 国产精品久久电影中文字幕| 久久中文看片网| 高潮久久久久久久久久久不卡| 久久久水蜜桃国产精品网| 丁香六月欧美| 欧美日韩亚洲综合一区二区三区_| 欧美绝顶高潮抽搐喷水| 99国产精品一区二区蜜桃av| 国产视频一区二区在线看| 91国产中文字幕| 999久久久精品免费观看国产| 欧美人与性动交α欧美精品济南到| 999精品在线视频| 精品熟女少妇八av免费久了| 神马国产精品三级电影在线观看 | 日韩有码中文字幕| 国产亚洲精品第一综合不卡| 亚洲精品国产一区二区精华液| 国产av一区二区精品久久| 制服人妻中文乱码| 亚洲精品国产精品久久久不卡| 午夜日韩欧美国产| 最近最新中文字幕大全免费视频| 少妇 在线观看| 桃红色精品国产亚洲av| 999久久久国产精品视频| 亚洲欧美激情综合另类| 精品国产一区二区三区四区第35| 欧美大码av| 亚洲熟妇熟女久久| 一本综合久久免费| a级毛片在线看网站| 亚洲 欧美 日韩 在线 免费| 中文字幕人妻熟女乱码| 亚洲成人精品中文字幕电影| 国产精品永久免费网站| 长腿黑丝高跟| 美国免费a级毛片| 99国产精品一区二区蜜桃av| 精品少妇一区二区三区视频日本电影| 99在线视频只有这里精品首页| 日韩精品青青久久久久久| 首页视频小说图片口味搜索| 黄色女人牲交| 欧美最黄视频在线播放免费| 久久人妻福利社区极品人妻图片| 国产精品一区二区精品视频观看| 国产av又大| 神马国产精品三级电影在线观看 | 欧美色欧美亚洲另类二区 | 两个人看的免费小视频| 亚洲无线在线观看| 欧美激情久久久久久爽电影 | 韩国精品一区二区三区| 99久久综合精品五月天人人| 天天躁夜夜躁狠狠躁躁| 欧美在线一区亚洲| 欧美日本视频| 成人免费观看视频高清| 看免费av毛片| 国产精品一区二区三区四区久久 | 亚洲av成人一区二区三| cao死你这个sao货| 亚洲 欧美 日韩 在线 免费| 国内精品久久久久久久电影| 亚洲精品久久成人aⅴ小说| 国产精品一区二区免费欧美| 嫩草影院精品99| 一边摸一边抽搐一进一出视频| 一本久久中文字幕| 久久婷婷成人综合色麻豆| 日本撒尿小便嘘嘘汇集6| 男女下面进入的视频免费午夜 | 精品久久久精品久久久| 日本 欧美在线| 身体一侧抽搐| 国产97色在线日韩免费| 久久天堂一区二区三区四区| 波多野结衣巨乳人妻| 精品免费久久久久久久清纯| 成人特级黄色片久久久久久久| 黄色片一级片一级黄色片| 日本欧美视频一区| 人人澡人人妻人| 亚洲性夜色夜夜综合| 欧美日韩精品网址| 别揉我奶头~嗯~啊~动态视频| 久久人人精品亚洲av| www.熟女人妻精品国产| av有码第一页| 曰老女人黄片| 午夜精品久久久久久毛片777| 精品久久久久久,| 级片在线观看| 亚洲人成77777在线视频| 9热在线视频观看99| 精品欧美国产一区二区三| 午夜免费成人在线视频| 国产片内射在线| 成人免费观看视频高清| 精品福利观看| av福利片在线| 国产又色又爽无遮挡免费看| 精品国产乱子伦一区二区三区| 丝袜美腿诱惑在线| 亚洲熟女毛片儿| 极品人妻少妇av视频| 婷婷六月久久综合丁香| 免费一级毛片在线播放高清视频 | 亚洲一区二区三区色噜噜| 女人精品久久久久毛片| 黄片小视频在线播放| 在线观看舔阴道视频| 中文字幕高清在线视频| 国产精品 国内视频| 午夜福利高清视频| 精品午夜福利视频在线观看一区| 88av欧美| 超碰成人久久| 男人操女人黄网站| 国产真人三级小视频在线观看| 亚洲在线自拍视频| 搡老妇女老女人老熟妇| 不卡一级毛片| 老司机午夜福利在线观看视频| 人人妻人人爽人人添夜夜欢视频| 亚洲专区字幕在线| 青草久久国产| 亚洲第一av免费看| 久久青草综合色| 亚洲精品中文字幕一二三四区| 欧美日韩黄片免| 成人国语在线视频| 免费在线观看完整版高清| 91九色精品人成在线观看| 丝袜人妻中文字幕| 午夜影院日韩av| 两个人看的免费小视频| 免费女性裸体啪啪无遮挡网站| 国产99白浆流出| 老司机午夜福利在线观看视频| 天堂动漫精品| 色老头精品视频在线观看| 99香蕉大伊视频| 国产精品香港三级国产av潘金莲| 国产成人免费无遮挡视频| 男女下面插进去视频免费观看| 亚洲精品中文字幕在线视频| 欧美成人免费av一区二区三区| 99精品欧美一区二区三区四区| 搞女人的毛片| 69精品国产乱码久久久| 精品一品国产午夜福利视频| 波多野结衣巨乳人妻| 久久天躁狠狠躁夜夜2o2o| 久久久久亚洲av毛片大全| 中文字幕最新亚洲高清| 国内毛片毛片毛片毛片毛片| 一区二区三区激情视频| 久久草成人影院| 1024香蕉在线观看| 国产精品亚洲美女久久久| 久久午夜亚洲精品久久| 亚洲国产看品久久| 99国产精品99久久久久| 又大又爽又粗| 人人妻人人澡欧美一区二区 | 久久中文字幕人妻熟女| 久久久久久免费高清国产稀缺| 精品欧美一区二区三区在线| 久久香蕉国产精品| 麻豆成人av在线观看| 十分钟在线观看高清视频www| 中文字幕人妻丝袜一区二区| 欧美不卡视频在线免费观看 | 麻豆国产av国片精品| 精品一品国产午夜福利视频| 午夜日韩欧美国产| 午夜福利视频1000在线观看 | 国产成人av激情在线播放| 日日干狠狠操夜夜爽| 亚洲 国产 在线| 久久午夜综合久久蜜桃| 51午夜福利影视在线观看| 国产精品综合久久久久久久免费 | 免费少妇av软件| 免费一级毛片在线播放高清视频 | 日韩精品青青久久久久久| av天堂久久9| 亚洲av第一区精品v没综合| 午夜精品国产一区二区电影| 亚洲精品国产色婷婷电影| 欧美激情高清一区二区三区| 少妇粗大呻吟视频| 久久久久精品国产欧美久久久| 午夜福利免费观看在线| 日韩中文字幕欧美一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 成人免费观看视频高清| 夜夜看夜夜爽夜夜摸| 午夜福利成人在线免费观看| 国产精品爽爽va在线观看网站 | 国产亚洲精品第一综合不卡| 中文字幕久久专区| 色精品久久人妻99蜜桃| 亚洲欧美日韩无卡精品| 国产成人免费无遮挡视频| 十八禁人妻一区二区| 日本欧美视频一区| 国产成人精品久久二区二区免费| 久久性视频一级片| 日本 av在线| 午夜福利影视在线免费观看| 久久香蕉激情| 午夜免费观看网址| 女性被躁到高潮视频| 美女 人体艺术 gogo| 日本精品一区二区三区蜜桃| 国语自产精品视频在线第100页| av片东京热男人的天堂|