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

    旋轉(zhuǎn)導(dǎo)向系統(tǒng)穩(wěn)定平臺自適應(yīng)動態(tài)面控制

    2024-01-01 00:00:00萬敏宋佳儒黃山山陳苗苗
    關(guān)鍵詞:RBF神經(jīng)網(wǎng)絡(luò)

    摘要:井下的多種干擾因素為旋轉(zhuǎn)導(dǎo)向系統(tǒng)穩(wěn)定平臺的控制器設(shè)計增加了復(fù)雜性。為了應(yīng)對未知摩擦力矩和建模誤差對穩(wěn)定平臺的不良影響,提出一種自適應(yīng)神經(jīng)網(wǎng)絡(luò)動態(tài)面控制方法,該方法使用RBF神經(jīng)網(wǎng)絡(luò)逼近摩擦及干擾力矩,設(shè)置狀態(tài)觀測器獲取由于相關(guān)參數(shù)不確定導(dǎo)致的建模誤差,并引入動態(tài)面方法避免傳統(tǒng)反步控制帶來的“微分爆炸”,最后使用李雅普諾夫法證明系統(tǒng)的穩(wěn)定性。結(jié)果表明,該控制方法在穩(wěn)定平臺模型存在摩擦力矩、未知干擾和建模誤差的情況下,仍能使工具面角準(zhǔn)確、快速地跟蹤輸入指令信號,具有較好的自適應(yīng)性和魯棒性。

    關(guān)鍵詞:旋轉(zhuǎn)導(dǎo)向系統(tǒng); 穩(wěn)定平臺; 動態(tài)面控制; RBF神經(jīng)網(wǎng)絡(luò); 狀態(tài)觀測器

    中圖分類號:TP 273"" 文獻(xiàn)標(biāo)志碼:A

    引用格式:萬敏,宋佳儒,黃山山,等.旋轉(zhuǎn)導(dǎo)向系統(tǒng)穩(wěn)定平臺自適應(yīng)動態(tài)面控制[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2024,48(5):190-199.

    WAN Min, SONG Jiaru, HUANG Shanshan, et al. Adaptive dynamic surface control of" stabilized platform in" rotary steerable drilling system [J]. Journal of China University of Petroleum (Edition of Natural Science),2024,48(5):190-199.

    Adaptive dynamic surface control of stabilized platform in

    rotary steerable drilling system

    WAN Min, SONG Jiaru, HUANG Shanshan, CHEN Miaomiao

    (School of Mechatronic" Engineering, Southwest Petroleum University, Chengdu 610500, China)

    Abstract:A variety of interference factors underground increase the complexity of the controller design of the stabilized platform of the rotary steerable drilling system. In order to deal with the adverse effects of unknown friction torque and modeling error on the stabilized platform, an adaptive neural network dynamic surface control method was proposed. The RBF neural network was used to approximate the friction and disturbance torque, and the state observer was set to obtain the modeling error caused by uncertainty of the relevant parameters. The dynamic surface method was introduced to avoid the “differential explosion” caused by the traditional backstepping control. Finally, the Lyapunov method was used to prove the stability of the system. The results show that the controller can make the toolface angle track the input command signal accurately and quickly under the condition of friction torque, unknown interference and modeling error in the stabilized platform model, and has good adaptability and robustness.

    Keywords: rotary steerable drilling system; stabilized platform; dynamic surface control; RBF neural network; state observer

    在導(dǎo)向鉆井中工具面角是井眼軌跡控制中的一個重要參數(shù)[1-4]。穩(wěn)定平臺作為井下鉆具組合中的控制部件,其作用是按要求調(diào)整工具面角使其快速平穩(wěn)地跟蹤給定的工具面角指令[5-7]。當(dāng)前對穩(wěn)定平臺的研究都集中于導(dǎo)向機(jī)構(gòu)及鉆頭所受的扭轉(zhuǎn)振動對其的影響[8-9],但實(shí)際上系統(tǒng)內(nèi)部結(jié)構(gòu)和所受到的干擾力矩對控制的影響也是需要關(guān)注的重點(diǎn)。湯楠等[10]建立了穩(wěn)定平臺的數(shù)學(xué)模型,指出平臺可以用成熟的控制理論來控制。湯楠等[11]建立了穩(wěn)定平臺的狀態(tài)空間模型,驗(yàn)證了系統(tǒng)能夠達(dá)到工具面角的控制要求。Huo等[12]對穩(wěn)定平臺中存在的軸承摩擦、脈動摩擦、黏滯摩擦等力矩進(jìn)行了理論推導(dǎo)及數(shù)值分析。薛啟龍等[5]考慮到渦輪電機(jī)的相互影響、鉆井液與部件的摩擦等井下復(fù)雜情況,提出一種采用偏差模糊控制和微調(diào)增量的控制算法。孫峰等[13]研究了捷聯(lián)式自動垂直鉆井系統(tǒng),對穩(wěn)定平臺模型進(jìn)行了分析建模,并通過捷聯(lián)式穩(wěn)定平臺地面聯(lián)調(diào)證明了工具的可行性。Wang等[8]開發(fā)了一種捷聯(lián)式穩(wěn)定平臺,其使用反饋控制來避免扭轉(zhuǎn)振動,讓墊塊推靠井壁產(chǎn)生側(cè)向力,從而達(dá)到預(yù)置角度?;魫矍宓龋?4]針對旋轉(zhuǎn)導(dǎo)向系統(tǒng)穩(wěn)定平臺提出了一種RBF神經(jīng)網(wǎng)絡(luò)滑模自適應(yīng)控制方法,并設(shè)計網(wǎng)絡(luò)節(jié)點(diǎn)的喚醒與激活闕值來減少網(wǎng)絡(luò)規(guī)模。Zhang等[15]建立了確定鉆頭移動所需的轉(zhuǎn)向力的力學(xué)和數(shù)學(xué)模型,提出一種為每個墊塊分配最優(yōu)力的控制算法。Wang等[16]提出一種基于模型的自抗擾控制算法的工具面控制方法,并基于系統(tǒng)模型設(shè)計了負(fù)載力矩估計器和外殼速度估計器來獲取外部干擾。Zhang等[17]針對存在系統(tǒng)不確定性和故障信號的情況下,采用自適應(yīng)容錯控制方法對旋轉(zhuǎn)導(dǎo)向系統(tǒng)的鉆井軌跡進(jìn)行跟蹤和校正。上述研究在穩(wěn)定平臺控制方面取得了豐富的研究成果,但是很少考慮實(shí)際摩擦和其他干擾因素影響,難以建立其精確模型。筆者提出一種由神經(jīng)網(wǎng)絡(luò)觀測器、自適應(yīng)神經(jīng)網(wǎng)絡(luò)控制和動態(tài)面控制相結(jié)合的旋轉(zhuǎn)導(dǎo)向系統(tǒng)穩(wěn)定平臺控制策略,使得穩(wěn)定平臺不僅在有干擾的情況下能保持穩(wěn)定,且參數(shù)攝動造成的誤差下也能有自適應(yīng)性。

    1 穩(wěn)定平臺模型

    1.1 動力學(xué)模型

    如圖1所示,穩(wěn)定平臺由上渦輪發(fā)電機(jī)、下渦輪力矩電機(jī)、電子控制單元、位置傳感器組成[18]。

    其中上渦輪電機(jī)負(fù)責(zé)供電[19],因其負(fù)載電阻不改變,其對穩(wěn)定平臺主軸產(chǎn)生的力矩可以視為恒定值。下渦輪力矩電機(jī)作為扭矩發(fā)生器,采用脈沖寬度調(diào)制技術(shù)(PWM),通過改變電流來調(diào)整輸出力矩 [20],從而改變工具面角度及速度。位置傳感器包含陀螺儀、線性加速度計、放大器等[16]。在此基礎(chǔ)上結(jié)合牛頓第二定律,則穩(wěn)定平臺動力學(xué)模型表示為

    Ldidt+iR+Ceω=u,

    Jdωdt=Cmi-Tf-e,

    didt=ω.(1)

    式中,u為電機(jī)電壓,V;i為電機(jī)電流,A;R為電阻,Ω;L為電感,H;Ce為反電動勢系數(shù);Cm為電機(jī)扭矩系數(shù);J為電機(jī)慣性力矩,kg·m2;θ為工具面角,rad;ω為角速度,rad/s; Tf為穩(wěn)定平臺中摩擦干擾力矩,N·M;e為建模誤差,包括了建模不確定性和參數(shù)偏差導(dǎo)致的誤差,為一個有界未知項(xiàng)。

    由此可以得到角度、角速度、加速度與力矩、電壓、摩擦及誤差之間的關(guān)系:

    dθdt=ω,

    LJCmd2ωdt2=u-

    RJCmdωdt-

    TfRCm-eRCm-

    LfCm-

    LCm-Ceω.

    (2)

    當(dāng)Tf和e為0時,對式(2)進(jìn)行拉普拉斯變換,可得到系統(tǒng)的輸入輸出關(guān)系為

    θ(s)=1sω(s),

    θ(s)U(s)=Cm

    (R+Ls)J2s+CeCms .(3)

    考慮加入陀螺儀轉(zhuǎn)換因子KW、脈寬調(diào)制的比例系數(shù)K

    PWM、電機(jī)電磁力矩與電流轉(zhuǎn)換因子KE和穩(wěn)定平臺內(nèi)部結(jié)構(gòu),很容易得到閉環(huán)控制系統(tǒng),如圖2所示。

    將狀態(tài)變量定義為x=[x1,x2]T=[θ,ω]T,且因電機(jī)中電感一般很小,將其忽略。

    得到穩(wěn)定平臺的狀態(tài)方程為

    1=x2,

    2=-(KPWMKWKECm+CeCm)RJx2-

    CmRJTf-CmRJe+

    KPWMKECmRJu,

    y=x1.(4)

    1.2 摩擦力矩

    圖3為穩(wěn)定平臺的結(jié)構(gòu)及摩擦力矩示意圖。

    將力矩電機(jī)產(chǎn)生的電磁力矩定義為Me,參考Huo等[12]對穩(wěn)定平臺內(nèi)部的摩擦力矩分析,穩(wěn)定平臺受到3個摩擦力矩如下。

    (1)鉆井液對穩(wěn)定平臺的黏滯摩擦力。穩(wěn)定平臺中的渦輪電機(jī)是通過鉆井液沖刷渦輪帶動發(fā)電機(jī)供電,鉆井液將充斥在鉆鋌內(nèi)部,當(dāng)鉆鋌旋轉(zhuǎn)時會帶動鉆鋌環(huán)空鉆井液旋轉(zhuǎn),從而產(chǎn)生一個對穩(wěn)定平臺的黏滯摩擦力矩M1。

    M1=4πμLr20r2r20-r2(ω0-ω). (5)

    式中,μ為鉆井液的黏度系數(shù);L為鉆井液長度,m;r0和r分別為鉆鋌的內(nèi)半徑和穩(wěn)定平臺的外半徑,m; w0 和w分別為鉆鋌和穩(wěn)定平臺的角速度,rad/s。

    (2)上下支撐軸承對穩(wěn)定平臺的摩擦力。穩(wěn)定平臺通過上下主支撐軸固定在鉆鋌中,在隨著鉆鋌同方向旋轉(zhuǎn)時,支撐軸與穩(wěn)定平臺之間保持的是一種相對旋轉(zhuǎn)的狀態(tài),而產(chǎn)生摩擦扭矩M2。

    M2=sign(n0-n)r1G(fu1sin α+fu2cos α).(6)

    式中,n0和 n分別為鉆鋌和穩(wěn)定平臺轉(zhuǎn)速,r/min;r1為軸承摩擦平均半徑,mm;G為穩(wěn)定平臺總重力,N; fu1和fu2分別為徑向和軸向軸承滾動摩擦系數(shù); α為傾角,rad。

    (3)上下渦輪電機(jī)軸承摩擦力。上下渦輪電機(jī)的核心軸都是電機(jī)的定子,鉆井液在驅(qū)動兩個渦輪電機(jī)時將會存在動摩擦力;由于上渦輪發(fā)電機(jī)與旋轉(zhuǎn)方向相同,力矩電機(jī)與旋轉(zhuǎn)方向相反,它們將提供相反的兩個摩擦力矩M3與M4。

    M3=r2G2(fu1sin α+fu2cos α),(7)

    M4=r3G3(fu1sin α+fu2cos α).(8)

    式中,r2 為上渦輪馬達(dá)軸承的平均摩擦半徑,mm; r3 為下渦輪馬達(dá)軸承的平均摩擦半徑, mm;G2為上渦輪電機(jī)的重力,N;G3為下渦輪電機(jī)的重力, N。

    平臺所受到的摩擦力為

    Tf=M1+M2+M3-M4.(9)

    假設(shè)狀態(tài)變量為x1=θ,x2=ω,且角速度與轉(zhuǎn)速的關(guān)系為2πn=ω,則最終的總摩擦力矩可表示為

    Tf=signn0-x22πr1G(fu1sin α+fu2cos α)+

    (r2G2-r3G3)(fu1sin α+fu2cos α)+

    4πμLr20r2r20-r2(2πn0-x2). (10)

    2 自適應(yīng)控制器設(shè)計

    針對穩(wěn)定平臺特性設(shè)計一個基于狀態(tài)觀測器的神經(jīng)網(wǎng)絡(luò)的自適應(yīng)控制律。利用RBF神經(jīng)網(wǎng)絡(luò)的逼近特性逼近無法測量的摩擦力矩,狀態(tài)觀測器用于控制因相關(guān)參數(shù)不確定而導(dǎo)致的建模誤差,然后利用加入一階濾波器的反演控制方法實(shí)現(xiàn)控制。

    將穩(wěn)定平臺方程改寫成狀態(tài)空間形式:

    =AX+B[u+f(X)+d],

    y=CTX.(11)

    其中

    X=[x1,x2]T,

    A=0100,

    B=0KPWMKECmRJ,

    C=10, f(X)=-KW+CeKPWMKEx2-TfKPWMKE, d=-eKPWMKE .

    且d≤dM,dM為一個未知正數(shù)。

    將基本狀態(tài)觀測器設(shè)計為

    ·=A+B

    [()+u]+K(y-CT),

    =CT. (12)

    其中

    =[1,2]T,

    K=k1k2,

    (),

    =y-CT.

    式中,為狀態(tài)觀測值; K為觀測器增益向量; ()為f(X)的狀態(tài)估計值;

    為y的估計誤差。

    設(shè)置觀測誤差項(xiàng)=X-,并定義δ=ε+d,則對其求導(dǎo)可得,觀測器誤差為

    ·=A+B[f(X)-()+d]-KCT(X-)=(A-KCT)+δ=

    Λ+δ.(13)

    通過選取合適的K可以保證Λ是一個嚴(yán)格的Hurwitz矩陣[21],其滿足

    ΛTP+PΛ=-2Q.(14)

    式中,P和Q為任意給定的正定矩陣。

    2.1 神經(jīng)網(wǎng)絡(luò)逼近

    徑向基函數(shù)(radial basis function)神經(jīng)網(wǎng)絡(luò)具有單隱層的三層前饋網(wǎng)絡(luò),針對非線性函數(shù)具有任意逼近特性[22-23]。采用RBF神經(jīng)網(wǎng)絡(luò)逼近f(x),網(wǎng)絡(luò)算法為

    σj(x)=exp-x-cj22bj j=1,…,m. (15)

    被逼近對象可表示為

    y=WTσj(x)+ε.(16)

    式中,x為網(wǎng)絡(luò)輸入;cj為隱含層的中心向量;bj為高斯函數(shù)寬度;W為神經(jīng)網(wǎng)絡(luò)權(quán)值;σj(x)為隱含層輸出;ε為網(wǎng)絡(luò)的逼近誤差,ε≤εN。

    f()神經(jīng)網(wǎng)絡(luò)估計值為

    ()=Tσj().(17)

    2.2 動態(tài)面控制器設(shè)計

    為了解決傳統(tǒng)backstepping控制算法所存在的“微分爆炸”問題,引入DSC技術(shù)[24]。

    取β1為α1的低通濾波器1τs+1的輸出,并滿足τ1+β1=α1,β1(0)=α1(0),τ為時間常數(shù),可得:

    1=α1-β1τ,所產(chǎn)生的濾波誤差為

    e1=β1-α1 .(18)

    定義第一個跟蹤誤差z1=y-yd和第二個速度誤差z2=2-β1,利用x2=2+2,得

    1=-d=x2-d=2+2-d=z2+α1+e1+2-d.(19)

    定義Lyapunov函數(shù)為

    V1=12XT

    P+12

    e21+12z21.(20)

    12XTP是為了評估觀測器的性質(zhì),且

    P=PTgt;0為一個正定矩陣。為保證

    1為正定,設(shè)置虛擬控制律α1=d-c1z1,其中c1為正常數(shù),將在之后進(jìn)行設(shè)計。從而可得其導(dǎo)數(shù)為

    1=12·TP+12TP·+e11+z11=

    12T[ΛTP+PΛ]+TPδ+e11+z11=

    -TQ+TPδ+e11+z1(z2+α1+e1+2-d)=

    -TQ+TPδ+e11+z1z2+z1e1+z12-c1z21.(21)

    對速度誤差z2=2-β求導(dǎo),得

    2=·2-1=u+()+K1-1=u+W*T

    σ()+ε′+K1-1=

    u+WTσ()+

    Tσ()+ε′+

    K1-1.(22)

    其中

    W*=W+,ε′=ε*-ε.

    定義Lyapunov函數(shù)為

    V2=V1+12z22+12γT. (23)

    式(23)中γ為待設(shè)計的正參數(shù),則對其求導(dǎo):

    2=1+z22+1γT·=1+z2(u+

    K1+WTσ()+Tσ()+ε′-1)+1γT·=1+z2(u+K1+WTσ()+ε′-1)+1γT[γz2σ()-].(24)

    得到控制輸入u與自適應(yīng)律W·分別為

    u=-c2z2-z1-K1-WTσ(2)+α1-β1τ ,

    =γz2σ(2)-φW.(25)

    此時:

    2=-TQ+TPδ+e11+z1e1+z12+z2ε′+φγTW-c1z21-c2z2.(26)

    為了顯示整體的控制器結(jié)構(gòu),上述穩(wěn)定平臺控制方案配置流程如圖4所示。

    2.3 穩(wěn)定性

    為了保證整個系統(tǒng)的穩(wěn)定性,易知選取的正定Lyapunov函數(shù)為

    V=V2=12XTP+12e21+12z21+12z22+12γT.(27)

    根據(jù)不等式2ab≤a2+b2可得

    TPδ≤12(2+Pδ2),(28)

    z12≤12(z21+22),(29)

    z2ε′≤12(z22+ε′2).(30)

    由Youngs不等式得

    z2e1≤z22+14e21.(31)

    對正定矩陣Q恒成立的不等式[25]為

    λmin(Q)2≤TQ≤λmax(Q)2.(32)

    由式(28) ~ (32)可得

    2≤-[λmin(Q)-1]

    2+122

    Pδ2+e11+14e21+12ε′2+φγTW-c1-32z21-c2-12z22.(33)

    根據(jù)Lyapunov穩(wěn)定定理,式(33)中如果滿足-λmin(Q)-1gt;0,可知狀態(tài)觀測器(12)在李雅普諾夫定義下是漸近穩(wěn)定的。

    通過式(18),對其求導(dǎo)可得

    1=1-1=-e1τ+B1(·).(34)

    根據(jù)文獻(xiàn)[26],B1(·)為一個連續(xù)的有界函數(shù),滿足B1(·)=B1(z1,e1,yd,d,d)=-1=c11-d,且有一個最大的絕對值M。由youngs不等式得

    e11≤-e21τ+e1B1(·)≤-e21τ+12e21+12M2.(35)

    φγTW≤-φ2γT+φ2γ*TW*≤-φ2γ

    2+φ2γ*T

    *.(36)

    將式(35)、(36)帶入式(33),可以進(jìn)一步推導(dǎo)得

    ≤-[λmin(Q)-1]

    2+

    12

    Pδ2

    -c1-32z21-

    c2-12z22-

    1τ-14e21+12M2+12ε′2+φ2γ*T*-

    φ2γ

    W*2≤-[λmin(Q)-1]

    2

    -φ2γ

    W*2

    -c1-32z21-c2-12z22-1τ-14e21+Π.(37)

    其中

    Π=12Pδ2+12ε′2+φ2γW*TW*+12M2.(38)

    定理1[27-28]:考慮非線性系統(tǒng)公式(4),狀態(tài)觀測器式(12)和實(shí)際控制器式(25)組成的閉環(huán)系統(tǒng),對任意給定的正常數(shù)p,如果V(0)≤p,則存在正設(shè)計參數(shù)φi,ci,τ,使得閉環(huán)系統(tǒng)的所有信號半全局一致最終有界,系統(tǒng)誤差可以收斂到任意小。

    根據(jù)定理1,可知要使系統(tǒng)收斂,則參數(shù)需滿足:

    [λmin(Q)-1]>0,c1-32>0,c2-12>0,14-1τ>0. (39)

    并定義正參數(shù):

    c=min2[λmin(Q)-1]λmax(P),2c1-32,2c2-12,

    21τ-14,φ.(40)

    則式(37)可以表示為

    ≤-cV+Π.(41)

    對式(41)進(jìn)行微分求解可得

    0≤V≤Πc+V(0)-Πce-ct.(42)

    由定理1可得,閉環(huán)系統(tǒng)中所有狀態(tài)都是半全局最終一致有界的,并需合理設(shè)計參數(shù),以保證跟蹤誤差收斂到0的很小的鄰域內(nèi)。不等式(42)表明,該穩(wěn)定系統(tǒng)中的所有信號都是半全局一致最終有界的(SGUUB),并且跟蹤誤差滿足不等式

    y-yd≤2V(0)e-ct+Πc .(43)

    不等式(43)意味著跟蹤誤差是有界的,可以選擇合適的參數(shù)設(shè)計相應(yīng)的控制律,使跟蹤誤差收斂到一個非常小的零鄰域。此外Πc的值與參數(shù)的選擇有關(guān),如果參數(shù)選擇合理,誤差可能足夠小,但會導(dǎo)致較大的輸入。因此在實(shí)際的井下設(shè)計中,需要權(quán)衡誤差與輸入信號之間的關(guān)系,選擇合適的控制參數(shù)。

    3 試驗(yàn)仿真驗(yàn)證

    3.1 數(shù)值仿真

    在進(jìn)行仿真之前,需要先對穩(wěn)定平臺各模型參數(shù)進(jìn)行設(shè)定,其中典型參數(shù)來自文獻(xiàn)[12]及文獻(xiàn)[14]。穩(wěn)定平臺模型參數(shù):

    J=0.025 kg·m2,KPWM=3.44,KW=5.74,KE=0.22,Cm=3.82, Ce=0.4,R=12.5 Ω。摩擦力矩特性參數(shù):α=0.25 rad , G=350 N,G2=20 N,G3=350 N,L=1 m, n0=20 rad/s,r=0.055 m, r0=0.0675 m,r1=0.0255 m,r2=0.015 m,r3=0.015 m,fu1=0.002,fu2=0.005,μ=0.5。

    對于RBF神經(jīng)網(wǎng)絡(luò)參數(shù),RBF神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)選取2-5-1,根據(jù)輸入信號的范圍,選取高斯函數(shù)中心均勻分布在[-1,1]內(nèi),即cj和bj分別?。?1,-0.5,0,0.5,1]和0.5,令權(quán)值的初始值W(0)=0,將觀測器增益設(shè)計為K=k1k2=51000來克服模型誤差。

    設(shè)計的控制律參數(shù)為γ=0.1,τ=0.05,φ=0.5,c1=40,c2=10。

    將穩(wěn)定平臺的建模誤差設(shè)置為e(t)=cos t,設(shè)置理想工具面角輸入信號yd=sin t,系統(tǒng)初始狀態(tài)[θ0,ω0]=[-0.1,0]?;谶@些給定參數(shù)進(jìn)行仿真,得到圖5~10的結(jié)果。由圖5可知,系統(tǒng)在0.2 s內(nèi)即可穩(wěn)定跟蹤輸入的理想工具面角信號;由圖6可知,系統(tǒng)穩(wěn)態(tài)跟蹤誤差在0.004 rad之內(nèi),可知在同時存在多個摩擦及建模誤差的情況下,所提出的控制律可以達(dá)到快速、準(zhǔn)確的跟蹤效果;圖7為控制輸入信號情況,在此輸入的控制下,系統(tǒng)可以快速穩(wěn)定地跟蹤所設(shè)定的工具面角;圖8為摩擦力矩的逼近情況,可知通過使用RBF神經(jīng)網(wǎng)絡(luò),所設(shè)計的控制律將根據(jù)輸出誤差自適應(yīng)調(diào)整網(wǎng)絡(luò)權(quán)值,使得摩擦力矩以較小的誤差被逼近補(bǔ)償;圖9為工具面角觀測值曲線,狀態(tài)觀測器通過狀態(tài)重構(gòu)可以實(shí)現(xiàn)對建模誤差的補(bǔ)償。

    在本次仿真試驗(yàn)中,由于設(shè)置的輸入工具面角是一個時變量,而不像已有文獻(xiàn)中使用固定的常量工具面角度,結(jié)果驗(yàn)證了所提出控制器的優(yōu)越性,達(dá)到了使輸出能夠良好跟蹤期望指令,逼近未知摩擦力矩并克服誤差擾動的目的,證明穩(wěn)定平臺在本控制律作用下能達(dá)到跟蹤時變的連續(xù)輸入工具面角的效果。

    3.2 自適應(yīng)性仿真驗(yàn)證

    考慮到實(shí)際井下作業(yè)時,由于溫度、鉆井液會造成參數(shù)如轉(zhuǎn)動慣量J、負(fù)載電阻R的變化,實(shí)際上本文中所考慮的摩擦項(xiàng)還不是全部的摩擦力矩,還存在其他的摩擦干擾未被考慮。將摩擦力矩擴(kuò)大為兩倍,電阻及轉(zhuǎn)動慣量減少20%,改變建模誤差項(xiàng)的頻率使e(t)=cos (5t+0.5),得到圖10、 11。

    可以看出,工具面角的跟蹤效果依然良好,仍在0.2 s內(nèi)便能跟蹤給定輸入信號,且穩(wěn)態(tài)誤差也保持在0.007 rad之內(nèi)。證明所提出的控制律在面對實(shí)際井下未知摩擦及干擾時能很好地克服,在參數(shù)攝動之下也能維持良好的控制效果,有較強(qiáng)的自適應(yīng)性和魯棒性,為實(shí)際井下工程提供新的參考依據(jù)。

    3.3 對比仿真驗(yàn)證

    為了更好地說明提出的控制律的有效性以及為實(shí)際穩(wěn)定平臺設(shè)計控制器的原因,將把所提出的自適應(yīng)神經(jīng)網(wǎng)絡(luò)動態(tài)表面控制(ANNDSC)算法與現(xiàn)階段廣泛應(yīng)用于井下控制的傳統(tǒng)PID控制器及無神經(jīng)網(wǎng)絡(luò)狀態(tài)觀測器的DSC控制器進(jìn)行對比試驗(yàn)。

    傳統(tǒng)的PID控制器由于其簡單可靠的特點(diǎn),在實(shí)際工程中得到了廣泛的應(yīng)用,在試驗(yàn)中使用試錯法選擇了3個增益:

    kp=800、ki=0.5、kd=500,達(dá)到較好的控制效果。對于DSC控制器,通過推導(dǎo)可以分別得到DSC控制器的虛擬控制輸入和實(shí)際控制輸入如下:

    α1=-c1z1+yd.(44)

    u=RJCm(KPWMKWKECm+CeCm)RJx2+

    CmRJTf-c2z2-z1+α1-β1τ. (45)

    為了公平性,對于仿真參數(shù),為DSC控制器選擇了與ANNDSC相同的參數(shù)進(jìn)行對比仿真試驗(yàn),并得到圖12。由圖12可以明顯看出,在相同的參數(shù)下DSC控制器沒有達(dá)到穩(wěn)定跟蹤效果,PID控制器的跟蹤誤差也在[0,0.02]的范圍波動,只有提出的ANNDSC控制器最終達(dá)到了穩(wěn)態(tài),誤差極小,控制性能明顯優(yōu)于其他兩個控制器。

    若考慮參數(shù)的改變及摩擦力矩增加,將摩擦力矩擴(kuò)大為兩倍,電阻及轉(zhuǎn)動慣量減少20%,并改變建模誤差項(xiàng)的頻率,使e(t)=cos(5t+0.5),得到圖13。由圖13可看出,PID控制器和普通DSC控制器都不同程度地受到影響,這樣的控制器在井下會極快地失效,造成不可逆損失。而新提出的ANNDSC控制器未受到參數(shù)及摩擦力矩改變的影響,跟蹤誤差仍穩(wěn)定在0.007 rad之內(nèi),并沒有因?yàn)楦蓴_增加而增加跟蹤誤差,有極強(qiáng)的自適應(yīng)能力。

    4 結(jié) 論

    (1)基于狀態(tài)觀測器的神經(jīng)網(wǎng)絡(luò)自適應(yīng)動態(tài)面控制律可以有效地克服未知的摩擦力矩和建模誤差,確保穩(wěn)定平臺能夠快速且穩(wěn)定地跟蹤工具面角。

    (2)神經(jīng)網(wǎng)絡(luò)自適應(yīng)動態(tài)面控制律能夠適應(yīng)實(shí)際井下的未知摩擦和干擾,驗(yàn)證所提出的控制器具有較好的自適應(yīng)性和魯棒性。

    (3)與PID控制和DSC控制相比,所提出的控制律展現(xiàn)出更好的抗干擾性和跟蹤精度,相對于現(xiàn)存的經(jīng)典控制策略具有顯著優(yōu)勢。

    參考文獻(xiàn):

    [1] LIUX S. Directional deflection equations for steerable drilling tools and the control mechanism of wellbore trajectory[J]. Petroleum Exploration and Development, 2017,44(5):834-839.

    [2] LIU Y, QIN X B, JIA J B, et al. Factors to influence the trajectory control ability of a reverse push-the-bit rotary steerable system[J]. Processes, 2022,10(8):1621.

    [3] SHENG L, NIU Y C, WANG W L, et al. Estimation of toolface for dynamic point-the-bit rotary steerable systems via nonlinear polynomial filtering [J]. IEEE Transactions on Industrial Electronics, 2022,69(7):7192-7201.

    [4] 韓來聚,王瑞和,劉新華,等.調(diào)制式旋轉(zhuǎn)導(dǎo)向鉆井系統(tǒng)穩(wěn)定平臺控制原理及性能分析[J].石油大學(xué)學(xué)報(自然科學(xué)版),2004,28(5):49-51,60.

    HAN Laiju, WANG Ruihe, LIU Xinhua, et al. Control principle and servo platform performance analysis of modulated rotary steerable drilling system [J]. Journal of the University of Petroleum, China (Edition of Natural Science), 2004,28(5):49-51,60.

    [5] 薛啟龍,韓來聚,楊錦舟,等.旋轉(zhuǎn)導(dǎo)向鉆井穩(wěn)定平臺控制系統(tǒng)仿真研究[J].石油鉆探技術(shù),2010,38(4):10-14.

    XUE Qilong, HAN Laiju, YANG Jinzhou, et al. Study on controlling simulation system for stabilizing platform in rotary steering drilling system[J]. Petroleum Drilling Techniques, 2010,38(4):10-14.

    [6] SPENCER R, HANSON J, HOFFMANN O, et al. New directional drilling simulation tool reveals link between dynamic stability and tool face control[R]. SPE-163467-MS, 2013.

    [7] 劉永旺,郭福貴,管志川,等.協(xié)同鉆頭推靠與指向作用控制井眼軌跡的旋轉(zhuǎn)導(dǎo)向系統(tǒng)[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2022,46(5):73-80.

    LIU Yongwang, GUO Fugui, GUAN Zhichuan, et al. A new rotary steering method for controlling wellbore trajectory by combining pushing and pointing forces on drill bit [J]. Journal of China University of Petroleum (Edition of Natural Science), 2022,46(5):73-80.

    [8] WANG R H, XUE Q L, HAN L J, et al. Torsional vibration analysis of push-the-bit rotary steerable drilling system[J]. Meccanica, 2014,49(7):1601-1615.

    [9] SONGX Y, VADALI M, XUE Y Z, et al. Tracking control of rotary steerable tool face in directional drilling//2016 IEEE International Conference on" Advanced Intelligent Mechatronics (AIM)[C]. Banff, AB, Canada: IEEE, 2016:1210-1215.

    [10] 湯楠,穆向陽.調(diào)制式旋轉(zhuǎn)導(dǎo)向鉆井工具穩(wěn)定平臺控制機(jī)構(gòu)研究[J].石油鉆采工藝,2003(3):9-12,81.

    TANG Nan, MU Xiangyang. Research on control mechanism of stabilized platform of modulated rotary steering drilling tool[J]. Oil Drilling amp; Production Technology, 2003(3):9-12,81.

    [11] 湯楠,霍愛清,崔琪琳.基于狀態(tài)空間法的旋轉(zhuǎn)導(dǎo)向鉆井工具控制系統(tǒng)研究[J].石油學(xué)報,2004,25(2):89-92.

    TANG Nan, HUO Aiqing, CUI Qilin. Control system in rotary steering drilling tool based on state space method [J]. Acta Petrolei Sinica, 2004,25(2):89-92.

    [12] HUO A Q, HE Y Y, WANG Y L, et al. Nonlinear analysis of friction torque of steady platform control system for rotary steering drilling//2009 IEEE International Conference on Intelligent Computing and Intelligent Systems[C]. Shanghai: IEEE, 2009:577-580.

    [13] 孫峰,呂官云,馬清明.捷聯(lián)式自動垂直鉆井系統(tǒng)[J].石油學(xué)報,2011,32(2):360-363.

    SUN Feng, L Guanyun, MA Qingming. A strap-down automatic vertical drilling system[J]. Acta Petrolei Sinica, 2011,32(2):360-363.

    [14] 霍愛清,邱龍,汪躍龍.旋轉(zhuǎn)導(dǎo)向鉆井穩(wěn)定平臺的RBF網(wǎng)絡(luò)滑模變結(jié)構(gòu)控制[J].西安石油大學(xué)學(xué)報(自然科學(xué)版),2016,31(4):103-108.

    HUO Aiqing, QIU Long, WANG Yuelong. Sliding mode variable structure control of stabilized platform in rotary steerable drilling system based on RBF neural network[J]. Journal of Xian Shiyou University (Natural Science Edition), 2016,31(4):103-108.

    [15] ZHANG L G, LIU G R, LI W, et al. Analysis and optimization of control algorithms for RSSTSP for horizontal well drilling[J]. Journal of Petroleum Exploration and Production Technology, 2018,8(4):1069-1078.

    [16] WANG W L, GENG Y F, WANG N, et al. Toolface control method for a dynamic point-the-bit rotary steerable drilling system[J]. Energies, 2019,12(10):1-20.

    [17] ZHANG C, ZOU W, CHENG N B, et al. Adaptive fault-tolerant control for trajectory tracking and rectification of directional drilling[J]. International Journal of Control, Automation and Systems, 2022,20(1):334-348.

    [18] ALTURBEHH, WHIDBRON J, LUK P, et al. Modelling and control of the roll-stabilized control unit of a rotary steerable system directional drilling tool [J]. The Journal of Engineering,2019(17):4555-4559.

    [19] 黃曉凡,林恩懷.石油鉆井用泥漿發(fā)電機(jī)綜述[J].石油儀器,2013,27(1):31-33,36,8.

    HUANG Xiaofan, LIN Enhuai. An overview of mud alternator used in oil drilling[J]. Petroleum Instruments, 2013,27(1):31-33,36,8.

    [20] BUJA G S, KAZMIERKOWSKI M P. Direct torque control of PWM inverter-fed AC motors: a survey[J]. IEEE Transactions on Industrial Electronics, 2004,51(4):744-757.

    [21] GHOLAMI A, MARKAZI A H D. A new adaptive fuzzy sliding mode observer for a class of MIMO nonlinear systems[J]. Nonlinear Dynamics, 2012,70(3):2095-2105.

    [22] HARTMAN E J, KEELER J D, KOWALSKI J M. Layeredneural networks with Gaussian hidden units as universal approximations [J]. Neural Computation, 1990,2:210-215.

    [23] SANNER R M, SLOTINE J E. Gaussian networks for direct adaptive control[J]. IEEE Trans Neural Netw, 1992,3(6):837-863.

    [24] CUI Y, ZHANG H G, WANG Y C, et al. Adaptive neural dynamic surface control for a class of uncertain nonlinear systems with disturbances[J]. Neurocomputing, 2015,165:152-158.

    [25] SHOJAEI F, AREFI M M, KHAYATIAN A, et al. Observer-based fuzzy adaptive dynamic surface control of uncertain nonstrict feedback systems with unknown control direction and unknown dead-zone[J]. IEEE Transactions on Systems, Man amp; Cybernetics, 2019,49(11):2340-2351.

    [26] HOU N, DONG H L, WANG Z D, et al. Non-fragile state estimation for discrete Markovian jumping neural networks[J]. Neurocomputing, 2016,179:238-245.

    [27] CHEN B, LIU X P, LIU K F, et al. Novel adaptive neural control design for nonlinear MIMO time-delay systems[J]. Automatica, 2009,45(6):1554-1560.

    [28] HASSAN K K.非線性系統(tǒng)[M].北京:電子工業(yè)出版社,2017:78-81.

    (編輯 沈玉英)

    猜你喜歡
    RBF神經(jīng)網(wǎng)絡(luò)
    基于BP和RBF的手機(jī)網(wǎng)民規(guī)模預(yù)測方法研究
    基于RBF神經(jīng)網(wǎng)絡(luò)的交通事故預(yù)測
    商情(2016年46期)2017-03-06 04:08:38
    無線Mesh網(wǎng)絡(luò)發(fā)展
    無線Mesh網(wǎng)絡(luò)安全性研究
    無線Mesh網(wǎng)絡(luò)基礎(chǔ)知識
    基于RBF神經(jīng)網(wǎng)絡(luò)的PID自校正控制研究
    基于RBF神經(jīng)網(wǎng)絡(luò)的一回路核動力裝置典型故障診斷
    科技視界(2016年18期)2016-11-03 23:59:24
    基于RBF神經(jīng)網(wǎng)絡(luò)的數(shù)據(jù)挖掘方法探究
    RBF神經(jīng)網(wǎng)絡(luò)在創(chuàng)業(yè)板新股定價問題上的研究
    基于遺傳優(yōu)化RBF神經(jīng)網(wǎng)絡(luò)的電動負(fù)載模擬器控制
    亚洲一级一片aⅴ在线观看| 国产高清有码在线观看视频| av女优亚洲男人天堂| 亚洲av综合色区一区| 国产日韩欧美亚洲二区| 九草在线视频观看| 在线观看免费日韩欧美大片 | 美女中出高潮动态图| 国产高清不卡午夜福利| 极品人妻少妇av视频| 国产高清国产精品国产三级| 亚洲精品成人av观看孕妇| 十八禁高潮呻吟视频| 岛国毛片在线播放| 一区在线观看完整版| 成人漫画全彩无遮挡| 国产 精品1| 国产 一区精品| 国产伦精品一区二区三区视频9| 桃花免费在线播放| 免费观看的影片在线观看| 熟妇人妻不卡中文字幕| 久久女婷五月综合色啪小说| 国产精品秋霞免费鲁丝片| 国产高清有码在线观看视频| 九九爱精品视频在线观看| 不卡视频在线观看欧美| 美女主播在线视频| 91精品国产国语对白视频| 91精品国产九色| 欧美激情 高清一区二区三区| 久久热精品热| 激情五月婷婷亚洲| 国产精品一二三区在线看| 欧美变态另类bdsm刘玥| 亚洲欧美一区二区三区黑人 | 如日韩欧美国产精品一区二区三区 | 欧美性感艳星| 久久99热6这里只有精品| 久久毛片免费看一区二区三区| 三级国产精品欧美在线观看| 97在线视频观看| 人妻少妇偷人精品九色| 亚洲丝袜综合中文字幕| 亚洲精品自拍成人| 久久久久人妻精品一区果冻| 天堂俺去俺来也www色官网| 18+在线观看网站| 日韩av不卡免费在线播放| 国国产精品蜜臀av免费| 人妻少妇偷人精品九色| 日韩精品免费视频一区二区三区 | 国产午夜精品久久久久久一区二区三区| 伦理电影免费视频| 如何舔出高潮| 欧美+日韩+精品| 欧美亚洲 丝袜 人妻 在线| 国产片内射在线| 97精品久久久久久久久久精品| 一级毛片aaaaaa免费看小| 人妻夜夜爽99麻豆av| 精品久久蜜臀av无| 免费日韩欧美在线观看| 久热这里只有精品99| 天天操日日干夜夜撸| 自拍欧美九色日韩亚洲蝌蚪91| 2021少妇久久久久久久久久久| 老女人水多毛片| 欧美日韩一区二区视频在线观看视频在线| 日本vs欧美在线观看视频| 91精品伊人久久大香线蕉| 欧美成人精品欧美一级黄| 精品少妇黑人巨大在线播放| 日韩伦理黄色片| 91在线精品国自产拍蜜月| 欧美成人精品欧美一级黄| 黄色一级大片看看| 亚洲国产av影院在线观看| av卡一久久| 三上悠亚av全集在线观看| www.av在线官网国产| 亚洲天堂av无毛| 色视频在线一区二区三区| 久久99热6这里只有精品| 亚洲国产精品999| freevideosex欧美| 免费av中文字幕在线| 欧美xxxx性猛交bbbb| 少妇被粗大的猛进出69影院 | 国产极品粉嫩免费观看在线 | 男女无遮挡免费网站观看| 国产免费福利视频在线观看| 久久久精品94久久精品| 欧美亚洲日本最大视频资源| 免费高清在线观看日韩| 国产成人精品婷婷| 国产日韩欧美视频二区| 日韩在线高清观看一区二区三区| 欧美日韩成人在线一区二区| 国产成人免费无遮挡视频| 黑人巨大精品欧美一区二区蜜桃 | 国产欧美日韩综合在线一区二区| 午夜91福利影院| 精品亚洲成国产av| 亚洲精品亚洲一区二区| av国产精品久久久久影院| 中文字幕亚洲精品专区| 亚洲第一区二区三区不卡| 国产一区亚洲一区在线观看| 人妻少妇偷人精品九色| 欧美日本中文国产一区发布| 七月丁香在线播放| 免费黄色在线免费观看| 欧美+日韩+精品| 久久久久精品久久久久真实原创| 欧美 亚洲 国产 日韩一| 国产色婷婷99| 国产精品成人在线| 午夜福利影视在线免费观看| 欧美老熟妇乱子伦牲交| 18在线观看网站| 26uuu在线亚洲综合色| 最黄视频免费看| 欧美精品国产亚洲| 天天躁夜夜躁狠狠久久av| 国产一区二区在线观看日韩| 大码成人一级视频| 黄色怎么调成土黄色| 99久久精品国产国产毛片| 欧美3d第一页| 午夜影院在线不卡| 人妻制服诱惑在线中文字幕| 久久久久久久久久久久大奶| av免费在线看不卡| 久久精品久久精品一区二区三区| 国产成人aa在线观看| 中文字幕人妻熟人妻熟丝袜美| 人妻一区二区av| 91久久精品电影网| 女人精品久久久久毛片| xxxhd国产人妻xxx| 精品亚洲成a人片在线观看| 亚洲第一区二区三区不卡| 国产精品国产三级国产av玫瑰| av电影中文网址| 精品视频人人做人人爽| 成人国产麻豆网| 热99国产精品久久久久久7| 久久这里有精品视频免费| 久久热精品热| 天美传媒精品一区二区| av黄色大香蕉| 精品人妻偷拍中文字幕| 日本爱情动作片www.在线观看| 亚洲av在线观看美女高潮| 天堂中文最新版在线下载| 黑人猛操日本美女一级片| 欧美亚洲日本最大视频资源| 三级国产精品欧美在线观看| 亚洲av电影在线观看一区二区三区| 特大巨黑吊av在线直播| 22中文网久久字幕| 午夜福利视频精品| 久久久国产欧美日韩av| 亚洲精品日韩av片在线观看| 99久久精品国产国产毛片| 国产av一区二区精品久久| 国产极品天堂在线| 国语对白做爰xxxⅹ性视频网站| 久久狼人影院| 亚洲美女黄色视频免费看| 少妇猛男粗大的猛烈进出视频| 久久久a久久爽久久v久久| 国产精品国产三级专区第一集| 日日爽夜夜爽网站| 一个人免费看片子| 国产精品无大码| 国产免费视频播放在线视频| 少妇被粗大的猛进出69影院 | 亚洲一级一片aⅴ在线观看| 少妇 在线观看| 欧美日本中文国产一区发布| 大香蕉久久成人网| 一级毛片 在线播放| av又黄又爽大尺度在线免费看| 久久青草综合色| 亚洲人成网站在线播| 亚洲中文av在线| 日韩一本色道免费dvd| 亚洲伊人久久精品综合| 国产综合精华液| 九九爱精品视频在线观看| 一级毛片黄色毛片免费观看视频| 国产成人免费观看mmmm| 亚洲色图综合在线观看| 亚洲精品中文字幕在线视频| 日韩视频在线欧美| 国产av精品麻豆| 女人精品久久久久毛片| 国产亚洲最大av| 久久99热这里只频精品6学生| 91aial.com中文字幕在线观看| 国产 一区精品| 视频中文字幕在线观看| 韩国av在线不卡| kizo精华| 一级,二级,三级黄色视频| 十分钟在线观看高清视频www| 国产免费视频播放在线视频| 欧美日韩亚洲高清精品| 久久毛片免费看一区二区三区| 成人亚洲精品一区在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲国产色片| 亚洲在久久综合| 又大又黄又爽视频免费| 精品少妇内射三级| 欧美日韩av久久| 久久婷婷青草| 日韩av在线免费看完整版不卡| 日本猛色少妇xxxxx猛交久久| 国精品久久久久久国模美| 简卡轻食公司| 国产欧美日韩一区二区三区在线 | 九九久久精品国产亚洲av麻豆| 美女内射精品一级片tv| 久久久久久久大尺度免费视频| 精品视频人人做人人爽| 夫妻性生交免费视频一级片| 美女大奶头黄色视频| 日韩av免费高清视频| 亚洲精品日本国产第一区| 久久狼人影院| 伊人亚洲综合成人网| av一本久久久久| 久久人人爽人人片av| 午夜免费男女啪啪视频观看| 大香蕉久久成人网| 校园人妻丝袜中文字幕| 晚上一个人看的免费电影| 国产成人精品福利久久| 免费播放大片免费观看视频在线观看| 亚洲成人手机| 亚洲少妇的诱惑av| 你懂的网址亚洲精品在线观看| 一本色道久久久久久精品综合| 午夜福利视频精品| 国产日韩欧美亚洲二区| 亚洲第一区二区三区不卡| 丰满乱子伦码专区| 欧美精品亚洲一区二区| 老熟女久久久| 国国产精品蜜臀av免费| 熟女人妻精品中文字幕| 少妇的逼水好多| 高清午夜精品一区二区三区| 国产极品粉嫩免费观看在线 | 少妇熟女欧美另类| 人妻 亚洲 视频| 制服丝袜香蕉在线| √禁漫天堂资源中文www| 免费日韩欧美在线观看| 夜夜骑夜夜射夜夜干| 乱码一卡2卡4卡精品| 嘟嘟电影网在线观看| 夜夜骑夜夜射夜夜干| 国产精品免费大片| 国产成人av激情在线播放 | 三上悠亚av全集在线观看| 欧美激情国产日韩精品一区| 亚洲美女黄色视频免费看| 日本爱情动作片www.在线观看| 日日啪夜夜爽| 永久网站在线| 91久久精品国产一区二区三区| 天堂俺去俺来也www色官网| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 大片免费播放器 马上看| 在线观看人妻少妇| 国产免费一级a男人的天堂| 欧美另类一区| 久久久精品区二区三区| 少妇高潮的动态图| 久久精品久久久久久噜噜老黄| 欧美精品亚洲一区二区| 在线看a的网站| 能在线免费看毛片的网站| 久久久久国产网址| 日韩成人av中文字幕在线观看| 男女无遮挡免费网站观看| 久久久久久久久大av| 欧美日韩亚洲高清精品| 亚洲天堂av无毛| 一级黄片播放器| 国产免费一级a男人的天堂| 亚洲婷婷狠狠爱综合网| 国产免费现黄频在线看| 久久韩国三级中文字幕| 久久鲁丝午夜福利片| 丰满乱子伦码专区| 赤兔流量卡办理| 亚洲精品色激情综合| 久久精品国产亚洲av天美| 亚洲av国产av综合av卡| 狂野欧美激情性bbbbbb| 天天操日日干夜夜撸| 日韩免费高清中文字幕av| 免费大片黄手机在线观看| 亚洲av电影在线观看一区二区三区| 精品视频人人做人人爽| 少妇的逼好多水| 免费黄频网站在线观看国产| 我的女老师完整版在线观看| 黑人巨大精品欧美一区二区蜜桃 | 啦啦啦中文免费视频观看日本| 免费看不卡的av| 精品久久国产蜜桃| 人妻 亚洲 视频| 国产成人精品婷婷| 精品一区二区三卡| 欧美精品国产亚洲| 日韩精品免费视频一区二区三区 | 欧美性感艳星| 97超碰精品成人国产| 国产精品免费大片| 精品一区二区三卡| 免费久久久久久久精品成人欧美视频 | 黄色视频在线播放观看不卡| 人妻夜夜爽99麻豆av| 欧美 亚洲 国产 日韩一| 美女脱内裤让男人舔精品视频| 男男h啪啪无遮挡| 国产一区二区三区综合在线观看 | 久久精品久久久久久久性| 亚洲av中文av极速乱| 黄色视频在线播放观看不卡| 日韩一本色道免费dvd| a级毛片黄视频| 肉色欧美久久久久久久蜜桃| 99久国产av精品国产电影| 老司机影院毛片| 久久久久久久精品精品| 不卡视频在线观看欧美| 午夜免费鲁丝| 国产探花极品一区二区| 亚洲欧美日韩卡通动漫| 九草在线视频观看| 热re99久久国产66热| 亚洲中文av在线| 亚洲,一卡二卡三卡| 高清视频免费观看一区二区| 青青草视频在线视频观看| 亚洲无线观看免费| 黄色怎么调成土黄色| 亚洲精华国产精华液的使用体验| 中国国产av一级| 久久毛片免费看一区二区三区| 国产精品欧美亚洲77777| 精品一区二区三卡| 多毛熟女@视频| 精品人妻在线不人妻| 少妇的逼好多水| 久久狼人影院| 激情五月婷婷亚洲| 色5月婷婷丁香| 中国国产av一级| 日韩在线高清观看一区二区三区| 久久人人爽av亚洲精品天堂| 精品少妇黑人巨大在线播放| 日日摸夜夜添夜夜爱| 精品酒店卫生间| av在线播放精品| 不卡视频在线观看欧美| 美女福利国产在线| 欧美最新免费一区二区三区| 中文字幕精品免费在线观看视频 | av又黄又爽大尺度在线免费看| av在线老鸭窝| 国产男人的电影天堂91| 午夜视频国产福利| 久久午夜综合久久蜜桃| 十八禁网站网址无遮挡| 国产成人a∨麻豆精品| 亚洲av二区三区四区| 18禁在线无遮挡免费观看视频| 女人久久www免费人成看片| 久久国产亚洲av麻豆专区| 看十八女毛片水多多多| 日韩欧美精品免费久久| 亚洲国产精品一区二区三区在线| 最近2019中文字幕mv第一页| 日韩一本色道免费dvd| 一级黄片播放器| 免费人妻精品一区二区三区视频| 免费看av在线观看网站| 国产国语露脸激情在线看| 99热这里只有精品一区| 最近最新中文字幕免费大全7| 91成人精品电影| 精品久久蜜臀av无| 中文字幕免费在线视频6| 一区二区三区精品91| 熟女人妻精品中文字幕| 欧美 日韩 精品 国产| 蜜桃久久精品国产亚洲av| 在线看a的网站| 男男h啪啪无遮挡| 人妻 亚洲 视频| 菩萨蛮人人尽说江南好唐韦庄| 三级国产精品片| 狠狠婷婷综合久久久久久88av| 人人澡人人妻人| 天天影视国产精品| 精品国产乱码久久久久久小说| 最近的中文字幕免费完整| 成人亚洲精品一区在线观看| 成人漫画全彩无遮挡| 国产色婷婷99| 国产国语露脸激情在线看| 天堂俺去俺来也www色官网| 日本黄大片高清| 国产精品久久久久久av不卡| 亚洲欧美成人精品一区二区| 精品久久蜜臀av无| 黑人巨大精品欧美一区二区蜜桃 | 成人毛片a级毛片在线播放| 亚洲欧美中文字幕日韩二区| 婷婷色麻豆天堂久久| 成人二区视频| www.色视频.com| 天天影视国产精品| 狂野欧美激情性xxxx在线观看| 熟女电影av网| 国产 精品1| 亚洲精品乱码久久久久久按摩| 天堂俺去俺来也www色官网| 日本午夜av视频| 高清毛片免费看| 国产精品国产三级国产av玫瑰| 美女中出高潮动态图| 亚洲欧洲精品一区二区精品久久久 | av播播在线观看一区| 妹子高潮喷水视频| 亚洲av成人精品一二三区| 亚洲人成网站在线播| 99热网站在线观看| 精品视频人人做人人爽| 9色porny在线观看| 日韩,欧美,国产一区二区三区| 一二三四中文在线观看免费高清| 精品国产露脸久久av麻豆| 亚洲国产色片| 如日韩欧美国产精品一区二区三区 | 亚洲久久久国产精品| 日日爽夜夜爽网站| 熟女电影av网| 在线免费观看不下载黄p国产| 街头女战士在线观看网站| 蜜桃久久精品国产亚洲av| 免费不卡的大黄色大毛片视频在线观看| 精品一品国产午夜福利视频| 综合色丁香网| 亚洲成色77777| 色婷婷久久久亚洲欧美| 91国产中文字幕| 高清毛片免费看| 亚洲av福利一区| 亚洲人成网站在线观看播放| 啦啦啦视频在线资源免费观看| 成人亚洲精品一区在线观看| 国产片特级美女逼逼视频| 久久99精品国语久久久| 亚洲av在线观看美女高潮| 一级片'在线观看视频| 18+在线观看网站| 亚洲,一卡二卡三卡| 在线天堂最新版资源| 国产 精品1| 亚洲精品日韩av片在线观看| 国产色爽女视频免费观看| 精品久久国产蜜桃| 成年人免费黄色播放视频| 午夜91福利影院| 高清黄色对白视频在线免费看| 亚洲精品第二区| 久久精品国产自在天天线| 久久久a久久爽久久v久久| 伦理电影免费视频| 国产免费现黄频在线看| 亚洲欧美中文字幕日韩二区| 日本午夜av视频| 免费黄频网站在线观看国产| 久久午夜综合久久蜜桃| 夫妻性生交免费视频一级片| 亚洲人与动物交配视频| 国产国语露脸激情在线看| 国产色婷婷99| 精品少妇久久久久久888优播| 九九在线视频观看精品| 亚洲人与动物交配视频| 久久精品国产自在天天线| 国产成人精品一,二区| 久久午夜综合久久蜜桃| 色网站视频免费| 亚洲人与动物交配视频| 日韩不卡一区二区三区视频在线| 丝瓜视频免费看黄片| 少妇被粗大猛烈的视频| 亚洲av在线观看美女高潮| 亚洲一级一片aⅴ在线观看| 国产高清三级在线| 精品午夜福利在线看| 午夜影院在线不卡| 亚洲综合色网址| 成人午夜精彩视频在线观看| 一二三四中文在线观看免费高清| 人成视频在线观看免费观看| 飞空精品影院首页| 国产一区亚洲一区在线观看| 亚洲久久久国产精品| 亚洲综合精品二区| av在线观看视频网站免费| 久久精品国产鲁丝片午夜精品| 亚洲内射少妇av| 欧美日本中文国产一区发布| 中文乱码字字幕精品一区二区三区| 国产精品无大码| 国产高清三级在线| 视频中文字幕在线观看| 制服丝袜香蕉在线| 飞空精品影院首页| 九色成人免费人妻av| 欧美日本中文国产一区发布| 国产片特级美女逼逼视频| 日本午夜av视频| 日韩视频在线欧美| 亚洲欧美日韩卡通动漫| 成人毛片a级毛片在线播放| 纯流量卡能插随身wifi吗| 国产男女内射视频| 国产成人精品福利久久| 国产成人aa在线观看| 亚洲av成人精品一二三区| 80岁老熟妇乱子伦牲交| 美女福利国产在线| 午夜日本视频在线| 涩涩av久久男人的天堂| 久久综合国产亚洲精品| 熟女电影av网| 婷婷色av中文字幕| 777米奇影视久久| 久久久久国产网址| 国产精品99久久99久久久不卡 | 精品视频人人做人人爽| 久久国产亚洲av麻豆专区| 色94色欧美一区二区| 中文字幕精品免费在线观看视频 | 国产 一区精品| 久久精品国产a三级三级三级| 各种免费的搞黄视频| 精品一品国产午夜福利视频| 少妇的逼好多水| av视频免费观看在线观看| 国产精品久久久久久av不卡| 视频在线观看一区二区三区| videos熟女内射| 成人18禁高潮啪啪吃奶动态图 | 黄色毛片三级朝国网站| 亚洲精品乱码久久久久久按摩| 黑人巨大精品欧美一区二区蜜桃 | 久久久精品区二区三区| 五月玫瑰六月丁香| 欧美日本中文国产一区发布| 午夜激情av网站| 亚洲国产av影院在线观看| 中国美白少妇内射xxxbb| 夜夜骑夜夜射夜夜干| 国产熟女午夜一区二区三区 | 久久青草综合色| 欧美精品人与动牲交sv欧美| 亚洲精品一区蜜桃| 亚洲国产精品一区三区| 亚洲国产精品一区二区三区在线| 成年av动漫网址| 久久ye,这里只有精品| 国产av精品麻豆| 蜜桃国产av成人99| 国产亚洲午夜精品一区二区久久| 亚洲成色77777| 久久久精品区二区三区| 国产国拍精品亚洲av在线观看| 国产精品一区二区三区四区免费观看| 国产精品一区www在线观看| 国产成人91sexporn| 麻豆成人av视频| 国产精品嫩草影院av在线观看| 视频区图区小说| 在线播放无遮挡| 亚洲av免费高清在线观看| 五月伊人婷婷丁香| 久久久欧美国产精品| 中文字幕精品免费在线观看视频 | 大片电影免费在线观看免费| 欧美激情 高清一区二区三区| 国精品久久久久久国模美| 99精国产麻豆久久婷婷| 精品卡一卡二卡四卡免费| 满18在线观看网站| 韩国av在线不卡| 日本色播在线视频| 内地一区二区视频在线| 成人国语在线视频| 一级毛片aaaaaa免费看小| 精品少妇久久久久久888优播| 99国产综合亚洲精品|