城市天然氣管道系統故障傳播影響評估方法及系統與流程

            文檔序號:11403556閱讀:301來源:國知局
            城市天然氣管道系統故障傳播影響評估方法及系統與流程
            本發明實施例涉及可靠性工程
            技術領域
            ,具體涉及城市天然氣管道系統故障傳播影響評估方法及系統。
            背景技術
            :隨著科技的發展,天然氣的使用越來越多,現在大部分分城市都使用天然氣作為能源,為人們的生活帶來了極大的方便。城市天然氣管道系統是接收外部氣源、向市區及郊區村鎮居民供應生活用氣、企業熱源或化工原料的重要樞紐。盡早的發現城市天然氣管道系統中的故障、定位故障、探究故障傳播的規律對于天然氣輸送有著重要的意義。現有技術中,對于城市天然氣管道系統故障方面的研究主要集中在如何分析和建立故障發生對整條管道中天然氣運行參數的影響。且多是針對一條管道進行研究,從局部角度解決問題,如將管道進行分段,分析其風險發展趨勢。但是,現有技術中大多是只考慮城市天然氣管道系統故障對局部的影響,但是因為天然氣管道是相互連通的,一處故障會對整個系統造成影響,如果僅僅考慮管道系統故障對局部的影響,這對城市天然氣管道系統的故障分析、故障傳播影響評估以及故障監測都會造成影響。因此,如何提出一種方案,能夠提高城市天然氣管道系統故障傳播影響評估的準確性,成為亟待解決的問題。技術實現要素:針對現有技術中的缺陷,本發明實施例提供一種城市天然氣管道系統故障傳播影響評估方法及系統。一方面,本發明實施例提出一種城市天然氣管道系統故障傳播影響評估方法,包括:根據城市天然氣管道系統建立所述城市天然氣管道系統的網絡結構拓撲圖,所述網絡結構拓撲圖包括:節點和表示天然氣管道的邊,其中所述節點和邊用來表示所述城市天然氣管道系統的連接關系,所述邊帶有表示天然氣流動的方向;根據所述網絡結構拓撲圖建立網絡模型,所述網絡模型包括節點信息表、邊信息表以及所述節點的連接矩陣;根據所述網絡模型中的所述節點信息表、所述邊信息表和所述連接矩陣以及預設的約束條件,獲取所述城市天然氣管道系統正常狀態對應的第一流量壓力分配信息;在所述網絡模型中模擬系統故障,并根據預設的故障傳播規則和所述約束條件,獲取所述城市天然氣管道系統在故障狀態下對應的第二流量壓力分配信息;比較所述第一流量壓力分配信息和所述第二流量壓力分配信息對所述城市天然氣管道系統進行故障傳播影響評估。另一方面,本發明實施例提供一種城市天然氣管道系統故障傳播影響評估系統,包括:拓撲圖建立模塊,用于根據城市天然氣管道系統建立所述城市天然氣管道系統的網絡結構拓撲圖,所述網絡結構拓撲圖包括:節點和表示天然氣管道的邊,其中所述節點和邊用來表示所述城市天然氣管道系統的連接關系,所述邊帶有表示天然氣流動的方向;網絡模型建立模塊,用于根據所述網絡結構拓撲圖建立網絡模型,所述網絡模型包括節點信息表、邊信息表以及所述節點的連接矩陣;第一流量壓力分配獲取模塊,用于根據所述網絡模型中的所述節點信息表、所述邊信息表和所述連接矩陣以及預設的約束條件,獲取所述城市天然氣管道系統正常狀態對應的第一流量壓力分配信息;第二流量壓力分配獲取模塊,用于在所述網絡模型中模擬系統故障,并根據預設的故障傳播規則和所述約束條件,獲取所述城市天然氣管道系統在故障狀態下對應的第二流量壓力分配信息;故障傳播影響評估模塊,用于比較所述第一流量壓力分配信息和所述第二流量壓力分配信息對所述城市天然氣管道系統進行故障傳播影響評估。本發明實施例提供的城市天然氣管道系統故障傳播影響評估方法及系統,通過建立城市天然氣管道系統的網絡模型,獲取城市天然氣管道系統正常狀態下的流量和壓力分配情況,再進行系統故障的模擬,獲取城市天然氣管道系統故障后的流量和壓力分配情況。將正常狀態的流量壓力的分配情況與故障狀態的流量壓力的分配情況進行比較,進一步進行城市天然氣管道系統的故障傳播影響評估。實現了天然氣管道局部故障對整個天然氣管道的影響的評估,提高了城市天然氣管道系統故障傳播影響評估的準確性。附圖說明圖1為本發明實施例中城市天然氣管道系統故障傳播影響評估方法的流程示意圖;圖2為本發明實施例中網絡結構拓撲圖的結構示意圖;圖3為本發明實施例中城市天然氣管道系統故障傳播影響評估方法的流程示意圖;圖4為本發明實施例中又一網絡結構拓撲圖的結構示意圖;圖5為本發明實施例中故障判據模型結構示意圖;圖6為本發明實施例中又一網絡結構拓撲圖的結構示意圖;圖7為本發明實施例中調壓室故障傳播流程示意圖;圖8為本發明實施例中城市天然氣管道系統故障傳播影響評估系統的結構示意圖。具體實施方式為使本發明實施例的目的、技術方案和優點更加清楚,下面將結合本發明實施例中的附圖,對本發明實施例中的技術方案進行清楚地描述,顯然,所描述的實施例是本發明一部分實施例,而不是全部的實施例。基于本發明中的實施例,本領域普通技術人員在沒有做出創造性勞動前提下所獲得的所有其他實施例,都屬于本發明保護的范圍。圖1為本發明實施例中城市天然氣管道系統故障傳播影響評估方法的流程示意圖,如圖1所示,本發明實施例提供的城市天然氣管道系統故障傳播影響評估方法包括:s1、根據城市天然氣管道系統建立所述城市天然氣管道系統的網絡結構拓撲圖,所述網絡結構拓撲圖包括:節點和表示天然氣管道的邊,其中所述節點和邊用來表示所述城市天然氣管道系統的連接關系,所述邊帶有表示天然氣流動的方向;具體地,本發明實施例根據城市天然氣管道系統建立出該城市的天然氣管道網絡結構拓撲圖,即用節點和邊來表示城市天然氣管道系統的連接關系以及天然氣的流動方向。其中網絡結構拓撲圖中的節點表示城市天然氣管道系統中的調壓室或閥門,邊表示城市天然氣管道系統中的管道,并且每條邊帶有表示城市天然氣管道系統中的天然氣的流動方向。具體可以將網絡結構拓撲圖中的節點和邊分別進行順序編號。s2、根據所述網絡結構拓撲圖建立網絡模型,所述網絡模型包括節點信息表、邊信息表以及所述節點的連接矩陣;具體地,建立好城市天然氣管道系統的網絡結構拓撲圖之后,根據網絡結構拓撲圖建立該城市的天然氣的網絡模型,其中網絡模型中包括節點信息表、邊信息表以及節點的連接矩陣。其中節點信息表中可以包括節點對應的類型、等級、消耗流量以及節點對應的供給需求關系;邊信息表可以包括各個邊對應的等級、起點、終點以及該邊代表的管道的長度和直徑。其中,節點的類型可以根據節點的具體作用進行劃分;節點和邊對應的等級可以根據管道壓力進行劃分,管道壓力越大,等級越高,管道壓力越小,等級越小;節點的消耗流量可以根據城市天然氣管道系統的具體情況進行設置;節點對應的供給需求關系是指該節點是向其他節點供給天然氣或需要其他節點向該節點供給天然氣,即天然氣在該節點處時是輸出還是輸入;邊的起點和終點根據該邊對應的天然氣的流動方向進行設置。其中,連接矩陣是利用矩陣的方式表示節點之間的連接關系,本發明實施例中將定義連接矩陣為a=[aij],其中,i表示網絡結構拓撲圖中的節點i,j表示網絡結構拓撲圖中的管道j。具體地,若節點i不在管道j上,則對應在連接矩陣中的元素aij表示為0;若節點i在管道j上,并且為管道j的起點,則對應在連接矩陣中的元素aij表示為-1;若節點i在管道j上,并且為管道j的終點,則對應在連接矩陣中的元素aij表示為1。如:圖2為本發明實施例中網絡結構拓撲圖的結構示意圖,如圖2所示,圖中圓圈表示節點,圓圈中的字母表示節點編號,帶箭頭的線段表示邊,即表示管道,箭頭表示天然氣的流動方向。本發明實施例中城市天然氣管道的網絡結構拓撲圖中有6個節點a、b、c、d、e、f,有6條邊ab、bc、cd、ae、df、ef,其中箭頭表示天然氣在管道中的流向,則其中連接矩陣可以表示為如下公式(1):s3、根據所述網絡模型中的所述節點信息表、所述邊信息表和所述連接矩陣以及預設的約束條件,獲取所述城市天然氣管道系統正常狀態對應的第一流量壓力分配信息;具體地,建立城市天然氣管道系統的網絡模型后,根據該網絡模型中的節點信息表、邊信息表以及預設的約束條件,獲取城市天然氣管道系統正常狀態對應的第一流量壓力分配信息,即計算城市天然氣管道系統正常狀態對應的流量和壓力分配情況,具體可以計算出各個節點對應的壓力和各個管道對應的流量。其中本發明實施例中的約束條件主要是根據質量守恒定律、節點流量平衡和環網壓力平衡進行建立的,具體可以表示為如下公式(2):式中:∑input——表示整個城市天然氣管道系統輸入的天然氣總量,即從上級調壓室輸出的天然氣總量,其數值根據需求設置,也可以參照根據《城鎮天然氣設計規范》(gb50028-93)設置;∑output——表示輸出的天然氣總量,即實際使用的需求量;∑storage——表示貯存在天然氣管道中的天然氣總量,在使用過程中忽略為零;a——表示節點的連接矩陣;q——表示管道流量矩陣;q——表示節點流量矩陣;p——表示節點壓力矩陣;δp——表示管壓降矩陣。其中q、q、p以及δp可以從網絡模型中的節點信息表以及邊信息表中獲得。具體獲取第一流量壓力分配信息可以根據建立的網絡模型和預設的約束條件,利用水利計算方程,計算出網絡模型中各個節點對應的壓力以及各個管道對應的流量,獲取到城市天然氣管道系統在正常狀態下的第一流量壓力分配信息。其中計算過程中所需的參數如天然氣壓縮因子、天然氣相對密度和天然氣溫度等,可以根據實際情況進行設置。s4、在所述網絡模型中模擬系統故障,并根據預設的故障傳播規則和所述約束條件,獲取所述城市天然氣管道系統在故障狀態下對應的第二流量壓力分配信息;具體地,在網絡模型建立好后,在網絡模型中模擬系統故障,如隨機選取其中某個節點或管道假設其出現故障,根據預設的故障傳播規則和約束條件,獲取城市天然氣管道系統在故障狀態下對應的第二流量壓力分配信息。可以根據實際管道系統出現故障的數據,通過將模型中某個節點或管道的壓力或流量改為非正常狀態,模擬該節點或管道出現故障。本發明實施例主要是通過在網絡模型中注入系統故障后,根據預設的故障傳播規則,修改網絡模型,如:修改網絡模型中的節點信息表中各節點的消耗流量或節點的連接關系等,或修改邊信息表中的相關信息。當網絡模型改變后,對應連接矩陣、節點的流量矩陣會相應的改變,進一步約束條件會相應的改變,可以計算出系統故障狀態對應的城市天然氣管道系統的流量和壓力分配情況,即獲取到第二流量壓力分配信息。s5、比較所述第一流量壓力分配信息和所述第二流量壓力分配信息對所述城市天然氣管道系統進行故障傳播影響評估。具體地,計算出正常狀態下整個城市天然氣管道系統對應的第一流量壓力分配信息,和系統故障狀態時整個城市天然氣管道系統對應的第二流量壓力分配信息后,將第一流量壓力分配信息和第二流量壓力分配信息進行比較,可以分析出系統故障對整個城市天然氣管道系統的影響,進一步進行城市天然氣管道系統的故障傳播影響評估。如:可以分別評估每個管道或節點出現故障對整個城市天然氣管道系統的影響,進一步評估該管道或節點的重要性,以及整個城市天然氣管道系統的可靠性。其中網絡模型可以通過計算機軟件掃描,將網絡模型數字化,后期的流量和壓力的分配情況的計算以及故障傳播影響的評估也可以通過計算機軟件進行。圖3為本發明實施例中城市天然氣管道系統故障傳播影響評估方法的流程示意圖,如圖3所示,本發明實施例中的城市天然氣管道系統故障傳播影響評估方法包括:t1、建立城市天然氣系統多層級網絡模型。即根據城市天然氣管道建立城市天然氣管道系統的網絡結構拓撲圖,再根據建立的網絡結構拓撲圖建立城市天然氣系統多層級網絡模型。t2、多層級網絡模型仿真。即根據建立的多層級網絡模型和約束條件,獲取正常狀態下城市天然氣管道系統的第一流量壓力分配信息。t3、故障注入。即在建立的多層級網絡模型中進行模擬故障。t4、故障傳播模型仿真。即根據故障傳播規則和約束條件獲取故障狀態下城市天然氣管道系統的第二流量壓力分配信息。t5、故障傳播影響分析。將第一流量壓力分配信息和第二流量壓力分配信息進行比較分析,對城市天然氣管道系統的故障傳播影響進行評估。如:可以分別評估每個管道或節點出現故障對整個城市天然氣管道系統的影響,進一步評估該管道或系統的重要性,以及整個城市天然氣管道系統的可靠性。本發明實施例提供的城市天然氣管道系統故障傳播影響評估方法,通過建立城市天然氣管道系統的網絡模型,獲取城市天然氣管道系統正常狀態下的流量和壓力分配情況,再進行系統故障的模擬,獲取城市天然氣管道系統故障后的流量和壓力分配情況。將正常狀態的流量壓力的分配情況與故障狀態的流量壓力的分配情況進行比較,進一步進行城市天然氣管道系統的故障傳播影響評估。實現了天然氣管道局部故障對整個天然氣管道的影響的評估,提高了城市天然氣管道系統故障傳播影響評估的準確性。在上述實施例的基礎上,所述故障傳播規則包括:故障類型和所述故障類型對應的網絡模型修改策略。具體地,本發明實施例根據城市天然氣管道系統的實際使用情況,預先設置故障傳播規則,如設置故障類型以及各故障類型對應的網絡模型修改策略。其中網絡模型修改策略可以根據城市天然氣管道系統實際出現故障后整個管道系統的變化進行設置如:網絡模型中的節點或邊是否需要移除、是否需要新增節點或邊、負荷如何進行再次分配等,當然還可以包括其他網絡模型的修改策略,本發明實施例不作具體限定。當在模型中進行模擬系統故障后,根據預設的故障傳播規則修改網絡模型,如:修改連接矩陣或節點信息表中的節點消耗流量,進一步改變約束條件,計算在故障狀態下城市天然氣管道系統放入流量和壓力分配情況。本發明實施例提供的城市天然氣管道系統故障傳播影響評估方法,通過預先設置故障傳播規則,在建立好的網絡模型中模擬系統故障,根據故障傳播規則修改網絡模型,獲取城市天然氣管道系統在故障狀態下對應的流量和壓力分配情況。實現了天然氣管道局部故障對整個天然氣管道的影響的評估,提高了城市天然氣管道系統故障傳播影響評估的準確性。在上述實施例的基礎上,所述故障類型包括:調壓室故障。具體地,在故障傳播規則中設置的故障類型可以包括城市天然氣管道系統各組成部件的不同故障模式,如調壓室故障,調壓室故障是指城市天然氣管道系統中放入調壓器出現故障;當然故障類型還可以包括其他情況,具體根據實際需要進行設置,本發明實施例不作具體限定。設置好故障類型后,根據實際管道系統出現故障時對應的情況,設置每個故障類型對應的網絡模型修改策略,具體網絡模型修改策略的同上述實施例一致,此處不再贅述。本發明實施例提供的城市天然氣管道系統故障傳播影響評估方法,根據城市天然氣管道系統在實際應用時出現的故障類型設置故障類型以及對應的網絡模型修改策略,進一步獲取故障狀態對應的第二流量壓力信息。實現了局部管道系統故障對整個管道系統的故障影響評估,提高了城市天然氣管道系統故障傳播影響評估的準確性,并可以進行整個城市天然氣管道系統的可靠性評估,為城市天然氣管道系統的設計提供依據。在上述實施例的基礎上,所述方法還包括:根據所述節點的作用將所述節點進行分類,所述節點的類型包括:源節點、第一類中轉節點、第二類中轉節點、第三類中轉節點以及目標節點。具體地,根據各個節點的作用,將節點進行分類,本發明實施例中節點的類型包括:源節點、第一類中轉節點、第二類中轉節點、第三類中轉節點以及目標節點。其中源節點是指整個城市天然氣管道系統中,長輸管線進入城市處的調壓器,源節點的作用是提供整個管道系統中天然氣的來源;第一類中轉節點在管道系統中代表的是不同壓力等級管網之間的調壓器,它的作用是降壓分配,把高壓天然氣的壓力調低;第二類中轉節點在管道系統中代表的是管道之間的閥門,它的作用是進行天然氣的分配,改變天然氣輸出路徑;第三類中轉節點在管道系統中代表的是某些用氣區的調壓器,它的作用相當于第一類中轉節點和第二類中轉節點作用的綜合,它除了會把上游來氣供給用戶需求外,也把剩余天然氣供給下游用戶;目標節點在管道系統中代表的是某些用氣區的調壓器,它作用是調壓供給用戶。每一類型的節點的壓力和流量的計算方法不同,相應的節點之間的管道的流量和壓力的計算方法也會不同,可以根據節點的類型設置其對應的壓力和流量的計算方法,進一步可以計算出節點之間管道的計算方法。此外,本發明實施例中節點的類型對建立預設條件的建立也有影響,例如:約束條件中的質量守恒定律中的輸入管道系統中的天然氣總量等于輸出的天然氣總量以及管道貯存的天然氣總量,其中管道系統中的天然氣的輸入和輸出可以根據節點類型進行計算。根據城市天然氣管道系統設置各個節點的連接關系,城市天然氣管道系統中的流向,各個節點的類型,以及根據管道系統中管道的壓力設置個各個節點以及邊的等級,即建立城市天然氣管道的網絡結構拓撲圖。圖4為本發明實施例中又一網絡結構拓撲圖的結構示意圖,如圖4所示,本發明實施例中包括3個等級,1個源節點,4個第一類中轉節點,7個第二類中轉節點,2個第三類中轉節點和1個目標節點,其中邊的箭頭方向表示天然氣在管道中的流向。表1為本發明實施例中根據圖2對應的網絡結構拓撲圖中建立的節點信息表,表2為本發明實施例中根據圖2對應的網絡結構拓撲圖中建立的邊信息表,如表1和表2所示,可以看出圖2中的網絡結構拓撲圖中的所有的節點和邊屬于同一等級n,當然節點信息表和邊信息表中還可以包括其他信息,如節點壓力等,本發明實施例不作具體限定。表1網絡層級節點編號節點類型消耗流量供給/需求na第一類中轉節點-200供給nb第二類中轉節點————nc目標節點+100需求nd第三類中轉節點+100需求ne第二類中轉節點————nf第二類中轉節點————-其中:因為在圖2中只能夠顯示一個等級的模型,所以表1中有些節點的供給和需求關系以及消耗流量和暫時不能確定,在節點信息表中用“-”表示。表2網絡層級邊編號起點終點長度直徑nabablabdabnbcbclbcdbcndcdcldcddcnadaeladdadnefeflefdefnfdfdlfddfd本發明實施例提供的城市天然氣管道系統故障傳播影響評估方法,根據各節點的作用將節點進行分類,每一類型的節點對應不同的消耗流量以及節點壓力,將各節點的類型以及對應的消耗流量等,存儲在節點信息表中,以便后期計算管道系統的流量和壓力分配情況,提高了城市天然氣管道系統故障傳播影響評估的準確性。在上述實施例的基礎上,所述方法還包括:設置故障判據,在所述網絡模型中模擬系統故障,獲取所述管道系統在故障狀態下的對應的流量壓力分配信息后,若根據所述故障判據和所述流量壓力分配信息判斷獲知所述管道系統處于穩定狀態,則將所述管道系統狀態信息作為所述第二流量壓力分配信息。具體地,在模擬系統故障時,在注入系統故障后,根據故障傳播規則以及約束條件,獲取管道系統在故障狀態下對應的流量壓力分配信息。由于注入故障,可能導致管道系統不穩定,本發明實施例通過設置故障判據,根據故障判據判斷當前狀態下對應的管道系統的流量和壓力的分配是否處于穩定狀態,若是,則將該管道系統狀態信息作為第二流量壓力分配信息。將第二流量壓力分配信息與第一流量壓力分配信息進行比較,進行管道系統的故障傳播影響評估。本發明實施例中的故障判據是指設置流量的下限閾值和上限閾值,以及壓力的下限閾值和上限閾值,只有當流量和壓力都在對應的上限值和下限值之間時,表示系統處于穩定狀態。其中流量和壓力的上下限閾值可以采用如下公式(3)進行定義:式中:c——表示容量;l——表示負荷,負荷是指網絡單元所承受的應力,即為天然氣流量和壓力,都作為觀測參數。α——表示余量參數,一般情況下α>0;j——表示節點編號。圖5為本發明實施例中故障判據模型結構示意圖,如圖5所示,橫坐標軸表示流量,縱坐標軸表示壓力,只有當流量和壓力都在矩形范圍內,表示該系統處于正常狀態,否則認為處于故障狀態。實際應用時,可以依次判斷網絡模型中的節點以及邊是否處于正常狀態,當所有的節點和邊都處于正常狀態時,則認為管道系統處于穩定狀態。若判斷網絡模型中有些節點或邊處于故障狀態,則根據故障傳播規則對相應的節點或邊進行處理,如移除或增加等,即修改網絡模型,重新進行天然氣的分配,計算天然氣管道的流量和壓力分配情況,再根據故障判據進行判斷,直至系統處于穩定狀態。本發明實施例提供的城市天然氣管道系統故障傳播影響評估方法,通過建立城市天然氣管道系統的網絡模型,獲取城市天然氣管道系統正常狀態下的流量和壓力分配情況,再進行系統故障的模擬,獲取城市天然氣管道系統故障后的流量和壓力分配情況。將正常狀態的流量壓力的分配情況與故障狀態的流量壓力的分配情況進行比較,進一步進行城市天然氣管道系統的故障傳播影響評估。實現了天然氣管道局部故障對整個天然氣管道的影響的評估,提高了城市天然氣管道系統故障傳播影響評估的準確性。下面結合故障類型中的調壓室故障介紹本發明實施例的技術方案,在建立網絡模型,獲取到管道系統的第一流量壓力信息后,進行模擬系統故障,本發明實施例中模擬的是調壓室故障。調壓室的出口壓力達不到指定要求,系統都處于降量運行的狀態,但不會因為壓力不足而關斷某些閥門。表3網絡層級節點編號節點類型消耗流量供給/需求11目標節點+150需求12第二類中轉節點0——13第三類中轉節點+120-需求14第三類中轉節點+100需求15第三類中轉節點+200需求16第二類中轉節點+150——27第二類中轉節點0——28目標節點+100需求29第三類中轉節點+400需求210第三類中轉節點+100需求211第二類中轉節點0——212源節點-1200供給圖6為本發明實施例中又一網絡結構拓撲圖的結構示意圖,如圖6所示,圖中圓圈表示節點,圓圈中的數字表示節點編號,帶箭頭的線段表示邊,即表示管道,箭頭表示天然氣的流動方向。如圖6所示,本發明實施例中包括14條邊和12個節點,其節點信息表如表3所示,邊信息表如表4所示。根據建立的網絡模型和約束條件獲取到正常狀態下,該城市天然氣管道系統的流量和壓力分配信息,其中正常狀態下節點的壓力分配表如表5所示,管道的流量分配表如表6所示:表4網絡層級邊編號起點終點長度直徑113680030012431200250136111003001423120025015548003001621900250175295025028751002502978950300210107120030021198800350212121190035021311912003502141210800350表5節點編號壓力(pa)11331.30821358.54731331.70741337.44651624.23261331.52971912.23281894.85791893.062101956.917111992.539122000表6邊編號流量(nm3/h)132.979267.444332.979465.5345267.4446117.0207182.55586009147.92010747.9201147.92012352.07913352.07914847.920圖7為本發明實施例中調壓室故障傳播流程示意圖,如圖7所示,本發明的提供的調壓室故障傳播包括以下步驟:r1、定位故障調壓室。即假設網絡模型中的一個表示調壓室的節點出現故障,可以通過輸入非正常的該節點的流量值或壓力值實現。r2、計算調壓室故障后調壓室對應的壓力降。其中調壓室故障后的壓力降可以根據如下公式(5)~(7)進行計算:δp=p(0)-p(1)(6)式中:q——故障后管道的流量初始值,可以根據水利計算方程進行修改;△p——調壓室故障后的壓力降;s——為管道的阻抗,該參數與流體的種類、流動狀態、已知環境條件和管道的基本參數有關;p(0)——正常情況下的節點基準壓力;p(1)——根據連接關系得到的上一層網絡連接處的節點壓力;λ——天然氣管道的摩擦阻力系數;ρ0——標準狀態下的天然氣密度(0.73kg/nm3);p0——標準大氣壓,p0=101325pa;t0——標準狀態下的絕對溫度(273.15k);t——燃氣絕對溫度(k);z0——標準狀態下的氣體壓縮因子;z——氣體壓縮因子;l——燃氣管道的長度(m)d——天然氣管道內經(mm);故障發生點的上一層網絡的基準壓力不變,改變的是該節點的流量,采用如下公式(8)可以計算調壓室故障后,各節點的流量值:q——節點的流量;qin——下一層網絡節點流量。r3、修改節點流量矩陣和管道流量矩陣。根據計算出的調壓室故障后的壓力降,以及各節點的流量,修改網絡模型中節點信息表和邊信息表,根據節點信息表和邊信息表,修改節點流量矩陣和管道流量矩陣。r4、根據約束條件計算本層流量壓力分配情況。根據修改后的節點流量矩陣和管道流量矩陣的約束條件,基于導納系數分配算法,計算該層的流量和壓力分配情況,并存儲。r5、判斷下層網絡是否計算完成。即判斷是否將網絡模型中的所有層級的系統的流量壓力分配情況是否全部計算出來,若是,執行r8,若否,執行r6。r6、根據連接關系計算下一層網絡的基準壓力。故障發生點的上一層網絡的基準壓力不變,改變的是該節點的流量。下一層網絡的基準壓力具體可以根據公式(5)~(7)進行計算:r7、計算該下一層網絡的壓力降。計算下一次級的壓力降后,返回r3,直至所有層級計算完。r8、輸出故障后整個系統的流量壓力分配信息。表7節點編號正常壓力(pa)降量運行壓力降量百分比%11331.308831.30837.55721358.547858.54736.80331331.707831.70737.54541337.446837.44637.38451624.2321124.23230.78361331.529831.52937.55071912.2321412.23226.14781894.8571394.85726.38791893.0621393.06226.412101956.9171456.91725.550111992.5391492.53925.093122000150025故障注入后,獲取到管道系統在調壓室故障后對應的第二流量壓力分配信息,將第一流量分配信息和第二流量壓力分配信息進行比較,表7為本發明實施例中節點壓力對比表,表8為本發明實施例中管道流量對比表,具體如表7和表8所示:表8邊編號流量(nm3/h)降量運行流量降量百分比%132.97929.9749.113267.44446.86430.51332.97918.72343.228465.53431.78451.5005267.44478.11070.7936117.02070.83339.4697182.555131.95127.7198600386.13235.6449147.920146.6770.51110747.920741.9690.7951147.92012.21074.52012352.079256.2427.22113352.079344.002.29414847.920329.09161.18由數據可知,節點的壓力降低較多,都在25%以上;管道的流量降低不均。在這種降量運行的條件下,管網無法滿足所有節點的流量、壓力需求,網絡幾乎癱瘓。管道系統的性能可靠度也大大降低。所以在進行管道維護和管理時,應當特別注意該故障的發生,降低該影響的波及范圍。本發明實施例提供的城市天然氣管道系統故障傳播影響評估方法,可以評估和分析局部的系統故障對整個管道系統的影響,可以得出不同節點或管道的重要度,進行系統可靠性設計和冗余設計等。同時可以根據大量的仿真數據在模型中定位故障。根據故障傳播評估的大量的故障數據,可以建立數據庫,作為故障字典。如果發生故障,通過對不同位置的監測信號比對和分析,搜索故障字典即可定位故障源,分析故障模式,采取相應的解決方案。圖8為本發明實施例中城市天然氣管道系統故障傳播影響評估系統的結構示意圖,如圖8所示,本發明實施例提供的城市天然氣管道系統故障傳播影響評估系統包括:拓撲圖建立模塊81、網絡模型建立模塊82、第一流量壓力分配獲取模塊83、第二流量壓力分配獲取模塊84和故障傳播影響評估模塊85,其中:拓撲圖建立模塊81用于根據城市天然氣管道系統建立所述城市天然氣管道系統的網絡結構拓撲圖,所述網絡結構拓撲圖包括:節點和表示天然氣管道的邊,其中所述節點和邊用來表示所述城市天然氣管道系統的連接關系,所述邊帶有表示天然氣流動的方向;網絡模型建立模塊82用于根據所述網絡結構拓撲圖建立網絡模型,所述網絡模型包括節點信息表、邊信息表以及所述節點的連接矩陣;第一流量壓力分配獲取模塊83用于根據所述網絡模型中的所述節點信息表、所述邊信息表和所述連接矩陣以及預設的約束條件,獲取所述城市天然氣管道系統正常狀態對應的第一流量壓力分配信息;第二流量壓力分配獲取模塊84用于在所述網絡模型中模擬系統故障,并根據預設的故障傳播規則和所述約束條件,獲取所述城市天然氣管道系統在所述系統故障對應的狀態下的第二流量壓力分配信息;故障傳播影響評估模塊85用于比較所述第一流量壓力分配信息和所述第二流量壓力分配信息對所述述城市天然氣管道系統進行故障傳播影響評估。具體地,拓撲圖建立模塊81根據城市天然氣管道系統建立出該城市的天然氣管道網絡結構拓撲圖,即用節點和邊來表示城市天然氣管道系統的連接關系以及天然氣的流動方向。其中網絡結構拓撲圖中的節點表示城市天然氣管道系統中的調壓室或閥門,邊表示城市天然氣管道系統中的管道,并且每條邊帶有表示城市天然氣管道系統中的天然氣的流動方向。網絡模型建立模塊82根據網絡結構拓撲圖建立該城市的天然氣的網絡模型,其中網絡模型中包括節點信息表、邊信息表以及節點的連接矩陣。其中節點信息表和邊信息表的內容,以及連接矩陣的建立方法同上述實施例一致,此處不再贅述。第一流量壓力分配獲取模塊83根據該網絡模型中的節點信息表、邊信息表以及預設的約束條件,獲取城市天然氣管道系統正常狀態對應的第一流量壓力分配信息,即計算城市天然氣管道系統正常狀態對應的流量和壓力分配情況,具體可以計算出各個節點對應的壓力和各個管道對應的流量。其中約束條件的建立同上述實時例一致,此處不再贅述。第二流量壓力分配獲取模塊84在網絡模型中模擬系統故障,如隨機選取其中某個節點或管道假設其出現故障,根據預設的故障傳播規則和約束條件,獲取城市天然氣管道系統在故障狀態下對應的第二流量壓力分配信息。故障傳播影響評估模塊85將第一流量壓力分配信息和第二流量壓力分配信息進行比較,可以分析出系統故障對整個城市天然氣管道系統的影響,進一步進行城市天然氣管道系統的故障傳播影響評估。具體的評估方法同上述實施例一致,此處不再贅述。本發明實施例提供的城市天然氣管道系統故障傳播影響評估系統,通過建立城市天然氣管道系統的網絡模型,獲取城市天然氣管道系統正常狀態下的流量和壓力分配情況,再進行系統故障的模擬,獲取城市天然氣管道系統故障后的流量和壓力分配情況。將正常狀態的流量壓力的分配情況與故障狀態的流量壓力的分配情況進行比較,進一步進行城市天然氣管道系統的故障傳播影響評估。實現了天然氣管道局部故障對整個天然氣管道的影響的評估,提高了城市天然氣管道系統故障傳播影響評估的準確性。在上述實施例的基礎上,所述故障傳播規則包括:故障類型和所述故障類型對應的網絡模型修改策略。具體地,本發明實施例根據城市天然氣管道系統的實際使用情況,預先設置故障傳播規則,如設置故障類型以及各故障類型對應的網絡模型修改策略。其中,網絡模型修改策略的設置同上述實施例一致,此處不再贅述。當在模型中進行模擬系統故障后,根據預設的故障傳播規則修改網絡模型,如:修改連接矩陣或節點信息表中的節點消耗流量,進一步改變約束條件,計算在故障狀態下城市天然氣管道系統放入流量和壓力分配情況。本發明實施例提供的城市天然氣管道系統故障傳播影響評估系統,通過預先設置故障傳播規則,在建立好的網絡模型中模擬系統故障,根據故障傳播規則修改網絡模型,獲取城市天然氣管道系統在故障狀態下對應的流量和壓力分配情況。實現了天然氣管道局部故障對整個天然氣管道的影響的評估,提高了城市天然氣管道系統故障傳播影響評估的準確性。在上述實施例的基礎上,所述故障類型包括:系統故障、調壓室故障和閥門故障。具體地,在故障傳播規則中設置的故障類型可以包括城市天然氣管道系統各組成部件的不同故障模式,如調壓室故障,調壓室故障是指城市天然氣管道系統中放入調壓器出現故障。當然故障類型還可以包括其他情況,具體根據實際需要進行設置,本發明實施例不作具體限定。設置好故障類型后,根據實際管道系統出現故障時對應的情況,設置每個故障類型對應的網絡模型修改策略,具體網絡模型修改策略的同上述實施例一致,此處不再贅述。本發明實施例提供的城市天然氣管道系統故障傳播影響評估系統,根據城市天然氣管道系統在實際應用時出現的故障類型設置故障類型以及對應的網絡模型修改策略,進一步獲取故障狀態對應的第二流量壓力信息。實現了局部管道系統故障對整個管道系統的故障影響評估,提高了城市天然氣管道系統故障傳播影響評估的準確性,并可以進行整個城市天然氣管道系統的可靠性評估,為城市天然氣管道系統的設計提供依據。在上述實施例的基礎上,所述拓撲圖建立模塊還用于:根據所述節點的作用將所述節點進行分類,所述節點的類型包括:源節點、第一類中轉節點、第二類中轉節點、第三類中轉節點以及目標節點。具體地,拓撲圖建立模塊在根據城市天然氣管道系統建立網絡結構拓撲圖時,根據各個節點的作用,將節點進行分類,本發明實施例中節點的類型包括:源節點、第一類中轉節點、第二類中轉節點、第三類中轉節點以及目標節點。其中各類型的節點的作用同上述實施例一致,此處不再贅述。本發明實施例提供的城市天然氣管道系統故障傳播影響評估系統,根據各節點的作用將節點進行分類,每一類型的節點對應不同的消耗流量以及節點壓力,將各節點的類型以及對應的消耗流量等,存儲在節點信息表中,以便后期計算管道系統的流量和壓力分配情況,提高了城市天然氣管道系統故障傳播影響評估的準確性。在上述實施例的基礎上,所述第二流量壓力分配獲取模塊還用于:設置故障判據,在所述網絡模型中模擬系統故障,獲取所述管道系統在所述故障狀態下的對應的流量壓力分配信息后,若根據所述故障判據和所述流量壓力分配信息判斷獲知所述管道系統處于穩定狀態,則將所述管道系統狀態信息作為所述第二流量壓力分配信息。具體地,第二流量壓力分配獲取在模塊模擬系統故障時,在注入系統故障后,根據故障傳播規則以及約束條件,獲取管道系統在故障狀態下對應的流量壓力分配信息。由于注入故障,可能導致管道系統不穩定,本發明實施例通過設置故障判據,根據故障判據判斷當前狀態下對應的管道系統的流量和壓力的分配是否處于穩定狀態,若是,則將該管道系統狀態信息作為第二流量壓力分配信息。將第二流量壓力分配信息與第一流量壓力分配信息進行比較,進行管道系統的故障傳播影響評估。其中故障判據的設置方法同上述實施例一致,此處不再贅述。本發明實施例中的系統用于執行上述實施例所述的方法,具體實施方式同上述實施例一致,此處不再贅述。本發明實施例提供的城市天然氣管道系統故障傳播影響評估系統,可以評估和分析局部的系統故障對整個管道系統的影響,可以得出不同節點或管道的重要度,進行系統可靠性設計和冗余設計等。同時可以根據大量的仿真數據在模型中定位故障。根據故障傳播評估的大量的故障數據,可以建立數據庫,作為故障字典。如果發生故障,通過對不同位置的監測信號比對和分析,搜索故障字典即可定位故障源,分析故障模式,采取相應的解決方案。最后應說明的是:以上實施例僅用以說明本發明的技術方案,而非對其限制;盡管參照前述實施例對本發明進行了詳細的說明,本領域的普通技術人員應當理解:其依然可以對前述各實施例所記載的技術方案進行修改,或者對其中部分技術特征進行等同替換;而這些修改或者替換,并不使相應技術方案的本質脫離本發明各實施例技術方案的精神和范圍。當前第1頁12
            當前第1頁1 2 
            網友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品