999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

深水鋼質懶波型立管技術進展及研究方向

2017-06-23 13:12:24王金龍范嘉堃中國特種設備檢測研究院壓力管道事業部北京0007中海石油氣電集團有限責任公司研發中心北京0009
石油礦場機械 2017年3期
關鍵詞:分析模型研究

王金龍,范嘉堃(.中國特種設備檢測研究院 壓力管道事業部,北京0007;.中海石油氣電集團有限責任公司 研發中心,北京0009)

深水鋼質懶波型立管技術進展及研究方向

王金龍1,范嘉堃2
(1.中國特種設備檢測研究院 壓力管道事業部,北京100017;2.中海石油氣電集團有限責任公司 研發中心,北京100029)

作為一種改良型的鋼懸鏈線立管(Steel Catenary Riser,SCR),鋼質懶波型立管(Steel Lazy-Wave Riser,SLWR)在深水油氣開發中的應用越來越廣泛。介紹了SLWR技術概況,調研了SLWR在海洋油氣開發中的應用現狀;分析了SLWR在力學分析、安裝分析等力學方面的研究進展;進行了SLWR分析的主要理論介紹,包括懸鏈線理論和大撓度梁理論;最后就SLWR技術作出了展望,可對中國SLWR技術發展起到一定的參考作用。

深水;立管;鋼質懶波型立管;大撓度梁理論;發展趨勢

鋼懸鏈線立管(Steel Catenary Riser,SCR)在使用過程中,隨著水深的增長,水中懸掛段長度越來越長,海洋環境載荷更加惡劣。特別是懸掛在大運動幅度的生產平臺上,SCR的頂部張力越來越大,對立管以及連接立管和平臺的接頭都帶來了重大的技術挑戰[1-2]。為了解決和適應此種工況下油氣的開采,SCR發展出了很多種形式:簡單懸鏈線式(Simple Catenary Riser)、懶波型立管(Lazy-wave Riser)、陡波型立管(Steep-wave Riser)、懶S型立管(Lazy S Riser)、陡S型立管(Steep S Riser)和順應波式立管(Pliant Wave Riser)[3]等形式。為了減小立管頂端的懸掛張力,其中一種方法就是在SCR中間段鋼管外側安裝多個浮筒模塊,抵抗一部分張力,成為深水中可行的選擇,這樣形成的立管就是鋼質懶波型立管(Steel Lazy-wave Riser,SLWR)。

本文介紹了SLWR技術概況,調研了SLWR在海洋油氣開發中的應用現狀;分析了SLWR在力學分析、安裝分析等力學方面的研究進展;進行了SLWR分析的主要理論介紹,包括懸鏈線理論和大撓度梁理論;最后就SLWR技術作出了展望,并提出了相應的研究建議。

1 應用現狀

鋼質懶波型立管(Steel Lazy-wave Riser,SLWR)形態如圖1所示。“懶”是指和海床接觸處的管道軸線近似平行于海床,“波”是波浪的簡稱,是用來描述由于增加了浮筒模塊而形成的立管形狀。雖然這種立管比簡單懸鏈線立管成本更高(需要額外制造及維護的浮筒和更高的安裝成本,浮筒一般為更輕質的聚乙烯材料),但能較大幅度降低立管的頂端張力,而成為在深水開發中的一種更加可行的SCR解決方案。由于浮筒的特性,增大了立管在海水中的浮力,頂端張力最大能下降50%左右。相對于簡單懸鏈線立管,SLWR在底部形成拱彎段和垂彎段,能夠將觸地點(Touchdown Point,TDP)運動和頂端運動隔離開,減小頂端載荷,提高疲勞壽命[4]。

圖1 簡單懸鏈線立管和鋼質懶波型立管結構形式

自1994年殼牌石油公司在墨西哥灣872 m水深的“Auger”張力腿平臺安裝了世界上第1條SCR[5]開始,經過20 a的發展,SCR在海上油氣田開發中得到了迅速發展和應用。目前已有上百條SCR投入使用[6],并且廣泛應用于各種水深和各種浮式平臺。

