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

    基于Smith-模糊PID解耦的固井水泥漿密度大滯后控制研究

    2025-08-07 00:00:00李勇樊恒冀宇程思達(dá)蒙健升
    石油機(jī)械 2025年7期
    關(guān)鍵詞:傳遞函數(shù)預(yù)估液位

    中圖分類(lèi)號(hào):TE256文獻(xiàn)標(biāo)識(shí)碼:ADOI:10.12473/CPM.202406064

    Long Time-delay Control ofCementSlurryDensityBased on Smith-fuzzy PID Decoupling

    LiYong1Fan Heng2JiYu2Cheng Sida'Meng Jiansheng2 (1.CNPCEngineeringTechnologyRamp;DCompanyLimited;2.ColegeofElectronicEngineeing,Xi’anShiyouUniesity)

    Abstract:The cement slurry mixing system has the characteristics of nonlinearity,time variation and long time-delay,which make it imposible to achieve stable and accurate control of automatic cementing devices by traditional techniques.This paper presents an advanced strategy of control based on Smith predictive compensation and fuzzy PID decoupling(Smith-fuzzy PID decoupling).Firstly,through the mathematical modeling and coupling analysis of the cement slurry mixing system,the feedforward compensator was used to compensate for the possible disturbances in the system inadvance,which effectively reduces the coupling effct.Then,the Smithfuzzy PID decoupling controler was designed,and the PID parameters were adjusted through fuzzy logic to adapt to the nonlinear and time-varying characteristics of the system,so as to achieve accurate control of density. Finally,simulation experiments and anti-interference analysis were caried out on the Matlab/Simulink platform.

    The results show that compared with Smith predictive compensation and fuzzy PID control,the proposed control based on Smith-fuzzy PID decoupling is significantly superior in dynamic performance,steady-state performance and anti-interference ability,and has a strong robustness.This control strategy can efectivelyreduce the overshoot,speed up the system adjustment,and improve the system stability.It provides areliable solution to the accuratecontrol ofcement slurry density,and will playatheoretical guiding role in the intellgent developmentof cement slurry mixing system.

    Keywords:cement slurry;Smith predictive control;density control; long time-delay;fuzzy control PI decoupling

    0 引言

    隨著石油勘探開(kāi)發(fā)領(lǐng)域的不斷拓展,大直徑井、深井及超深井的應(yīng)用日趨廣泛,這無(wú)疑對(duì)固井技術(shù)的智能化水平提出了更高的要求[1-2]。固井作業(yè)技術(shù)是油氣資源開(kāi)發(fā)過(guò)程中的關(guān)鍵技術(shù),直接關(guān)系到油氣開(kāi)采的效率和生產(chǎn)壽命[3]。水泥漿密度的精確控制是固井作業(yè)中尤為關(guān)鍵的一環(huán),直接影響著固井的質(zhì)量和效果,對(duì)油氣井的整體穩(wěn)定性和安全性產(chǎn)生深遠(yuǎn)影響。然而,由于現(xiàn)場(chǎng)條件的限制和測(cè)量技術(shù)的局限性,特別是存在大滯后問(wèn)題,一直是影響固井過(guò)程中的水泥漿密度智能控制的主要瓶頸。在工業(yè)過(guò)程控制中,大滯后過(guò)程往往包含顯著的純滯后環(huán)節(jié),對(duì)控制系統(tǒng)的穩(wěn)定性構(gòu)成嚴(yán)峻挑戰(zhàn)。鑒于時(shí)滯問(wèn)題的復(fù)雜性和挑戰(zhàn)性,眾多學(xué)者致力于尋找切實(shí)可行的解決策略。

    韓曹悅等4提出了一種改進(jìn)的二自由度內(nèi)??刂?,能有效改善大時(shí)滯對(duì)象的控制效果,提高了系統(tǒng)的魯棒性。苗河泉等基于改進(jìn)型Smith預(yù)估控制對(duì)拖拉機(jī)自動(dòng)轉(zhuǎn)向系統(tǒng)存在的時(shí)滯、外部干擾等進(jìn)行了補(bǔ)償設(shè)計(jì),結(jié)果表明,該系統(tǒng)能較好地跟蹤參考信號(hào),且能在受到干擾后較快地恢復(fù)到穩(wěn)定狀態(tài)。任長(zhǎng)輝等采用SmithPSO-BP-PID反饋控制解決熱軋帶鋼層流冷卻精冷段的滯后問(wèn)題,顯著提高了控制效果。張磊等[考慮到煤炭地下氣化反應(yīng)具有時(shí)變性和大滯后特點(diǎn),設(shè)計(jì)了內(nèi)??刂扑惴ㄒ詰?yīng)對(duì)地下氣化過(guò)程中模型不確定性的挑戰(zhàn)。張皓等[8]針對(duì)溫度模型存在的大滯后現(xiàn)象,為提升控制精度和響應(yīng)速度,提出了一種創(chuàng)新的控制策略—基于T-S型模型模糊神經(jīng)網(wǎng)絡(luò)的內(nèi)模自適應(yīng)PID控制方法;該方法結(jié)合了模糊神經(jīng)網(wǎng)絡(luò)的強(qiáng)大逼近能力和內(nèi)模控制的預(yù)測(cè)補(bǔ)償特性,以及自適應(yīng)PID的靈活調(diào)整能力,從而有效應(yīng)對(duì)溫度模型的大滯后問(wèn)題。

    盡管現(xiàn)有控制方法在處理時(shí)滯方面取得了進(jìn)展,但在復(fù)雜非線性系統(tǒng)中仍面臨適應(yīng)性和參數(shù)調(diào)整的挑戰(zhàn)。為此,筆者提出了一種融合Smith預(yù)估補(bǔ)償和模糊PID控制的策略,旨在有效解決大滯后和多變量耦合問(wèn)題,實(shí)現(xiàn)水泥漿密度的精準(zhǔn)穩(wěn)定控制。

    1固井水泥漿混漿系統(tǒng)分析

    1.1固井水泥漿混漿系統(tǒng)工作原理

    干灰通過(guò)計(jì)量閥送入高能混合器,清水由離心泵輸送至混合器形成高速水流,吸人干灰。干灰與水在混合器初步混合后,經(jīng)攪拌器進(jìn)一步攪拌,水泥漿循環(huán)泵加壓輸送至混漿槽,進(jìn)行多次混合。傳感器實(shí)時(shí)監(jiān)測(cè)水泥漿密度等信息,PLC根據(jù)數(shù)據(jù)調(diào)整執(zhí)行器,維持密度穩(wěn)定。符合要求的水泥漿由灌注泵輸送至固井泵,注入井筒,完成固井作業(yè)。具體流程部分如圖1所示。

    圖1固井水泥漿混漿系統(tǒng)自動(dòng)控制原理示意圖Fig.1Schematic diagram of automatic control oncement slurry mixing

    1.2數(shù)學(xué)模型的推導(dǎo)過(guò)程

    采用機(jī)理建模法以及結(jié)合工程經(jīng)驗(yàn)對(duì)固井水泥漿系統(tǒng)建立數(shù)學(xué)模型。在固井水泥漿連續(xù)混拌的過(guò)程中,根據(jù)流量守恒關(guān)系,在穩(wěn)態(tài)工作點(diǎn)[9],其穩(wěn)態(tài)方程式可表示為:

    Qw0+Qc0-Qs0=0

    式中: Qw0 為輸入穩(wěn)態(tài)水流量, m3/min ; Qc0 為輸入 穩(wěn)態(tài)干灰流量, m3/min ; Qs0 為輸出穩(wěn)態(tài)水泥漿流

    量, m3/min 。

    動(dòng)態(tài)條件下的動(dòng)態(tài)方程式為:

    根據(jù)質(zhì)量守恒定律,建立水泥漿連續(xù)混拌系統(tǒng)質(zhì)量守恒方程:

    式中: Qw 為輸入水流量, m3/min ; Qc 為輸入干灰流量, m3/min ; Qs 為輸出水泥漿流量, m3/min : ρw 為清水密度, kg/m3 : ρc 為干灰密度, kg/m3 : ρs 為水泥漿密度, kg/m3 ; A 為混漿槽底面積, m2 h 為混漿槽的液位高度, m : χt 為時(shí)間變量,s。

    水泥槳從混漿槽流出的流量與混漿槽的體積和底面積相關(guān)。

    在系統(tǒng)工作點(diǎn)狀態(tài)附近小范圍內(nèi)將系統(tǒng)進(jìn)行線性化處理,可以得到:

    式中: Δh(t) 為 Φt 時(shí)刻液位高度差, m ; h0 為穩(wěn)態(tài)液 位高度, m 。

    那么式(2)和式(3)就可以寫(xiě)為:

    式中: ΔQw(t) 為 χt 時(shí)刻水流量的增量; ΔQc(t) 為干灰流量的增量; 為液位增量的一階導(dǎo)數(shù); ρs0 為初始穩(wěn)態(tài)水泥漿密度; Δρs(t) 為 Φt 時(shí)刻水泥漿密度的增量; 為 Δρs(t) 的一階導(dǎo)數(shù)。

    固井水泥混漿過(guò)程是一個(gè)具有非線性、強(qiáng)耦合且存在大滯后的過(guò)程,其數(shù)學(xué)模型可以視為一階慣性延遲系統(tǒng),其傳遞函數(shù)描述如下:

    式中: K 為被控對(duì)象的放大系數(shù); T 為慣性時(shí)間常數(shù); τ 為滯后時(shí)間常數(shù); s 為頻率復(fù)變量。

    經(jīng)推導(dǎo)可知固井水泥漿混漿系統(tǒng)傳遞函數(shù)如下:

    式中: H(s) 為混漿槽的液位高度, m . ρs(s) 為水泥漿混漿密度, kg/m3 ; Qw(s) 為清水流量, m3/min 為干灰流量, m3/min 。

    已知系統(tǒng)傳遞函數(shù),就可以依據(jù)其特性對(duì)系統(tǒng)進(jìn)行精確的仿真和控制;在不影響實(shí)際生產(chǎn)的前提下評(píng)估各種控制算法的性能,可以識(shí)別潛在的穩(wěn)定性問(wèn)題;并對(duì)控制系統(tǒng)進(jìn)行魯棒性分析,以確定其對(duì)過(guò)程變化、外部干擾和模型不確定性的適應(yīng)能力,進(jìn)而為實(shí)際操作提供指導(dǎo)。

    1.3解耦控制分析與設(shè)計(jì)

    在實(shí)際工業(yè)控制過(guò)程中,許多過(guò)程參數(shù)都隨生產(chǎn)流程的持續(xù)進(jìn)行而不斷變化。同時(shí),這些參數(shù)之間存在復(fù)雜的多變量耦合作用[10],這些耦合作用的關(guān)系往往不可避免,所以導(dǎo)致過(guò)程控制達(dá)不到理想的控制效果。

    固井水泥漿混漿系統(tǒng)中的液位控制和密度控制之間存在著顯著的耦合作用。液位控制通過(guò)實(shí)時(shí)采集混漿罐液位、排量等數(shù)據(jù),自動(dòng)調(diào)節(jié)清水和干灰的計(jì)量閥,將混漿罐內(nèi)液位控制在設(shè)定范圍內(nèi)。液位發(fā)生變化時(shí),混合罐內(nèi)的水泥槳體積與密度會(huì)隨之變化,從而影響水泥槳的密度。而密度控制則調(diào)節(jié)混漿中干灰和水的比例,以確?;鞚{的密度符合預(yù)定要求。同樣,水泥漿的密度波動(dòng)也會(huì)影響液位的穩(wěn)定性。密度和液位耦合控制原理圖如2所示。

    圖2水泥漿密度和液位耦合控制原理圖Fig.2 Schematic diagram of couplingcontrol on cement slurry density and liquid level

    圖2中, Gcl(s) 為液位控制回路的PID控制器傳遞函數(shù), Gc2(s) 為密度控制回路的PID控制器傳遞函數(shù), G11(s) 為液位控制的傳遞函數(shù), G22(s) 為密度控制的傳遞函數(shù), G12(s) 為密度控制對(duì)液位控制產(chǎn)生影響的傳遞函數(shù), G21(s) 為液位控制對(duì)密度控制產(chǎn)生影響的傳遞函數(shù)。

    通過(guò)計(jì)算系統(tǒng)的相對(duì)增益,判斷通道間的相互影響程度,評(píng)估不同通道的性能表現(xiàn),進(jìn)而選擇最適合特定應(yīng)用的控制策略。

    假設(shè)在多變量解耦系統(tǒng)中有 q 個(gè)輸入和 p 個(gè)輸出,則其傳遞函數(shù)陣可表示為:

    式中: Y(s) 為輸出信號(hào); U(s) 為輸人信號(hào); Gpq(s)

    為第 q 個(gè)輸入與第 p 個(gè)輸出間的傳遞函數(shù)。

    然后引入相對(duì)增益來(lái)討論各個(gè)輸入輸出之間的關(guān)聯(lián)程度。相對(duì)增益 λij 表示第 j 個(gè)輸入 u 對(duì)第 i 個(gè)輸出 y 的影響力,其定義如下所示:

    器如下:

    式中: GD1(s) 為液位控制回路的前饋補(bǔ)償器;GD2(s) 為密度控制回路的前饋補(bǔ)償器

    通過(guò)前饋補(bǔ)償解耦,可以預(yù)先對(duì)系統(tǒng)中可能出現(xiàn)的擾動(dòng)進(jìn)行補(bǔ)償,有效降低或消除這些耦合效應(yīng)[13],使得系統(tǒng)更加接近于理想的無(wú)耦合狀態(tài),從而改善系統(tǒng)的穩(wěn)定性。

    2 控制算法原理

    式中:分子為第一放大系數(shù),表示其他輸入 um 保持不變(開(kāi)環(huán))的情況下 uj 的改變對(duì) yi 產(chǎn)生的增益;分母為第二放大系數(shù),表示在其他輸出 ym 不變(閉環(huán))的情況下 uj 的改變對(duì) yi 產(chǎn)生的增益; ?m 為其他變量; c 為第 m 個(gè)物理量的實(shí)際值。

    固井水泥漿混漿系統(tǒng)是一個(gè)雙輸人雙輸出的耦合過(guò)程,因此采用解析法來(lái)求解。系統(tǒng)輸入輸出關(guān)系可表示為如下形式:

    式中: y1 ! y2 分別為輸出1、2; Kij 為第 i 個(gè)被控量Δui 相對(duì)于第 j 個(gè)控制變量 Δyj 的靜態(tài)增益; 分別為輸入1、2。

    則可以求得系統(tǒng)的相對(duì)增益矩陣:

    相對(duì)增益矩陣反映了系統(tǒng)中各輸入與輸出之間的相對(duì)影響程度[]。相對(duì)增益矩陣的每個(gè)元素表示了系統(tǒng)中某一輸人對(duì)某一輸出的相對(duì)影響程度。當(dāng)矩陣中的元素接近1時(shí),表示該輸入對(duì)該輸出的影響較為直接和顯著;而當(dāng)元素遠(yuǎn)離1時(shí),則表示該輸人與輸出之間存在較強(qiáng)的耦合關(guān)系[12]。由此可知,液位控制和密度控制之間存在較強(qiáng)的耦合性,應(yīng)該采用多變量解耦控制方法對(duì)水泥混漿系統(tǒng)進(jìn)行分析。

    采用需要更少裝置數(shù)的前饋補(bǔ)償解耦方法對(duì)液位控制和密度控制進(jìn)行解耦,得到的前饋補(bǔ)償

    2.1 模糊PID自適應(yīng)原理

    模糊PID控制是一種結(jié)合了模糊邏輯和傳統(tǒng)PID控制的方法,旨在通過(guò)實(shí)時(shí)優(yōu)化PID參數(shù)來(lái)提高控制系統(tǒng)的性能[14]。其基本原理是:將控制系統(tǒng)中的誤差、誤差變化率以及誤差累積量作為輸人,經(jīng)過(guò)模糊化處理后輸人到模糊控制器中[15];應(yīng)用模糊推理機(jī)制,結(jié)合模糊輸入變量和規(guī)則,推導(dǎo)出模糊輸出變量的集合;最后,通過(guò)去模糊化的方法將模糊輸出轉(zhuǎn)化為具體的控制命令,作為PID控制器的輸出。模糊PID控制原理圖如圖3所示。

    圖3模糊PID控制原理圖Fig.3Schematic diagram offuzzyPID control

    當(dāng)系統(tǒng)外界環(huán)境發(fā)生變化時(shí),系統(tǒng)根據(jù)2個(gè)輸人量動(dòng)態(tài)調(diào)節(jié)控制參數(shù) KP 、 KI 、 Kp ,使控制系統(tǒng)仍能保證良好控制精度。其中:

    式中: KP0 : KI0 、 KD0 分別為初始PID整定參數(shù);

    ΔKP 、 ΔKI : ΔKp 分別為模糊控制調(diào)節(jié)參數(shù)的增量。

    2.2Smith預(yù)估控制基本原理

    Smith預(yù)估控制是一種經(jīng)典的補(bǔ)償控制策略,用于解決具有滯后環(huán)節(jié)的系統(tǒng)的控制問(wèn)題。Smith預(yù)估控制的基本思想是:在控制器中引入一個(gè)模型,用于預(yù)測(cè)系統(tǒng)的輸出,并根據(jù)預(yù)測(cè)值進(jìn)行控制信號(hào)的補(bǔ)償[16-17]。

    含有純滯后的被控對(duì)象閉環(huán)控制系統(tǒng)如圖4所示。圖4中, R(s) 為系統(tǒng)的輸入, E(s) 為誤差信號(hào),U(s) 為控制器輸出的控制信號(hào), Y(s) 為系統(tǒng)的輸出,Gc(s) 為系統(tǒng)的控制器, GP(s) 為被控對(duì)象去除純滯后時(shí)間部分的傳遞函數(shù), e-τs 為系統(tǒng)的時(shí)滯環(huán)節(jié)。

    圖4含有純滯后的被控對(duì)象閉環(huán)控制系統(tǒng) Fig.4Closed-loop control system for controlled objectwithpuredelay

    那么系統(tǒng)的傳遞函數(shù)可表示為:

    由于系統(tǒng)傳遞函數(shù)的分母中包含純滯后環(huán)節(jié),而純滯后環(huán)節(jié)會(huì)使系統(tǒng)對(duì)輸入信號(hào)的反應(yīng)速度變慢,影響系統(tǒng)的動(dòng)態(tài)性能和跟蹤精度。所以引入一個(gè)并聯(lián)于被控對(duì)象的補(bǔ)償器,即Smith預(yù)估補(bǔ)償器 Gm(s) [18]。通過(guò)加入Smith預(yù)估補(bǔ)償器,可以對(duì)被控對(duì)象的傳遞函數(shù)進(jìn)行適當(dāng)?shù)恼{(diào)整,從而減小純滯后環(huán)節(jié)對(duì)系統(tǒng)性能的影響。帶有Smith預(yù)估器的系統(tǒng)如圖5所示。

    圖5純滯后補(bǔ)償閉環(huán)控制系統(tǒng) Fig.5Closed-loop control system with pure delay compensation

    控制量 U(s) 與反饋量 Y(s) 間的傳遞函數(shù)為:

    因此可以得到Smith預(yù)估補(bǔ)償器為:

    Gm(s)=GP(s)(1-e-τs

    然而實(shí)際上,Smith預(yù)估補(bǔ)償器雖然在邏輯上與被控對(duì)象相關(guān)聯(lián),但在物理連接上并非直接并聯(lián)。相反,它是反向連接在控制器 Gc(s) 上的[19-20]。實(shí)際的純滯后補(bǔ)償控制系統(tǒng)如圖6所示。

    圖6實(shí)際純滯后補(bǔ)償閉環(huán)控制系統(tǒng) Fig. 6Actual closed-loop control system with pure delay compensation

    通過(guò)引入Smith預(yù)估補(bǔ)償器,可以有效消除時(shí)滯因子,使得控制器能夠更快地響應(yīng)系統(tǒng)的變化??刂破髂軌蛱崆邦A(yù)測(cè)系統(tǒng)輸出的變化,以適應(yīng)被控對(duì)象的動(dòng)態(tài)變化并減小誤差,提升控制系統(tǒng)的穩(wěn)定性和準(zhǔn)確性。

    3控制算法設(shè)計(jì)與仿真

    3.1模糊PID解耦控制設(shè)計(jì)與仿真

    設(shè)計(jì)采用二維模糊控制器,能夠同時(shí)處理輸入變量偏差和偏差的變化率。設(shè)輸人變量 e ! ec 與輸出變量 ΔKr 、 ΔKI : ΔKp 的模糊集為{NB,NM,NS,ZO,PS,PM, 。根據(jù)工況,液位控制的基本論域?yàn)? , , 0.4], ΔKp (204號(hào) =[-0.1 0.1]。密度控制的基本論域?yàn)?e= [-1 8001800], , ΔKp= [-0.0080.008], ΔKD=[-0.0010.001] 。對(duì)于水泥漿液位和密度這類(lèi)具有一定波動(dòng)范圍的控制對(duì)象,三角形隸屬度函數(shù)能夠較好地描述輸人變量與輸出變量之間的模糊關(guān)系,因此隸屬度函數(shù)采用三角形函數(shù)[21-23]。

    針對(duì)固井水泥混漿系統(tǒng)工況和專(zhuān)家經(jīng)驗(yàn),通過(guò)計(jì)算當(dāng)前系統(tǒng)誤差和誤差變化量,總結(jié)參數(shù)ΔKP 、 ΔKI 1 ΔKp 的整定規(guī)律,如表1所示。

    表1 ΔKp.ΔKI.ΔKp"模糊控制規(guī)則表

    最后經(jīng)解模糊化處理將模糊控制器輸出的模糊值轉(zhuǎn)換為精確的控制信號(hào),以便于驅(qū)動(dòng)清水比例閥和干灰比例閥,進(jìn)而保持混漿槽液位和水泥漿密度的穩(wěn)定。

    基于模糊PID解耦控制的固井混漿液位密度控制仿真結(jié)構(gòu)圖如圖7所示。

    其中液位控制的模糊PID控制器封裝模塊,即控制器的封裝子系統(tǒng)如圖8所示。

    圖7模糊PID解耦控制仿真結(jié)構(gòu)圖
    Fig.7Simulation structure of fuzzy PID decoupling control圖8模糊PID控制器的封裝子系統(tǒng)

    系統(tǒng)運(yùn)行后得到的液位控制和密度控制仿真結(jié)果如圖9和圖10所示。對(duì)于液位控制曲線,其上升時(shí)間為 5.762s ,峰值時(shí)間是圖9中1號(hào)標(biāo)為8.034s ,超調(diào)量為 12.625% 。這意味著系統(tǒng)在達(dá)到穩(wěn)態(tài)時(shí)波動(dòng)較小,控制更加穩(wěn)定。同時(shí),調(diào)節(jié)時(shí)間是圖9中2號(hào)標(biāo)為33.951s,保證了系統(tǒng)能夠快速恢復(fù)到穩(wěn)定狀態(tài)。

    對(duì)于密度控制曲線,其上升時(shí)間為8.921s,峰值時(shí)間是圖10中1號(hào)標(biāo)為 12.838s ,超調(diào)量?jī)H為 1.125% 。這表明密度控制在達(dá)到穩(wěn)態(tài)時(shí)幾乎沒(méi)有波動(dòng),控制精度極高。此外,調(diào)節(jié)時(shí)間是圖10中2號(hào)標(biāo)為26.539s,系統(tǒng)的響應(yīng)速度較快。

    Fig.8Encapsulated subsystem of a fuzzyPIDcontroller圖9液位模糊PID解耦控制仿真曲線 Fig. 9 Simulation curve of fuzzy PID decoupling control onliquid level
    圖10密度模糊PID解耦控制仿真曲線Fig.10Simulation curve of fuzzy PID decouplingcontrol on density

    模糊PID解耦控制方法在固井水泥漿控制系統(tǒng)中確實(shí)取得了顯著的效果,不僅提高了控制精度,還增強(qiáng)了系統(tǒng)的穩(wěn)定性。雖然該方法已經(jīng)大幅降低了超調(diào)量,但超調(diào)現(xiàn)象仍未完全消除,意味著控制性能仍有提升空間。

    3.2Smith-模糊PID解耦控制設(shè)計(jì)與仿真

    采用Smith-模糊PID解耦控制方法,旨在進(jìn)一步降低超調(diào)量并改善系統(tǒng)的動(dòng)態(tài)響應(yīng)。Smith預(yù)估器預(yù)測(cè)水泥漿密度的變化趨勢(shì),為模糊PID控制提供先驗(yàn)信息;而模糊PID控制則根據(jù)預(yù)測(cè)結(jié)果和實(shí)時(shí)狀態(tài)動(dòng)態(tài)調(diào)整控制策略,其控制仿真結(jié)構(gòu)圖如圖11所示。得到的液位控制和密度控制仿真曲線如圖12和圖13所示。

    圖11Smith-模糊PID解耦控制仿真結(jié)構(gòu)圖
    Fig.11Simulation structure of Smith-fuzzy PID decoupling control

    由圖12和圖13可知,液位控制曲線的上升時(shí)間是圖12中1號(hào)標(biāo)為 6.783s ,雖然相較于之前的控制方法略有增加,但系統(tǒng)無(wú)超調(diào)現(xiàn)象。此外,圖12中2號(hào)標(biāo)調(diào)節(jié)時(shí)間也縮短至30.4s,這顯示出系統(tǒng)迅速恢復(fù)穩(wěn)定狀態(tài)的能力。密度控制曲線的上升時(shí)間是圖13中1號(hào)標(biāo)為12.107s,無(wú)超調(diào)現(xiàn)象。圖13中2號(hào)標(biāo)調(diào)節(jié)時(shí)間為 32.058s ,進(jìn)一步確保了系統(tǒng)能夠快速響應(yīng)并達(dá)到穩(wěn)定狀態(tài)。

    圖12液位Smith-模糊PID解耦控制仿真曲線
    圖13密度Smith-模糊PID解耦控制仿真曲線Fig.13Simulation curves of Smith-fuzzy PID decoupling control on density

    3.3 仿真結(jié)果分析

    將PID控制、PID解耦控制、Smith-PID解耦控制和Smith-模糊PID解耦控制這4種控制方法進(jìn)行比較,由文獻(xiàn)[10]所得數(shù)據(jù),得到系統(tǒng)仿真結(jié)果性能指標(biāo)對(duì)比如表2所示。由表2可知,以上控制方法均能夠使曲線達(dá)到預(yù)設(shè)值,傳統(tǒng)的PID控制的超調(diào)量較大,為 46.1% ,調(diào)節(jié)時(shí)間長(zhǎng),為123.503s,這表明在沒(méi)有解耦的情況下,系統(tǒng)受到液位和密度之間耦合效應(yīng)的影響較大。在引入前饋解耦方法后,消除了液位和密度之間的耦合效應(yīng),性能改善,但超調(diào)仍存在。

    表2仿真結(jié)果性能指標(biāo)對(duì)比Table 2 Comparison of performance indicators of simulation results

    為了進(jìn)一步優(yōu)化系統(tǒng)性能,采用了模糊PID控制方法。該方法融合了模糊邏輯與PID控制的優(yōu)勢(shì),能夠基于系統(tǒng)實(shí)時(shí)狀態(tài)自動(dòng)調(diào)整控制參數(shù)。模糊PID控制策略實(shí)施后,系統(tǒng)超調(diào)量有所減小,但未徹底消除,調(diào)節(jié)時(shí)間大約穩(wěn)定在 30s 。最后結(jié)合Smith預(yù)估補(bǔ)償器和模糊PID控制,即Smith-模糊PID解耦控制,通過(guò)這些改進(jìn),控制系統(tǒng)的性能已經(jīng)優(yōu)化至最佳狀態(tài)。Smith預(yù)估器能夠提前預(yù)測(cè)系統(tǒng)輸出的變化趨勢(shì),從而允許模糊PID控制器根據(jù)預(yù)測(cè)結(jié)果和實(shí)時(shí)反饋信號(hào)進(jìn)行更加精準(zhǔn)的控制。這種結(jié)合使得控制系統(tǒng)能夠更好地應(yīng)對(duì)不確定性,減少超調(diào)現(xiàn)象,并加快調(diào)節(jié)速度。盡管系統(tǒng)上升時(shí)間有所延長(zhǎng),但超調(diào)量已基本消除,并且調(diào)節(jié)時(shí)間變得更短,這有助于系統(tǒng)保持良好的穩(wěn)定性和可靠性。

    鑒于前饋解耦控制在消除液位和密度控制耦合效應(yīng)方面的有效性,將前饋解耦控制策略添加到以下3種方法:PID控制、Smith-PID控制和Smith-模糊PID控制,并進(jìn)行比較,得到的液位控制和密度控制的仿真曲線如圖14和圖15所示。從圖14和圖15可以看出,Smith-模糊PID解耦控制能使控制系統(tǒng)的性能狀態(tài)達(dá)到最優(yōu)。

    圖14液位控制仿真對(duì)比曲線
    Fig.14Comparison curve of liquid level control simulation圖15密度控制仿真對(duì)比曲線Fig.15Comparison curve of density control simulation

    為分析系統(tǒng)的抗干擾能力,設(shè)置仿真時(shí)間為300s ,待系統(tǒng)按照設(shè)定值達(dá)到穩(wěn)定狀態(tài)后,在仿真時(shí)間到達(dá)100s時(shí)給系統(tǒng)添加一個(gè)階躍輸入信號(hào)來(lái)分析系統(tǒng)的抗干擾能力,結(jié)果如圖16和圖17所示。根據(jù)圖16和圖17,干擾信號(hào)的引人導(dǎo)致控制系統(tǒng)出現(xiàn)了不同程度的震蕩,但最終所有系統(tǒng)都逐漸恢復(fù)到了穩(wěn)定狀態(tài)。

    相較于PID解耦和模糊PID解耦控制,Smith-模糊PID解耦控制曲線在受到干擾后,受影響程度最小,且恢復(fù)速度更快,顯示出更強(qiáng)的魯棒性。Smith-模糊PID解耦控制方法有效地減少了超調(diào)量,其中系統(tǒng)液位的超調(diào)量分別下降了 13.95% 和10.7% ,幾乎實(shí)現(xiàn)了零超調(diào);其系統(tǒng)密度超調(diào)量分別降低了 6.1875% 和 0.25% ,同樣接近于零超調(diào)。在這3種控制策略中,Smith-模糊PID控制系統(tǒng)展現(xiàn)出了更優(yōu)的動(dòng)態(tài)性能和更佳的穩(wěn)態(tài)性能,同時(shí)具備了更強(qiáng)的抗干擾能力,從而確保了系統(tǒng)的高效穩(wěn)定性。

    圖16液位控制加入干擾仿真對(duì)比曲線
    圖17密度控制加入干擾仿真對(duì)比曲線Fig.17Comparison curve of interference simulation with density control

    通過(guò)對(duì)比這4種控制方法可以得出:Smith-模糊PID解耦控制在減小超調(diào)、加快調(diào)節(jié)速度以及提高系統(tǒng)穩(wěn)定性方面表現(xiàn)最優(yōu)。因此,在固井水泥漿控制系統(tǒng)中,Smith-模糊PID解耦控制策略能夠取得更加理想的控制效果。

    4結(jié)束語(yǔ)

    針對(duì)固井水泥漿密度控制大滯后問(wèn)題設(shè)計(jì)了控制算法,提出了一種結(jié)合Smith預(yù)估補(bǔ)償和模糊PID控制的解耦控制策略,通過(guò)Smith預(yù)估器補(bǔ)償系統(tǒng)滯后,模糊PID算法動(dòng)態(tài)調(diào)整參數(shù),并解耦液位與密度的耦合關(guān)系,解決了傳統(tǒng)方法中超調(diào)嚴(yán)重和控制精度不足的難題。

    結(jié)果表明,相較于PID解耦和模糊PID解耦控制方法,采用Smith-模糊PID解耦控制能夠顯著減少超調(diào)量,系統(tǒng)液位超調(diào)量分別降低了 13.95% 和10.7% ,幾乎實(shí)現(xiàn)了零超調(diào);密度超調(diào)量分別降低了 6.1875% 和 0.25% ,同樣接近于零超調(diào)。

    綜合仿真結(jié)果,Smith-模糊PID解耦控制算法在動(dòng)態(tài)性能、穩(wěn)態(tài)性能及抗干擾能力方面表現(xiàn)出顯著的優(yōu)勢(shì),能夠?qū)崿F(xiàn)密度的穩(wěn)定精確控制,且無(wú)超調(diào)現(xiàn)象發(fā)生,為水泥漿密度控制提供了可靠的解決方案。

    參考文獻(xiàn)

    [1]于永金,夏修建,王治國(guó),等.深井、超深井固井關(guān)鍵技術(shù)進(jìn)展及實(shí)踐[J].新疆石油天然氣,2023,19(2):24-33.YUYJ,XIA XJ,WANG ZG, et al.Progress andapplication of the key technologies of deep and ultra-deepwell cementing[J].Xinjiang Oil amp; Gas,2023,19(2):24-33.

    [2]廖權(quán)文,胡建均,史懷忠,等.文23儲(chǔ)氣庫(kù)鉆井工程關(guān)鍵技術(shù)[J].石油鉆采工藝,2023,45(2):160-166.LIAOQW,HUJJ,SHIHZ,etal.Key technolo-gies in drilling engineering of Wen 23 underground gasstorage[J].Oil Drillingamp; Production Technology,2023,45(2):160-166.

    [3]張懷文,馬立國(guó),王琦,等.固井自動(dòng)監(jiān)控水泥頭及閘閥系統(tǒng)研制與應(yīng)用[J].石油機(jī)械,2022,50(6):16-21.ZHANG H W, MA LG,WANGQ,et al.Develop-ment and application of cement head and gate valvesystem for automatic monitoring cementing [J]. ChinaPetroleumMachinery,2022,50(6):16-21.

    [4]韓曹悅,汪思源,王文標(biāo).大時(shí)滯過(guò)程下的改進(jìn)二自由度內(nèi)模控制[J/OL].控制工程,2025:1-8(2024-04-28)[2024-06-01].https://doi.org/10.14107/j.cnki.kzgc.20210320.HAN CY,WANG SY,WANG WB. An improved two-degree-of-freed-om internal model control for processeswith large time delay[J/OL].Control Engineering ofChina, 2025: 1-8 (2024-04-28)[2024-06-01].https://doi.org/10.14107/j.cnki. kzgc.20210320.

    [5]苗河泉,刁培松,徐廣飛,等.基于改進(jìn)史密斯預(yù)估控制的電液轉(zhuǎn)向時(shí)滯補(bǔ)償研究[J].農(nóng)機(jī)化研究,2023,45(7):232-237.MIAO HQ,DIAO PS,XUGF,et al. Research on time-delaycompensation of electro-hydraulic steeringbased onimproved smith predictor control [J]. Journal of Agricul-tural Mechanization Research,2023,45(7):232-237.

    [6]任長(zhǎng)輝,田海,陶震,等.熱軋帶鋼層流冷卻的精冷段SmithPSO-BP-PID反饋控制[J].鍛壓技術(shù),2023,48(12):177-181.RENCH,TIANH,TAOZ,et al.SmithPSO-BP-PID feedback control on fine cold section for laminar cool-ing of hot rolled strip steel [J].Forgingamp; StampingTechnology,2023,48(12):177-181.

    [7]張磊,于建林,周明軍,等.一種模型不確定的煤炭地下氣化控制算法研究[J].天然氣與石油,2023,41 (1):141-146.ZHANG L, YU JL, ZHOU MJ,et al. Study on con-trol algorithm of undergroundcoal gasification processwith model uncertainty[J].Natural Gas and Oil,2023,41 (1):141-146.

    [8]張皓,高瑜翔,唐軍,等.基于T-S型模糊神經(jīng)網(wǎng)絡(luò)的內(nèi)模自適應(yīng)PID算法研究[J].工業(yè)儀表與自動(dòng)化裝置,2021 (6):118-124.ZHANGH,GAO YX,TANGJ,et al.Research oninternal model adaptive PID algorithm based on T-S fuzzyneural network [J]. Industrial Instrumentation amp; Auto-mation,2021 (6):118-124.

    [9]樊恒,冀宇,程思達(dá),等.基于前饋解耦的固井水泥漿密度控制仿真研究[J].石油管材與儀器,2024,10 (2): 47-53.FAN H,JIY,CHENGSD,et al. Simulation on den-sitycontrol of cementing slurry based on feed - forwarddecoupling [J].Petroleum Instruments,2024,10(2):47-53.

    [10]劉啟蒙,賈濤,崔俊國(guó),等.大滯后固井系統(tǒng)預(yù)估補(bǔ)償控制算法研究[J].自動(dòng)化與儀表,2023,38(10):100-103,114.LIUQM,JIA T,CUI JG,et al.Research on autono-mous control algorithm of cementing system with largedelay[J].Automationamp; Instrumentation,2023,38(10):100-103,114.

    [11]謝震,楊曙昕,代鵬程,等.構(gòu)網(wǎng)型全功率風(fēng)電機(jī)組網(wǎng)側(cè)變流器耦合分析及抑制策略[J].電工技術(shù)學(xué)報(bào),2023,38(14):3745-3758,3768.XIE Z,YANG SX,DAI P C,et al. Grid-Side cou-pling analysis and suppression strategy of Grid-FormingFull-Power wind turbines [J]. Transactions of ChinaElectrotechnical Society,2023,38(14):3745-3758,3768.

    [12]武浩,陳海勝,黃克謹(jǐn).內(nèi)部熱耦合中間隔壁塔的設(shè)計(jì)與控制[J].北京化工大學(xué)學(xué)報(bào)(自然科學(xué)版),2023,50(4):10-20.WU H,CHEN HS,HUANG K J. Design and controlofaninternally heat integratedmiddle dividing wallcolumn[J]. Journal of Beijing University of ChemicalTechnology(Natural Science Edition),2023,50(4):10-20.

    [13]吳瓊,蘇建徽,解寶,等.基于最大凈功率輸出的PEMFC陰極供氣系統(tǒng)優(yōu)化控制研究[J].太陽(yáng)能學(xué)報(bào),2024,45(2):283-290.WUQ,SUJH,XIE B,et al. Study on optimal con-trol of cathode gas supply system of PEMFC based onmaximum net power output [J].Acta Energiae SolarisSinica,2024,45(2):283-290.

    [14]韓愛(ài)國(guó),宋福豪,田韶鵬,等.燃料電池系統(tǒng)建模與供氣系統(tǒng)控制方法[J].江蘇大學(xué)學(xué)報(bào)(自然科學(xué)版),2024,45(2):147-153.HAN AG,SONGFH,TIAN SP,et al.Modelingof fuel cell system and control method of air supplysystem[J]. Journal of Jiangsu University (Natural Sci-enceEdition),2024,45(2):147-153.

    [15]陳歲繁,楊松,李其朋.基于BP-NSGA-Ⅱ優(yōu)化的高速電梯轎廂水平振動(dòng)變論域模糊PID控制[J].噪聲與振動(dòng)控制,2024,44(2):63-69,81.CHEN S F,YANG S,LI Q P.Variable universe fuzzyPID control of horizontal vibration of high-speed elevatorCars based on BP-NSGA-II optimization [J].Noiseand Vibration Control,2024,44(2):63-69,81.

    [16]冷紅菱,袁致濤,盧冀偉,等.高壓輥磨過(guò)程的料倉(cāng)質(zhì)量自動(dòng)控制研究[J].金屬礦山,2023(6):178-183.LENG HL,YUAN ZT,LUJW,et al.Research onautomatic control of silo weight in high-pressure grindingrolls process [J].Metal Mine,2023(6):178-183.

    [17]莫岳平,徐迪.基于自抗擾控制的異步電機(jī)泵水系統(tǒng)控制研究[J].計(jì)算機(jī)仿真,2023,40(4):263-266,271.MO YP, XU D. Research on water pumping systemcontrol of induction motor based on active disturbancerejection control [J].Computer Simulation,2023,40(4):263-266,271.

    [18]張學(xué)恒,趙麒,毋天峰,等.基于磁粉離合器的電解電容卷繞張力系統(tǒng)研究[J].微處理機(jī),2021,42 (3):39-43.ZHANG XH,ZHAO Q,WUTF,et al.Researchon electrolytic capacitor winding tension system basedon magnetic powder clutch[J].Microprocessors,2021,42(3):39-43.

    [19]潘聰,葉宇,顧伯忠,等. 2.5m 大視場(chǎng)高分辨率望遠(yuǎn)鏡消光筒溫度控制[J].紅外與激光工程,2023,52 (9):183-190.PANC,YEY,GUBZ,et al.Tem-perature con-trolofextinctiontubeforthe2.5-meterlarge-field andhigh-resolution telescope [J]. Infrared and Laser Engi-neering,2023,52(9):183-190.

    [20]解英杰,尤洋,謝慕君.Smith-Fuzzy-PID在集中供熱控制系統(tǒng)中的應(yīng)用研究[J].計(jì)算機(jī)測(cè)量與控制,2014,22(9):2823-2825.XIE YJ,YOUY,XIE MJ.Applications study ofSmith-Fuzzy-PID in central heating system [J].Com-puterMeasurementamp;Control,2014,22(9):2823-2825.

    [21]魏秋月,王文慶.大滯后系統(tǒng)模糊自適應(yīng)PI-Smith控制[J].自動(dòng)化與儀表,2011,26(4):27-30.WEIQY,WANGWQ.Fuzzy self-adaptivePI-Smithcontrol for long time-delay system [J].Automationamp;Instrumentation,2011,26(4):27-30.

    [22]董偉烽,李鐘慎.加入滯后時(shí)間削弱器的大滯后系統(tǒng)的最優(yōu)模糊PID控制[J].福州大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,36(增刊1):36-38.DONG WF,LI Z S. The optimal fuzzy PID control forlargetimedelaysystemwith decreasingtimedelaycon-troller[J]. Journal of Fuzhou University(Natural Sci-enceEdition),2008,36(Sup.1):36-38.

    [23]謝濤,周邵萍,王佳碩,等.基于BP神經(jīng)網(wǎng)絡(luò)的Smith-Fuzzy-PID算法在閥門(mén)定位中的應(yīng)用研究[J].華東理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2024,50(7):5-7.XIET,ZHOUSP,WANGJS,et al.Applicationresearch of Smith -Fuzzy-PIDalgorithm based on BPneural network in valve positioning [J]. Journal of EastChina University of Science and Technology,2024,50 (7):5-7.

    第一

    作者簡(jiǎn)介:李勇,高級(jí)工程師,生于1984年,2020年畢業(yè)于中國(guó)石油勘探開(kāi)發(fā)研究院油氣井工程專(zhuān)業(yè),獲博士學(xué)位,現(xiàn)從事固井工藝技術(shù)、智能化固井研究工作。地址:(102206)北京市昌平區(qū)。電話(huà):(010)80162257。email: liyongdri@cnpc.com.cn。通信作者:程思達(dá),工程師。email:chengsddr@cnpc.com.cn。

    收稿日期:2024-06-20 修改稿收到日期:2024-10-03(本文編輯劉鋒)

    猜你喜歡
    傳遞函數(shù)預(yù)估液位
    中考命題難度控制體系的構(gòu)建
    基于頻域前饋?zhàn)赃m應(yīng)算法的車(chē)內(nèi)噪聲主動(dòng)控制仿真與試驗(yàn)分析
    基于STM32燃?xì)鉄崴仩t液位控制系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)
    石油儲(chǔ)罐液位開(kāi)關(guān)的應(yīng)用分析
    寶馬530車(chē)?yán)鋮s液液位過(guò)低報(bào)警
    智能液位控制器在排水系統(tǒng)中的應(yīng)用
    科技傳播(2012年10期)2012-06-06 09:46:22
    婷婷色综合www| 一级二级三级毛片免费看| 制服丝袜香蕉在线| 在线精品无人区一区二区三 | 男女边摸边吃奶| 三级国产精品欧美在线观看| 久久久久久久久久人人人人人人| 成人午夜精彩视频在线观看| 亚洲国产毛片av蜜桃av| 国产 一区 欧美 日韩| 熟妇人妻不卡中文字幕| 国产无遮挡羞羞视频在线观看| 成年免费大片在线观看| av在线老鸭窝| 国产亚洲精品久久久com| 亚洲色图综合在线观看| 最近中文字幕高清免费大全6| 久久综合国产亚洲精品| 深爱激情五月婷婷| 久久97久久精品| 人人妻人人添人人爽欧美一区卜 | 国语对白做爰xxxⅹ性视频网站| 亚洲高清免费不卡视频| 天天躁日日操中文字幕| 啦啦啦视频在线资源免费观看| 精品久久久久久久末码| 婷婷色综合www| av不卡在线播放| 成年人午夜在线观看视频| 麻豆成人午夜福利视频| 国产精品久久久久久精品电影小说 | 久久国产亚洲av麻豆专区| 久久97久久精品| 好男人视频免费观看在线| 18禁在线无遮挡免费观看视频| 一本久久精品| 在现免费观看毛片| 老司机影院成人| 亚洲欧美中文字幕日韩二区| 97热精品久久久久久| 一级毛片我不卡| 插阴视频在线观看视频| 夜夜骑夜夜射夜夜干| 热99国产精品久久久久久7| 免费观看的影片在线观看| 丰满少妇做爰视频| 午夜老司机福利剧场| 成年人午夜在线观看视频| 国产 精品1| 久久精品久久久久久久性| 国产精品久久久久久久久免| 欧美成人精品欧美一级黄| 如何舔出高潮| 丝袜喷水一区| 久久这里有精品视频免费| 人人妻人人看人人澡| 免费少妇av软件| 大又大粗又爽又黄少妇毛片口| 日韩欧美一区视频在线观看 | 秋霞伦理黄片| 欧美+日韩+精品| 国产精品.久久久| 欧美最新免费一区二区三区| 你懂的网址亚洲精品在线观看| 两个人的视频大全免费| 亚洲精品自拍成人| 日韩av在线免费看完整版不卡| 91久久精品电影网| 久久热精品热| 亚洲熟女精品中文字幕| 国产成人91sexporn| 国产男女内射视频| 免费av中文字幕在线| 丰满迷人的少妇在线观看| 91精品伊人久久大香线蕉| 日韩亚洲欧美综合| 51国产日韩欧美| 国产精品久久久久久久久免| 乱码一卡2卡4卡精品| 午夜免费观看性视频| 亚洲无线观看免费| 国产在线免费精品| 你懂的网址亚洲精品在线观看| 久久精品人妻少妇| 搡女人真爽免费视频火全软件| 亚洲av.av天堂| 男女边吃奶边做爰视频| 国产av国产精品国产| 夜夜骑夜夜射夜夜干| 日韩一区二区三区影片| 国产精品蜜桃在线观看| 久久久久久久亚洲中文字幕| 亚洲av中文字字幕乱码综合| 十八禁网站网址无遮挡 | 欧美高清成人免费视频www| 国产一区二区三区综合在线观看 | 亚洲av中文av极速乱| 自拍偷自拍亚洲精品老妇| 国产精品偷伦视频观看了| 亚洲精品,欧美精品| 草草在线视频免费看| 一级片'在线观看视频| 成人漫画全彩无遮挡| 91在线精品国自产拍蜜月| 一本久久精品| 亚洲av综合色区一区| videos熟女内射| av视频免费观看在线观看| 伊人久久精品亚洲午夜| www.av在线官网国产| 身体一侧抽搐| 一级毛片黄色毛片免费观看视频| 中文在线观看免费www的网站| 亚洲第一av免费看| 亚洲久久久国产精品| 一级毛片黄色毛片免费观看视频| 男的添女的下面高潮视频| 亚洲欧美精品专区久久| 亚洲美女视频黄频| 精品久久久久久久久av| 黄色怎么调成土黄色| 精品人妻熟女av久视频| 午夜精品国产一区二区电影| 免费少妇av软件| 亚洲最大成人中文| 永久网站在线| 国产亚洲5aaaaa淫片| 爱豆传媒免费全集在线观看| 99九九线精品视频在线观看视频| 久久国产亚洲av麻豆专区| 丝袜脚勾引网站| 纯流量卡能插随身wifi吗| 秋霞伦理黄片| 色网站视频免费| 亚洲欧美精品专区久久| 黑人高潮一二区| 嫩草影院新地址| 男人爽女人下面视频在线观看| 欧美 日韩 精品 国产| 日日啪夜夜爽| 亚洲无线观看免费| 两个人的视频大全免费| 日韩制服骚丝袜av| 亚洲欧洲日产国产| 91精品国产国语对白视频| 免费黄色在线免费观看| tube8黄色片| 不卡视频在线观看欧美| 日本-黄色视频高清免费观看| 街头女战士在线观看网站| 又大又黄又爽视频免费| 又粗又硬又长又爽又黄的视频| 亚洲精品久久午夜乱码| 女性被躁到高潮视频| 精品99又大又爽又粗少妇毛片| av专区在线播放| 人体艺术视频欧美日本| 婷婷色麻豆天堂久久| 亚洲国产av新网站| 我要看日韩黄色一级片| 毛片女人毛片| 中文字幕制服av| 免费看不卡的av| 99视频精品全部免费 在线| 亚洲一区二区三区欧美精品| 国产一区二区在线观看日韩| 亚洲国产欧美人成| 亚洲高清免费不卡视频| 国产无遮挡羞羞视频在线观看| 日韩视频在线欧美| 成人亚洲欧美一区二区av| 国产精品福利在线免费观看| xxx大片免费视频| 91精品国产国语对白视频| 麻豆乱淫一区二区| 久久国产精品男人的天堂亚洲 | 高清黄色对白视频在线免费看 | 精品国产乱码久久久久久小说| 国产精品精品国产色婷婷| 日韩亚洲欧美综合| 免费黄网站久久成人精品| 欧美变态另类bdsm刘玥| 亚洲电影在线观看av| 91aial.com中文字幕在线观看| 丰满乱子伦码专区| 插阴视频在线观看视频| 晚上一个人看的免费电影| 日韩成人伦理影院| 寂寞人妻少妇视频99o| 日本与韩国留学比较| 国产男人的电影天堂91| 好男人视频免费观看在线| 香蕉精品网在线| 久久国产精品男人的天堂亚洲 | 欧美精品人与动牲交sv欧美| 最近最新中文字幕免费大全7| 亚洲,一卡二卡三卡| 美女高潮的动态| 欧美人与善性xxx| 男人舔奶头视频| 精品一区在线观看国产| 最近的中文字幕免费完整| 久久久久久久亚洲中文字幕| 精品酒店卫生间| 婷婷色综合大香蕉| 人人妻人人添人人爽欧美一区卜 | 日本欧美视频一区| av天堂中文字幕网| 视频区图区小说| 天美传媒精品一区二区| 三级国产精品片| 青春草视频在线免费观看| 国产有黄有色有爽视频| 国产精品秋霞免费鲁丝片| 一二三四中文在线观看免费高清| 99热全是精品| 亚洲精品一二三| 国产爱豆传媒在线观看| 一本色道久久久久久精品综合| 天天躁日日操中文字幕| 亚洲精品456在线播放app| kizo精华| 少妇被粗大猛烈的视频| 国产成人91sexporn| 日日摸夜夜添夜夜爱| 建设人人有责人人尽责人人享有的 | 成年免费大片在线观看| 国产成人午夜福利电影在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲精品色激情综合| 毛片一级片免费看久久久久| 一级黄片播放器| 日韩欧美一区视频在线观看 | 身体一侧抽搐| 99久久中文字幕三级久久日本| 久久97久久精品| 久久精品久久久久久久性| 国产精品不卡视频一区二区| 这个男人来自地球电影免费观看 | 日日摸夜夜添夜夜爱| 热re99久久精品国产66热6| 国产伦理片在线播放av一区| 男女边摸边吃奶| 男女边吃奶边做爰视频| 国产美女午夜福利| 日韩一区二区视频免费看| 亚洲欧美精品自产自拍| 少妇人妻 视频| 日日摸夜夜添夜夜爱| 干丝袜人妻中文字幕| 国产午夜精品一二区理论片| 人妻夜夜爽99麻豆av| 亚洲一级一片aⅴ在线观看| 好男人视频免费观看在线| 内地一区二区视频在线| 亚洲中文av在线| 亚洲精品久久久久久婷婷小说| 久久久精品94久久精品| 51国产日韩欧美| 成人午夜精彩视频在线观看| 国产女主播在线喷水免费视频网站| 一级毛片 在线播放| 日韩欧美精品免费久久| 国产色爽女视频免费观看| 大陆偷拍与自拍| 日韩精品有码人妻一区| 嫩草影院新地址| 在线亚洲精品国产二区图片欧美 | 亚洲欧美日韩卡通动漫| 亚洲av综合色区一区| 日韩国内少妇激情av| 少妇人妻一区二区三区视频| 99九九线精品视频在线观看视频| 亚洲精品456在线播放app| 日本黄色日本黄色录像| 26uuu在线亚洲综合色| 国产精品久久久久久av不卡| 老熟女久久久| av天堂中文字幕网| 国产成人91sexporn| 丰满少妇做爰视频| 久久国产精品男人的天堂亚洲 | 国产精品成人在线| 久久久精品94久久精品| 最近的中文字幕免费完整| 熟女av电影| 亚洲精品日韩av片在线观看| 国产免费视频播放在线视频| videos熟女内射| 久久久久久久久久人人人人人人| 色婷婷av一区二区三区视频| 午夜福利影视在线免费观看| 在线亚洲精品国产二区图片欧美 | 王馨瑶露胸无遮挡在线观看| 精品亚洲乱码少妇综合久久| 最近中文字幕2019免费版| 午夜福利网站1000一区二区三区| 性高湖久久久久久久久免费观看| 一本一本综合久久| 性高湖久久久久久久久免费观看| 啦啦啦中文免费视频观看日本| 99视频精品全部免费 在线| 人妻夜夜爽99麻豆av| 成人国产麻豆网| 日韩欧美精品免费久久| 联通29元200g的流量卡| 欧美精品一区二区大全| 成人影院久久| 如何舔出高潮| 一个人看的www免费观看视频| 久久久久久久久大av| 99久久精品国产国产毛片| 日韩欧美一区视频在线观看 | 只有这里有精品99| 嫩草影院新地址| 欧美成人午夜免费资源| 高清午夜精品一区二区三区| 在线观看三级黄色| 一区二区三区乱码不卡18| 99热这里只有精品一区| 免费看光身美女| 视频区图区小说| 免费观看的影片在线观看| a级一级毛片免费在线观看| 只有这里有精品99| 国产综合精华液| 午夜精品国产一区二区电影| 欧美区成人在线视频| 2018国产大陆天天弄谢| 最近2019中文字幕mv第一页| 午夜视频国产福利| av免费观看日本| 六月丁香七月| 黄色欧美视频在线观看| 永久免费av网站大全| 极品少妇高潮喷水抽搐| 人人妻人人澡人人爽人人夜夜| 成人国产麻豆网| 亚洲va在线va天堂va国产| av专区在线播放| 亚洲美女黄色视频免费看| 卡戴珊不雅视频在线播放| 国产大屁股一区二区在线视频| 国产爱豆传媒在线观看| 美女xxoo啪啪120秒动态图| 久久久久久久国产电影| 亚洲电影在线观看av| 丰满乱子伦码专区| 国产一级毛片在线| 99视频精品全部免费 在线| 黑人猛操日本美女一级片| 国产无遮挡羞羞视频在线观看| 黄色日韩在线| 一边亲一边摸免费视频| 欧美成人精品欧美一级黄| 欧美区成人在线视频| 一级毛片电影观看| 国产精品久久久久久久久免| 欧美日韩一区二区视频在线观看视频在线| 2018国产大陆天天弄谢| 99热这里只有是精品在线观看| 亚洲av二区三区四区| 欧美亚洲 丝袜 人妻 在线| 国产精品人妻久久久影院| 啦啦啦视频在线资源免费观看| 国产黄色视频一区二区在线观看| 极品少妇高潮喷水抽搐| 不卡视频在线观看欧美| 久久热精品热| av又黄又爽大尺度在线免费看| 高清不卡的av网站| 99热全是精品| 联通29元200g的流量卡| 成人亚洲欧美一区二区av| 最新中文字幕久久久久| 舔av片在线| 精品一区二区三区视频在线| 欧美高清性xxxxhd video| 欧美+日韩+精品| 午夜福利在线在线| av又黄又爽大尺度在线免费看| 午夜免费鲁丝| 亚洲欧美一区二区三区国产| 久久久久久人妻| 尾随美女入室| 亚洲一区二区三区欧美精品| 我的女老师完整版在线观看| 国产免费视频播放在线视频| 美女高潮的动态| 蜜桃在线观看..| 日韩欧美 国产精品| 在线观看一区二区三区| 大片免费播放器 马上看| 男女无遮挡免费网站观看| 亚洲成人中文字幕在线播放| 大香蕉97超碰在线| 亚洲无线观看免费| 国产精品一及| 少妇熟女欧美另类| 成人亚洲精品一区在线观看 | 91精品伊人久久大香线蕉| 国产高清三级在线| 久久99精品国语久久久| 好男人视频免费观看在线| 97超碰精品成人国产| 亚洲美女视频黄频| 久久久久国产精品人妻一区二区| 欧美激情极品国产一区二区三区 | 国产精品精品国产色婷婷| 亚洲精品一区蜜桃| 亚洲av中文av极速乱| 好男人视频免费观看在线| 国产成人a区在线观看| 亚洲美女黄色视频免费看| 国产亚洲精品久久久com| 一个人看视频在线观看www免费| 97超视频在线观看视频| 美女视频免费永久观看网站| 国产熟女欧美一区二区| 激情 狠狠 欧美| 久久久久久久久久久丰满| 成人毛片a级毛片在线播放| 美女内射精品一级片tv| 蜜桃在线观看..| 中国三级夫妇交换| 亚洲精品456在线播放app| 寂寞人妻少妇视频99o| 午夜福利视频精品| 国产精品人妻久久久久久| 久久久久视频综合| 成人二区视频| av.在线天堂| 视频中文字幕在线观看| 大又大粗又爽又黄少妇毛片口| 精华霜和精华液先用哪个| 久久久久精品性色| 三级经典国产精品| 97在线视频观看| 中文在线观看免费www的网站| 国产一区二区三区av在线| 国产人妻一区二区三区在| 精品国产一区二区三区久久久樱花 | 精品一区二区免费观看| 热99国产精品久久久久久7| 亚洲中文av在线| 少妇的逼好多水| 91精品国产九色| 亚洲一区二区三区欧美精品| 国产成人免费无遮挡视频| 成年女人在线观看亚洲视频| 一区二区av电影网| 国产精品久久久久久精品古装| 午夜福利网站1000一区二区三区| 嫩草影院入口| 网址你懂的国产日韩在线| 黄片无遮挡物在线观看| 亚洲国产最新在线播放| 伦精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 成人国产av品久久久| 午夜激情久久久久久久| 欧美极品一区二区三区四区| 男人舔奶头视频| 一级毛片我不卡| 十八禁网站网址无遮挡 | 国产午夜精品久久久久久一区二区三区| 欧美性感艳星| 欧美日韩亚洲高清精品| 亚洲成色77777| av免费在线看不卡| 在线 av 中文字幕| 国产色婷婷99| 亚洲av国产av综合av卡| 校园人妻丝袜中文字幕| 高清毛片免费看| 寂寞人妻少妇视频99o| 国产精品爽爽va在线观看网站| 色视频www国产| 性高湖久久久久久久久免费观看| 亚洲av男天堂| 国产片特级美女逼逼视频| 欧美高清成人免费视频www| 哪个播放器可以免费观看大片| 少妇猛男粗大的猛烈进出视频| 涩涩av久久男人的天堂| 国产免费一级a男人的天堂| 18禁在线播放成人免费| av视频免费观看在线观看| 国产成人精品婷婷| 国产精品麻豆人妻色哟哟久久| 欧美bdsm另类| 性色av一级| 亚洲一区二区三区欧美精品| 在线 av 中文字幕| 夜夜看夜夜爽夜夜摸| 丰满少妇做爰视频| 久久久久久久久久久免费av| 欧美性感艳星| 国产又色又爽无遮挡免| 一级a做视频免费观看| 中文字幕久久专区| 一本—道久久a久久精品蜜桃钙片| 精品99又大又爽又粗少妇毛片| 国产精品福利在线免费观看| 亚洲欧美日韩另类电影网站 | 精品久久久久久久久av| 免费少妇av软件| 97超视频在线观看视频| 最近最新中文字幕大全电影3| www.色视频.com| 久久久久久久久久久丰满| 在线观看一区二区三区| 久久久久性生活片| 日韩av在线免费看完整版不卡| 亚洲精品日本国产第一区| 亚洲欧美精品自产自拍| 国模一区二区三区四区视频| 久久久久久久久久久免费av| 亚洲色图综合在线观看| 高清视频免费观看一区二区| 一级二级三级毛片免费看| 尾随美女入室| 中国国产av一级| 欧美日本视频| 精品一区在线观看国产| 精品一品国产午夜福利视频| 一级毛片 在线播放| 又爽又黄a免费视频| 美女国产视频在线观看| 国产在线一区二区三区精| 少妇丰满av| 女的被弄到高潮叫床怎么办| 日本av免费视频播放| 久热久热在线精品观看| 欧美日韩精品成人综合77777| 亚洲精品中文字幕在线视频 | 看免费成人av毛片| 精品少妇久久久久久888优播| 久久精品国产亚洲av涩爱| 一级毛片久久久久久久久女| 亚洲,欧美,日韩| 精品少妇久久久久久888优播| 99re6热这里在线精品视频| 欧美性感艳星| 91狼人影院| 人人妻人人添人人爽欧美一区卜 | 精品久久久久久久久亚洲| 国产欧美日韩一区二区三区在线 | 日韩电影二区| 久久人人爽人人爽人人片va| 精品一品国产午夜福利视频| 多毛熟女@视频| 久久毛片免费看一区二区三区| a级毛片免费高清观看在线播放| 国产又色又爽无遮挡免| 狂野欧美激情性xxxx在线观看| 97超碰精品成人国产| 丰满迷人的少妇在线观看| 国产毛片在线视频| 久久青草综合色| 国产精品一区二区三区四区免费观看| 卡戴珊不雅视频在线播放| 黄色日韩在线| 日本一二三区视频观看| 永久网站在线| 乱系列少妇在线播放| 亚洲真实伦在线观看| 性色av一级| 免费观看a级毛片全部| 欧美成人午夜免费资源| 久久久久久久久大av| 免费黄频网站在线观看国产| 又粗又硬又长又爽又黄的视频| 80岁老熟妇乱子伦牲交| 成人影院久久| 亚洲久久久国产精品| 日日摸夜夜添夜夜爱| 免费av不卡在线播放| 国产精品久久久久成人av| 国产爱豆传媒在线观看| 午夜激情福利司机影院| 亚洲激情五月婷婷啪啪| 欧美丝袜亚洲另类| 久久国产乱子免费精品| 精品99又大又爽又粗少妇毛片| 亚洲伊人久久精品综合| 丝袜喷水一区| 国产精品三级大全| 18+在线观看网站| 最近中文字幕高清免费大全6| 女性生殖器流出的白浆| 精品久久久精品久久久| 身体一侧抽搐| 欧美xxxx性猛交bbbb| 免费人妻精品一区二区三区视频| 久久久久久久大尺度免费视频| 亚洲av.av天堂| 国产淫语在线视频| 青春草亚洲视频在线观看| 日韩国内少妇激情av| 99久久综合免费| 成人毛片60女人毛片免费| 国产视频内射| 最新中文字幕久久久久| 最近最新中文字幕大全电影3| 夫妻午夜视频| 久久精品熟女亚洲av麻豆精品| 欧美国产精品一级二级三级 | 一个人看视频在线观看www免费| 男男h啪啪无遮挡| 成年美女黄网站色视频大全免费 | 精品久久久久久久久亚洲| av在线app专区| 九色成人免费人妻av| 欧美人与善性xxx|