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

    徐深氣田新增產(chǎn)能管網(wǎng)障礙拓?fù)鋬?yōu)化

    2016-10-13 06:28:32陳雙慶魏立新王玉普
    關(guān)鍵詞:集氣集氣站多邊形

    陳雙慶, 劉 揚(yáng), 魏立新, 官 兵, 徐 哲, 張 爽, 王玉普

    ( 1.東北石油大學(xué) 石油工程學(xué)院,黑龍江 大慶 163318; 2. 中國石油昆侖燃?xì)庥邢薰?燃?xì)饧夹g(shù)研究院,黑龍江 哈爾濱 150000; 3. 中鐵哈爾濱鐵路局 綏化工務(wù)段,黑龍江 綏化 152000 )

    ?

    徐深氣田新增產(chǎn)能管網(wǎng)障礙拓?fù)鋬?yōu)化

    陳雙慶1, 劉揚(yáng)1, 魏立新1, 官兵1, 徐哲2, 張爽3, 王玉普1

    ( 1.東北石油大學(xué) 石油工程學(xué)院,黑龍江 大慶163318;2. 中國石油昆侖燃?xì)庥邢薰?燃?xì)饧夹g(shù)研究院,黑龍江 哈爾濱150000;3. 中鐵哈爾濱鐵路局 綏化工務(wù)段,黑龍江 綏化152000 )

    為降低徐深氣田新建產(chǎn)能管網(wǎng)系統(tǒng)投資,考慮集氣支線的連接方式和障礙對整體建設(shè)費(fèi)用的影響,將集氣支線連接方式作為優(yōu)化變量之一,采用R函數(shù)法和分層求凸包法,推導(dǎo)并建立任意障礙多邊形的數(shù)學(xué)表達(dá)式和可行布局約束條件,以管網(wǎng)建設(shè)費(fèi)用最小為目標(biāo),建立拓?fù)洳季謨?yōu)化數(shù)學(xué)模型。根據(jù)數(shù)學(xué)模型的結(jié)構(gòu)層次,構(gòu)建改進(jìn)的混合遺傳算法求解策略,設(shè)計(jì)幾何位置實(shí)數(shù)編碼和拓?fù)潢P(guān)系整數(shù)編碼的多參數(shù)級聯(lián)編碼方式,調(diào)整選擇復(fù)制遺傳算子的操作方式,建立考慮集氣支線氣量均勻性的自適應(yīng)種群進(jìn)化的適應(yīng)度函數(shù);結(jié)合 K—中心聚類法、叉積法、貪心算法和Prim算法給出初始種群的建立方法,驗(yàn)證模型和算法的有效性。結(jié)果表明:集氣支線的連接方式和障礙是管網(wǎng)布局優(yōu)化的重要影響因素,基于合理初值的多參數(shù)級聯(lián)編碼遺傳算法比常規(guī)遺傳算法的尋優(yōu)效果和速度更好。

    集氣支線; 拓?fù)鋬?yōu)化; R函數(shù); 障礙; 混合遺傳算法; 徐深氣田

    0 引言

    徐深氣田是大慶油田的主力天然氣產(chǎn)區(qū),已經(jīng)逐步發(fā)展并形成“高壓采氣,集中換熱、分離、脫水”的多井集氣工藝。徐深氣田的集輸流程:天然氣經(jīng)由集氣站處理后,通過集氣干線閥室掛接到集氣干線;再由集氣干線匯總到集氣總站。隨著徐深氣田區(qū)塊產(chǎn)能井和加密井等新增產(chǎn)能試井投產(chǎn),新建管網(wǎng)系統(tǒng)的規(guī)劃成為徐深氣田開發(fā)建設(shè)的關(guān)鍵。

    為保障集氣干線管道安全運(yùn)行,新建集氣支線管道通常掛接到已建集氣干線閥室,由于集氣站位置和管網(wǎng)連接形式耦合變化,集氣站與集氣干線閥室的不同連接方式導(dǎo)致整體的管線長度和建設(shè)費(fèi)用發(fā)生擾動。對于輻射—枝狀組合式管網(wǎng),通常描述為低級別站廠經(jīng)串接后直接與高級別站廠連接[1-3],忽略集氣干線閥室作為管網(wǎng)節(jié)點(diǎn)的作用,因此將集氣支線連接方式作為優(yōu)化變量、開展管網(wǎng)拓?fù)鋬?yōu)化研究具有重要意義。村屯、湖泊和生態(tài)保護(hù)區(qū)等障礙是油氣田集輸管網(wǎng)布局優(yōu)化的關(guān)鍵,甄寶軍[4]、魏立新[5]、雷楊[6]、左敏[7]等將障礙作為多邊形,確定管線和節(jié)點(diǎn)與障礙多邊形的位置關(guān)系和越障路徑,為管網(wǎng)障礙拓?fù)洳季謨?yōu)化研究奠定一定的基礎(chǔ)。這些研究沒有給出障礙多邊形的數(shù)學(xué)表達(dá)式,在判斷站址是否位于障礙內(nèi)部時,需要采用射線法和面積判斷法進(jìn)行繁瑣的計(jì)算與判斷。射線法沒有考慮障礙多邊形過于復(fù)雜引起的一些臨界特殊情況,如射線與障礙的一角相交的同時重合于障礙的一條邊,或者射線穿越障礙的多個頂點(diǎn)等不能準(zhǔn)確做出判斷[8]。

    筆者采用R函數(shù)方法和分層求凸包法,推導(dǎo)并構(gòu)建障礙多邊形的數(shù)學(xué)表征函數(shù),通過計(jì)算函數(shù)值的大小和解的個數(shù)確定目標(biāo)節(jié)點(diǎn)和管線與障礙的位置關(guān)系。對于徐深氣田的管網(wǎng)結(jié)構(gòu),考慮集氣站與集氣干線閥室的連接方式對拓?fù)洳季值挠绊?,建立含障礙的集輸管網(wǎng)新增產(chǎn)能拓?fù)鋬?yōu)化數(shù)學(xué)模型;設(shè)計(jì)染色體的編碼方式,改進(jìn)遺傳算法的優(yōu)化算子,優(yōu)化初始染色體形成方法,確定新增產(chǎn)能管網(wǎng)的最優(yōu)拓?fù)洳季帧?/p>

    1 障礙表示及越障路徑

    障礙的一般數(shù)學(xué)描述包括初等函數(shù)疊加和多邊形逼近的方式,文中采用多邊形逼近的方式。引入R函數(shù)和分層凸包法[9],將障礙多邊形各邊的隱式函數(shù)進(jìn)行集合操作,得到障礙多邊形的隱式函數(shù)?;谡系K多邊形函數(shù),判斷目標(biāo)節(jié)點(diǎn)與障礙的位置關(guān)系,通過最短路徑算法求得越障的最優(yōu)路徑。

    1.1多邊形表示

    對于n邊形障礙,逆時針遍歷多邊形,分別計(jì)算相鄰兩條邊的向量積,確定各頂點(diǎn)的凹凸性,進(jìn)而應(yīng)用分層凸包法計(jì)算多邊形的各層凸包集合,其中每一條邊都是某一層凸包的一個葉子節(jié)點(diǎn)。多邊形凸包劃分示意見圖1,凸包樹三級分層示意見圖2。對于分層葉子節(jié)點(diǎn),在奇數(shù)層中葉子結(jié)點(diǎn)表示的邊是凸的,在偶數(shù)層中葉子結(jié)點(diǎn)表示的邊是凹的。根據(jù)R函數(shù)的定義,兩個凸葉子節(jié)點(diǎn)的運(yùn)算為交運(yùn)算(∧),兩個凹葉子節(jié)點(diǎn)的運(yùn)算為并運(yùn)算(∨):

    (1)

    (2)

    圖1 多邊形凸包劃分示意

    圖2 多邊形凸包樹分層示意

    圖1的多邊形可以表述為

    (3)

    式中:fi(i=1,2,…i)為各直線方程的隱函數(shù)。

    推廣并構(gòu)建統(tǒng)一的函數(shù)表征方法。由式(1)和(2)可知,交、并運(yùn)算滿足交換律,即

    (4)

    (5)

    因此,式(1)和(2)等同于式(4)和(5),式(3)可以改寫為

    (6)

    障礙多邊形的表征方法:采用分層凸包法,將每一層凸包的葉子節(jié)點(diǎn)集合定義為φi={φi+1,φi+2,…fj…fk…},i=1,2,…m-1,j

    (7)

    按照從高層凸包向低層凸包求交、并運(yùn)算的法則,障礙多邊形i的隱式數(shù)學(xué)函數(shù)可以歸納為

    (8)

    1.2最優(yōu)越障路徑

    對于管線沿直線方向是否與障礙相交,主要采用射線法[8]和叉積法[4,6]等判斷。文中采用交點(diǎn)法與叉積法結(jié)合進(jìn)行判斷,對于端點(diǎn)位于障礙外的管線判斷步驟為

    (1)將目標(biāo)管線的直線方程l(x,y)與障礙多邊形函數(shù)聯(lián)立,采用Newton法求得直線方程的根的個數(shù),如果交點(diǎn)個數(shù)不小于2,則轉(zhuǎn)步驟2;否則,轉(zhuǎn)步驟4。

    (2)采用叉積法,遍歷判斷交點(diǎn)所在多邊形的邊的端點(diǎn)是否位于直線方程的兩側(cè),對于兩條邊具有公共頂點(diǎn)的情況,則判斷兩條邊是否位于l(x,y)的兩側(cè),如果是,則轉(zhuǎn)步驟3;否則,轉(zhuǎn)步驟4。

    (3)目標(biāo)管線與障礙多邊形相交,輸出計(jì)算結(jié)果。

    (4)目標(biāo)管線不與障礙多邊形相交,輸出計(jì)算結(jié)果。

    根據(jù)計(jì)算結(jié)果分別得到兩節(jié)點(diǎn)間的距離,由于徐深氣田地處松嫩平原腹地,地面高程變化在100 m內(nèi),只考慮兩者間的平面距離。當(dāng)兩節(jié)點(diǎn)間線段不與障礙相交時,兩節(jié)點(diǎn)間為歐式距離;當(dāng)兩節(jié)點(diǎn)間線段與障礙相交時,兩節(jié)點(diǎn)間距離為采用Dijkstra算法求得的最短路徑距離,最短路徑采用一條折線段表示。

    2 數(shù)學(xué)模型

    在輻射—枝狀集氣系統(tǒng)中,氣井與集氣站連接形式為輻射狀,集氣站與集氣干線和集氣總站之間連接形式為枝狀。為保證管網(wǎng)系統(tǒng)的供氣安全和方便管理,集氣站的數(shù)量和處理量應(yīng)該在一定范圍內(nèi),集氣站與氣井之間的連接關(guān)系應(yīng)該滿足集輸半徑約束和隸屬唯一性約束,集氣支線的輸氣量和連接關(guān)系應(yīng)該滿足一定約束條件,集氣站的幾何位置應(yīng)該處在障礙外,建立拓?fù)鋬?yōu)化數(shù)學(xué)模型:

    (9)

    (10)

    (11)

    (12)

    (13)

    (14)

    (15)

    (16)

    (17)

    (18)

    (19)

    (20)

    (21)

    (22)

    式(10)為集輸半徑約束;式(11)和(12)為集氣站處理量約束,其中式(11)表示集氣站的處理量應(yīng)該與其所轄各氣井氣量相等;式(13)和(14)為集氣支線輸量約束,其中式(13)表示各集氣站的氣量總和應(yīng)該等于各個集氣支線的輸量之和,式(14)表示集氣支線的輸量應(yīng)該在一定范圍內(nèi);式(15)表示新建集氣站節(jié)點(diǎn)數(shù)量約束;式(16)表示集氣站節(jié)點(diǎn)與集氣干線閥室節(jié)點(diǎn)之間的新建連接關(guān)系約束,即新建的集氣管線和新建集氣支線數(shù)量之和應(yīng)該等于新建集氣站的數(shù)量;式(17)表示集氣站節(jié)點(diǎn)和集氣干線閥室節(jié)點(diǎn)連接關(guān)系約束,即每一個集氣站只能直接或串接后與一個干線閥室相連;式(18)表示井站隸屬關(guān)系唯一性約束;式(19)為可行布局約束,即集氣站不能位于障礙范圍內(nèi);式(20-22)為取值范圍約束。

    3 模型求解

    數(shù)學(xué)模型中集氣站址與拓?fù)溥B接關(guān)系存在耦合關(guān)系,通常采用遺傳算法[2,5,10]、粒子群算法[11]和蟻群算法[12]等進(jìn)行求解。文中采用操作簡單和易獲取全局最優(yōu)解的遺傳算法進(jìn)行求解,將管網(wǎng)結(jié)構(gòu)轉(zhuǎn)化為有向連通圖,應(yīng)用遺傳算法整體求解集氣站的幾何位置和各級節(jié)點(diǎn)之間的連接關(guān)系,建立一種高效的混合遺傳求解策略。

    首先,根據(jù)集氣站處理能力給出集氣站的數(shù)量,結(jié)合Prim算法、K—中心點(diǎn)聚類算法和貪心算法,給出初始種群,應(yīng)用幾何位置實(shí)數(shù)編碼和拓?fù)潢P(guān)系整數(shù)編碼的多參數(shù)級聯(lián)編碼方式。其次,應(yīng)用改進(jìn)的選擇復(fù)制算子和交叉、變異算子執(zhí)行遺傳操作,并判斷染色體表征的管網(wǎng)形態(tài)、集輸半徑、障礙布局和集氣站處理量等是否滿足約束條件,進(jìn)而對染色體進(jìn)行評估。最后,通過調(diào)整集氣站的數(shù)量實(shí)現(xiàn)多種布局的對比優(yōu)選,求得管網(wǎng)的最優(yōu)布局,求解算法流程見圖3。

    遺傳算法過程主要包括染色體編碼方式和適應(yīng)度函數(shù)的確定,以及交叉、變異、選擇復(fù)制等遺傳算子的設(shè)計(jì)等??紤]數(shù)學(xué)模型的結(jié)構(gòu)特點(diǎn),將幾何位置求解和拓?fù)溥B接關(guān)系求解結(jié)合,給出遺傳算法操作步驟。

    3.1染色體編碼

    將集氣站的幾何位置優(yōu)化為連續(xù)變量,拓?fù)溥B接關(guān)系優(yōu)化為離散變量。采用多參數(shù)級聯(lián)編碼,首先以所有新增集氣站的坐標(biāo)序列作為染色體基因,采用實(shí)數(shù)編碼方式,構(gòu)建染色體的幾何位置并優(yōu)化子串;其次,為充分利用已建集氣系統(tǒng),新建氣井可以與已建集氣站連接,對集氣站節(jié)點(diǎn)進(jìn)行排序,將氣井隸屬的集氣站的序號作為染色體基因,構(gòu)建“井—站”連接關(guān)系染色體子串;再次,考慮新增集氣站可以與已建集氣站、新建集氣站、集氣干線閥室三類節(jié)點(diǎn)連接,在集氣站節(jié)點(diǎn)排序的基礎(chǔ)上,對集氣干線閥室節(jié)點(diǎn)進(jìn)行排序,將與新建集氣站連接的節(jié)點(diǎn)序號作為染色體基因,建立“站—站、站—閥”連接關(guān)系子串;最后,將各染色體子串按照一定順序連接成為一個完整染色體串:

    (23)

    3.2初始群體產(chǎn)生

    遺傳算法結(jié)果受初始群體的影響較大,良好的初始染色體可以有效降低遺傳算法收斂時間。遺傳操作初始群體產(chǎn)生步驟:

    (2)獲取集氣干線的拐點(diǎn)坐標(biāo),應(yīng)用叉積法,判斷距離各個集氣站(已建和新建)最近的集氣干線管段與集氣站的關(guān)系,將位于同一側(cè)的集氣站進(jìn)行標(biāo)記劃分。

    圖3 數(shù)學(xué)模型求解算法流程Fig.3 Flow chart of algorithm for solving mathematical model

    (4)將染色體子串按順序排列,構(gòu)成一個完整染色體。

    3.3適應(yīng)度函數(shù)確定

    為保證集氣管網(wǎng)的運(yùn)行安全,在考慮距離因素的基礎(chǔ)上,每條集氣支線串接的集氣站輸氣量應(yīng)該較為均勻。適應(yīng)度函數(shù)遞增的方向應(yīng)該與種群更新的方向相同,目標(biāo)函數(shù)式(9)為極小化問題。采用混合遺傳模擬退火算法的思路,為了優(yōu)選每條集氣支線輸氣量較為均勻的染色體,采用罰函數(shù)法對適應(yīng)度函數(shù)進(jìn)行設(shè)計(jì),通過構(gòu)建罰函數(shù)使集氣支線之間輸氣量差別較大的染色體在進(jìn)化過程中被淘汰,得到自適應(yīng)種群進(jìn)化的適應(yīng)度函數(shù)表達(dá)式:

    (24)

    式中:Fi為第i個染色體對應(yīng)的目標(biāo)函數(shù);Fmin、Fmax分別為當(dāng)代所有染色體的最小、最大目標(biāo)函數(shù);MT為較大的正實(shí)數(shù),可以憑經(jīng)驗(yàn)給定,也可取為當(dāng)代染色體集氣支線的最大輸氣量的T倍;T為集氣支線的條數(shù);Qi、Qj分別為第i條和第j條集氣支線的輸氣量,由串接的所有集氣站氣量之和求得。種群進(jìn)化初期Fmax較大,各染色體適應(yīng)度差距較小,優(yōu)秀個體數(shù)量較多,能夠避免出現(xiàn)早熟現(xiàn)象;隨著種群繼續(xù)進(jìn)化,F(xiàn)max的數(shù)值減小,集氣支線輸氣量逐漸均勻,優(yōu)秀個體優(yōu)勢逐漸增大,劣質(zhì)個體被逐漸淘汰,可以促進(jìn)算法的收斂。

    3.4交叉操作

    交叉操作是促使染色體基因充分重組、加速收斂的有效操作,隨機(jī)選擇種群中的k(k為小于P的偶數(shù))個染色體進(jìn)行配對交叉,對于每一代交叉操作,隨機(jī)選取算術(shù)交叉和兩點(diǎn)交叉的方法,將父代染色體中的基因以一定的概率遺傳到子代。在實(shí)際交叉過程中,兩父代染色體的相同類別子串基因?qū)?yīng)交叉。對于染色體的整數(shù)基因,交叉產(chǎn)生的實(shí)數(shù)基因在向下取整后進(jìn)入下一代染色體。

    3.5變異操作

    為提升算法局部搜索能力,輔助交叉操作完成全局最優(yōu)解的獲得,采用隨機(jī)選擇單點(diǎn)基因變異的策略,其整數(shù)基因的單點(diǎn)變異改變量為隨機(jī)整數(shù),其實(shí)數(shù)基因的單點(diǎn)變異改變量為隨機(jī)實(shí)數(shù)。

    3.6染色體可行性判別及調(diào)整步驟

    (1)由于存在障礙,產(chǎn)生新個體的部分基因可能位于障礙范圍,為避免重復(fù)生成新個體,采用距離不可行基因最近的障礙頂點(diǎn)坐標(biāo)代替。

    (2)為保證集氣站與集氣干線閥室之間的連接關(guān)系為樹狀結(jié)構(gòu),應(yīng)檢查集氣站節(jié)點(diǎn)之間是否成環(huán),如果成環(huán),則應(yīng)用破圈法隨機(jī)選擇環(huán)中一條管線,將它連接到環(huán)外的其他集氣站或集氣干線閥室上,并更改相應(yīng)的基因數(shù),直到該管網(wǎng)形態(tài)滿足樹狀結(jié)構(gòu)為止。

    (3)產(chǎn)生的染色體需要進(jìn)行可行性判斷,根據(jù)染色體的表現(xiàn)形式計(jì)算采氣管線的長度、集氣站處理量、集氣支線輸氣量,如果滿足約束條件,則產(chǎn)生的染色體可行;否則,重新執(zhí)行交叉和變異操作,重新生成染色體。

    3.7改進(jìn)的選擇復(fù)制操作

    3.8終止準(zhǔn)則

    當(dāng)連續(xù)進(jìn)化q代后最優(yōu)染色體沒有變化時,可以認(rèn)為當(dāng)前最優(yōu)染色體即為最優(yōu)解。并設(shè)定最大進(jìn)化代數(shù)nmax,以保證算法在一定時間內(nèi)終止計(jì)算。

    4 管材造價

    一般在生產(chǎn)中計(jì)算標(biāo)準(zhǔn)大氣壓下的氣體流量,需要將標(biāo)況下的體積轉(zhuǎn)換成工況下的體積,采用文獻(xiàn)[16]的方法進(jìn)行轉(zhuǎn)換。其中流速的選取應(yīng)該滿足經(jīng)濟(jì)流速的要求,徐深氣田采氣管道的經(jīng)濟(jì)流速為4~6 m/s,集氣管道和集氣支線管道的經(jīng)濟(jì)流速為10~15 m/s。在選取管徑時,應(yīng)該按大于計(jì)算管徑的最小管道規(guī)格進(jìn)行選擇。各種規(guī)格管材造價見表1。

    表1 各種規(guī)格管材造價

    5 計(jì)算實(shí)例

    徐深氣田某區(qū)塊有已建集氣總站1座、集氣站6座,已建氣井38口、集氣干線1條、集氣干線閥室4個,新建投產(chǎn)共42口井,新增產(chǎn)能226.8×104m3/d,集氣支線輸氣量范圍為(20~120)×104m3/d,集氣站處理量范圍為(20~80)×104m3,集氣半徑為3 000 m,初始井位布局見圖4?;贑#語言編制求解程序,分別采用基本混合遺傳算法和文中方法對某區(qū)塊進(jìn)行管網(wǎng)拓?fù)鋬?yōu)化。初始參數(shù)設(shè)置:種群規(guī)模為60,進(jìn)化代數(shù)參數(shù)為nmax,q=20,初始交叉概率為0.95,初始變異概率為0.01,MT=5×106。對于集氣支線的處理,基本混合遺傳算法是在求得管網(wǎng)其他拓?fù)浣Y(jié)構(gòu)的基礎(chǔ)上,按照距離最小原則將集氣站串接到集氣干線閥室上,采用基本混合遺傳算法求得的管網(wǎng)拓?fù)洳季忠妶D5,采用文中方法求得的管網(wǎng)拓?fù)洳季忠妶D6,兩種方法的優(yōu)化結(jié)果見表2-3,求解性能結(jié)果見圖7。由表2-3和圖5-7可以看出,采用文中方法可以有效降低投資,獲得更優(yōu)的數(shù)值解和求解效率,并且所得布局的集氣支線氣量更加均勻。

    圖4 徐深氣田某區(qū)塊初始井位布局Fig.4 The initial wells layout of a block in Xushen gas field

    圖5 徐深氣田某區(qū)塊基本混合遺傳算法優(yōu)化布局

    圖6 徐深氣田某區(qū)塊文中方法優(yōu)化布局

    圖7 基本混合遺傳算法和文中方法優(yōu)化結(jié)果

    管道規(guī)格/(?×b:mm×mm)基本混合遺傳算法文中方法管道長度/m費(fèi)用/萬元管道長度/m費(fèi)用/萬元50×421708.152172.6721014.862103.3060×46114.78703.815540.88637.7476×512627.021642.9011484.261494.2276×61751.38231.372660.82351.5289×53378.48476.081534.85230.43108×71579.13260.204744.34783.48114×91727.75311.25127×101493.26298.91

    表3 文中方法與基本混合遺傳算法新建管道結(jié)果

    6 結(jié)論

    (1)基于R函數(shù)法和分層求凸包法,推導(dǎo)并建立任意障礙多邊形的數(shù)學(xué)函數(shù)表達(dá)式,簡化判別目標(biāo)節(jié)點(diǎn)和目標(biāo)管線與障礙多邊形關(guān)系的計(jì)算方法,為障礙描述提供一種可行有效的方法。

    (2)考慮集氣干線閥室和障礙的客觀存在對拓?fù)洳季值挠绊?,以管網(wǎng)建設(shè)費(fèi)用最小為目標(biāo),建立徐深氣田輻射—枝狀管網(wǎng)形態(tài)的新增產(chǎn)能拓?fù)洳季謨?yōu)化數(shù)學(xué)模型。該數(shù)學(xué)模型也適用于其他輻射—枝狀管網(wǎng)新增產(chǎn)能布局優(yōu)化。

    (3)將遺傳算法、貪心算法、K-中心聚類法和Prim算法結(jié)合,優(yōu)化初始種群的產(chǎn)生方法,設(shè)計(jì)整數(shù)和實(shí)數(shù)多參數(shù)級聯(lián)編碼方式,改進(jìn)選擇復(fù)制遺傳算子和染色體評價函數(shù),形成一種改進(jìn)的混合遺傳算法求解策略。該算法可以實(shí)現(xiàn)集氣站位置和拓?fù)溥B接關(guān)系的整體優(yōu)化求解,具有更好的尋優(yōu)效果和效率。

    (4)基于文中方法,開發(fā)徐深氣田集氣系統(tǒng)優(yōu)化設(shè)計(jì)軟件對徐深氣田采氣區(qū)塊進(jìn)行試算,優(yōu)化后可降低8.61%系統(tǒng)投資,驗(yàn)證模型及算法的有效性。

    [1]Sanayea S, Mahmoudimehra J. Optimal design of a natural gas transmission network layout [J]. Chemical Engineering Research and Design, 2013,91(12):2465-2476.

    [2]李自力,孫云峰,張子波,等.基于遺傳算法的氣田集輸管網(wǎng)整體優(yōu)化方法[J].天然氣工業(yè),2011,31(8):86-89.

    Li Zili, Sun Yunfeng, Zhang Zibo, et al. A global optimization method based on genetic algorithms for gas gathering pipeline network in a gas field [J]. Natural Gas Industry, 2011,31(8):86-89.

    [3]商冠琪.枝狀管網(wǎng)布局優(yōu)化中的3類問題[J].化工機(jī)械,2014,41(3):342-344.

    Shang Guanqi. Three problems in optimal design of branched pipeline network [J]. Chemical Machinery, 2014,41(3):342-344.

    [4]甄寶軍.考慮凸形障礙物的天然氣集輸管網(wǎng)系統(tǒng)規(guī)劃研究[D].成都:西南石油大學(xué),2005.

    Zhen Baojun. The study of gas gathering pipe network system considering the convex obstacles [D]. Chengdu: Southwest Petroleum University, 2005.

    [5]魏立新,劉揚(yáng).油氣集輸系統(tǒng)障礙拓?fù)洳季謨?yōu)化設(shè)計(jì)方法 [J].石油學(xué)報(bào),2006,27(6):120-124.

    Wei Lixin, Liu Yang. Obstacle topological layout optimization design of oil gas gathering and transferring system [J]. Acta Petrolei Sinica, 2006,27(6):120-124.

    [6]雷楊.基于凸形障礙物的注水管網(wǎng)優(yōu)化研究[D].北京:中國石油大學(xué),2011.

    Lei Yang. Study on optimization for oilfield water-injection pipe network under complicated terrain [D]. Beijing: China University of Petroleum, 2011.

    [7]左敏.磨溪?dú)馓镌嚥杉敼芫W(wǎng)設(shè)計(jì)優(yōu)化研究[D].成都:西南石油大學(xué),2014

    Zuo Min. Study on design optimization of gas gathering and transportation pipeline network in the gas field of Moxi [D]. Chengdu: Southwest Petroleum University, 2014.

    [8]劉潤濤,劉玉珍.點(diǎn)在多邊形內(nèi)測試的新算法[J].工程圖學(xué)學(xué)報(bào),2008,28(2):89-93.

    Liu Runtao, Liu Yuzhen. A new algorithm for determining whether a point is inside of a polygon [J]. Journal of Engineering Graphics, 2008,28(2):89-93.

    [9]吳堅(jiān),王小椿,姜虹,等.多邊形的隱函數(shù)表示法[J].計(jì)算機(jī)工程與應(yīng)用,2003,39(32):87-89.

    Wu Jian, Wang Xiaochun, Jiang Hong, et al. Representation for polygons with implicit function [J]. Computer Engineering and Applications, 2003,39(32):87-89.

    [10]劉揚(yáng),鞠志忠,鮑云波.一類多級星式網(wǎng)絡(luò)的拓?fù)鋬?yōu)化設(shè)計(jì)方法[J].大慶石油學(xué)院學(xué)報(bào),2009,33(2):68-73.

    Liu Yang, Ju Zhizhong, Bao Yunbo. Topological optimization design of a multilevel star network [J]. Journal of Daqing Petroleum Institute, 2009,33(2):68-73.

    [11]Ibrahim Mtolera, Li Haibin, Liu Ye, et al. Optimization of tree pipe networks layout and size, using particle swarm optimization [J]. WSEAS Transactions on Computers, 2014,13(1):219-230.

    [12]Afshar M H. Layout and size optimization of tree-like pipe network by incremental solution building ants [J]. Canadian Journal of Civil Engineering, 2008,35(2):129-139.

    [13]楊建軍,戰(zhàn)紅,劉揚(yáng),等.星狀原油集輸管網(wǎng)拓?fù)鋬?yōu)化的混合遺傳算法[J].西南石油大學(xué)學(xué)報(bào):自然科學(xué)版,2008,30(4):166-169.

    Yang Jianjun, Zhan Hong, Liu Yang, et al. Hybrid genetic algorithm for topology optimization of stellated oil gathering and transportation pipeline network [J]. Journal of Southwest Petroleum University: Science & Technology Edition, 2008,30(4):166-169.

    [14]周軍,李曉平,鄧濤,等.集輸系統(tǒng)優(yōu)化設(shè)計(jì)研究的體系結(jié)構(gòu)與發(fā)展方向[J].油氣儲運(yùn),2014,33(7):707-713.

    Zhou Jun, Li Xiaoping, Deng Tao, et al. Structure and trend of optimal design of gathering system [J]. Oil & Gas Storage and Transportation, 2014,33(7):707-713.

    [15]黃維和.大型天然氣管網(wǎng)系統(tǒng)可靠性[J].石油學(xué)報(bào),2013,34(2):401-404.

    Huang Weihe. Reliability of large-scale natural gas pipeline network [J]. Acta Petrolei Sinica, 2013,34(2):401-404.

    [16]中國國家標(biāo)準(zhǔn)化管理委員會,SY/T 6143—1996,天然氣流量的標(biāo)準(zhǔn)孔板計(jì)量方法[S].北京:中國標(biāo)準(zhǔn)出版社,1996.

    China National Standardization Management Committee. SY/T 6143—1996, Standard orifice plate metering method for natural gas flow [S]. Beijing: China Standard Press, 1996.

    2016-03-04;編輯:任志平

    國家科技支撐計(jì)劃項(xiàng)目(2012BAH28F03);國家自然科學(xué)基金面上項(xiàng)目(51674086);東北石油大學(xué)創(chuàng)新科研項(xiàng)目(YJSCX2015-012NEPU)

    陳雙慶(1990-), 男,博士研究生,主要從事油氣集輸系統(tǒng)優(yōu)化與節(jié)能降耗技術(shù)方面的研究。

    10.3969/j.issn.2095-4107.2016.04.012

    TE863

    A

    2095-4107(2016)04-0096-10

    猜你喜歡
    集氣集氣站多邊形
    多邊形中的“一個角”問題
    鋁電解槽集氣結(jié)構(gòu)的數(shù)值模擬研究
    一種用于橡膠擠出機(jī)的環(huán)保伸縮集氣罩
    天然氣凈化廠和集氣站在消防監(jiān)督檢查中的常見問題分析
    延安氣田集氣站的無人值守改造
    化工管理(2021年7期)2021-05-13 00:46:32
    集氣站管道典型缺陷及失效案例分析
    多邊形的藝術(shù)
    解多邊形題的轉(zhuǎn)化思想
    普光氣田集氣總站水洗脫氯工程項(xiàng)目管理難點(diǎn)與對策探析
    多邊形的鑲嵌
    少妇裸体淫交视频免费看高清| 国产精品一区二区性色av| 精品无人区乱码1区二区| 久久精品国产亚洲av香蕉五月| 中文资源天堂在线| 非洲黑人性xxxx精品又粗又长| 99在线人妻在线中文字幕| 女同久久另类99精品国产91| 国产久久久一区二区三区| 两个人视频免费观看高清| 亚洲精品日韩av片在线观看| 国产午夜精品一二区理论片| 国产中年淑女户外野战色| 听说在线观看完整版免费高清| 欧美+亚洲+日韩+国产| 亚洲成av人片在线播放无| 国产69精品久久久久777片| 日韩成人av中文字幕在线观看| 99在线人妻在线中文字幕| 亚洲av.av天堂| 成人毛片60女人毛片免费| 一本精品99久久精品77| 插阴视频在线观看视频| 国产成人精品婷婷| 欧美+亚洲+日韩+国产| 亚洲,欧美,日韩| 校园春色视频在线观看| 成年女人永久免费观看视频| 亚洲精品亚洲一区二区| 国产精品不卡视频一区二区| 亚洲精品456在线播放app| av卡一久久| 日韩三级伦理在线观看| 亚洲精品亚洲一区二区| 国产探花极品一区二区| 级片在线观看| 国产一区二区亚洲精品在线观看| 亚洲精品自拍成人| 国产成人freesex在线| 国内久久婷婷六月综合欲色啪| 男女视频在线观看网站免费| 久久久久久久久久黄片| 久久久午夜欧美精品| 看免费成人av毛片| 乱码一卡2卡4卡精品| 国产人妻一区二区三区在| 免费大片18禁| 婷婷六月久久综合丁香| 美女脱内裤让男人舔精品视频 | 中文字幕制服av| 高清毛片免费观看视频网站| 五月玫瑰六月丁香| 亚洲自拍偷在线| 国产毛片a区久久久久| 久久人妻av系列| 国产精品美女特级片免费视频播放器| 噜噜噜噜噜久久久久久91| 欧美bdsm另类| 国产成人精品久久久久久| 国产伦理片在线播放av一区 | 亚洲性久久影院| 尤物成人国产欧美一区二区三区| 给我免费播放毛片高清在线观看| 亚洲内射少妇av| 激情 狠狠 欧美| 国产蜜桃级精品一区二区三区| 亚洲经典国产精华液单| 欧美潮喷喷水| 特大巨黑吊av在线直播| 免费观看的影片在线观看| 我要看日韩黄色一级片| 草草在线视频免费看| 一区二区三区高清视频在线| 欧美高清成人免费视频www| 女的被弄到高潮叫床怎么办| 久久午夜亚洲精品久久| 亚洲最大成人av| 99久久中文字幕三级久久日本| 久久热精品热| 亚洲色图av天堂| 国产成人精品一,二区 | 网址你懂的国产日韩在线| 国产欧美日韩精品一区二区| 人妻久久中文字幕网| 国产精品久久久久久精品电影小说 | 深爱激情五月婷婷| 国产精品国产三级国产av玫瑰| 日韩av不卡免费在线播放| 成年版毛片免费区| 国内精品久久久久精免费| 日本黄色片子视频| 禁无遮挡网站| 久久久精品欧美日韩精品| avwww免费| 此物有八面人人有两片| 亚洲国产精品合色在线| 男女那种视频在线观看| 如何舔出高潮| 人妻少妇偷人精品九色| 亚洲婷婷狠狠爱综合网| 噜噜噜噜噜久久久久久91| av又黄又爽大尺度在线免费看 | 91狼人影院| 可以在线观看毛片的网站| 午夜精品国产一区二区电影 | 亚洲精品粉嫩美女一区| 麻豆成人午夜福利视频| 亚洲aⅴ乱码一区二区在线播放| 国产黄a三级三级三级人| 五月伊人婷婷丁香| 此物有八面人人有两片| 小说图片视频综合网站| 看免费成人av毛片| 尾随美女入室| 久久这里有精品视频免费| 级片在线观看| 中文字幕av在线有码专区| 乱码一卡2卡4卡精品| 欧美极品一区二区三区四区| 国产麻豆成人av免费视频| 卡戴珊不雅视频在线播放| 噜噜噜噜噜久久久久久91| 又粗又硬又长又爽又黄的视频 | 欧美日韩乱码在线| 国产男人的电影天堂91| 日韩中字成人| 国产av在哪里看| 亚洲中文字幕一区二区三区有码在线看| 精品少妇黑人巨大在线播放 | av黄色大香蕉| 蜜桃亚洲精品一区二区三区| 国产精品美女特级片免费视频播放器| 亚洲精品影视一区二区三区av| 丰满乱子伦码专区| 国产成人精品婷婷| 我的女老师完整版在线观看| 少妇丰满av| 91aial.com中文字幕在线观看| 永久网站在线| 成人亚洲欧美一区二区av| 欧美不卡视频在线免费观看| 亚洲性久久影院| 91精品一卡2卡3卡4卡| 国产激情偷乱视频一区二区| 欧美一级a爱片免费观看看| 色尼玛亚洲综合影院| 欧洲精品卡2卡3卡4卡5卡区| 乱码一卡2卡4卡精品| 中出人妻视频一区二区| 麻豆国产av国片精品| 18禁裸乳无遮挡免费网站照片| 狂野欧美激情性xxxx在线观看| 久久久久久久久大av| 在线观看一区二区三区| 99国产极品粉嫩在线观看| 日韩一本色道免费dvd| 91久久精品电影网| 我的老师免费观看完整版| 国产亚洲精品久久久久久毛片| 国产精品99久久久久久久久| 波多野结衣巨乳人妻| 日本一本二区三区精品| 级片在线观看| 国产91av在线免费观看| 国产精品伦人一区二区| 亚洲一区高清亚洲精品| 乱系列少妇在线播放| 99在线视频只有这里精品首页| 毛片一级片免费看久久久久| 国产黄片视频在线免费观看| 青青草视频在线视频观看| 最近视频中文字幕2019在线8| 变态另类丝袜制服| 久久久久久伊人网av| 亚洲av男天堂| 丰满的人妻完整版| 黄片wwwwww| 3wmmmm亚洲av在线观看| 免费观看a级毛片全部| 黄片wwwwww| 亚洲美女视频黄频| 热99在线观看视频| 久久亚洲精品不卡| 成人国产麻豆网| 草草在线视频免费看| 亚洲综合色惰| 国产午夜福利久久久久久| 可以在线观看毛片的网站| 国内久久婷婷六月综合欲色啪| 成人亚洲精品av一区二区| 久久精品久久久久久久性| 成人二区视频| 亚洲av电影不卡..在线观看| 欧美色欧美亚洲另类二区| 国产69精品久久久久777片| 又爽又黄无遮挡网站| 女的被弄到高潮叫床怎么办| 亚洲成人久久性| 色综合亚洲欧美另类图片| 18禁在线播放成人免费| 中文字幕精品亚洲无线码一区| 国产精品麻豆人妻色哟哟久久 | 3wmmmm亚洲av在线观看| 国产不卡一卡二| 哪个播放器可以免费观看大片| 亚洲成人精品中文字幕电影| 亚洲欧美精品专区久久| 五月伊人婷婷丁香| 色综合亚洲欧美另类图片| 久久99热这里只有精品18| 99热只有精品国产| 爱豆传媒免费全集在线观看| 国产高清视频在线观看网站| 亚洲欧洲日产国产| 免费大片18禁| 天天躁日日操中文字幕| 亚洲欧美清纯卡通| 能在线免费观看的黄片| 夜夜夜夜夜久久久久| 中出人妻视频一区二区| 中出人妻视频一区二区| 女的被弄到高潮叫床怎么办| 又爽又黄a免费视频| 亚洲av成人精品一区久久| 亚洲av中文字字幕乱码综合| 亚洲性久久影院| h日本视频在线播放| 成人特级黄色片久久久久久久| 精品熟女少妇av免费看| 亚洲精品乱码久久久久久按摩| 日韩人妻高清精品专区| 国产亚洲5aaaaa淫片| 色吧在线观看| 91久久精品国产一区二区三区| 嫩草影院精品99| 国产黄色小视频在线观看| 一边摸一边抽搐一进一小说| 最后的刺客免费高清国语| 简卡轻食公司| 成年免费大片在线观看| 亚洲国产精品国产精品| av天堂中文字幕网| 天天躁夜夜躁狠狠久久av| 国产三级在线视频| 色哟哟哟哟哟哟| 白带黄色成豆腐渣| 久久久久久久久久久免费av| 国内精品美女久久久久久| 国产亚洲5aaaaa淫片| 亚洲美女视频黄频| 国产av不卡久久| 综合色丁香网| 在线观看午夜福利视频| 免费观看的影片在线观看| 99热6这里只有精品| 高清毛片免费看| 99久久精品一区二区三区| 日韩,欧美,国产一区二区三区 | 国产精品国产高清国产av| 日本成人三级电影网站| ponron亚洲| 国产精品福利在线免费观看| 丰满的人妻完整版| 亚洲内射少妇av| 久久6这里有精品| 精品一区二区三区视频在线| 久久精品国产99精品国产亚洲性色| 一卡2卡三卡四卡精品乱码亚洲| 一本精品99久久精品77| kizo精华| 亚洲人成网站在线播放欧美日韩| 边亲边吃奶的免费视频| 国产 一区精品| 亚州av有码| 精品免费久久久久久久清纯| 嫩草影院入口| 一进一出抽搐动态| 内地一区二区视频在线| 欧美潮喷喷水| 亚洲色图av天堂| 悠悠久久av| 男女下面进入的视频免费午夜| 亚州av有码| 国产高清有码在线观看视频| 精品一区二区三区视频在线| 中文字幕免费在线视频6| 亚洲三级黄色毛片| 亚洲色图av天堂| 一进一出抽搐gif免费好疼| av免费观看日本| 国国产精品蜜臀av免费| 卡戴珊不雅视频在线播放| 日韩 亚洲 欧美在线| 尤物成人国产欧美一区二区三区| 成人午夜精彩视频在线观看| 变态另类成人亚洲欧美熟女| 日韩精品青青久久久久久| 国内揄拍国产精品人妻在线| 好男人视频免费观看在线| av卡一久久| 久久99热6这里只有精品| 成人美女网站在线观看视频| 大香蕉久久网| 久久久午夜欧美精品| 99热只有精品国产| 免费看美女性在线毛片视频| 欧美xxxx黑人xx丫x性爽| 一级黄片播放器| 丝袜喷水一区| 人妻制服诱惑在线中文字幕| 男女啪啪激烈高潮av片| 免费电影在线观看免费观看| 亚洲最大成人av| 一本一本综合久久| 校园春色视频在线观看| 午夜视频国产福利| 国产极品精品免费视频能看的| av在线观看视频网站免费| 国产精品一区二区三区四区久久| 校园人妻丝袜中文字幕| 简卡轻食公司| 亚洲欧美日韩东京热| 亚洲国产精品sss在线观看| 国产男人的电影天堂91| 欧美xxxx黑人xx丫x性爽| 91精品国产九色| 久久久久久久久中文| 变态另类成人亚洲欧美熟女| 美女大奶头视频| 久久人人爽人人爽人人片va| 中国美女看黄片| 国产精品伦人一区二区| 波野结衣二区三区在线| 一级二级三级毛片免费看| 精品人妻一区二区三区麻豆| 久久久久久九九精品二区国产| 久久国产乱子免费精品| 97在线视频观看| 99久久精品国产国产毛片| 99热这里只有精品一区| 国产午夜精品论理片| 久久精品人妻少妇| АⅤ资源中文在线天堂| 亚洲av成人精品一区久久| 麻豆av噜噜一区二区三区| 一级毛片aaaaaa免费看小| 亚洲第一区二区三区不卡| 99热网站在线观看| 亚洲aⅴ乱码一区二区在线播放| 少妇人妻精品综合一区二区 | 国产精品伦人一区二区| 中文资源天堂在线| 日韩人妻高清精品专区| 99在线人妻在线中文字幕| 精品久久久久久久末码| 亚洲av不卡在线观看| 国产成人一区二区在线| av天堂在线播放| 淫秽高清视频在线观看| 亚洲欧美日韩高清专用| 老女人水多毛片| 插逼视频在线观看| 91麻豆精品激情在线观看国产| 久久6这里有精品| 成熟少妇高潮喷水视频| 日韩中字成人| 国产成人精品婷婷| 亚洲欧美日韩高清专用| 黄色欧美视频在线观看| 国产精品一区二区性色av| 国产精品一区二区三区四区久久| 一区二区三区免费毛片| 久久久久久久久久久丰满| 在线观看66精品国产| a级一级毛片免费在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲自偷自拍三级| 在现免费观看毛片| 99热精品在线国产| 欧美一区二区精品小视频在线| 91麻豆精品激情在线观看国产| 国内精品宾馆在线| 国语自产精品视频在线第100页| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 3wmmmm亚洲av在线观看| 久久精品夜色国产| 一进一出抽搐gif免费好疼| 毛片一级片免费看久久久久| 不卡视频在线观看欧美| 美女xxoo啪啪120秒动态图| 国产成人a区在线观看| 成人国产麻豆网| 国产日本99.免费观看| 欧美激情国产日韩精品一区| 久久精品综合一区二区三区| 99热这里只有是精品50| 看非洲黑人一级黄片| 一级黄片播放器| .国产精品久久| 国产精品爽爽va在线观看网站| 精品99又大又爽又粗少妇毛片| 亚洲精品自拍成人| 久99久视频精品免费| 亚洲av二区三区四区| 国产精品女同一区二区软件| 欧美最黄视频在线播放免费| av国产免费在线观看| 国产黄色视频一区二区在线观看 | 国产伦一二天堂av在线观看| 精品人妻视频免费看| 亚洲最大成人手机在线| kizo精华| 亚洲精品日韩av片在线观看| 能在线免费观看的黄片| 成年版毛片免费区| 久久久久久伊人网av| 在线免费观看不下载黄p国产| 天天一区二区日本电影三级| 久久精品影院6| 免费人成视频x8x8入口观看| 国产高清三级在线| 亚洲乱码一区二区免费版| 国产片特级美女逼逼视频| 成人高潮视频无遮挡免费网站| 少妇被粗大猛烈的视频| 桃色一区二区三区在线观看| 欧美日韩一区二区视频在线观看视频在线 | 少妇高潮的动态图| 欧美性猛交╳xxx乱大交人| 亚洲在线观看片| 日韩精品有码人妻一区| 桃色一区二区三区在线观看| 国产极品精品免费视频能看的| 一个人免费在线观看电影| 国内久久婷婷六月综合欲色啪| 成年女人看的毛片在线观看| 最近视频中文字幕2019在线8| 日韩一区二区三区影片| 人妻夜夜爽99麻豆av| 国产精品女同一区二区软件| 少妇丰满av| 欧美色视频一区免费| 国产精品国产高清国产av| 夜夜爽天天搞| 插逼视频在线观看| 1024手机看黄色片| 亚洲欧美精品专区久久| 在线播放国产精品三级| 亚洲av免费在线观看| 看片在线看免费视频| 国产一区二区三区在线臀色熟女| 赤兔流量卡办理| 超碰av人人做人人爽久久| 看免费成人av毛片| 国产亚洲精品久久久com| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲网站| 色噜噜av男人的天堂激情| 国产综合懂色| 淫秽高清视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 男女视频在线观看网站免费| 国产成人精品久久久久久| 非洲黑人性xxxx精品又粗又长| 国产精品福利在线免费观看| 色尼玛亚洲综合影院| 亚洲一级一片aⅴ在线观看| 天堂av国产一区二区熟女人妻| 欧美性猛交黑人性爽| 久久鲁丝午夜福利片| 久久久久九九精品影院| 国内少妇人妻偷人精品xxx网站| 99热全是精品| 成人国产麻豆网| 国内精品久久久久精免费| 久久久精品94久久精品| www.av在线官网国产| 99久久精品国产国产毛片| 日韩精品青青久久久久久| 成人鲁丝片一二三区免费| 国产蜜桃级精品一区二区三区| 日本一本二区三区精品| 国产精品精品国产色婷婷| 成人高潮视频无遮挡免费网站| 精品久久久久久久末码| 中文字幕av成人在线电影| 精品日产1卡2卡| 美女国产视频在线观看| 免费一级毛片在线播放高清视频| 99九九线精品视频在线观看视频| 少妇高潮的动态图| 麻豆成人午夜福利视频| 五月伊人婷婷丁香| 成人综合一区亚洲| 免费看美女性在线毛片视频| 一个人看的www免费观看视频| 免费观看的影片在线观看| 日日啪夜夜撸| 色尼玛亚洲综合影院| 波多野结衣高清作品| 久久亚洲国产成人精品v| 免费人成视频x8x8入口观看| 天堂中文最新版在线下载 | 亚洲av二区三区四区| 在线免费十八禁| 18禁黄网站禁片免费观看直播| 久久精品国产鲁丝片午夜精品| 精品午夜福利在线看| 国产亚洲精品av在线| 欧美在线一区亚洲| 午夜福利高清视频| 国产精品久久久久久久久免| 精品午夜福利在线看| 伦精品一区二区三区| 一边摸一边抽搐一进一小说| 久久热精品热| 亚洲精品乱码久久久v下载方式| 国产高清视频在线观看网站| 日韩国内少妇激情av| 能在线免费看毛片的网站| 神马国产精品三级电影在线观看| 午夜福利成人在线免费观看| 免费不卡的大黄色大毛片视频在线观看 | 99久国产av精品国产电影| 村上凉子中文字幕在线| 日本在线视频免费播放| 久久精品影院6| 中文字幕av在线有码专区| 在线观看美女被高潮喷水网站| 在线播放无遮挡| 亚洲av一区综合| 麻豆成人av视频| 亚洲国产精品成人综合色| 国产伦理片在线播放av一区 | 欧美人与善性xxx| 久久国内精品自在自线图片| 成人三级黄色视频| 乱码一卡2卡4卡精品| 一级二级三级毛片免费看| 你懂的网址亚洲精品在线观看 | 成人永久免费在线观看视频| 亚洲五月天丁香| 两个人的视频大全免费| 特大巨黑吊av在线直播| 麻豆国产97在线/欧美| 国产黄片美女视频| 韩国av在线不卡| 久久精品夜夜夜夜夜久久蜜豆| 精品国内亚洲2022精品成人| 亚洲av第一区精品v没综合| 你懂的网址亚洲精品在线观看 | 国产精品99久久久久久久久| 久久久色成人| 26uuu在线亚洲综合色| 天堂av国产一区二区熟女人妻| av视频在线观看入口| 亚洲一区二区三区色噜噜| 女人被狂操c到高潮| 国产精品一及| 搞女人的毛片| 久久午夜亚洲精品久久| 村上凉子中文字幕在线| 观看免费一级毛片| 国产黄片视频在线免费观看| 欧美一区二区精品小视频在线| 特大巨黑吊av在线直播| 成年女人永久免费观看视频| 三级男女做爰猛烈吃奶摸视频| 狂野欧美白嫩少妇大欣赏| 日本色播在线视频| 亚洲不卡免费看| 久久精品夜夜夜夜夜久久蜜豆| 校园春色视频在线观看| 精品国产三级普通话版| 国产又黄又爽又无遮挡在线| 亚洲精品粉嫩美女一区| 久久久久久久久久成人| 午夜福利在线观看吧| 精品熟女少妇av免费看| 少妇人妻精品综合一区二区 | 久久久精品94久久精品| 天天一区二区日本电影三级| 亚洲五月天丁香| 女同久久另类99精品国产91| 毛片女人毛片| 99国产极品粉嫩在线观看| 校园春色视频在线观看| 日韩欧美一区二区三区在线观看| 国产大屁股一区二区在线视频| 国产一区亚洲一区在线观看| 天堂av国产一区二区熟女人妻| 乱人视频在线观看| 白带黄色成豆腐渣| 国产片特级美女逼逼视频| 国产午夜精品久久久久久一区二区三区| 日韩欧美一区二区三区在线观看| 国产一级毛片在线| 高清日韩中文字幕在线| 国产精品电影一区二区三区| 狂野欧美激情性xxxx在线观看| 在线播放无遮挡| 一级毛片久久久久久久久女| 日本成人三级电影网站| 在线免费观看不下载黄p国产| 2021天堂中文幕一二区在线观| 麻豆av噜噜一区二区三区| 男女下面进入的视频免费午夜| 日本三级黄在线观看| 久久精品国产清高在天天线| 国内精品久久久久精免费| 欧美性猛交黑人性爽| 九草在线视频观看| 午夜福利视频1000在线观看| 久久精品人妻少妇| 中文字幕免费在线视频6|