作為一種新型的、更先進的SCR,SLWR的第1次應用是在2009年位于巴西Campos Basin區域的Parque das conchas (BC-10)深水項目(如圖2所示)。在該區塊的1期工程Ostra、Abalone 和Argonauta B-west中,由殼牌石油公司在1 800 m水深處安裝了7根SLWR在轉塔系泊式FPSO(Floating Production Storage and Offloading,浮式生產儲油卸油裝置)上[7-8]。此舉是為了在深水、超深水油氣田開發中,提高立管的疲勞壽命,降低立管頂端載荷,解決SCR在這方面的不足。之后在2012年,安納達科石油公司(Anadarko)在墨西哥灣約1 515 m水深的Caesar Tonga項目(如圖3所示)中的Green Canyon 680區塊安裝了2根SLWR,并于2012-03-12投產。2根由管中管組成的SLWR回接到約1 515 m水深處的“Constitution”Spar(單柱式平臺)上,它們的內部設計壓力87.56 MPa(12 700 psi),能耐93.33 ℃高溫[9]。

圖2 Parque das conchas (BC-10)油田布局

圖3 Caesar Tonga油田布局

我國于2014-04-23投產第1個1 500 m水深的深水油氣田——荔灣3-1(LW3-1)氣田,采用水下生產系統+淺水導管架平臺的生產模式,未使用深水立管。國內還無深水SCR和SLWR的實際應用記錄,缺少相關技術積累,并且在深水SCR和SLWR的技術儲備方面,與國外深水油氣開發先進國家美國和巴西等差距明顯。隨著我國南海深水油氣開發,相信SLWR在我國也會很快得到應用。SLWR的設計技術、分析技術和安裝技術,需要進行深入詳細的研究,以便提高深水立管技術儲備水平,更好地服務于我國的深水油氣田開發戰略。

2 研究進展

SLWR分析主要有數值模擬、模型試驗及實際監測分析3種方法。模型試驗及實際監測分析成本高,而SLWR的數值模擬具有成本低、有效性高的特點。對于SLWR的動力響應、靜力學模型和安裝力學分析等研究還非常少,因此對于SLWR的一些關鍵力學問題進行詳細而全面研究是很必要的,具有重要的學術價值和工程意義。

2.1 深水SLWR力學分析

深水SLWR在海水中除了受到風浪流等環境載荷作用外,還受到底端海床土壤的約束以及頂端生產平臺運動的影響。SLWR在海水中管道長度非常長,立管頂端接近豎直,底端接近水平,而且由于安裝了浮筒,呈現出拱彎段和垂彎段2個傾角和曲率變化都非常大的彎曲形態,具有非常強的幾何非線性特征。因此,如何準確計算SLWR的形態以及受力狀態成為各國研究人員關注的重點。

國內外很多學者對SCR的靜力學形態、動力響應、渦激振動(Vortex Induced Vibration,VIV)、疲勞、管土作用、可靠性、管道鋪設安裝等特性進行了大量全面深入的研究,也取得了非常重要的成果。由于SLWR是一種特殊的SCR,SLWR和SCR在底端海床土壤約束、力學模型等諸多方面有很多相似之處,SCR的相關研究成果對于SLWR的研究有著重要的參考價值。很多研究人員參考SCR研究成果,對SLWR的力學模型提出過很多理論和方法,進行過一些非常重要的研究。

