亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        新疆呼圖壁地下儲氣庫的InSAR形變監(jiān)測與模擬

        2016-09-21 02:10:23余鵬飛馮光財(cái)喬學(xué)軍
        大地測量與地球動力學(xué) 2016年9期
        關(guān)鍵詞:呼圖壁儲氣庫氣井

        陳 威 余鵬飛 熊 維 李 杰 馮光財(cái) 喬學(xué)軍

        1 中國地震局地震研究所(地震大地測量重點(diǎn)實(shí)驗(yàn)室),武漢市洪山側(cè)路40號,430071 2 新疆維吾爾自治區(qū)地震局,烏魯木齊市北京南路42號,830011 3 中南大學(xué)地球科學(xué)與信息物理學(xué)院,長沙市麓山南路932號,410083

        ?

        新疆呼圖壁地下儲氣庫的InSAR形變監(jiān)測與模擬

        陳威1余鵬飛1熊維1李杰2馮光財(cái)3喬學(xué)軍1

        1中國地震局地震研究所(地震大地測量重點(diǎn)實(shí)驗(yàn)室),武漢市洪山側(cè)路40號,430071 2新疆維吾爾自治區(qū)地震局,烏魯木齊市北京南路42號,830011 3中南大學(xué)地球科學(xué)與信息物理學(xué)院,長沙市麓山南路932號,410083

        利用TerraSAR-X衛(wèi)星2013-08~2014-08的17景雷達(dá)影像,采用小基線集(small baseline subset,SBAS)InSAR技術(shù)獲取呼圖壁(HTB)地下儲氣庫(underground gas storage,UGS)運(yùn)行期間的地表形變序列,并結(jié)合UGS注(采)氣井口的壓力數(shù)據(jù),采用多點(diǎn)源Mogi模型,對HTB UGS的形變場進(jìn)行模擬。結(jié)果表明,整個UGS區(qū)域的形變特征為非連續(xù)分布,形變與注(采)氣壓力變化具有較好的相關(guān)性;注(采)氣期間沿衛(wèi)星視線向(LOS)的形變峰值分別為10 mm和-8 mm;采用自適應(yīng)前向搜索法,基于多點(diǎn)源Mogi模型初步模擬注(采)氣期間的形變過程,當(dāng)UGS的注(采)氣平均氣壓為18 MPa和15 MPa時(shí),LOS的形變可達(dá)7 mm和-4 mm,地表形變的大小與注(采)氣井口密度有關(guān);UGS的儲氣分布呈非均勻狀態(tài),即地下氣庫結(jié)構(gòu)復(fù)雜多變。

        地下儲氣庫;呼圖壁;小基線集InSAR;Mogi模型;地殼形變

        UGS在注(采)氣期間的應(yīng)力交變不但會引起儲氣庫蓋層地表周期性的呼吸狀變形, 而且可能引起儲層巖體裂紋擴(kuò)張,造成斷層活化并誘發(fā)地震[1]。向地殼深部地層注入流體的過程中,誘發(fā)地震的頻次及最大震級一般與注入量、注入速率、地層滲透率及地層應(yīng)力大小正相關(guān),且大部分地震發(fā)生在注入過程中,并集中在注入井井底附近的一定范圍內(nèi)[2]。因此,加強(qiáng)UGS的形變監(jiān)測,將有利于深入了解儲氣庫內(nèi)部的構(gòu)造運(yùn)動特征、儲層巖石應(yīng)力應(yīng)變狀態(tài)及地震活動性隨注(采)氣壓力變化的規(guī)律等,對確保儲氣庫的穩(wěn)定、高效運(yùn)行具有重要的理論和實(shí)用意義。

        HTB UGS屬枯竭氣田儲氣庫,是中國目前最大的地下儲氣庫,位于準(zhǔn)噶爾盆地南緣(圖1),屬北天山山前坳陷,蓋層為第四紀(jì)超厚沉積物,厚度可達(dá)3.1~3.7 km。該區(qū)經(jīng)歷了多期構(gòu)造運(yùn)動,特別是喜山期,受北天山強(qiáng)烈活動影響,山前區(qū)強(qiáng)烈褶皺并伴生一系列大型逆掩斷裂,造成深淺層構(gòu)造差異很大[3]。

        虛線粉色框內(nèi)為地下儲氣庫范圍,紅色正方形為呼圖壁注采氣井位,黃色圓點(diǎn)為2013-06~08間的地震分布,藍(lán)色四邊形為SAR圖像覆蓋區(qū)域圖1 呼圖壁儲氣庫DEM地形、斷層、井口分布、地震活動及SAR影像概況Fig.1 The general situation of Hutubi gas storage DEM topography,wells,faults,seismic activities distribution and SAR images

        2013-06 HTB UGS投入運(yùn)行。兩個月內(nèi),其周邊發(fā)生多次小地震(圖1),打破了該地區(qū)地震的平靜期。本文基于SBAS-InSAR技術(shù),利用TerraSAR-X衛(wèi)星2013-08~2014-08的觀測數(shù)據(jù),對HTB UGS及周邊區(qū)域開展地表形變的時(shí)序分析。在獲取LOS形變場的基礎(chǔ)上,結(jié)合UGS運(yùn)行期間注(采)氣井口的壓力變化(表1),利用多點(diǎn)源的Mogi模型對注(采)氣期間的地表形變進(jìn)行相應(yīng)的數(shù)值模擬,揭示儲氣庫的形變與注(采)氣壓力變化的關(guān)系。

        表1 HTB UGS運(yùn)行期的注(采)氣時(shí)間及井口壓力統(tǒng)計(jì)

        1 SBAS-InSAR 技術(shù)

        假定在時(shí)間(t1,t2,…,tN)之間獲取N+1幅SAR影像,且所有影像均在同一坐標(biāo)系下。根據(jù)設(shè)定的時(shí)空基線閾值生成M幅干涉圖,其中M應(yīng)滿足:

        (1)

        對滿足要求的時(shí)空基線閾值的干涉圖進(jìn)行濾波及解纏處理,去除平地相位和地形相位影響,根據(jù)影像的相干信息選取高相干點(diǎn),得到高相干點(diǎn)的相位為:

        (2)

        為獲取有物理意義的形變時(shí)序,可用SAR影像兩個獲取時(shí)間之間的平均速度來表示式(2)中的相位:

        (3)

        則式(2)可表示為:

        (4)

        將式(4)寫成矩陣形式:

        (5)

        式(5)是M×N的矩陣。應(yīng)用奇異值分解法(singular value decomposition,SVD)可獲平均速度時(shí)間序列,最后通過各個時(shí)間段內(nèi)速度的積分可獲取形變時(shí)間序列。

        圖2 SBAS空間及時(shí)間基線分布Fig.2 Temporal and spatial baseline distribution of SBAS-InSAR pairs

        2  SAR數(shù)據(jù)選取、處理及時(shí)序分析

        2.1SAR數(shù)據(jù)選取與處理

        選取2013-08-29~2014-08-27期間的17景TerraSAR-X衛(wèi)星降軌數(shù)據(jù),影像數(shù)據(jù)的具體參數(shù)見表2。利用ISCE(InSAR scientific computing environment)軟件處理生成48景垂直基線小于400 m、時(shí)間間隔小于3個月的多視差分干涉圖。在差分干涉處理過程中,利用3″的SRTM數(shù)字高程模型去除地形相位,并用多項(xiàng)式對軌道參數(shù)進(jìn)行擬合以減少軌道誤差。為抑制相位噪聲,對距離向和方位向分別作4視和12視的多視處理,并采用加權(quán)功率譜法對生成的干涉圖進(jìn)行濾波,最后利用SNAPHU程序進(jìn)行相位解纏。

        表2 TerraSAR-X影像基本參數(shù)

        2.2InSAR時(shí)序計(jì)算及分析

        生成干涉圖后,進(jìn)行SBAS計(jì)算和時(shí)序分析,對大氣延遲引起的干涉相位進(jìn)行改正。采用Doin等[4]的方法,結(jié)合歐洲中期預(yù)報(bào)中心ECMWF(European Centre for Medium-Range Weather Forecasts)發(fā)布的氣象數(shù)據(jù),對HTB UGS地區(qū)的大氣延遲誤差進(jìn)行改正。采用網(wǎng)絡(luò)Deramping方法進(jìn)行軌道誤差估計(jì)。首先利用最小二乘法對每幅干涉圖單獨(dú)進(jìn)行軌道誤差估算;然后通過線性反演重新估算整個干涉圖像網(wǎng)絡(luò)的軌道參數(shù),以確保一致性;最后將重新估算的軌道參數(shù)合并,生成與干涉圖像網(wǎng)絡(luò)一致的軌道改正及每幅干涉圖的軌道改正[5]。

        為減小因相干性較低引起的誤差,選取相干系數(shù)大于0.5的點(diǎn)進(jìn)行時(shí)序計(jì)算。以2013-08-29為參考時(shí)間,最終獲取了其他圖像2013-10-01~2014-08-27的HTB UGS地表LOS形變時(shí)間序列(圖3,黑色虛線為UGS,五角星與圓點(diǎn)是兩個井口位置)??梢钥闯?,HTB UGS地表形變的隆升或者下沉?xí)r間與儲氣庫注(采)氣周期基本一致,在1~2月份儲氣庫的采氣運(yùn)行期間,由采氣過程引起的地表下沉量可達(dá)-8 mm;在夏、秋季儲氣庫注氣運(yùn)行期間,由注氣引起的地表隆升可達(dá)到10 mm。

        提取儲氣庫氣井HUK9和HUK19附近點(diǎn)位的形變時(shí)序(圖4),并與井壓的時(shí)間序列進(jìn)行比對。結(jié)果表明,UGS的地表LOS形變與井壓變化正相關(guān)。在UGS運(yùn)行期間,從6月左右開始,隨著井壓增大,地表開始隆升,并隨氣壓的增大達(dá)到極值;從冬季開始,隨著井壓減小,地表下沉,并在2月左右達(dá)到極值。

        圖3 HTB UGS形變時(shí)間序列圖Fig.3 The deformation time series of HTB UGS

        圖4 所選點(diǎn)位的形變與井壓時(shí)間序列圖Fig.4 Deformation time series of the chosen ground points

        3 形變數(shù)值模擬

        HTB UGS的注(采)氣主要通過26個氣井進(jìn)行控制,因此,將儲氣庫的注(采)氣井作為壓力點(diǎn)源,注(采)氣井口壓力數(shù)據(jù)作為輸入,采用多點(diǎn)源的Mogi模型(表3)和自適應(yīng)前向搜索法[6],對HTB UGS地區(qū)的地表形變進(jìn)行模擬(圖5,圓大小代表半徑,顏色代表埋深)。模擬獲取的HTB UGS地區(qū)的地表形變及殘差結(jié)果見圖6。結(jié)果顯示,地表形變的大小與注(采)氣井口密度有關(guān),在儲氣庫注氣期間,儲氣庫地表形變最大可達(dá)到7 mm;在采氣期間,儲氣庫地表形變最大可達(dá)到-4 mm。

        圖6顯示,儲氣庫注(采)氣井密度較大地區(qū)的殘差較小,范圍在-2~1 mm,其他地區(qū)誤差值較大,范圍在-4~2 mm,表明模擬研究雖然提取了HTB UGS地區(qū)由注(采)氣引起的形變,但是一些由局部形變和噪聲影響的點(diǎn)無法模擬。

        圖5 Mogi模型空間分布Fig.5 The spatial distribution of Mogi model

        井號埋深/km半徑/m井號埋深/km半徑/mHUK13.1125HUK143.3150HUK23.2150HUK153.4125HUK33.1175HUK163.4175HUK43.0150HUK173.5175HUK53.1125HUK183.6175HUK63.1175HUK193.8150HUK73.2125HUK203.6125HUK83.3175HUK213.5175HUK93.5175HUK223.3150HUK103.4150HUK233.5125HUK113.2175HUK243.6125HUK123.3150HUK253.4150HUK133.3175HUK263.2150

        圖6 Mogi模型模擬結(jié)果與殘差(紅色正方形為注(采)氣井位)Fig.6 Residual result and simulation surface displacement based on Mogi model (red squares are gas ingection wells)

        4 討 論

        在SBAS-InSAR數(shù)據(jù)處理中,考慮軌道誤差、DEM誤差和大氣延遲誤差等因素對HTB UGS地表形變的影響。首先采用多項(xiàng)式模型在SAR影像干涉處理過程中初步改正軌道誤差,并在SBAS時(shí)序分析過程中利用網(wǎng)絡(luò)Deramping方法消除殘余的軌道誤差。DEM誤差對地表形變量的影響與空間基線、衛(wèi)星入射角和斜距有關(guān),當(dāng)空間基線小于400 m時(shí),利用TerraSAR衛(wèi)星數(shù)據(jù)和SRTM3 DEM 可以獲得的精度優(yōu)于2 mm。Emardson等[7]通過Stacking方法處理長時(shí)間序列的SAR數(shù)據(jù),可以獲取1 mm/a的形變。Doin等[8]在Emardson等的基礎(chǔ)上定量給出了大氣對流層延遲誤差的影響,獲得改正后的大氣延遲誤差殘差精度優(yōu)于2 mm/km。本文利用Doin等的方法,結(jié)合ECMWF的氣象數(shù)據(jù),對HTB UGS地區(qū)的大氣延遲誤差進(jìn)行改正。

        InSAR時(shí)序結(jié)果表明,UGS運(yùn)行初期,雖然氣壓與形變存在正相關(guān),但形變量較小,這可能是因?yàn)樗ソ邭馓锔慕ǖ腢GS在第一個注采氣周期處于殘余油氣、水汽和注入氣體等混合動態(tài)平衡過程,而這種混合程度影響儲氣庫產(chǎn)水量和體積變化速率。同時(shí),在UGS運(yùn)行初期,其蓋層巖性(如孔隙度、滲透率、抗壓強(qiáng)度、抗拉張強(qiáng)度等)均處于良好狀態(tài),不易造成應(yīng)力積累和體積變化[9]。因此,相同的氣壓條件下,在UGS投入運(yùn)行的初期,形變量較小(圖4)。整個UGS的形變幅度從注(采)氣井密集的儲氣庫中心區(qū)域向四周逐漸變小,且儲氣庫周邊區(qū)域大部分表現(xiàn)為下沉,尤其是儲氣庫東側(cè)區(qū)域下沉幅度較大。這可能是由于HTB氣田開發(fā)鉆井大部分集中在此區(qū)域,且該區(qū)域巖性以粉砂巖、粉砂泥質(zhì)巖和泥質(zhì)粉砂巖為主,儲氣層層間非均勻性強(qiáng),地層富含的敏感礦物容易與水、酸等流體產(chǎn)生化學(xué)反應(yīng)生成沉積物[10]。 另外,HTB UGS地處天山北麓沖洪積平原中下游農(nóng)灌區(qū),屬于地下水嚴(yán)重超采區(qū),位于地下水“降落漏斗”地帶,其漏斗中心水位下降6~23.25 m[11]。

        Mogi模型獲取的注(采)氣井埋深范圍為3.1~3.7 km,與地球物理方法測井技術(shù)所得井口埋深范圍一致,模擬的儲氣庫總?cè)萘繛?5.35億m3,與生產(chǎn)容量45.6億m3相符[12]。Mogi模型是基于理想彈性半無限空間的結(jié)果,而HTB UGS地區(qū)地質(zhì)構(gòu)造復(fù)雜,多為層狀不均勻介質(zhì),模擬結(jié)果必然存在一定誤差。形變殘差圖(圖6)表明,儲氣庫埋深較大的注(采)氣井集中區(qū)域的形變模擬具有較高精度,而儲氣庫外圍區(qū)域的殘差較大,驗(yàn)證了Mogi模型適用于埋深比較大的腔體區(qū)域[13],也說明了基于簡單假設(shè)的Mogi模型難以模擬整個區(qū)域的形變,因此,需要建立以地震地質(zhì)、地球物理及儲氣井的壓力與氣量等參數(shù)為約束的地質(zhì)力學(xué)模型。

        5 結(jié) 語

        本文利用2013-08~2014-08觀測的17景TerraSAR-X衛(wèi)星雷達(dá)影像,采用SBAS-InSAR技術(shù)獲取了HTB UGS的地表LOS形變序列。在UGS運(yùn)行初期,由于儲氣庫蓋層巖性(孔隙度、滲透率、抗壓強(qiáng)度、抗拉張強(qiáng)度等)尚處于“靜穩(wěn)”狀態(tài),不易造成應(yīng)力積累和體積變化,隨著UGS注(采)氣的正常運(yùn)行,地表形變與UGS注(采)氣氣壓密切相關(guān)。UGS的儲氣分布呈非均勻狀態(tài),即地下氣庫結(jié)構(gòu)復(fù)雜多變,需要進(jìn)一步融入GPS、精密水準(zhǔn)和重力測量等結(jié)果,并建立精細(xì)的地質(zhì)力學(xué)模型,以更好地揭示UGS內(nèi)部結(jié)構(gòu)并研究其形變機(jī)理。

        [1]Shapiro S A, Dinske C, Kummerow J. Probability of a Given-Magnitude Earthquake Induced by a Fluid Injection[J]. Geophys Res Lett, 2007,34(L22 314)

        [2]McGarr A, Spottiswoode S M, Ortlepp W D. Relationship of Mine Tremors to Induced Stresses and to Rock Properties in Focal Region[J].Bull Seism Soc Am, 1979,65(9):81-93

        [3]陳立春. 北天山烏魯木齊轉(zhuǎn)換區(qū)構(gòu)造系晚第四紀(jì)活動性[D].北京:中國地震局地質(zhì)研究所, 2011(Chen Lichun. Late Quaternary Behavior of the Active Tectonic System in the Urumqi Transform Region of the North Tianshan[D].Beijing: Institute of Geology, CEA, 2011)

        [4]Doin M P, Lasserre C, Peltzer G, et al. Corrections of Stratified Tropospheric Delays in SAR Interferometry: Validation with Global Atmospheric Models[J]. J App Geophys, 2009,69(1):35-50

        [5]Biggs J,Wright T, Lu Z, et al. Multi-Interferogram Method for Measuring Interseismic Deformation: Denali Fault, Alaska[J].Geophysical Journal International, 2007, 170(3): 1 165-1 179

        [6]周暉, 徐晨, 邵世煌,等.自適應(yīng)搜索優(yōu)化算法[J].計(jì)算機(jī)科學(xué), 2008, 35(10): 188-191(Zhou Hui, Xu Chen, Shao Shihuang, et al. Adaptive Free Search Algorithm[J].Computer Science, 2008, 35(10): 188-191 )

        [7]Emardson T R, Simons M, Webb F H. Neutral Atmospheric Delay in Interferometric Synthetic Aperture Radar Applications: Statistical Description and Mitigation [J]. J Geophys Res, 2003, 108(B5)

        [8]Doin M P, Grandin R, Lasserre C, et al. Demcor-Rections before Unwrapping in a Small Baseline Strategy for InSAR Time Series Analysis[J]. Geoscience and Remote Sensing, 2011, 69(10): 1 353-1 356

        [9]Jha B, Bottazzi F, Wojcik R, el at. Reservoir Characterization in an Underground Gas Storage Field Using Joint Inversion of Flow and Geodetic Data[J]. Int J Numer Anal Meth Geomech, 2015, 39(14): 1 619-1 638

        [10] 張明, 趙金省. 呼圖壁氣田儲層地質(zhì)特征分析及潛在傷害因素研究[J].石油天然氣學(xué)報(bào), 2008, 30(5): 189-193(Zhang Ming, Zhao Jinsheng. Analysis of Geological Characteristics and Potential Damage Factors of Reservoir in Hutubi Gas Field[J]. Journal of Oil and Gas Technology, 2008, 30(5): 189-193)

        [11] 楊宏偉, 朱瑾. 準(zhǔn)噶爾盆地南緣平原區(qū)地下水動態(tài)及地下水資源可持續(xù)利用對策[J].新疆地質(zhì), 2012,30(3):350-354(Yang Hongwei, Zhu Jin. The Groundwater Dynamic and Sustainable Utilization Countermeasures of Groundwater Resources in Southern Plain of Junggar Basin in Recent 30 Years[J]. Xinjiang Geology, 2012,30(3):350-354)

        [12] 曹錫秋. 新疆某地衰竭氣藏地下儲氣庫地應(yīng)力特征研究[D].北京:中國地質(zhì)大學(xué)(北京), 2013(Cao Xiqiu. Study on In-Situ Stress Characteristics of Underground Gas Storage in a Depleted Gas Reservoir in Xinjiang[D]. Beijing: China University of Geosciences(Beijing), 2013)

        [13] 胡亞軒, 施行覺, 王慶良,等. 騰沖火山區(qū)地表垂直形變分析[J]. 大地測量與地球動力學(xué), 2003, 23(2): 37-41(Hu Yaxuan, Shi Xingjue, Wang Qingliang, et al. Analysis of Vertical Deformation in Tengchong Volvanic Area [J]. Journal of Geodesy and Geodynamics, 2003, 23(2): 37-41)

        Foundation support:National Natural Science Foundation of China,No.41474097, 41274027,41474016,41574005; Director Fund of Institute of Seismology, CEA, No.IS201326129,IS201116013.

        About the first author:CHEN Wei, postgraduate, majors in InSAR deformation monitoring and geodynamics, E-mail: chenweicug@126.com.Corresponding author:QIAO Xuejun, PhD, researcher, majors in InSAR and GNSS deformation monitoring and geodynamics, E-mail: xuejunq@sohu.com.

        InSAR Deformation Monitoring and Simulation of Underground Gas Storage in Hutubi, Xinjiang

        CHENWei1YUPengfei1XIONGWei1LIJie2FENGGuangcai3QiaoXuejun1

        1Key Laboratory of Earthquake Geodesy, Institute of Seismology, CEA, 40 Hongshance Road, Wuhan 430071, China 2Earthquake Administration of Xinjiang Uygur Autonomous Region, 42 South-Beijing Road, Urumqi 830011, China 3School of Geoscience and Info-Physics, Central South University, 932 South-Lushan Road, Changsha 410083, China

        In this paper, 17 scenes TerraSAR-X radar images acquired from August 2013 to August 2014 are used by the small baseline subset (SBAS) InSAR method to obtain the surface deformation series

        during the operation of underground gas storage (UGS) in Hutubi (HTB). We combine these data with data from injection/production pressure wells, using multi-point source Mogi model to simulate the deformation field of UGS in Hutubi. The results show that the deformation characteristics of the whole UGS area is a discontinuous distribution, and the deformation peak value along satellite line of sight (LOS) is 10 mm and -8 mm during gas injection and gas production respectively. The retrieved deformation sequences conform to injection/production pressure changes very well. Based on the multi-point source Mogi model, we simulate the deformation process of HTB UGS and apply an adaptive forward search method to obtain the radius and depth of point source. The simulated results indicate that when the injection/production average pressure of HTB UGS is 18 MPa and 15 MPa, LOS deformation is up to 7 mm and -4 mm respectively; surface deformation is related to the density of gas injection/production wells. The UGS gas distribution is not uniform, indicating that the structure of underground gas storage is complex.

        underground gas storage (UGS); Hutubi; small baseline subset (SBAS) InSAR; Mogi model; crustal deformation

        2016-04-29

        喬學(xué)軍,博士,研究員,主要研究方向?yàn)镚NSS與InSAR形變監(jiān)測與地球動力學(xué),E-mail:xuejunq@sohu.com。

        10.14075/j.jgg.2016.09.011

        1671-5942(2016)09-0803-05

        P237

        A

        項(xiàng)目來源:國家自然科學(xué)基金(41474097,41274027,41474016,41574005);中國地震局地震研究所所長基金(IS201326129,IS201116013)。第一作者簡介:陳威,碩士生,主要研究方向?yàn)镮nSAR形變監(jiān)測與地球動力學(xué),E-mail:chenweicug@126.com。

        猜你喜歡
        呼圖壁儲氣庫氣井
        全國最大壓縮機(jī)在呼圖壁儲氣庫成功安裝
        港華鹽穴儲氣庫的運(yùn)營特點(diǎn)及其工藝改進(jìn)
        煤氣與熱力(2022年4期)2022-05-23 12:44:52
        氣井用水合物自生熱解堵劑解堵效果數(shù)值模擬
        基于鉆孔應(yīng)變觀測約束的2016年新疆呼圖壁M6.2地震的發(fā)震斷層研究
        地震研究(2021年1期)2021-04-13 01:05:08
        呼圖壁氣槍數(shù)據(jù)中面波振幅變化研究
        地震研究(2021年1期)2021-04-13 01:04:52
        “花兒昌吉、美麗呼圖壁”環(huán)境保護(hù)計(jì)劃分析
        基于STM32F207的便攜式氣井出砂監(jiān)測儀設(shè)計(jì)
        電子制作(2018年9期)2018-08-04 03:31:16
        氣井出砂動態(tài)監(jiān)測技術(shù)研究
        電子制作(2018年2期)2018-04-18 07:13:33
        氣井節(jié)點(diǎn)分析技術(shù)在橋口井區(qū)的應(yīng)用
        鹽穴儲氣庫注采集輸系統(tǒng)優(yōu)化
        寂寞人妻渴望被中出中文字幕| 91精品国产91| 国产亚洲精品成人av在线| 久久久亚洲免费视频网| 国产精品国产三级国产av品爱网| 免费观看黄网站在线播放| 中文字幕久久久久久精| 久久精品成人欧美大片| 亚洲网站免费看| 亚洲免费一区二区av| 18禁在线永久免费观看| 亚洲欧美在线观看| 粗大挺进尤物人妻一区二区| 国产视频一区二区三区久久亚洲| 国产午夜av秒播在线观看| 国产一线二线三线女| 久久久久久久久久免免费精品| 日本频道一区二区三区| 无码爆乳护士让我爽| 国内精品久久久久久久影视麻豆| 无码AV无码免费一区二区| 亚洲一区二区三区免费网站| 中文字幕在线日亚洲9| 亚洲白白色无码在线观看| 久久精品这里就是精品| 日本激情网站中文字幕 | 欧美 变态 另类 人妖| 99久久国产综合精品女乱人伦| 亚洲av本道一本二本三区| 久久亚洲精品成人av无码网站 | 日韩丝袜人妻中文字幕| 亚洲人成综合第一网站| 另类内射国产在线| 欧美精品久久久久久三级| 久久99精品综合国产女同| 欧美另类人妖| 国产成人拍精品免费视频| 国产自拍精品在线视频| 国产精品区一区二区三在线播放| 天堂网在线最新版www中文网| 国产精品自在在线午夜出白浆|