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

    天津市腦卒中高住院風險的氣象原因探析*

    2022-07-08 03:06:26郝天依韓素芹姚青王曉佳樊文雁辛江龍
    氣象學報 2022年3期
    關鍵詞:變溫氣壓天津市

    郝天依 韓素芹 姚青 王曉佳 樊文雁 辛江龍

    1.天津市環(huán)境氣象中心,天津,300074

    2.中國氣象局-南開大學大氣環(huán)境與健康研究聯(lián)合實驗室,天津,300074

    3.天津市健康氣象交叉創(chuàng)新中心,天津,300074

    4.天津市健康醫(yī)療大數(shù)據(jù)有限公司,天津,301700

    1 引言

    腦卒中是由急性腦循環(huán)障礙引起腦功能損傷的腦血管病,已成為中國第一位的致死及致殘疾病(《中國卒中中心報告2019》編寫組,2021;《中國腦卒中防治報告》編寫組,2020)。對于已存在發(fā)病基礎且內(nèi)環(huán)境處于易發(fā)病狀態(tài)的個體,外界環(huán)境的劇烈變化可實時誘發(fā)、啟動腦卒中發(fā)生(Wang,et a1,2013;黃如訓,2005)。所以通過研究氣象因素的健康效應,減少可干預危險因素對腦卒中發(fā)病的不利影響,對做好腦卒中的預防工作具有重大意義。

    氣溫、氣壓和濕度與腦卒中發(fā)病均存在較大關聯(lián)。Takahashi 等(1957)便觀察到腦出血的死亡率與地理分布有關,并通過實驗證實了其相關是由氣溫的差異造成的。中國腦卒中發(fā)病和死亡也存在明顯的地域差別,腦卒中發(fā)病率和死亡率最高的是東北地區(qū),中部地區(qū)次之,最低的是南部地區(qū)(Wang,et a1,2017b);東北、華北、西南的西藏和貴州地區(qū)發(fā)病率明顯高于其他地區(qū),東南沿海地區(qū)相對較低(Zhou,et a1,2016)。溫度過低或過高均會增加腦卒中發(fā)病風險,但各個地區(qū)存在不同的最適溫度(Gasparrini,et a1,2015)。中國北方城市的最適日均溫度為20.1—24.4℃,南方城市則為25.2—30.8℃(Yang,et a1,2016)。氣溫對機體的影響機制研究表明,一旦遇到強冷、暖氣流的刺激易引起血流動力學改變,促成血壓波動,是腦卒中風險增高的主要因素之一(Wang,et a1,2017a;Zhao,et a1,2019)。其次,氣溫驟變導致心率異常和血液動力學變化是導致腦卒中的誘因之一,尤其是在氣溫較低的冬季和較高的夏季(Simpson,et a1,2015;廉慧,2016)。另外,氣溫的變化可影響血液內(nèi)成分的變化,如血小板、紅細胞、纖維蛋白原、凝血因子及脂蛋白等改變是誘發(fā)腦卒中的病因之一(Braune,et a1,2016;Cambras,et a1,2017)。

    除了氣溫,普遍還認為腦卒中發(fā)病率與氣壓呈正相關(Jimenez-Conde,et a1,2008;Larrieu,et a1,2007),但俄羅斯學者Shaposhnikov 等(2014)發(fā)現(xiàn)腦卒中隨日平均氣壓升高而穩(wěn)定下降,與氣壓的變化呈負相關;Honig 等(2016)和Hori 等(2012)均發(fā)現(xiàn)氣壓降低可增加腦深部出血的風險。Jimenez-Conde 等(2008)研究表明腦卒中的發(fā)病與氣壓、氣溫不存在相關,但其變化可導致腦卒中的發(fā)病增多。劉博等(2014)的研究也表明腦卒中對負變溫和正變壓反應敏感??諝庀鄬穸扰c腦卒中關系的探討各地結果差別較大,且缺乏明確的機制研究。歐洲的研究表明相對濕度越低腦血管病死亡率越高(Díaz,et a1,2002),而中國的研究則表明腦卒中病例數(shù)與相對濕度呈顯著正相關(Cao,et a1,2016),且當相對濕度超過74%時病例數(shù)明顯增多(黃明北等,1999)。Lim 等(2017)的研究發(fā)現(xiàn)韓國北部地區(qū)腦卒中發(fā)病與濕度變化不存在相關??傊?,氣象因素誘發(fā)腦卒中的機制復雜,不是單一機制而是多種機制疊加的共同作用。

    腦卒中發(fā)病與氣象因素關系的研究結果因地域差異而不盡相同,所以進行特定區(qū)域的細致研究顯得至關重要。近年來天津市腦卒中疾病負擔越來越重,相關研究多是基于其流行特征、危險因素、疾病負擔等方面開展(胡亞會,2020;王志梅,2018;魏常松等,2018),其中危險因素主要涉及生活方式、病史、家族史等,關于氣象環(huán)境因素的研究相對較少。文中利用2016—2020 年天津市腦卒中日住院數(shù)據(jù)和氣象數(shù)據(jù),開展高影響氣象因子不同時間尺度變化、滯后以及累積效應對腦卒中住院風險的影響研究,旨在為腦卒中的綜合防控提供參考,最終達到降低其發(fā)病率和死亡率的目的。

    2 資料與方法

    2.1 資 料

    腦卒中醫(yī)療數(shù)據(jù)來源于天津市健康醫(yī)療大數(shù)據(jù)有限公司,數(shù)據(jù)來源的醫(yī)院級別及數(shù)量見表1,基本能夠代表天津市腦卒中住院的整體情況。數(shù)據(jù)包括2016—2020 年的腦卒中日住院人數(shù)、患者年齡及性別等信息。同期氣象數(shù)據(jù)來源于天津市氣象局西青觀測站,包括氣壓、溫度、相對濕度、風速等氣象要素的小時數(shù)據(jù)及日均數(shù)據(jù)。

    表1 醫(yī)療數(shù)據(jù)來源Table 1 Number of hospitals based on medical data source

    2.2 方法

    2.2.1 累積概率分布與斯皮爾曼(Spearman)相關系數(shù)

    采用腦卒中日住院人數(shù)累積概率分布對其住院風險進行等級劃分。根據(jù)日住院人數(shù)累積比例的25%、50%、75%和95%分位數(shù),確定低風險、中低風險、中風險、中高風險、高風險5 個風險等級。采用Spearman 相關系數(shù)進行氣象因子與腦卒中各風險等級住院人數(shù)的相關分析。分析各風險等級腦卒中住院人數(shù)與當天及滯后6 d 內(nèi)氣象因子的相關,確定高影響氣象因子。

    2.2.2 分布滯后非線性模型與半?yún)?shù)廣義相加模型

    為研究氣象因素對腦卒中影響的滯后效應和累積效應,本研究利用R 軟件4.3.2 版本中的dlnm程序包進行分布滯后非線性模型(Distributed Lag Non-linear Model,DLNM)(Gasparrini,et a1,2010)建模,擬合高影響氣象因子的“交叉基”函數(shù),確定影響的最大滯后時間?;谝褬嫿ê玫摹敖徊婊焙瘮?shù),利用R 軟件4.3.2 版本中的mgcv 程序包進行半?yún)?shù)廣義相加模型(Generalized Additive Models,GAM)(Hastie,et a1,1995)建模,明確高影響氣象因子與腦卒中住院人數(shù)的暴露-反應關系。為驗證模型的穩(wěn)定性,本研究通過改變模型中高影響氣象因子的時間序列自由度(3—6),對其進行敏感性分析。

    3 腦卒中住院人群與住院風險分布特征

    3.1 腦卒中住院人數(shù)整體情況

    從天津市2016—2020 年腦卒中住院人數(shù)月變化(圖1)可見,2016—2020 年腦卒中天津市全年住院人數(shù)依次為139779、179408、201488、227533 和133387。其中2018 和2019 年全年住院人數(shù)超過20 萬,2019 年全年住院人數(shù)最多,2020 年可能受新型冠狀病毒肺炎疫情的影響住院人數(shù)明顯減少。結合2016—2020 年數(shù)據(jù)來源的醫(yī)院數(shù)量并不能得出天津市近年腦卒中住院人數(shù)呈現(xiàn)增多的趨勢,但2019 與2020 年數(shù)據(jù)來源的醫(yī)院數(shù)基本已囊括天津市的二級以上醫(yī)院,所以2019 年數(shù)據(jù)在一定程度上可代表天津市腦卒中住院人數(shù)的基本情況。腦卒中住院人數(shù)在全年中存在兩個峰值,分別出現(xiàn)在3—4 月和11—12 月,2 月腦卒中的住院人數(shù)遠少于其他月份,這可能與傳統(tǒng)節(jié)日春節(jié)的習俗有關。2019 年腦卒中日平均住院人數(shù)為623;2 月的住院人數(shù)最少(14063),4 月的住院人數(shù)最多(21135),其次為12 月(21101)。

    圖1 2016—2020 年腦卒中住院人數(shù)月變化Fig.1 Changes in stroke inpatients from 2016 to 2020

    3.2 腦卒中住院人群特征

    天津市腦卒中男、女住院人數(shù)比為8:5,男性住院患者明顯多于女性(圖2a)。男、女住院人數(shù)在全年中分布基本一致。50 歲以上尤其是60 歲以上人群為易患腦卒中疾病的高危人群(圖2b)。50 歲以上男、女住院患者分別占各自總人數(shù)的92.4%、97.2%。男、女住院患者均為70 歲以上的比例最高,占比分別為37.7%和51.4%。50 歲以上人群為腦卒中的危險人群的主要原因是老年人身體機能衰退,對外界和體內(nèi)環(huán)境改變的適應能力降低。當氣溫等發(fā)生劇烈變化時,老年人的機體調(diào)節(jié)能力往往不足以應對氣象條件的驟變對機體的影響。

    圖2 2016—2020 年月均腦卒中住院患者男、女(a)人數(shù)和(b)年齡分布Fig.2 Male and female(a)number of inpatients and(b)age distribution of stroke inpatients from 2016 to 2020

    3.3 腦卒中住院風險分布特征

    基于腦卒中日住院人數(shù)累積概率分布的25%、50%、75%和95%分位數(shù)確定低風險、中低風險、中風險、中高風險、高風險等級(圖3)。日住院人數(shù)累積概率達到25%、50%、75%和95%分位數(shù)時對應的日住院人數(shù)分別為517、619、705 和812 人。各風險等級對應的日住院人數(shù)范圍及日數(shù)占比見表2。全年住院中高風險等級以上的日數(shù)占比為18.3%,其中高風險占比3.3%。中高風險等級對應的日住院人數(shù)范圍為706 至812。高風險等級的日住院人數(shù)在812 人以上。

    圖3 腦卒中日住院人數(shù)的累積概率分布Fig.3 Cumulative probability distribution of daily stroke hospitalization

    表2 腦卒中5 種住院風險等級劃分統(tǒng)計Table 2 Statistics for the classification of five types of stroke hospitalization risk level

    5 種住院風險等級在各月出現(xiàn)的概率如圖4 所示,住院中高以上風險等級出現(xiàn)概率在全年中存在2 個峰值,分別為11 月和4 月,概率為41.1%和36.7%。12 月中高以上風險等級出現(xiàn)概率也較高,為33.3%。可見深秋、初冬和春季為天津腦卒中住院高風險季節(jié),且深秋的風險高于春季。同時住院中高以上風險等級出現(xiàn)概率在全年中存在2 個谷值,分別為2 月和8 月,這一結果與中國其他城市冬季腦卒中高發(fā)不同,但與夏季低發(fā)一致(林亞楠等,2018;劉娜等,2018)。

    圖4 各月腦卒中住院風險等級出現(xiàn)概率(綠線為中高風險與高風險出現(xiàn)概率之和)Fig.4 Occurrence probabilities of five stroke hospitalization risk levels in each month(green line is the sum of the occurrence probabilities of medium to high risk and high risk)

    4 腦卒中高住院風險的氣象原因探析

    4.1 氣象因素與腦卒中住院人數(shù)的相關分析

    住院人數(shù)與氣象因子呈較好的相關關系,且氣象因子與高風險等級的住院人數(shù)相關最強。通過相關分析得到高影響氣象因子為氣壓、氣溫、氣溫日較差、24 h 變壓和24 h 變溫。表3 為高風險等級日住院人數(shù)與氣象因子累積影響的相關系數(shù)。日住院人數(shù)與氣壓累積2 d 的影響相關系數(shù)最大,且為正相關,其中與最高氣壓的相關系數(shù)最大(0.424)。與氣溫累積3 d 的影響相關系數(shù)最大,呈負相關,與最高氣溫的相關系數(shù)為-0.413。腦卒中住院人數(shù)與氣壓呈正相關,與氣溫呈負相關,這一結果與太原和上海等地的研究結果一致(王瑞,2015)。同時腦卒中住院人數(shù)還與氣壓和氣溫日較差存在相關,與24 h 變壓呈正相關,與24 h 變溫呈負相關。

    表3 當天及前7 d 氣象因子累積與高風險等級住院人數(shù)的相關系數(shù)Table 3 Correlation coefficients of the accumulation of meteorological factors on the day and the previous 7 d with the number of high-risk hospitalizations

    4.2 高影響氣象因子與腦卒中住院人數(shù)的關系

    為得到腦卒中住院人數(shù)在高影響氣象因子取值范圍內(nèi)的分布特征,將日最高氣壓、日最低氣溫、24 h 變壓和24 h 變溫的取值范圍按其0.5 個單位劃分為等間隔的小區(qū)間,并在每個區(qū)間內(nèi)求腦卒中住院人數(shù)的均值,繪制氣泡圖,氣泡的大小和顏色深淺反映了所在區(qū)間內(nèi)數(shù)據(jù)點的數(shù)量。當日最高氣壓低于1042 hPa 時,日住院人數(shù)基本維持在700 人以下的水平,當日最高氣壓高于1042 hPa時,日住院人數(shù)急劇增多,均在850 人以上(圖5a)。24 h 變壓在正、負4 hPa 的區(qū)間內(nèi),腦卒中住院人數(shù)變化不大,在500 至600 人之間波動,且該區(qū)域樣本數(shù)量較大。隨24 h 變壓絕對值增大,住院人數(shù)呈增長趨勢。24 h 負變壓10—15 hPa 和24 h 正變壓大于18 hPa 的情況下,腦卒中住院人數(shù)相對較多(圖5b)。在日最低氣溫低于-12℃時,住院人數(shù)相對較多(圖5c)。當24 h 變溫為-5℃—7℃時,住院人數(shù)無明顯波動,隨變溫增大,住院人數(shù)呈增多趨勢。且負變溫情況的住院人數(shù)多于正變溫(圖5d)。整體來看在氣溫較低、氣壓較高、24 h 變壓和變溫較大的情況下對應的腦卒中住院人數(shù)較多。

    圖5 (a)日最高氣壓、(b)24 h 變壓、(c)日最低氣溫、(d)24 h 變溫與腦卒中住院人數(shù)的關系Fig.5 Relationship between(a)daily maximum barometric pressure,(b)24 h pressure variation,(c)daily minimum temperature,(d)24 h temperature variation and the number of stroke hospitalizations

    4.3 不同時間尺度氣溫變化對腦卒中住院人數(shù)的影響

    4.3.1 月變溫對腦卒中住院人數(shù)的影響

    天津全年月均溫度呈倒“U”型分布,7 月最高,1 月最低。圖6 中月變溫表示本月平均氣溫與上月平均氣溫的差,2 月到7 月月均氣溫呈上升趨勢,其中3—5 月升溫趨勢明顯,升溫幅度均在6℃以上。8 月至次年1 月月均氣溫呈下降趨勢,10—12 月降溫趨勢明顯,降溫幅度均在7℃以上。3 月和10 月為正、負月變溫最大的月份,均超過了8.5℃,中高以上住院風險出現(xiàn)概率在4 月和11 月出現(xiàn)了峰值,其中負月變溫顯著的深秋和初冬腦卒中住院風險最高,正月變溫顯著的春季腦卒中住院風險次高。腦卒中中高以上住院風險出現(xiàn)概率峰值滯后月變溫峰值1 個月,即出現(xiàn)最大正、負月變溫的下一月,腦卒中高住院風險出現(xiàn)的概率最大。而月變溫較小的1、2、7 和8 月為腦卒中住院的低風險期??梢?,顯著的月變溫為影響腦卒中住院的重要因素,且月變溫為負值的深秋腦卒中住院風險最高。

    圖6 月變溫與腦卒中住院風險的關系Fig.6 Relationship between monthly temperature change and hospitalization risk of stroke

    4.3.2 24 h 負變溫對腦卒中住院人數(shù)的影響

    由24 h 負變溫與腦卒中住院人數(shù)的關系(圖7)可知,24 h 負變溫對3、4、5、10、11、12 月的腦卒中住院人數(shù)影響較大。3—5 月24 h 負變溫超過6℃時住院人數(shù)明顯增多,超過720 人。24 h 負變溫在11 和12 月對腦卒中住院人數(shù)影響最大,24 h負變溫超過4℃時住院人數(shù)可超過740 人,由于11 和12 月負月變溫顯著,再疊加上24 h 負變溫,可增強不同尺度負變溫對腦卒中住院人數(shù)的影響。所以24 h 負變溫在深秋、初冬和春季對腦卒中住院人數(shù)影響較大,對夏季腦卒中的影響相反,24 h 降溫幅度越大,腦卒中住院人數(shù)反而越少。所以顯著的月際溫度變化與顯著的日尺度的24 h 變溫相疊加,可對腦卒中的危險人群造成較大影響。

    圖7 24 h 負變溫與腦卒中住院人數(shù)的關系Fig.7 Relationship between 24 h negative temperature change and the number of stroke hospitalizations

    5 高影響氣象因子對腦卒中住院風險的影響

    5.1 高影響氣象因子的滯后效應

    圖8 為高影響氣象因子與腦卒中住院人數(shù)的暴露-滯后-反應關系。30℃以上高溫和-5℃以下低溫均可引起腦卒中住院高風險,且對腦卒中住院人數(shù)產(chǎn)生即時影響,但低溫的風險高于高溫(圖8a)。隨滯后天數(shù)增加,高溫效應逐漸消失,低溫效應維持,低溫效應在滯后2 至3 d 內(nèi)達到最強。24 h 變溫在波動較小時腦卒中住院風險較低,隨變溫絕對值增大住院風險升高,且負變溫的風險高于正變溫(圖8b)。當24 h 正變溫超過8℃時,住院風險明顯升高,且強的正變溫對住院人數(shù)的影響表現(xiàn)為“即時效應”,主要影響當天和滯后1 d。7℃以上24 h負變溫存在明顯的“滯后效應”,在滯后3—5 d 的住院風險最高。氣溫日較差對腦卒中住院風險的影響表現(xiàn)為明顯的“滯后效應”(圖8c)。當氣溫日較差為14—20℃,滯后6—8 d 的住院風險最高。當氣溫日較差高于20℃時,在滯后2—4 d 的住院風險最高。所以隨氣溫日較差增大,住院高風險的出現(xiàn)時間提前。

    圖8 (a)日平均氣溫、(b)24 h 變溫、(c)氣溫日較差與腦卒中住院人數(shù)的暴露-滯后-反應關系(色階:相對危險度)Fig.8 Exposure-lag-response relationship between(a)daily average temperature,(b)24 h temperature variation,(c)daily temperature range and stroke hospitalizations(shaded:relative risk)

    5.2 高影響氣象因子的累積效應

    通過高影響氣象因子的滯后效應可知產(chǎn)生影響的最長滯后天數(shù)可達10 d,但均在滯后3 d 左右對相對危險度產(chǎn)生最大影響,所以選取最大影響的滯后天數(shù)為3 d,研究當天和滯后3 d 的高影響因子對腦卒中住院人數(shù)影響的累積效應。圖9 為累積3 d 的高影響氣象因子與腦卒中住院人數(shù)的暴露-反應關系。氣溫對腦卒中住院風險的影響表現(xiàn)為“低溫效應”,平均氣溫低于0℃時腦卒中住院風險較高,隨溫度進一步降低,風險呈近似線性的增長趨勢。24 h 負變溫的風險大于正變溫,隨24 h 負變溫增大住院相對風險逐漸增大。氣溫日較差與腦卒中住院人數(shù)的暴露-反應關系均近似呈“J”型分布。氣溫日較差高于17℃時腦卒中住院風險較高,且隨氣溫日較差增大而增大。

    圖9 (a)平均氣溫、(b)24 h 變溫、(c)氣溫日較差與腦卒中住院人數(shù)的暴露-反應關系(陰影:相對危險度95%置信區(qū)間)Fig.9 Exposure-response relationship between(a)mean temperature,(b)24 h temperature variation,(c)daily temperature range and stroke hospitalizations(shaded:95% confidence interval of relative risk )

    5.3 敏感性分析結果

    為驗證模型的穩(wěn)定性進行了敏感性分析試驗。結果(圖10)表明,通過改變平均氣溫和24 h變溫模型中時間序列自由度(3—6),本研究所得結果亦較為穩(wěn)定,相對危險度值變化幅度較小,如改變平均氣溫的時間序列自由度,腦卒中住院人數(shù)的相對危險度值變化范圍為1.00104—1.00129。改變24 h 變溫的時間序列自由度,腦卒中住院人數(shù)的相對危險度值變化范圍為1.00653—1.00723。將交互作用項添加到半?yún)?shù)廣義相加模型后對結果有影響,但與調(diào)整前變化不大,模型結果基本穩(wěn)定。

    圖10 (a)平均氣溫和(b)24 h 變溫在不同時間自由度下對住院風險影響的敏感性分析Fig.10 Sensitivity analysis of the effects of(a)mean temperature and(b)24 h temperature variation on the risk of hospitalization at different time degrees of freedom

    6 結論

    基于2016—2020 年天津市腦卒中日住院數(shù)據(jù)和氣象數(shù)據(jù),采用腦卒中日住院人數(shù)累積概率分布進行住院風險等級劃分,通過Spearman 相關分析確定高影響氣象因子,分析高影響氣象因子的變化對腦卒中住院的影響。在此基礎上利用分布滯后非線性模型與半?yún)?shù)廣義相加模型得到高影響氣象因子對腦卒中住院的滯后以及累積效應。主要結論如下:

    (1)天津市腦卒中年住院人數(shù)為22.5 萬左右,最高月住院人數(shù)可達2.1 萬,日均住院人數(shù)為623。男、女住院人數(shù)比為8:5,50 歲以上人群為易患腦卒中疾病的高危人群,50 歲以上男、女住院患者分別占各自總人數(shù)的92.4%、97.2%。

    (2)氣象因子與住院人數(shù)呈較好的相關關系,尤其是與高風險等級的住院人數(shù)相關系數(shù)最大。高相關氣象因子為氣壓、氣溫、氣溫日較差、24 h變壓和24 h 變溫。且與氣壓呈正相關,與氣溫呈負相關,與氣溫日較差呈負相關,與24 h 變壓呈正相關,與24 h 變溫呈負相關。

    (3)深秋、初冬和春季為天津市腦卒中住院的高風險季節(jié),深秋的住院風險最高。腦卒中住院高風險出現(xiàn)概率峰值滯后月變溫峰值1 個月。6℃以上顯著的月變溫和6℃以上的24 h 負變溫相疊加可增大腦卒中住院風險,且月變溫為負變溫時,腦卒中住院風險最高。

    (4)30℃以上高溫和-5℃以下低溫均可引起腦卒中住院高風險,但低溫的風險高于高溫,低溫效應在滯后2—3 d 內(nèi)達到最強。24 h 負變溫的住院風險高于正變溫,在深秋、初冬和春季對腦卒中住院影響最大,且變溫強度越大住院風險越高,7℃以上24 h 負變溫在滯后3—5 d 的住院風險最高。

    天津市腦卒中高住院風險的氣象原因主要是月尺度和日尺度溫度劇烈變化的疊加,其次為極端氣溫及其24 h 劇烈變化的影響。本研究結果存在一定的局限性,一是研究中氣象數(shù)據(jù)僅來自1 個氣象站,無法代表個體暴露水平,會導致結果存在一定偏差;二是受限于醫(yī)療數(shù)據(jù)收集,對腦卒中病人住院的具體醫(yī)學原因不得而知,沒有進行具體病因下氣象原因的探討,需要進一步的深入研究。

    致 謝:感謝成都信息工程大學王式功教授在研究過程中給予的指導與幫助!

    猜你喜歡
    變溫氣壓天津市
    看不見的氣壓
    幼兒畫刊(2021年5期)2021-12-02 04:24:04
    有趣的云
    如果畫筆會說話?
    藝術啟蒙(2020年6期)2020-07-16 03:48:04
    壓力容器氣壓端蓋注射模設計
    模具制造(2019年4期)2019-06-24 03:36:46
    氯乙烯生產(chǎn)中變溫吸附脫水工藝的使用及改進
    天津市第三屆“未來之星”初中數(shù)學邀請賽
    凍融處理對甘薯變溫壓差膨化干燥動力學的影響
    非共面四頻激光陀螺變溫零偏周期性波動
    Mn摻雜ZnSe量子點變溫發(fā)光性質(zhì)研究
    中國光學(2015年5期)2015-12-09 09:00:43
    電滲—堆載聯(lián)合氣壓劈烈的室內(nèi)模型試驗
    欧美另类亚洲清纯唯美| 久久精品91无色码中文字幕| 老熟妇仑乱视频hdxx| 女同久久另类99精品国产91| 久久久久国产精品人妻aⅴ院 | av不卡在线播放| 99久久99久久久精品蜜桃| 国产单亲对白刺激| 波多野结衣一区麻豆| 午夜福利在线观看吧| 91av网站免费观看| 99国产综合亚洲精品| 少妇裸体淫交视频免费看高清 | 人成视频在线观看免费观看| av网站免费在线观看视频| a级毛片在线看网站| 在线观看午夜福利视频| 一级片免费观看大全| 电影成人av| 国产黄色免费在线视频| 久久久久国内视频| 久久久久国内视频| 后天国语完整版免费观看| 成人手机av| 日本一区二区免费在线视频| 精品一区二区三区四区五区乱码| 亚洲成a人片在线一区二区| 久久久久久亚洲精品国产蜜桃av| 亚洲成国产人片在线观看| www.精华液| 国产精品欧美亚洲77777| 高清在线国产一区| 丰满迷人的少妇在线观看| 中文字幕av电影在线播放| 欧美最黄视频在线播放免费 | 十分钟在线观看高清视频www| 18禁裸乳无遮挡动漫免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 可以免费在线观看a视频的电影网站| 老司机福利观看| 欧美成人免费av一区二区三区 | 亚洲国产精品合色在线| 国产亚洲一区二区精品| 亚洲全国av大片| 国产高清激情床上av| 日本黄色日本黄色录像| 91av网站免费观看| 免费看十八禁软件| 欧美人与性动交α欧美精品济南到| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩亚洲国产一区二区在线观看 | 啦啦啦 在线观看视频| 中文字幕人妻丝袜一区二区| 日韩视频一区二区在线观看| 俄罗斯特黄特色一大片| 久久亚洲真实| 免费观看a级毛片全部| 久久久久久亚洲精品国产蜜桃av| 久久久久精品国产欧美久久久| 我的亚洲天堂| 99久久99久久久精品蜜桃| 少妇 在线观看| 欧美日韩亚洲高清精品| 免费看a级黄色片| 国产主播在线观看一区二区| 大香蕉久久网| 欧美日韩乱码在线| 丝袜人妻中文字幕| 18禁黄网站禁片午夜丰满| 免费不卡黄色视频| 一级片'在线观看视频| 亚洲精品国产色婷婷电影| 男女高潮啪啪啪动态图| 高潮久久久久久久久久久不卡| 日日爽夜夜爽网站| 淫妇啪啪啪对白视频| 91在线观看av| 变态另类成人亚洲欧美熟女 | 午夜老司机福利片| 久久久久国产一级毛片高清牌| 男女高潮啪啪啪动态图| 99精品在免费线老司机午夜| 欧美日韩国产mv在线观看视频| 激情在线观看视频在线高清 | 久久性视频一级片| 亚洲一区高清亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产中文字幕在线视频| 精品少妇久久久久久888优播| av电影中文网址| 手机成人av网站| 日韩中文字幕欧美一区二区| 亚洲少妇的诱惑av| 另类亚洲欧美激情| 久久久久久久午夜电影 | 叶爱在线成人免费视频播放| 最近最新免费中文字幕在线| 久99久视频精品免费| 最新在线观看一区二区三区| 色尼玛亚洲综合影院| 日本wwww免费看| svipshipincom国产片| 精品国内亚洲2022精品成人 | 亚洲av成人不卡在线观看播放网| 国产无遮挡羞羞视频在线观看| 村上凉子中文字幕在线| a级片在线免费高清观看视频| 深夜精品福利| 国产精品久久久av美女十八| 国产不卡av网站在线观看| av有码第一页| 精品福利永久在线观看| 在线观看免费午夜福利视频| 国产成人欧美在线观看 | 欧美午夜高清在线| 岛国毛片在线播放| 夜夜爽天天搞| 波多野结衣一区麻豆| 欧美成人免费av一区二区三区 | 久久精品国产a三级三级三级| 久久精品国产亚洲av香蕉五月 | 真人做人爱边吃奶动态| 欧美日韩瑟瑟在线播放| 大型黄色视频在线免费观看| 亚洲片人在线观看| 国产一区二区三区综合在线观看| 免费女性裸体啪啪无遮挡网站| 欧美日韩视频精品一区| 精品欧美一区二区三区在线| 怎么达到女性高潮| 777久久人妻少妇嫩草av网站| 中文字幕制服av| 国产主播在线观看一区二区| 亚洲性夜色夜夜综合| a在线观看视频网站| 欧美丝袜亚洲另类 | 成人三级做爰电影| 高潮久久久久久久久久久不卡| 久久精品亚洲熟妇少妇任你| 大型黄色视频在线免费观看| av视频免费观看在线观看| 人妻一区二区av| 亚洲国产精品一区二区三区在线| av福利片在线| 久久ye,这里只有精品| 午夜福利影视在线免费观看| 婷婷成人精品国产| 午夜福利欧美成人| 国产欧美亚洲国产| 精品免费久久久久久久清纯 | 国产精品99久久99久久久不卡| 91字幕亚洲| 免费在线观看完整版高清| 51午夜福利影视在线观看| 国产熟女午夜一区二区三区| 久久久国产一区二区| 9191精品国产免费久久| 黄色成人免费大全| 十八禁高潮呻吟视频| 91精品三级在线观看| 久久精品成人免费网站| 久久久久视频综合| 色综合婷婷激情| 又紧又爽又黄一区二区| 亚洲欧美色中文字幕在线| 成在线人永久免费视频| 制服诱惑二区| 自拍欧美九色日韩亚洲蝌蚪91| 91精品国产国语对白视频| 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区三| 国产伦人伦偷精品视频| 久久这里只有精品19| 怎么达到女性高潮| 下体分泌物呈黄色| 男女午夜视频在线观看| 亚洲av美国av| 啪啪无遮挡十八禁网站| 日本一区二区免费在线视频| 黄色a级毛片大全视频| 成人特级黄色片久久久久久久| 国产又爽黄色视频| av有码第一页| 69精品国产乱码久久久| 我的亚洲天堂| 变态另类成人亚洲欧美熟女 | 高清av免费在线| 国产成+人综合+亚洲专区| 欧美日韩亚洲高清精品| 欧美黄色片欧美黄色片| 午夜免费鲁丝| 国产精品成人在线| 欧美另类亚洲清纯唯美| www日本在线高清视频| 91成人精品电影| 国产亚洲欧美在线一区二区| 午夜精品久久久久久毛片777| 欧美国产精品一级二级三级| 国产乱人伦免费视频| 久久久国产欧美日韩av| 老司机影院毛片| 久久久国产成人免费| 久久久国产一区二区| 欧美 亚洲 国产 日韩一| 三上悠亚av全集在线观看| 国产欧美日韩一区二区三| 一级片免费观看大全| 99久久精品国产亚洲精品| 夜夜躁狠狠躁天天躁| 少妇 在线观看| 欧美在线黄色| 天天添夜夜摸| 香蕉久久夜色| 色老头精品视频在线观看| 午夜免费成人在线视频| 国产午夜精品久久久久久| 69精品国产乱码久久久| √禁漫天堂资源中文www| 老熟妇乱子伦视频在线观看| 在线观看一区二区三区激情| 成人亚洲精品一区在线观看| 超色免费av| 国产男女内射视频| 99国产精品一区二区三区| 精品国产超薄肉色丝袜足j| 1024香蕉在线观看| 99riav亚洲国产免费| 免费久久久久久久精品成人欧美视频| 欧美精品啪啪一区二区三区| tube8黄色片| 香蕉丝袜av| 国产国语露脸激情在线看| 国产精品一区二区在线不卡| 欧美日韩av久久| 巨乳人妻的诱惑在线观看| 欧美激情高清一区二区三区| 人人妻人人澡人人看| 国产精品一区二区免费欧美| 午夜两性在线视频| 在线永久观看黄色视频| 天天添夜夜摸| 女警被强在线播放| 制服人妻中文乱码| 亚洲精品成人av观看孕妇| 亚洲精品久久成人aⅴ小说| 精品久久久久久久毛片微露脸| 久99久视频精品免费| 国产成人精品久久二区二区免费| 999久久久精品免费观看国产| 久久天堂一区二区三区四区| 国产亚洲精品久久久久5区| 欧美丝袜亚洲另类 | 91麻豆av在线| 老司机深夜福利视频在线观看| 国产男靠女视频免费网站| 露出奶头的视频| 操出白浆在线播放| 国产片内射在线| 麻豆成人av在线观看| 国产精品二区激情视频| 一级片免费观看大全| 老司机福利观看| 曰老女人黄片| 村上凉子中文字幕在线| 午夜亚洲福利在线播放| 国产成人免费观看mmmm| 亚洲av片天天在线观看| 人人妻,人人澡人人爽秒播| 免费高清在线观看日韩| 国产精品秋霞免费鲁丝片| 咕卡用的链子| 韩国精品一区二区三区| 成人免费观看视频高清| 精品无人区乱码1区二区| 久久 成人 亚洲| 免费在线观看影片大全网站| 国产一区二区三区综合在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 成人精品一区二区免费| 国产男靠女视频免费网站| 99国产综合亚洲精品| 在线av久久热| 欧美不卡视频在线免费观看 | 高清欧美精品videossex| 日韩成人在线观看一区二区三区| 日韩视频一区二区在线观看| 亚洲免费av在线视频| 精品久久久久久久久久免费视频 | 91大片在线观看| 又紧又爽又黄一区二区| 国产一卡二卡三卡精品| 两性夫妻黄色片| 脱女人内裤的视频| 身体一侧抽搐| 国产午夜精品久久久久久| 性少妇av在线| 成年人午夜在线观看视频| 国产成人精品在线电影| 亚洲欧美一区二区三区黑人| 欧美另类亚洲清纯唯美| 丝袜美足系列| 亚洲av成人不卡在线观看播放网| 日韩视频一区二区在线观看| 99久久精品国产亚洲精品| 又黄又粗又硬又大视频| 精品国产一区二区三区四区第35| 69av精品久久久久久| 国产精品永久免费网站| 麻豆乱淫一区二区| 国产不卡av网站在线观看| 一区二区三区精品91| 久久ye,这里只有精品| 一区二区三区激情视频| 老汉色av国产亚洲站长工具| 一进一出抽搐动态| 一本综合久久免费| 美女午夜性视频免费| 欧美日韩成人在线一区二区| 欧美久久黑人一区二区| 日韩免费av在线播放| 国产亚洲精品久久久久5区| 中文字幕色久视频| 777久久人妻少妇嫩草av网站| 亚洲国产看品久久| 国产精品久久视频播放| 人人澡人人妻人| 在线观看www视频免费| 成年动漫av网址| 天天操日日干夜夜撸| 久久久久久久久免费视频了| 欧美亚洲日本最大视频资源| 成在线人永久免费视频| 午夜久久久在线观看| 侵犯人妻中文字幕一二三四区| 一本一本久久a久久精品综合妖精| 亚洲人成电影免费在线| 美国免费a级毛片| 亚洲片人在线观看| 成人黄色视频免费在线看| 99久久99久久久精品蜜桃| 国产蜜桃级精品一区二区三区 | www日本在线高清视频| 国产精品久久久人人做人人爽| 国产欧美日韩综合在线一区二区| 国产精品久久久久成人av| 黄色女人牲交| 99精品久久久久人妻精品| 成人永久免费在线观看视频| 国产高清激情床上av| 久久九九热精品免费| 亚洲国产欧美网| 亚洲国产精品sss在线观看 | 久久久国产精品麻豆| 熟女少妇亚洲综合色aaa.| 久久婷婷成人综合色麻豆| 久久人妻福利社区极品人妻图片| 伦理电影免费视频| 老司机靠b影院| 欧美激情高清一区二区三区| 国内久久婷婷六月综合欲色啪| 国产成人欧美| 狠狠狠狠99中文字幕| 夜夜躁狠狠躁天天躁| 欧美精品高潮呻吟av久久| 国产黄色免费在线视频| 久久精品国产综合久久久| 亚洲国产精品sss在线观看 | 亚洲熟妇中文字幕五十中出 | 国产精品免费大片| 不卡一级毛片| 性色av乱码一区二区三区2| 丰满迷人的少妇在线观看| 国产男靠女视频免费网站| 大型av网站在线播放| av天堂久久9| 99久久人妻综合| 亚洲五月天丁香| 免费观看精品视频网站| avwww免费| 亚洲一卡2卡3卡4卡5卡精品中文| 不卡av一区二区三区| 午夜福利,免费看| 操美女的视频在线观看| 亚洲精品成人av观看孕妇| 成在线人永久免费视频| www.熟女人妻精品国产| 天天躁夜夜躁狠狠躁躁| 一夜夜www| 亚洲色图综合在线观看| 国产在线观看jvid| 久久性视频一级片| 91在线观看av| 另类亚洲欧美激情| 很黄的视频免费| 精品少妇一区二区三区视频日本电影| 大香蕉久久网| 一级片'在线观看视频| 在线观看www视频免费| 午夜亚洲福利在线播放| 亚洲精品av麻豆狂野| 久久精品熟女亚洲av麻豆精品| 欧美日韩一级在线毛片| 国产精品免费视频内射| 久久精品人人爽人人爽视色| 正在播放国产对白刺激| 精品亚洲成a人片在线观看| 中亚洲国语对白在线视频| 热re99久久国产66热| 丰满的人妻完整版| 久久久精品国产亚洲av高清涩受| 又紧又爽又黄一区二区| 免费女性裸体啪啪无遮挡网站| 久久青草综合色| 国产蜜桃级精品一区二区三区 | 精品人妻1区二区| 两性夫妻黄色片| 国产精品99久久99久久久不卡| 欧美日韩中文字幕国产精品一区二区三区 | 精品福利永久在线观看| 国产精品一区二区在线不卡| 亚洲国产毛片av蜜桃av| 成人手机av| 激情在线观看视频在线高清 | 看片在线看免费视频| 日韩欧美一区视频在线观看| 女人精品久久久久毛片| 大型黄色视频在线免费观看| 一区二区三区精品91| 亚洲午夜精品一区,二区,三区| 亚洲人成77777在线视频| 两个人看的免费小视频| 精品视频人人做人人爽| 国产极品粉嫩免费观看在线| 成年版毛片免费区| 啦啦啦 在线观看视频| svipshipincom国产片| 91成人精品电影| 久久香蕉激情| 一边摸一边做爽爽视频免费| 国产成人精品无人区| 老司机深夜福利视频在线观看| 女人久久www免费人成看片| 国产人伦9x9x在线观看| tube8黄色片| 18禁黄网站禁片午夜丰满| 亚洲一区中文字幕在线| av有码第一页| 9色porny在线观看| 国产又色又爽无遮挡免费看| 免费在线观看亚洲国产| 12—13女人毛片做爰片一| 一级,二级,三级黄色视频| 又黄又粗又硬又大视频| 国产亚洲欧美在线一区二区| 国产在线精品亚洲第一网站| 亚洲第一av免费看| 国产欧美日韩综合在线一区二区| 精品少妇一区二区三区视频日本电影| 天堂中文最新版在线下载| 成人精品一区二区免费| 捣出白浆h1v1| 黑人欧美特级aaaaaa片| 在线观看免费午夜福利视频| 水蜜桃什么品种好| 香蕉丝袜av| 啦啦啦视频在线资源免费观看| tube8黄色片| 国产麻豆69| 成人av一区二区三区在线看| 丝袜人妻中文字幕| 国产精品久久久久成人av| 色94色欧美一区二区| 国产在线精品亚洲第一网站| 国产亚洲av高清不卡| 亚洲av日韩精品久久久久久密| 午夜老司机福利片| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利免费观看在线| 人人澡人人妻人| 亚洲视频免费观看视频| 精品无人区乱码1区二区| 日韩中文字幕欧美一区二区| 亚洲av熟女| 搡老熟女国产l中国老女人| 精品久久久久久久毛片微露脸| 丁香六月欧美| 人人妻人人爽人人添夜夜欢视频| 夫妻午夜视频| 国产xxxxx性猛交| 女警被强在线播放| 亚洲欧美激情在线| 电影成人av| 曰老女人黄片| 亚洲av熟女| 女性被躁到高潮视频| 国产男女超爽视频在线观看| 日日摸夜夜添夜夜添小说| 捣出白浆h1v1| 免费在线观看日本一区| 男女床上黄色一级片免费看| 在线播放国产精品三级| 美女午夜性视频免费| 中文字幕色久视频| 国产一区有黄有色的免费视频| 精品久久久久久电影网| 伊人久久大香线蕉亚洲五| 一级作爱视频免费观看| 久久婷婷成人综合色麻豆| 波多野结衣一区麻豆| 在线观看www视频免费| 中文欧美无线码| 久热爱精品视频在线9| 国产精品久久视频播放| 97人妻天天添夜夜摸| 色精品久久人妻99蜜桃| 国产精品99久久99久久久不卡| 亚洲精华国产精华精| 国产精品1区2区在线观看. | 午夜精品在线福利| 久久久精品区二区三区| 中文欧美无线码| 在线观看www视频免费| 久久性视频一级片| 国产99白浆流出| 欧美日韩成人在线一区二区| 又紧又爽又黄一区二区| 嫁个100分男人电影在线观看| 欧美日韩一级在线毛片| 国产野战对白在线观看| 日韩欧美一区视频在线观看| 免费高清在线观看日韩| 校园春色视频在线观看| 一本大道久久a久久精品| 亚洲成人免费av在线播放| 最近最新中文字幕大全免费视频| 可以免费在线观看a视频的电影网站| av不卡在线播放| 亚洲第一青青草原| 亚洲欧美激情综合另类| 国产精品久久久久成人av| 久久精品aⅴ一区二区三区四区| 69精品国产乱码久久久| 欧美久久黑人一区二区| 新久久久久国产一级毛片| 国产成人精品久久二区二区91| 丰满迷人的少妇在线观看| 国产欧美日韩一区二区三区在线| 老鸭窝网址在线观看| 成人18禁高潮啪啪吃奶动态图| ponron亚洲| 日韩欧美免费精品| 欧美成人免费av一区二区三区 | 女同久久另类99精品国产91| 亚洲精品国产精品久久久不卡| 夜夜夜夜夜久久久久| 99国产精品一区二区三区| 人妻一区二区av| 十分钟在线观看高清视频www| 国产乱人伦免费视频| 久久中文看片网| 很黄的视频免费| 国产精品1区2区在线观看. | 欧美精品人与动牲交sv欧美| 欧美最黄视频在线播放免费 | 夫妻午夜视频| 丝瓜视频免费看黄片| 在线观看午夜福利视频| 99久久99久久久精品蜜桃| 久久天堂一区二区三区四区| 欧美日韩一级在线毛片| 手机成人av网站| 日本黄色日本黄色录像| 免费久久久久久久精品成人欧美视频| 婷婷成人精品国产| 91av网站免费观看| 精品少妇久久久久久888优播| 亚洲国产毛片av蜜桃av| 午夜亚洲福利在线播放| 最新美女视频免费是黄的| 欧美+亚洲+日韩+国产| 国产成人精品无人区| 黄色毛片三级朝国网站| 久久狼人影院| 在线观看免费午夜福利视频| 精品福利观看| 精品久久久久久久毛片微露脸| 亚洲精品国产色婷婷电影| 免费黄频网站在线观看国产| 国产av又大| a级毛片黄视频| 露出奶头的视频| 黄色a级毛片大全视频| 国产成人系列免费观看| 久久香蕉激情| 少妇 在线观看| 久久中文字幕人妻熟女| 真人做人爱边吃奶动态| 人人妻人人爽人人添夜夜欢视频| 午夜成年电影在线免费观看| 亚洲精品在线观看二区| av电影中文网址| 免费黄频网站在线观看国产| 99香蕉大伊视频| 亚洲精品一卡2卡三卡4卡5卡| 一个人免费在线观看的高清视频| 国产人伦9x9x在线观看| 嫩草影视91久久| 99精国产麻豆久久婷婷| 午夜免费鲁丝| 久久精品aⅴ一区二区三区四区| 亚洲综合色网址| ponron亚洲| 19禁男女啪啪无遮挡网站| 国产精品秋霞免费鲁丝片| 欧美人与性动交α欧美软件|