Li[10]將 SLWR管道模擬成3段簡單的懸鏈線,并且研究了管道的動力響應,但懸鏈線理論忽略了管道的彎曲剛度,不能完整反映立管的力學特性。Quéau[11-12]使用懸鏈線理論、邊界層、Winkler地基梁理論對SLWR進行了理論分析,但該文分析局限于靜態,并且海水中大部分管道采用理想懸鏈線理論模擬,無法考慮海流和管道內流載荷,有很大局限性。da Silveira[13]提出了基于理想纜繩理論的SLWR的靜態常微分方程的數學模型,對SLWR進行大量的參數分析,得到最優形態的SLWR。Campbell[14]詳敘了最優最合適懶波形態的SLWR選擇方法,提出了SLWR靜態分析、動態分析、疲勞分析的一整套方法。Yue[15]使用商業軟件OrcaFlex分析了幾種不同拱彎形式的SLWR,并進行了比較。Torres[16]對注水用的SLWR進行了技術可行性的詳細研究以及一根外輸原油的SLWR進行了完整性和可行性的研究。Petruska[17]研究了在同一平臺選擇SCR、SLWR等不同立管的相關可能性。Balena[18]詳細介紹了立管浮筒段的相關細節及浮筒段模擬成一等效單層管,并且對SLWR進行了疲勞分析。Santillan[19-20]系統研究了懶S型和陡S型立管的力學模型,并通過實驗驗證了在不同形態和邊界條件下的數值結果的正確性。懶S型與懶波型立管形態幾乎類似,不同之處僅在于懶S型立管的拱彎段浮力是通過浮筒和懸繩將大浮力加載在管道的一個點上,而SLWR是一段特定長度的浮筒段立管提供很長一段長度的豎直向上的浮力。Sun[21]采用集中質量法對陡波形立管進行了靜態和動態分析,并針對浮筒長度對立管的動態性能的影響進行了參數敏感性分析。

深水臍帶纜具有高柔性特征,但若臍帶纜以懶波形態的形式安裝布置在海水中時,其形態與SLWR形態相似,計算模型具備極大的借鑒意義。

Sun[22]使用有限元軟件ABAQUS對懶波型臍帶纜進行了參數分析。李清泉[23]對海洋環境載荷作用下懶波型臍帶纜的方位角、拖曳力系數及浮筒半徑等影響參數進行分析,分析結果表明臍帶纜位置構型與這些參數之間有很大相關性。宋磊建[24]在給出懸鏈線和懶波型理論公式的基礎上,研究了某一深海臍帶纜進行懸鏈線布局和懶波型布局2種不同布局形式的選型設計,計算得到臍帶纜在2種布局下的靜態和動力響應,及其有效張力和彎曲曲率的分布特性。Ruan[25]提出了一種考慮環境載荷、彈性海床、邊界層效應的深水懶波型臍帶纜靜力學理論模型,將懶波型臍帶纜分為懸掛在海水中的纜繩段、鄰近觸地點的邊界層段、海床段3段結構。

2.2 深水SLWR安裝分析

深水SLWR在安裝時與SCR有很多相類似的點,在安裝過程中遇到安裝工期、突發狀況等情況時,立管會從安裝船舶棄置到海底,或者從海底回收至安裝船舶的棄置回收過程,從安裝船舶轉移到生產平臺的轉移安裝過程等。SLWR由于浮筒段的存在,這些棄置回收和轉移安裝等過程非常復雜,對工程安裝提出了相當大的挑戰,而且立管安裝對于立管的服役壽命也存在影響。因此在立管安裝前,需要對管道進行大量的計算分析以保證安裝的安全性[26],選擇最佳的安裝方法與路徑。一些工程師和學者對SLWR安裝問題進行了大量的深入分析和研究。

Hoffman[27]敘述了SLWR的相關設計和安裝內容,特別研究了SLWR與平臺接頭的細節。Lahey[9]敘述了在Caesar-Tonga 油田的Green Canyon 680區塊1 514.856 m水深處安裝的2根SLWR,并且詳細描述了SLWR浮筒的安裝過程(如圖4所示)。

圖4 浮筒安裝現場

Senra[28]提出了SLWR設計面臨的挑戰,并且詳細分析了SLWR的安裝過程。Yue[29]對比研究了不同形態的SLWR,并且對SLWR使用J-lay方法安裝的可行性進行了研究。Thomas[4]介紹了SLWR的棄置回收過程以及它在安裝過程中遇到的挑戰(如圖5所示),但該文沒有具體描述棄置回收的力學模型,僅介紹了使用商業軟件計算得到的結果。

(立管由上部的牽引纜繩懸掛。從左至右為棄置過程,反向則為回收過程)圖5 SLWR棄置回收過程中立管形態

對于SLWR從安裝船舶轉移到生產平臺的轉移安裝問題,目前研究非常少。Thomas[4]介紹了SLWR在轉移安裝過程中遇到的挑戰,但沒有具體描述轉移安裝時SLWR的數學模型,僅介紹了使用商業軟件計算得到的結果。

