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

    渤海環(huán)境條件對FPSO運動及系泊響應(yīng)的影響研究

    2017-05-10 08:11:49李牧
    資源節(jié)約與環(huán)保 2017年2期
    關(guān)鍵詞:波高系泊監(jiān)測數(shù)據(jù)

    李牧

    (中海油能源發(fā)展采油服務(wù)公司天津300451)

    渤海環(huán)境條件對FPSO運動及系泊響應(yīng)的影響研究

    李牧

    (中海油能源發(fā)展采油服務(wù)公司天津300451)

    環(huán)境因素是影響FPSO船體運動的主導(dǎo)因素,因此掌握環(huán)境因素對船體運動影響的基本規(guī)律,有利于作業(yè)者根據(jù)環(huán)境條件對船體運動情況做出大致判斷,合理選擇不同環(huán)境條件下的裝配載措施。本文主要運用統(tǒng)計分析方法對渤海曹妃甸作業(yè)區(qū)112FPSO現(xiàn)場監(jiān)測數(shù)據(jù)中的環(huán)境監(jiān)測數(shù)據(jù)和船體運動監(jiān)測數(shù)據(jù)進行研究,總結(jié)得到環(huán)境數(shù)據(jù)對船體運動、系泊響應(yīng)之間的統(tǒng)計回歸關(guān)系,從而為海上進行FPSO運動及系泊力快速預(yù)報提供理論支持依據(jù)。

    FPSO運動;系泊響應(yīng);監(jiān)測數(shù)據(jù);統(tǒng)計分析;回歸關(guān)系

    1 引言

    目前國內(nèi)所應(yīng)用的浮式生產(chǎn)儲油裝置(FPSO)均采用單點系泊系統(tǒng),此種系泊系統(tǒng)可使得FPSO在海上作業(yè)時具備風(fēng)向標(biāo)效應(yīng),使其受到最小的系泊力。由于FPSO處于海上風(fēng)浪多變的復(fù)雜環(huán)境之中,自身存儲原油載重變化不斷,加上FPSO連接外輸油船對船體本身產(chǎn)生拉力,因此FPSO軟剛臂系泊裝置受力異常復(fù)雜[1-2]。近年來,海洋環(huán)境條件愈加惡劣,超強臺風(fēng)頻頻出現(xiàn),頻發(fā)的單點事故導(dǎo)致了嚴(yán)重的經(jīng)濟損失及社會影響。對FPSO軟剛臂系泊系統(tǒng)進行實時、全面、長期的現(xiàn)場監(jiān)測,實時掌握單點系泊系統(tǒng)的狀態(tài),是對系泊結(jié)構(gòu)進行安全評估最有效的途徑[3]。能夠為作業(yè)者提供及時準(zhǔn)確的系泊安全性信息,及時對單點潛在危險進行預(yù)警,確保單點系泊系統(tǒng)的安全生產(chǎn)。經(jīng)過監(jiān)測系統(tǒng)一段時間的運行,已經(jīng)記錄了一定的實測數(shù)據(jù),這些數(shù)據(jù)的分析處理作為系泊在線監(jiān)測系統(tǒng)的延續(xù)成為亟待研究的問題。

    本文目標(biāo)監(jiān)測數(shù)據(jù)來源于位于渤海曹妃甸11-1/2油田的海洋石油112現(xiàn)場監(jiān)測系統(tǒng)。監(jiān)測內(nèi)容主要包括:單點受力測量、環(huán)境監(jiān)測和FPSO運動監(jiān)測三個部分。通過對系泊監(jiān)測數(shù)據(jù)的處理,研究并掌握實測數(shù)據(jù)的分析方法,研究環(huán)境對FPSO船體運動情況的影響,挖掘數(shù)據(jù)潛在的意義,并對日常生產(chǎn)及應(yīng)急管理提供合理化建議,也可為上層決策提供科學(xué)依據(jù)。

    2 理論基礎(chǔ)

    2.1 數(shù)據(jù)分析方法

    海洋工程結(jié)構(gòu)物在服役過程中,將不間斷地遭受到各種條件海洋環(huán)境的作用。因此,海洋環(huán)境是影響海洋工程結(jié)構(gòu)的安全服役的主要因素之一,沒有對相應(yīng)位置的特殊環(huán)境條件的清晰認(rèn)識,海洋工程結(jié)構(gòu)物設(shè)計將無從談起。通常來講影響海洋工程結(jié)構(gòu)物設(shè)計的環(huán)境因素包括:風(fēng)、浪、流等,在某些特定的區(qū)域,還應(yīng)考慮內(nèi)波、地震和海冰的影響。

    2.1.1 風(fēng)設(shè)計值

    風(fēng)設(shè)計值的獲取與波浪設(shè)計值的獲取方法相同。依靠足夠次數(shù)的觀察,建立經(jīng)驗分布函數(shù)?;谶@些函數(shù),修正韋布爾律或其他漸近律。通過外推,得出十年一遇、百年一遇等設(shè)計值。

    風(fēng)速的長期分布一般采用通用的分布函數(shù)或散布圖來表示。除非實測數(shù)據(jù)顯示嚴(yán)重不符,否則在參考高度的處10分鐘平均風(fēng)速的變化情況用韋布爾分布來表示[4]:

    式中,系數(shù)和由實測數(shù)據(jù)確定。

    對應(yīng)于一定回歸周期內(nèi)的最大分速,其概率可近似表示為:

    2.1.2 浪

    海洋結(jié)構(gòu)物在生存期間將經(jīng)歷無數(shù)次波浪的作用,對于海上的油田設(shè)施而言,一般認(rèn)為三小時內(nèi)的海況是不變的,即波浪的平均統(tǒng)計特性不變。每個海況包括一系列波高和周期可變化的波浪。通過對這些波浪應(yīng)用極其精確的方式加以描述,從而可將它與水動力載荷聯(lián)系在一起。目前對于波浪測量數(shù)據(jù)的分析方法主要有:逐波分析和譜分析[4]。

    2.1.3 流

    在某一給定近海位置,流是洋流、潮流和風(fēng)共同作用下表面水團的驅(qū)動效應(yīng)的疊加結(jié)果。流的強度和方向在不同深度的水層中會發(fā)生變化,一般來說,流剖面圖分為十年一遇的速度“剖面圖”、百年一遇的速度“剖面圖”等。

    典型洋流可分為:風(fēng)致流、潮汐流、循環(huán)流、渦流、孤立流、沿岸流。

    2.1.4 相關(guān)性及回歸分析

    測度相關(guān)關(guān)系的方式有三種,即相關(guān)表、相關(guān)圖和相關(guān)系數(shù)。相關(guān)表是一種反映變量之間相關(guān)關(guān)系的統(tǒng)計表。是將某一變量按其取值的大小排列,然后在將與其相關(guān)的另一變量的對應(yīng)值平行排列,便可得到簡單的相關(guān)表。相關(guān)表和相關(guān)圖是研究相關(guān)關(guān)系的直觀工具,它們只能對現(xiàn)象之間存在的相關(guān)關(guān)系的方向、形式和密切程度作大致的判斷,而不能說明其密切程度的大小。因此,需要計算相關(guān)系數(shù)[5]。

    根據(jù)線性相關(guān)變量的多少、分析問題的角度不同,相關(guān)系數(shù)可以分為簡單相關(guān)系數(shù)、偏相關(guān)系數(shù)和復(fù)相關(guān)系數(shù)。

    2.2 基礎(chǔ)數(shù)據(jù)

    本次研究主要針對海洋工程設(shè)施現(xiàn)場監(jiān)測數(shù)據(jù)的后處理方法展開,目標(biāo)監(jiān)測數(shù)據(jù)來源于位于渤海曹妃甸11-1/2油田的海洋石油112現(xiàn)場監(jiān)測系統(tǒng)。監(jiān)測內(nèi)容主要包括:SPM受力測量、環(huán)境監(jiān)測和FPSO運動監(jiān)測三個部分,本報告中用到的監(jiān)測數(shù)據(jù)如下:

    1)該系統(tǒng)環(huán)境數(shù)據(jù)的監(jiān)測內(nèi)容主要包括:

    -風(fēng)速、風(fēng)向-氣溫、濕度和氣壓

    -浪高、周期和波向-剖面流速和流向測量

    本報告研究內(nèi)容主要針對風(fēng)、浪、流三種環(huán)境要素。每種環(huán)境要素均采用不同的采樣頻率。具體監(jiān)測內(nèi)容詳細(xì)說明如下表:

    FPSO系泊力的監(jiān)測采用間接測量的手段,通過監(jiān)測YOKE和系泊鏈角度,利用力學(xué)原理反算出系泊系統(tǒng)受力。各個監(jiān)測要素的數(shù)據(jù)說明如下:

    表2 ~3系泊力監(jiān)測數(shù)據(jù)說明

    表2 ~2船體運動監(jiān)測數(shù)據(jù)說明

    表2 ~1海洋環(huán)境監(jiān)測數(shù)據(jù)說明

    2)FPSO運動和位置監(jiān)測

    -FPSO艏向

    -FPSO運動姿態(tài)和位置

    3)SPM系泊力監(jiān)測

    -系泊錨鏈角度

    -系泊支撐結(jié)構(gòu)應(yīng)變

    針對FPSO運動和位置的監(jiān)測主要包括:船體位置(緯度、經(jīng)度)、艏向、橫搖、縱搖、東向速度、北向速度、高度(垂蕩),各個監(jiān)測要素的數(shù)據(jù)說明如下:

    上述數(shù)據(jù)來自2015年1月至2015年12月系統(tǒng)自動保存的原始數(shù)據(jù),因此在數(shù)據(jù)分析過程中根據(jù)技術(shù)要求進行必要的甄別和篩查工作。

    3環(huán)境條件對船體運動影響規(guī)律

    環(huán)境因素是影響船體運動的主導(dǎo)因素,因此掌握環(huán)境因素對船體運動影響的基本規(guī)律,有利于作業(yè)者根據(jù)環(huán)境條件對船體運動情況做出大致判斷,合理選擇不同環(huán)境條件下的裝、配載措施。

    3.1 艏向

    由于單點系泊系統(tǒng)本身具有風(fēng)向標(biāo)效應(yīng)的特點,隨著環(huán)境條件的變化艏向也隨之發(fā)生變化。本節(jié)主要針對這一變化規(guī)律進行了總結(jié),分別就風(fēng)、浪、流與船體艏向之間的影響規(guī)律進行逐個分析,得出了在不同強度的風(fēng)、浪、流條件下,兩者之間相對夾角之間的關(guān)系。

    圖3 .1 1風(fēng)速&相對風(fēng)向聯(lián)合分布圖

    圖3.1~2風(fēng)速VS相對風(fēng)向變化規(guī)律

    圖3 .1-1直觀展示了不同風(fēng)速下對應(yīng)的相對方向變化情況,由圖可知,大部分工況出現(xiàn)在風(fēng)速10m/s以下,10m/s以上的風(fēng)速出現(xiàn)概率較小。除此之外,不同風(fēng)速下不同相對風(fēng)向都有可能出現(xiàn),但出現(xiàn)概率不同。(注:相對風(fēng)向表示風(fēng)向與艏向之間的夾角)

    圖3 .1~3波高VS相對浪向聯(lián)合分布圖

    圖3 .1~4波高VS相對浪向變化規(guī)律

    艏向的變化受多種環(huán)境因素的影響,因此不同風(fēng)速與相對風(fēng)向的變化并非一一對應(yīng),但兩者之間存在一定的總體變化規(guī)律,如圖3.1~2。隨著風(fēng)速的增大相對風(fēng)向逐步減小,兩者之間存在近似線性的關(guān)系。由于目前現(xiàn)有數(shù)據(jù)中大于12m/s風(fēng)速的工況較少,個別工況不具有代表性,因此在回歸分析過程中作為錯誤數(shù)據(jù)進行處理。圖3.1~3直觀展示了波高對應(yīng)的相對浪向分布情況,由圖可知,在波高較低的工況下,相對浪向的分布較為平均,隨之波浪的增加相對浪向在小夾角出現(xiàn)的概率逐步增大。由圖3.1~4可知,隨著波高的增加相對浪向的角度逐步減小,并大致呈現(xiàn)線性關(guān)系。由擬合曲線第一點工況數(shù)量足夠大但不符合變化規(guī)律可知,波高并不能唯一確定艏向的位置,艏向位置應(yīng)由環(huán)境中較強的要素予以確定。圖3.1~5直觀展示了流速對應(yīng)相對流向的分布情況,由圖可知,相對流向主要集中在180°附近,即流向與艏向為反向。且兩者之間的變化規(guī)律相對明顯,說明流對艏向的影響相對強烈。由圖3.1~6可知,隨著流速的增加相對流向逐步向180°接近并趨于平穩(wěn),整體變化趨勢呈二次變化關(guān)系。結(jié)合圖3.1~5及3.1~6可知,流速對艏向的影響系數(shù)較大,但艏向的最終確定依然是風(fēng)、浪、流等環(huán)境要素聯(lián)合作用的結(jié)果。

    圖3 .1~5流速VS相對流向聯(lián)合分布圖

    圖3 .1~6流速VS相對流向變化規(guī)律

    3.2 橫搖運動

    橫搖是衡量FPSO運動性能的一項重要指標(biāo),橫搖運動過大會直接影響海上的作業(yè)安全,從而造成停產(chǎn),因此有必要分析環(huán)境要素與FPSO橫搖運動之間的影響規(guī)律,以便在生產(chǎn)設(shè)計中作為參考。由于橫搖運動的均值經(jīng)常受FPSO裝載情況影響,因此本小節(jié)不對橫搖運動的均值做重點研究,研究過程中僅考慮波浪與橫搖運動之間的影響規(guī)律。由于波浪相對船體的入射角各不相同,為了考慮入射角對波高的影響,引入了相對波高的概念,即將波高在船體艏向的垂直方向上進行了分解。

    圖3 .2~1波高、周期與橫搖角標(biāo)準(zhǔn)差散布圖(0°~45°)

    在不考慮波浪入射角度影響的情況下,浪向與艏向夾角小于45°時,波高及周期對船體橫搖運動的影響如圖3.2~1所示。

    圖3 .2~2相對波高、周期與橫搖角標(biāo)準(zhǔn)差散布圖(0°~45°)

    考慮波浪入射角的影響,將波高在船體艏向垂直方向上進行分解,得到相對波高、周期對船體橫搖運動的影響圖3.2~2,總體變化規(guī)律與圖3.2~1相似。

    圖3 .2~3波高、周期與橫搖角標(biāo)準(zhǔn)差散布圖(45°~135°)

    在不考慮波浪入射角度影響的情況下,浪向與艏向夾角為45°~135°之間時,波高及周期對船體橫搖運動的影響如圖3.2~4所示。由圖可知,整體規(guī)律及標(biāo)準(zhǔn)差分布情況不具有特別明顯的規(guī)律性。

    圖3 .2~4相對波高、周期與橫搖角標(biāo)準(zhǔn)差散布圖(45°~135°)

    綜上分析,波浪對船體橫搖運動的影響關(guān)系相對復(fù)雜,涉及到波高、周期和船體裝載情況等因素。因此上述分析中規(guī)律不顯著的主要原因有:1.橫搖運動涉及影響因素眾多,目前監(jiān)測數(shù)據(jù)有限無法完成各分情況的討論;2.目前監(jiān)測數(shù)據(jù)中環(huán)境條件較為溫和,橫搖運動相對較小,各個影響因素的影響效果難以區(qū)分。

    3.3 縱搖運動

    縱搖運動的大小直接關(guān)系到系泊鏈張力的變化,且縱搖運動的固有周期處于一般波浪頻率范圍內(nèi),且普遍認(rèn)為是波浪誘導(dǎo)的波頻運動響應(yīng),本小節(jié)主要討論波浪與船體縱搖運動之間的關(guān)系。

    圖3 .3~1波高、周期縱搖角標(biāo)準(zhǔn)差散布圖(0°~45°)

    在不考慮波浪入射角度影響的情況下,浪向與艏向夾角小于45°時,波高和周期對船體縱搖運動的影響如圖3.3~1所示??v搖角隨波高增加而逐步增大的規(guī)律較為明顯,但波周期與縱搖角之間的變化規(guī)律較為模糊,這與波浪周期主要集中在6s以內(nèi),而波周期大于6s的工況數(shù)量有限。

    圖3 .3~2相對波高、周期與縱搖角標(biāo)準(zhǔn)差散布圖(0°~45°)

    在考慮波浪入射角度的影響后,浪向與艏向夾角小于45°時,波高和周期對船體縱搖運動的影響規(guī)律與上圖相同,但與目前監(jiān)測數(shù)據(jù)中環(huán)境條件較為溫和有關(guān),因此入射角對縱搖運動的影響,需要進一步進行討論。

    在不考慮波浪入射角度影響的情況下,浪向與艏向夾角在45°~135°之間時,波高和周期對縱搖角的影響如圖3.3~3。波高和周期對縱搖角的影響有限,但也與環(huán)境條件較為溫和有關(guān),需要進一步進行分析。

    圖3 .3~4相對波高、周期與縱搖角標(biāo)準(zhǔn)差散布圖(45°~135°)

    在考慮波浪入射角度影響的情況下,波向與艏向夾角在45°~135°之間時,波高和周期對縱搖角的影響規(guī)律與上圖相似,說明在溫和的環(huán)境條件下,波浪入射角度對船體縱搖運動的影響有限。

    綜上分析,在溫和環(huán)境條件下,波浪對船體縱搖運動的影響主要發(fā)生在迎浪條件下,且在周期較小時,波浪周期對船體縱搖運動的影響有限。對于波浪入射角的影響還需要進一步進行討論。

    3.4 縱蕩運動

    縱蕩運動是影響系泊系統(tǒng)回復(fù)力的直接因素,掌握環(huán)境條件對船體縱蕩運動的影響規(guī)律即可對照剛度曲線估算系泊系統(tǒng)回復(fù)力。故本小節(jié)分別針對風(fēng)、浪、流對船體縱蕩運動的影響規(guī)律進行了總結(jié)。

    圖3 .4~1風(fēng)速VS縱蕩運動變化規(guī)律

    由圖3.4~1所示,隨著風(fēng)速的增大船體縱蕩運動逐步增大,且縱蕩運動的大小與風(fēng)速的平方大致呈線性關(guān)系。

    圖3 .4~2波高VS縱蕩運動變化規(guī)律

    由圖3.4~2所示,波高與船體縱蕩運動呈現(xiàn)正相關(guān)關(guān)系,船體縱蕩運動隨波高的增大而增大。

    圖3 .4~3流速VS縱蕩運動變化規(guī)律

    由圖3.4~3所示,縱蕩運動流速之間的關(guān)系呈線性關(guān)系,但線性關(guān)系斜率較小,因此流速對船體縱蕩運動的影響較小。綜上分析,船體縱蕩運動主要由風(fēng)速和波高決定,流速對船體縱蕩運動的影響較小。

    4船體運動對系泊響應(yīng)影響規(guī)律

    船體運動與系泊響應(yīng)兩者之間的變化關(guān)系即系泊系統(tǒng)剛度,根據(jù)現(xiàn)場監(jiān)測得到的船體運動數(shù)據(jù)及系泊響應(yīng)數(shù)據(jù),對比兩者之間的數(shù)據(jù),對照剛度曲線即可檢驗監(jiān)測數(shù)據(jù)的可靠性及設(shè)計系泊剛度曲線的合理性。船體距單點中的縱蕩距離是影響系泊力變化最主要的因素,船體其他運動變化只會引起系泊力在小幅范圍內(nèi)變化,因此系泊系統(tǒng)剛度的分析以縱蕩運動為主。

    圖4 ~1系泊鏈張力VS縱蕩散布圖

    圖4 ~2水平恢復(fù)力VS縱蕩運動散布圖

    由圖4~1可知,系泊鏈張力和船體縱蕩運動之間僅存在大致的趨勢關(guān)系,但并不存在明顯的比例關(guān)系,這主要由于系泊鏈張力是隨船體運動變化最為劇烈的要素,且受系泊鏈垂直方向夾角的影響,因此系泊鏈張力的變化在趨勢線周圍分布較為離散,在實際設(shè)計過程常常以水平恢復(fù)力作為標(biāo)準(zhǔn)進行剛度曲線的繪制,由圖4~2可知,目前監(jiān)測得到的縱蕩運動范圍較小,對照剛度曲線的變化趨勢處于線性段的范圍內(nèi),因此利用線性關(guān)系進行擬合,擬合結(jié)果較為理想。綜上所述,船體運動與系泊響應(yīng)之間的關(guān)系主要取決于剛度曲線的變化,在實際校核過程中可利用水平回復(fù)力與縱蕩運動之間的關(guān)系進行校核。

    5 結(jié)論

    根據(jù)現(xiàn)有數(shù)據(jù)的分析結(jié)果,可用于工程實際及可進一步推廣到設(shè)施延壽及可靠性評估等領(lǐng)域的結(jié)論如下:

    1.船體艏向及周向位置長期以兩個方向為主,這就造成單點主軸承各部分的疲勞、磨損等壽命不均,通過對船體位置的詳細(xì)跟蹤記錄,可對單點主軸承各部分的實際疲勞損傷進行評估,達到對單點主軸承的全壽命周期跟蹤管理;

    2.根據(jù)監(jiān)測數(shù)據(jù)得到不同裝載條件的下船體六自由度運動固有周期,通過實測固有周期可以校準(zhǔn)數(shù)值模型,為日后數(shù)值預(yù)報提供可能,及時發(fā)現(xiàn)危險工況,并及時做出應(yīng)急響應(yīng);

    3.運用系泊鏈張力的現(xiàn)場實測數(shù)據(jù)可以真實記錄其受力變化情況,結(jié)合統(tǒng)計分析方法,不僅可以用于系泊鏈實際疲勞累積損傷的評估,而且可以利用系泊鏈張力在不同受力區(qū)間內(nèi)的發(fā)生概率,直接推算系泊鏈張力的受力極值,避免其他間接計算方法造成的誤差累積;

    4.利用環(huán)境條件對船體運動及系泊響應(yīng)之間的變化規(guī)律,可以檢驗監(jiān)測數(shù)據(jù)的一致性,細(xì)化錯誤數(shù)據(jù)的篩查工作,及時發(fā)現(xiàn)監(jiān)測設(shè)備的故障,保證珍貴數(shù)據(jù)的真實有效。

    [1]劉莉峰、趙玉梁.FPSO單點軟剛臂系泊監(jiān)測系統(tǒng)[J].船海工程,2015,44(3):83-86.

    [2]梁修鋒,楊建民,李欣,于洋.FPSO甲板上浪的數(shù)值模擬[J].水動力學(xué)研究與進展,2007,22(2):229-236.

    [3]劉濤.FPSO軟剛臂單點系泊系統(tǒng)監(jiān)測技術(shù)研究[D].大連:大連理工大學(xué),2013.

    [4]劉水庚.海洋工程水動力學(xué)[M].北京:國防工業(yè)出版社,2012.

    [5]章舜仲、王樹梅.相關(guān)系數(shù)矩陣與多元線性相關(guān)分析

    [J].大學(xué)數(shù)學(xué),2011,27(1):195-198.

    李牧1986年天津大學(xué)碩士研究生,工程師。浮體結(jié)構(gòu)專業(yè)方向依托項目泊監(jiān)測數(shù)據(jù)分析系統(tǒng)開發(fā)及應(yīng)用研究。依托項目:系泊監(jiān)測數(shù)據(jù)分析系統(tǒng)開發(fā)及應(yīng)用研究,項目編號:HFXMLZ-CY201403

    猜你喜歡
    波高系泊監(jiān)測數(shù)據(jù)
    基于FHDI-GNWM 數(shù)據(jù)的全球超越概率波高宏觀分布特征分析
    基于CFD方法的系泊狀態(tài)FSO舵載荷計算
    高韌性海洋用三級系泊鏈鋼的研制
    山東冶金(2022年2期)2022-08-08 01:50:50
    基于漂流浮標(biāo)的南大洋衛(wèi)星高度計有效波高研究
    海洋通報(2021年3期)2021-08-14 02:20:46
    非平整港池的多向不規(guī)則波試驗研究
    GSM-R接口監(jiān)測數(shù)據(jù)精確地理化方法及應(yīng)用
    基于懸鏈線方程的系泊系統(tǒng)分析
    GPS異常監(jiān)測數(shù)據(jù)的關(guān)聯(lián)負(fù)選擇分步識別算法
    基于小波函數(shù)對GNSS監(jiān)測數(shù)據(jù)降噪的應(yīng)用研究
    飽和秋色
    女報seaside(2014年10期)2014-04-29 21:02:51
    av视频在线观看入口| 国产一级毛片七仙女欲春2| 黄色丝袜av网址大全| 欧美黑人巨大hd| 婷婷六月久久综合丁香| 欧美日韩瑟瑟在线播放| 女的被弄到高潮叫床怎么办 | 日本 欧美在线| 美女黄网站色视频| 国产亚洲精品久久久久久毛片| 久久国产精品人妻蜜桃| 欧美在线一区亚洲| 少妇被粗大猛烈的视频| 免费av毛片视频| 国产精品久久久久久久电影| 大又大粗又爽又黄少妇毛片口| 亚洲人与动物交配视频| 久久午夜亚洲精品久久| 99热只有精品国产| 久久久精品大字幕| 日韩欧美一区二区三区在线观看| 亚洲精品乱码久久久v下载方式| 97人妻精品一区二区三区麻豆| 女人被狂操c到高潮| 3wmmmm亚洲av在线观看| 成人亚洲精品av一区二区| 搡女人真爽免费视频火全软件 | av.在线天堂| 99riav亚洲国产免费| 女同久久另类99精品国产91| 成人无遮挡网站| 全区人妻精品视频| 久久九九热精品免费| 成年女人看的毛片在线观看| 99久久精品国产国产毛片| 国产成人福利小说| 午夜影院日韩av| .国产精品久久| 国产在线精品亚洲第一网站| 亚洲av免费高清在线观看| 狂野欧美激情性xxxx在线观看| 久久国产乱子免费精品| 亚洲av五月六月丁香网| 国产 一区 欧美 日韩| 午夜精品久久久久久毛片777| 久久久精品欧美日韩精品| 日本色播在线视频| 三级男女做爰猛烈吃奶摸视频| 国产探花在线观看一区二区| 麻豆精品久久久久久蜜桃| 99久久无色码亚洲精品果冻| 亚洲精品色激情综合| 亚洲精品在线观看二区| 又黄又爽又免费观看的视频| 级片在线观看| 在线免费观看不下载黄p国产 | 欧美最黄视频在线播放免费| xxxwww97欧美| 久久午夜福利片| 麻豆国产97在线/欧美| 日韩欧美精品v在线| 国产老妇女一区| 中文字幕高清在线视频| 成人亚洲精品av一区二区| 给我免费播放毛片高清在线观看| 欧美一区二区国产精品久久精品| 日韩欧美精品v在线| 中文亚洲av片在线观看爽| 国产伦精品一区二区三区四那| 欧美成人免费av一区二区三区| 嫁个100分男人电影在线观看| 综合色av麻豆| 久久亚洲精品不卡| av在线蜜桃| 人妻制服诱惑在线中文字幕| 老女人水多毛片| aaaaa片日本免费| 最近最新免费中文字幕在线| 色哟哟·www| 国产爱豆传媒在线观看| 啦啦啦啦在线视频资源| 亚洲真实伦在线观看| 亚洲久久久久久中文字幕| 看黄色毛片网站| 亚洲av一区综合| 高清日韩中文字幕在线| 日本成人三级电影网站| 国产探花在线观看一区二区| 中出人妻视频一区二区| 国产淫片久久久久久久久| 久久6这里有精品| 女生性感内裤真人,穿戴方法视频| 婷婷色综合大香蕉| 中文字幕av成人在线电影| 真人一进一出gif抽搐免费| 91在线精品国自产拍蜜月| 亚洲熟妇熟女久久| 日韩精品有码人妻一区| 人人妻,人人澡人人爽秒播| 精品一区二区三区视频在线观看免费| 国产成人一区二区在线| 啦啦啦观看免费观看视频高清| 成熟少妇高潮喷水视频| bbb黄色大片| 99riav亚洲国产免费| 一夜夜www| 色哟哟·www| 99热只有精品国产| 一区福利在线观看| 自拍偷自拍亚洲精品老妇| 国产一级毛片七仙女欲春2| 看十八女毛片水多多多| 久久久久免费精品人妻一区二区| 淫妇啪啪啪对白视频| 一卡2卡三卡四卡精品乱码亚洲| aaaaa片日本免费| 在线免费十八禁| 一边摸一边抽搐一进一小说| 最近中文字幕高清免费大全6 | 高清毛片免费观看视频网站| 两个人视频免费观看高清| 国产一区二区激情短视频| 99久久中文字幕三级久久日本| 色噜噜av男人的天堂激情| av中文乱码字幕在线| 91久久精品电影网| 女同久久另类99精品国产91| 人妻制服诱惑在线中文字幕| 国产亚洲精品综合一区在线观看| 亚洲av成人精品一区久久| 91在线精品国自产拍蜜月| 十八禁网站免费在线| 久久这里只有精品中国| 国产三级在线视频| 国产高潮美女av| 九九在线视频观看精品| 国产亚洲精品av在线| 亚洲在线观看片| 高清在线国产一区| 女的被弄到高潮叫床怎么办 | 不卡一级毛片| 国产精品久久久久久久电影| 国产乱人伦免费视频| 伊人久久精品亚洲午夜| netflix在线观看网站| 亚洲av中文字字幕乱码综合| 国产色爽女视频免费观看| 国产av一区在线观看免费| 午夜福利成人在线免费观看| 99热这里只有是精品在线观看| 国产麻豆成人av免费视频| 精品一区二区三区视频在线| 男人的好看免费观看在线视频| 亚洲精品影视一区二区三区av| 日本a在线网址| 亚洲av一区综合| 99精品在免费线老司机午夜| 中文在线观看免费www的网站| 欧美另类亚洲清纯唯美| 久久草成人影院| 久久精品国产99精品国产亚洲性色| 久久精品国产亚洲网站| 亚洲成人久久爱视频| 亚洲av免费高清在线观看| 禁无遮挡网站| 亚洲va在线va天堂va国产| 一区二区三区免费毛片| 欧美zozozo另类| 最近中文字幕高清免费大全6 | 日韩 亚洲 欧美在线| 可以在线观看的亚洲视频| 99精品久久久久人妻精品| 国产人妻一区二区三区在| 欧美色视频一区免费| 99精品久久久久人妻精品| 亚洲最大成人手机在线| 一本久久中文字幕| 久久久久九九精品影院| 日韩高清综合在线| 日本一二三区视频观看| 久久国内精品自在自线图片| 午夜福利在线观看吧| 在线观看午夜福利视频| 亚洲精华国产精华精| 精品久久久久久久末码| 嫩草影院入口| 最后的刺客免费高清国语| 桃红色精品国产亚洲av| 久久久久九九精品影院| 免费观看的影片在线观看| 亚洲四区av| 99精品久久久久人妻精品| 亚洲成人久久爱视频| 非洲黑人性xxxx精品又粗又长| 亚洲欧美激情综合另类| 国产69精品久久久久777片| 麻豆av噜噜一区二区三区| 啦啦啦啦在线视频资源| 免费看日本二区| 精品99又大又爽又粗少妇毛片 | 国产一区二区三区av在线 | 成年免费大片在线观看| 免费看av在线观看网站| 精品午夜福利视频在线观看一区| 精品久久久久久成人av| 国产在视频线在精品| 天堂网av新在线| 国产在线男女| 国产精品亚洲美女久久久| 欧美日韩国产亚洲二区| 成人高潮视频无遮挡免费网站| 成人国产综合亚洲| 亚洲人成网站在线播放欧美日韩| 精品一区二区三区视频在线观看免费| 亚洲一区二区三区色噜噜| 在线观看av片永久免费下载| 精品人妻1区二区| 91av网一区二区| 久久婷婷人人爽人人干人人爱| 国产私拍福利视频在线观看| 国产高清视频在线观看网站| 国产一区二区在线av高清观看| 丰满乱子伦码专区| 免费观看人在逋| 中亚洲国语对白在线视频| 给我免费播放毛片高清在线观看| 亚洲人与动物交配视频| 成人鲁丝片一二三区免费| 校园春色视频在线观看| 亚洲第一电影网av| 亚洲在线观看片| 亚洲成av人片在线播放无| 日本熟妇午夜| 亚洲中文字幕日韩| 男女啪啪激烈高潮av片| 午夜影院日韩av| 九九久久精品国产亚洲av麻豆| 亚洲av免费在线观看| 露出奶头的视频| 女人十人毛片免费观看3o分钟| 一边摸一边抽搐一进一小说| 欧美极品一区二区三区四区| 91在线观看av| 日本-黄色视频高清免费观看| 香蕉av资源在线| 国产国拍精品亚洲av在线观看| 一个人观看的视频www高清免费观看| 午夜福利成人在线免费观看| 国产精品1区2区在线观看.| 国产91精品成人一区二区三区| 国产成人影院久久av| 99热这里只有是精品50| 精品一区二区三区人妻视频| 免费高清视频大片| 美女cb高潮喷水在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产日韩欧美精品在线观看| 亚洲美女黄片视频| 欧美性感艳星| 麻豆av噜噜一区二区三区| 久久久久久久久中文| 亚洲专区中文字幕在线| 国产一区二区在线观看日韩| 欧美激情在线99| 成人亚洲精品av一区二区| 黄色丝袜av网址大全| 中文字幕高清在线视频| 男人舔奶头视频| 黄片wwwwww| 欧美zozozo另类| 国产大屁股一区二区在线视频| 免费在线观看影片大全网站| 哪里可以看免费的av片| 给我免费播放毛片高清在线观看| 亚洲人成网站高清观看| 中文字幕高清在线视频| 国产一级毛片七仙女欲春2| 香蕉av资源在线| 色播亚洲综合网| 欧美国产日韩亚洲一区| av黄色大香蕉| 国产久久久一区二区三区| 97碰自拍视频| 国产一区二区三区av在线 | 久久99热6这里只有精品| 18+在线观看网站| 无遮挡黄片免费观看| 国产一区二区激情短视频| 国产美女午夜福利| 国产亚洲精品综合一区在线观看| 91av网一区二区| 高清毛片免费观看视频网站| 成年免费大片在线观看| 精品人妻一区二区三区麻豆 | 国产成人aa在线观看| 最新在线观看一区二区三区| av在线老鸭窝| 国产毛片a区久久久久| 亚洲精品亚洲一区二区| 一个人免费在线观看电影| 少妇猛男粗大的猛烈进出视频 | 欧美zozozo另类| 日本欧美国产在线视频| 欧美高清性xxxxhd video| 搡老熟女国产l中国老女人| 1000部很黄的大片| 丰满的人妻完整版| 日日干狠狠操夜夜爽| 国内精品久久久久久久电影| 人妻丰满熟妇av一区二区三区| 看十八女毛片水多多多| 亚洲美女黄片视频| 琪琪午夜伦伦电影理论片6080| 国模一区二区三区四区视频| 免费看日本二区| 亚洲成av人片在线播放无| 1000部很黄的大片| 我的女老师完整版在线观看| 岛国在线免费视频观看| 日本撒尿小便嘘嘘汇集6| 成人一区二区视频在线观看| 中文亚洲av片在线观看爽| 国内揄拍国产精品人妻在线| 午夜亚洲福利在线播放| 在线免费观看不下载黄p国产 | 国产精品,欧美在线| 亚洲国产精品成人综合色| 99热这里只有是精品在线观看| 国产精品野战在线观看| 亚洲电影在线观看av| 久久99热这里只有精品18| 久久6这里有精品| 免费看日本二区| 亚洲性久久影院| 午夜免费成人在线视频| 久久精品综合一区二区三区| 亚洲四区av| 国产高清三级在线| 国产精品亚洲一级av第二区| 少妇丰满av| 日本a在线网址| 一级av片app| 两性午夜刺激爽爽歪歪视频在线观看| 午夜a级毛片| 免费人成视频x8x8入口观看| 91精品国产九色| 免费看光身美女| 美女免费视频网站| 美女xxoo啪啪120秒动态图| 国产视频内射| 亚洲性久久影院| 亚洲欧美日韩无卡精品| 欧美性猛交黑人性爽| 三级毛片av免费| 国产伦精品一区二区三区四那| 一区二区三区高清视频在线| АⅤ资源中文在线天堂| 欧美色视频一区免费| 免费在线观看日本一区| 国产高潮美女av| 国产精品日韩av在线免费观看| 一进一出抽搐动态| 午夜视频国产福利| 悠悠久久av| 精品久久久久久久久av| 18禁黄网站禁片午夜丰满| 在线观看舔阴道视频| 日本与韩国留学比较| 男女之事视频高清在线观看| 极品教师在线视频| 人妻制服诱惑在线中文字幕| 男人舔女人下体高潮全视频| 国产精品久久久久久av不卡| 国内久久婷婷六月综合欲色啪| 欧美区成人在线视频| 亚洲av免费在线观看| 精品一区二区三区视频在线| 午夜视频国产福利| 亚洲真实伦在线观看| 亚洲av一区综合| 老司机福利观看| 国产精品1区2区在线观看.| 久久久久久久久大av| 给我免费播放毛片高清在线观看| 成人综合一区亚洲| 国产三级中文精品| 超碰av人人做人人爽久久| 老司机福利观看| 在线观看66精品国产| 免费电影在线观看免费观看| 1024手机看黄色片| 午夜视频国产福利| 国产精品福利在线免费观看| 狠狠狠狠99中文字幕| 高清在线国产一区| 成年女人看的毛片在线观看| 少妇的逼好多水| 日本黄色片子视频| 欧美不卡视频在线免费观看| 国产伦精品一区二区三区视频9| 日韩欧美精品免费久久| 精品一区二区三区av网在线观看| 22中文网久久字幕| 国内毛片毛片毛片毛片毛片| 校园春色视频在线观看| 欧美日韩乱码在线| 国产精品久久电影中文字幕| 久久香蕉精品热| 欧美精品国产亚洲| 国内精品久久久久久久电影| 欧美绝顶高潮抽搐喷水| 少妇被粗大猛烈的视频| 午夜老司机福利剧场| 亚洲av成人精品一区久久| 亚洲专区中文字幕在线| 亚洲av一区综合| 日本在线视频免费播放| 三级国产精品欧美在线观看| 亚洲精品国产成人久久av| 国产v大片淫在线免费观看| 小说图片视频综合网站| 国产高清视频在线观看网站| 亚洲人成伊人成综合网2020| 少妇高潮的动态图| 精品午夜福利在线看| 韩国av在线不卡| h日本视频在线播放| 亚洲国产欧洲综合997久久,| aaaaa片日本免费| 免费无遮挡裸体视频| 国产欧美日韩精品亚洲av| 国产日本99.免费观看| 国产白丝娇喘喷水9色精品| 免费高清视频大片| 成人鲁丝片一二三区免费| 小说图片视频综合网站| 欧美日韩中文字幕国产精品一区二区三区| 91精品国产九色| 免费观看在线日韩| 超碰av人人做人人爽久久| 日韩欧美国产在线观看| 欧美+亚洲+日韩+国产| 久久久久久久精品吃奶| 亚洲在线自拍视频| 夜夜看夜夜爽夜夜摸| 十八禁网站免费在线| 一级av片app| 国产极品精品免费视频能看的| 欧美性猛交╳xxx乱大交人| 别揉我奶头 嗯啊视频| 久久热精品热| 日韩av在线大香蕉| 亚洲欧美日韩高清在线视频| 日本与韩国留学比较| 一本久久中文字幕| av天堂中文字幕网| 日日干狠狠操夜夜爽| 免费看光身美女| 琪琪午夜伦伦电影理论片6080| 日日夜夜操网爽| 国产精品一区二区三区四区久久| 国产在视频线在精品| 草草在线视频免费看| 国产淫片久久久久久久久| 国产午夜精品论理片| 欧美又色又爽又黄视频| 国产黄a三级三级三级人| 欧美zozozo另类| 成人亚洲精品av一区二区| 中文资源天堂在线| 亚洲精品粉嫩美女一区| 久久久久久久亚洲中文字幕| 高清日韩中文字幕在线| 久久6这里有精品| 亚洲五月天丁香| 一区福利在线观看| 亚洲七黄色美女视频| 欧美另类亚洲清纯唯美| 一级黄色大片毛片| 99久久精品热视频| 亚洲av免费高清在线观看| 国产伦一二天堂av在线观看| 免费看a级黄色片| 在线观看一区二区三区| 久久6这里有精品| 午夜视频国产福利| 久久精品国产亚洲av香蕉五月| 亚洲人与动物交配视频| 日本成人三级电影网站| 亚洲最大成人av| 国产又黄又爽又无遮挡在线| 亚洲精品一区av在线观看| 99riav亚洲国产免费| 欧美日韩国产亚洲二区| 天堂动漫精品| 校园春色视频在线观看| 久久欧美精品欧美久久欧美| 春色校园在线视频观看| 天天一区二区日本电影三级| 变态另类丝袜制服| 亚洲第一电影网av| 午夜影院日韩av| 国产人妻一区二区三区在| av天堂中文字幕网| 观看美女的网站| 看片在线看免费视频| 色吧在线观看| 成人鲁丝片一二三区免费| 一进一出抽搐gif免费好疼| 国模一区二区三区四区视频| 欧美绝顶高潮抽搐喷水| 成年女人永久免费观看视频| 级片在线观看| www日本黄色视频网| 一级黄片播放器| 三级国产精品欧美在线观看| 久久午夜福利片| 伊人久久精品亚洲午夜| 亚洲精品国产成人久久av| 成人美女网站在线观看视频| 亚洲美女视频黄频| 成人二区视频| 22中文网久久字幕| 毛片女人毛片| 在线观看66精品国产| 日日干狠狠操夜夜爽| 在线观看舔阴道视频| 亚洲精品粉嫩美女一区| 国产 一区 欧美 日韩| 成人国产麻豆网| 成人二区视频| av天堂在线播放| 国产真实乱freesex| 在线观看免费视频日本深夜| 两个人视频免费观看高清| 亚洲无线观看免费| 久久久久久久午夜电影| 亚洲av成人精品一区久久| 欧美日韩瑟瑟在线播放| 舔av片在线| 国产在线男女| 少妇的逼好多水| 亚洲av成人精品一区久久| www.色视频.com| 亚洲av不卡在线观看| 精品不卡国产一区二区三区| 欧美zozozo另类| 不卡视频在线观看欧美| 国产精品女同一区二区软件 | 国产蜜桃级精品一区二区三区| 免费看日本二区| 99国产极品粉嫩在线观看| 国产美女午夜福利| 日本成人三级电影网站| 成年女人永久免费观看视频| 蜜桃亚洲精品一区二区三区| 全区人妻精品视频| 国产一区二区三区在线臀色熟女| 午夜福利成人在线免费观看| 国产午夜福利久久久久久| 国产一区二区三区在线臀色熟女| 欧美另类亚洲清纯唯美| 国产免费一级a男人的天堂| 国产精品女同一区二区软件 | 人妻夜夜爽99麻豆av| 国产伦人伦偷精品视频| 久久国产精品人妻蜜桃| 免费看a级黄色片| 日韩欧美精品免费久久| 在线国产一区二区在线| 搡女人真爽免费视频火全软件 | 亚洲国产精品合色在线| 亚洲av免费在线观看| 欧美色欧美亚洲另类二区| 国产乱人伦免费视频| 亚洲熟妇熟女久久| 一级av片app| 成年女人看的毛片在线观看| 俄罗斯特黄特色一大片| 久久久久精品国产欧美久久久| 老熟妇仑乱视频hdxx| 一区二区三区四区激情视频 | 欧美日本视频| 久久久久久久久大av| 一级黄片播放器| 此物有八面人人有两片| 中文字幕精品亚洲无线码一区| 成人国产一区最新在线观看| 久99久视频精品免费| 午夜福利高清视频| 在线免费十八禁| 日韩人妻高清精品专区| bbb黄色大片| 亚洲av中文av极速乱 | 不卡一级毛片| 午夜精品在线福利| 两个人的视频大全免费| 成人性生交大片免费视频hd| ponron亚洲| 美女 人体艺术 gogo| 久久欧美精品欧美久久欧美| 亚洲欧美激情综合另类| 国产精品福利在线免费观看| 女人被狂操c到高潮| 18禁黄网站禁片免费观看直播| 亚洲av不卡在线观看| 久久亚洲精品不卡| 99久久成人亚洲精品观看| 亚洲一级一片aⅴ在线观看| 亚洲在线自拍视频| 亚洲内射少妇av| 成人精品一区二区免费| 国产精品精品国产色婷婷| 久久久久久久久久成人| 99国产精品一区二区蜜桃av|