以上是關于SLWR的安裝方面的研究,主要是一些工程上的計算探討,均未提出全面詳細的模擬SLWR鋪設、棄置回收、轉移安裝等不同安裝過程的力學模型。SLWR的靜力學、動力學、棄置回收、轉移安裝等安裝過程還需要進行深入研究,提出全面詳細的力學模型是非常有必要和有意義的。

3 主要力學模型

由于實際情況中立管的理論分析很復雜,具有復雜的受力狀況,因此需采取一定的簡化措施以及適當的狀態方程來盡可能地真實反映立管實際的工作狀況。運用純理論的方法對立管進行分析很困難。本文主要敘述目前SLWR力學計算中用到的2種主要的力學模型:懸鏈線理論模型和大撓度梁模型。

3.1 懸鏈線理論模型

SLWR(如圖6所示)的某一位置管道外側安裝上浮筒,增加此段的浮力,使得此管段浮力大于重力,因而向上拱起,形成了拱彎段。本模型為2維模型。SLWR為O-A-B-C-D-E之間的管段,以管道和海床剛接觸的觸地點TDP處標記為全局坐標系(x,y)的零點O。立管底端O點與海底管道相接,頂端E點通過柔性接頭連接在海洋平臺上,整段立管暴露在海中,受到海洋載荷及海洋平臺傳遞過來的載荷的作用。SLWR中間安裝有浮筒,由此分成3部分:下懸鏈線段O-A(0≤S≤Sl,S為從O點開始的管道弧長,Sl是下懸鏈線立管弧長),浮筒段A-B-C(Sl≤S≤Sl+Sb,Sb為浮筒段立管弧長),上懸鏈線段C-D-E(Sl+Sb≤S≤SL+Sb+Su,Su為上懸鏈線段立管弧長)。點B是拱彎段的最高點,稱為拱彎點,其高度被定義為拱彎點高度Yarch,AB段稱為拖曳段,BC段稱為舉升段。上懸鏈線段C-D-E呈懸鏈線形態,在D點附近處形成了懸垂段。D點是本段的最低點,定義為垂彎點,垂彎點高度Ysag也是SLWR的一個重要參數,CD段稱為跨接段,DE段稱為懸掛段。拱彎點與垂彎點之間高度差稱為拱彎高度,高度越大,彎曲程度越大。E點為懸掛點,它的懸掛角為立管安裝在海洋平臺上的重要參數。懸掛角βhang-off是懸掛點切線與豎直y向的夾角,懸掛點傾斜角θE是懸掛點切線與水平x向的夾角,二者成互余角,βhang-off=90°-θE。水平跨度x是TDP到懸掛點E之間的水平距離。

根據懶波型立管布局及形態特征,下懸鏈線段OA、拖曳段AB、舉升段BC、跨接段CD、懸掛段DE這5段可以看作為自由懸掛的懸鏈線。根據懸鏈線理論可以得到懶波型立管的形態理論公式[10,30]如下:

圖6 SLWR模型

(1)

(2)

(3)

(4)

(5)

S4+S5=Su

(6)

S2+S3=Sb

(7)

S=Sl

(8)

y1+y2-y3-y4+y5=V

(9)

(10)

(11)

(12)

(13)

(14)

(15)

(16)

(17)

式中:x1-5、y1-5、S1-5分別為下懸鏈線段OA、拖曳段AB、舉升段BC、跨接段CD、懸掛段DE的水平跨度、垂直高度和管道弧長;T0是管道O點處水平張力;wr為單位長度的管道在海水中的沉沒重力;we為單位長度的浮筒段在海水中的等效沉沒重力,wb=-we。

利用上述公式,即可求解懶波型立管形態。

3.2 大撓度梁理論模型

根據非線性大變形梁理論[31],將管道模擬為不可伸長的梁單元,不考慮扭轉,采用局部坐標系(θ,S)(傾斜角θ為梁與水平x向傾斜角,管道長度S為從O點開始的管道弧長),在法向和切向的微分方程[32]如下:

(18)

(19)

式中:T為軸向張力;EIj為彎曲剛度;wj為單位長度的管道在海水中的沉沒重力。

在全局坐標系(x,y)中SLWR水平位移xj、豎直位移yj,沿著立管分布的彎矩Mj和剪力Fj可通過如下公式得到:

(20)

(21)

(22)

(23)

4 結語

各國的科學研究人員和工程技術人員在SLWR設計、分析、安裝、運營等應用開發研究中取得了較為豐富的研究成果。隨著海洋油氣開發不斷深入和科技不斷發展,SLWR需要不斷更新理念和技術來滿足應用中的各種要求。在針對SLWR的動力學、安裝力學等研究方面,還需要大量的深入研究。

1) 完善觸地區管土作用模型研究。實際中海床土壤特性非常復雜,觸地點的疲勞效應、管土作用機理需要深入研究,需要采用更復雜、更符合海床土壤實際特性的管土作用模型。

2) 開展SLWR安裝過程中的動力學理論研究。SLWR的棄置回收和轉移安裝等安裝過程的各個不同階段的動態載荷復雜,分析難度大,應該考慮船體運動、錨鏈系統和波浪、海流等環境條件等因素對立管和纜繩受力的影響。

3) 開展SLWR模型試驗。應該通過準確的SLWR模型試驗,發現尺度效應的規律,不斷修正模型試驗結果,提高模型試驗精度,為今后SWLR應用提供參考。

[1] Zhang H,Li J.Effects of volumetric allocation on heave response of semi-submersible in deep sea[J].Science China Technology Sciences,2009,52(3):651-657.

[2] Borges F C L,Roitman N,Magluta C,et al.A concept to reduce vibrations in steel catenary risers by the use of viscoelastic materials[J].Ocean Engineering,2014(77):1-11.

[3] Bai Y,Bai Q.Subsea pipelines and risers[M].Amesterdam,Netherlands:Elsevier Science Ltd,2005:401-403.

[4] Thomas B,Benirschke A,Sarkar T.Parque das conchas (BC-10) steel lazy wave riser installation pre abandonment,recovery and transfer challenges[C]//Offshore Technology Conference.Houston,USA,2010:1-12.

[5] Phifer E H,Kopp F,Swanson R C,et al.Design and installation of Auger steel catenary risers[C]//Offshore Technology Conference.Houston,USA,1994:1-10.

[6] 陳景皓.深水鋼懸鏈線立管結構強度及可靠性研究[D].北京:中國石油大學(北京),2013:1-5.

[7] Stanton P N,Caldwell J,Yu C.Lazy wave SCR versus conventional and bottom wetted SCRs for tying back to a spar “I” tube[C]// Offshore Technology Conference.Houston,USA,2014:1-15.

[8] Newport A,Haheim S,Martineau E.Espirito Santo:Operational feedback on the use of steel rises on a turret moored FPSO[C]// Offshore Technology Conference.Houston,USA,2014:1-13.

[9] Lahey E,Burke R W,Balch E R,et al.Caesar-Tonga project steel lazy-wave riser design[C]// Offshore Technology Conference.Houston,USA,2013:1-16.

[10] Li S,Chau N.Dynamic response of deepwater lazy-wave catenary riser[C]// Deep Offshore Technology International.Amsterdam,Netherlands,2010:1-20.

[11] Quéau L M,Kimiaei M,Randolph M F.Analytical estimation of static stress range in oscillating steel catenary risers at touchdown areas and its application with dynamic amplification factors[J].Ocean Enginering,2014( 88):63-80.

[12] Quéau L M,Kimiaei M,Randolph M F.Lazy wave catenary risers:scaling factors and analytical approximation of the static stress range in the touchdown zone[C]//The 32nd International Conference on Ocean,Offshore and Arctic Engineering.Nantes,France,2013:1-10.

[13] da Silveira L M Y,de Arruda Martins C,Balena R,et al.Static analysis of lazy-wave steel risers[C]// The 18th International Congress of Mechanical Engineering.Ouro Preto,Brazil,2005:1-8.

[14] Campbell M.The complexities of fatigue analysis for deepwater risers[C]//The 2nd International Deepwater Pipeline Technology Conference.New Orleans,USA,1999:1-8.

[15] Yue B,Walters D,Yu W,et al.Lazy wave SCR on turret moored FPSO[C]// The 21st International Offshore and Polar Engineering Conference.Hawaii,USA,2011:1-8.

[16] Torres A L F L,Mourelle M M,Senra S F,et al.Influence of fatigue issues on the design of SCRs for deepwater offshore Brazil[C]// Offshore Technology Conference.Texas,USA,2008:1-9.

[17] Petruska D J,Zimmermann C A,Krafft K M,et al.Riser system selection and design for a deepwater FSO in the gulf of mexico[C]// Offshore Technology Conference.Houston,USA,2002:1-12.

[18] Balena R,Martins C A,Silveira L M Y,et al.Fatigue analysis of lazy-wave steel risers[C]// The 18th International Congress of Mechanical Engineering.Ouro Preto,Brazil,2005:1-8.

[19] Santillan S T,Virgin L N,Plaut R H.Static and dynamic behavior of highly deformed risers and piplines[J].Journal of Offshore Mechanics and Arctic Engineering,2010( 132):1-6.

[20] Santillan S T,Virgin L N.Numerical and experimental analysis of the static behavior of highly deformed risers[J].Ocean Engineering,2011( 38):1397-1402.

[21] Sun L,Bo Q.Global analysis of a flexible riser[J].Journal of Marine Science and Application,2011( 10):478-484.

[22] Sun H,Wang D.Sensitivity analysis of buoyancy modules parameters of lazy-wave flexible riser[C]// The 32nd International Conference on Ocean,Offshore and Arctic Engineering.Nantes,France,2013:1-7.

[23] 李清泉,楊和振.深海緩坡型臍帶纜干涉分析研究[J].振動與沖擊,2012,31(15):180-184.

[24] 宋磊建,付世曉,陳希恰,等.深海臍帶纜總體響應特性比較研究[J].振動與沖擊,2014,33(1):119-124.

[25] Ruan W,Bai Y,Cheng P.Static analysis of deepwater lazy-wave umbilical on elastic seabed[J].Ocean Engineering,2014( 91):73-83.

[26] Wang X,Zhou X,Ginnard A J,et al.Installation method evaluation of export SCRs[C]// Offshore Technology Conference.Houston,USA,2001:1-9.

[27] Hoffman J,Yun H,Modi A,et al.Parque das conchas (BC-10) pipeline,flowline and riser system design,installation and challenges[C]// Offshore Technology Conference.Houston,USA,2010:1-17.

[28] Senra S F,Mourelle M M,de Andrade E Q,et al.Challenges faced in the design of SLWR configuration for the pre-salt area[C]// The 30th International Conference on Ocean,Offshore and Arctic Engineering.Rotterdam,Netherlands,2011:1-9.

[29] Yue B,Campbell M,Walters D,et al.Improved SCR design for dynamic vessel application[C]// The 29th International Conference on Ocean,Offshore and Arctic Engineering.Shanghai,China,2010:1-10.

[30] 李艷,李欣,羅勇,等.深水鋼懸鏈線立管(SCR)的設計與研究進展[J].中國海洋平臺,2013,28(2):6-13.

[31] Chucheepsakul S,Monprapussorn T,Huang T.Large strain formulations of extensible flexible marine pipes transporting fluid[J].Journal of Fluids and Structues,2003,17(2):185-224.

[32] 王金龍,段夢蘭,田凱.海流作用下的深水懶波型立管形態研究[J].應用數學和力學,2014,35(9):959-968.

State-of-the-Art Advancement and Research Direction of Deepwater Steel Lazy-wave Riser

WANG Jinlong1,FAN Jiakun2
(1.PressurePipelineDivision,ChinaSpecialEquipmentInspectionandResearchInstitute,Beijing100017,China;2.ResearchandDevelopmentCenter,CNOOCGas&PowerGroup,Beijing100029,China)

As a kind of advanced SCR,Steel lazy-wave riser (SLWR) has gained more application in the development of deepwater oil and gas field.The technology overview of SLWR and its application status in offshore oil and gas development are introduced.The advancement of mechanic analysis and installation analysis is presented.The two main theories of SLWR analysis,including the catenary theory and large deflection beam theory is proposed.Finally,the outlook of SLWR in the analysis and the development direction in China is discussed,and it can provide certain guidance on the development of SLWR.

deep water;riser;steel lazy-wave riser;large deflection beam theory;developing trend

2016-12-19

國家重點基礎研究發展計劃(973 計劃)“深海工程結構的極端環境作用與壽命服役安全”(2011CB13702)

王金龍(1988-),男,工程師,博士,2015年畢業于復旦大學流體力學專業,現從事海底管道分析、檢測技術工作,E-mail:wangjinlong132@126.com。

1001-3482(2017)03-0001-07

TE95

A

10.3969/j.issn.1001-3482.2017.03.001

猜你喜歡
分析模型研究
一半模型
FMS與YBT相關性的實證研究
遼代千人邑研究述論
隱蔽失效適航要求符合性驗證分析
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
視錯覺在平面設計中的應用與研究
科技傳播(2019年22期)2020-01-14 03:06:54
EMA伺服控制系統研究
電力系統不平衡分析
電子制作(2018年18期)2018-11-14 01:48:24
電力系統及其自動化發展趨勢分析
主站蜘蛛池模板: 99视频有精品视频免费观看| 亚洲天堂视频在线观看| av天堂最新版在线| 亚洲AV无码乱码在线观看裸奔| 精品免费在线视频| 乱色熟女综合一区二区| 日本精品视频一区二区| 亚洲国产精品一区二区第一页免 | 亚洲日本韩在线观看| 精品亚洲麻豆1区2区3区 | 久久久久久久久亚洲精品| 国产高清无码第一十页在线观看| 本亚洲精品网站| 欧美午夜理伦三级在线观看| 亚洲天天更新| 精品撒尿视频一区二区三区| 色噜噜久久| 国产成人禁片在线观看| 亚洲av无码专区久久蜜芽| 精品欧美日韩国产日漫一区不卡| 国产精品永久久久久| 91福利国产成人精品导航| 国产精品国产三级国产专业不| 91美女视频在线观看| 国产农村妇女精品一二区| 国产美女91呻吟求| JIZZ亚洲国产| 成人字幕网视频在线观看| 亚洲精品欧美日韩在线| 国产日韩欧美精品区性色| 在线观看亚洲人成网站| 亚洲午夜国产精品无卡| 国产精品视频白浆免费视频| 香蕉久久国产超碰青草| 波多野结衣在线se| 亚洲综合第一页| 精品天海翼一区二区| 一本大道东京热无码av| 日本一区二区三区精品AⅤ| 亚洲一区二区精品无码久久久| 国产91在线|日本| 亚洲无码熟妇人妻AV在线| 色视频国产| 影音先锋丝袜制服| 国产国产人成免费视频77777| 自拍欧美亚洲| 超碰91免费人妻| 日韩一级二级三级| AV网站中文| 亚洲无码免费黄色网址| 免费激情网站| 国产精品2| 丁香婷婷激情网| 日韩精品无码一级毛片免费| 国产AV无码专区亚洲精品网站| 日韩精品无码一级毛片免费| 九色国产在线| 国产人成在线观看| 欧美一区福利| 污网站在线观看视频| 国产精品福利一区二区久久| 永久免费无码成人网站| 免费国产黄线在线观看| 国产精品.com| 亚洲无码视频喷水| 国产色图在线观看| 久久这里只精品国产99热8| 在线观看91香蕉国产免费| 亚洲精品第一页不卡| av无码一区二区三区在线| 精品福利国产| 国产精品无码AV中文| 成人av专区精品无码国产| 亚洲色无码专线精品观看| 性欧美在线| 永久成人无码激情视频免费| 久久国产精品麻豆系列| 国产一区二区色淫影院| 高潮毛片免费观看| 国产美女自慰在线观看| 国产成人AV男人的天堂| 中国毛片网|