• <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)模型試驗
    少妇高潮的动态图| 亚洲精品影视一区二区三区av| 欧美日韩精品成人综合77777| 国产免费一级a男人的天堂| 国产精品亚洲一级av第二区| 天天一区二区日本电影三级| 美女被艹到高潮喷水动态| av在线老鸭窝| 又爽又黄a免费视频| 精品日产1卡2卡| 一边摸一边抽搐一进一小说| 亚洲性夜色夜夜综合| 天堂动漫精品| 3wmmmm亚洲av在线观看| 日产精品乱码卡一卡2卡三| 精品一区二区三区人妻视频| 小说图片视频综合网站| 久久久久国产网址| 联通29元200g的流量卡| 国产成人aa在线观看| 久久九九热精品免费| 六月丁香七月| 最近中文字幕高清免费大全6| 97在线视频观看| 成人三级黄色视频| 97碰自拍视频| 亚洲国产高清在线一区二区三| 色5月婷婷丁香| 亚洲在线自拍视频| 亚洲在线观看片| 最近最新中文字幕大全电影3| 别揉我奶头 嗯啊视频| 如何舔出高潮| 久久精品国产99精品国产亚洲性色| 亚洲欧美中文字幕日韩二区| 成人av在线播放网站| 老司机影院成人| 国产真实乱freesex| 可以在线观看的亚洲视频| 国产欧美日韩一区二区精品| 给我免费播放毛片高清在线观看| 日韩av在线大香蕉| 嫩草影院新地址| 91久久精品国产一区二区成人| 男女视频在线观看网站免费| 国产精品精品国产色婷婷| 最近视频中文字幕2019在线8| 一级毛片久久久久久久久女| 日韩欧美免费精品| 免费大片18禁| 亚洲av第一区精品v没综合| 久久精品综合一区二区三区| 深夜精品福利| 天堂动漫精品| 床上黄色一级片| 日韩欧美三级三区| 国产亚洲91精品色在线| 国产一区二区三区av在线 | 又爽又黄无遮挡网站| 一级黄片播放器| 欧美极品一区二区三区四区| 久久精品国产鲁丝片午夜精品| 精品无人区乱码1区二区| 国产中年淑女户外野战色| 亚洲欧美中文字幕日韩二区| 毛片女人毛片| 免费人成视频x8x8入口观看| 日韩成人av中文字幕在线观看 | 成年av动漫网址| 性欧美人与动物交配| 久久精品国产亚洲网站| 免费在线观看成人毛片| 中文字幕精品亚洲无线码一区| 亚洲成人久久爱视频| 狂野欧美白嫩少妇大欣赏| 秋霞在线观看毛片| 色av中文字幕| 日本三级黄在线观看| 亚洲国产精品国产精品| 国产高清三级在线| 免费观看精品视频网站| 国产欧美日韩精品一区二区| 亚洲精华国产精华液的使用体验 | 亚洲中文日韩欧美视频| 日日摸夜夜添夜夜添av毛片| 久久久久免费精品人妻一区二区| 午夜免费男女啪啪视频观看 | 高清日韩中文字幕在线| 亚洲熟妇熟女久久| 久久久久久国产a免费观看| 一级毛片电影观看 | 美女高潮的动态| 天堂网av新在线| 桃色一区二区三区在线观看| aaaaa片日本免费| 成人二区视频| 国产精品永久免费网站| av黄色大香蕉| 亚洲三级黄色毛片| 国产亚洲精品av在线| 国产精品99久久久久久久久| 精品不卡国产一区二区三区| 12—13女人毛片做爰片一| 国产男靠女视频免费网站| 免费看a级黄色片| 国产三级中文精品| 一级a爱片免费观看的视频| 久久精品国产清高在天天线| 十八禁网站免费在线| 国产精品日韩av在线免费观看| 欧美性感艳星| 亚洲欧美清纯卡通| 日韩亚洲欧美综合| 51国产日韩欧美| 日本黄色片子视频| 精品久久久久久久久久久久久| 国产精品一区二区免费欧美| 国产色婷婷99| 国产亚洲精品综合一区在线观看| 久久鲁丝午夜福利片| 亚洲电影在线观看av| 日产精品乱码卡一卡2卡三| 日本 av在线| 我的女老师完整版在线观看| 亚洲国产欧洲综合997久久,| 俄罗斯特黄特色一大片| 色播亚洲综合网| 日韩欧美精品免费久久| 男人的好看免费观看在线视频| 国产成人aa在线观看| 如何舔出高潮| 中文在线观看免费www的网站| 天堂av国产一区二区熟女人妻| 国产精品人妻久久久久久| 亚洲欧美日韩东京热| av.在线天堂| 欧美又色又爽又黄视频| www日本黄色视频网| 菩萨蛮人人尽说江南好唐韦庄 | 国产一区二区亚洲精品在线观看| 久久久精品94久久精品| 日本免费a在线| 午夜福利视频1000在线观看| 国产欧美日韩一区二区精品| 男人狂女人下面高潮的视频| 直男gayav资源| 国产在线精品亚洲第一网站| 波多野结衣巨乳人妻| 成年版毛片免费区| 久久精品91蜜桃| 久久国产乱子免费精品| 五月玫瑰六月丁香| 男女边吃奶边做爰视频| 听说在线观看完整版免费高清| 欧美激情在线99| 在线免费十八禁| 插阴视频在线观看视频| av卡一久久| 国产精品嫩草影院av在线观看| 女人被狂操c到高潮| 日韩欧美精品免费久久| 日韩高清综合在线| 我要搜黄色片| 99久国产av精品国产电影| 麻豆精品久久久久久蜜桃| 成年av动漫网址| 看片在线看免费视频| 亚洲经典国产精华液单| 久久久久久九九精品二区国产| 午夜爱爱视频在线播放| 国产单亲对白刺激| 久久久久精品国产欧美久久久| 国产亚洲精品av在线| 九九在线视频观看精品| 久久久久久九九精品二区国产| 少妇丰满av| 淫秽高清视频在线观看| 亚洲四区av| 极品教师在线视频| 国产又黄又爽又无遮挡在线| 亚洲无线观看免费| 国产免费一级a男人的天堂| 久久热精品热| 免费电影在线观看免费观看| 天美传媒精品一区二区| 人妻少妇偷人精品九色| 欧美国产日韩亚洲一区| 色5月婷婷丁香| 国产成人精品久久久久久| 91av网一区二区| 亚洲av美国av| 欧美激情在线99| 99热全是精品| 久久久精品大字幕| 久久鲁丝午夜福利片| 久久精品人妻少妇| 亚洲综合色惰| 91在线精品国自产拍蜜月| 十八禁网站免费在线| 色吧在线观看| 亚洲性夜色夜夜综合| 99久国产av精品国产电影| 精品久久久久久成人av| 尾随美女入室| 一区二区三区高清视频在线| 国产伦精品一区二区三区视频9| 床上黄色一级片| 亚洲精品久久国产高清桃花| 最新在线观看一区二区三区| 亚洲在线观看片| 国产白丝娇喘喷水9色精品| 免费看a级黄色片| 国产精品一区二区三区四区久久| 嫩草影院新地址| 久久精品国产亚洲网站| 精品久久久噜噜| 日韩精品中文字幕看吧| 亚洲真实伦在线观看| avwww免费| 久久综合国产亚洲精品| 亚洲成人久久爱视频| 亚洲在线观看片| 国语自产精品视频在线第100页| 成人毛片a级毛片在线播放| 免费av毛片视频| 91久久精品国产一区二区成人| 非洲黑人性xxxx精品又粗又长| 高清毛片免费看| 国内精品一区二区在线观看| 亚洲乱码一区二区免费版| 精品午夜福利在线看| 亚洲精品影视一区二区三区av| 欧美成人精品欧美一级黄| 亚洲欧美成人精品一区二区| 在线观看av片永久免费下载| 久久精品国产清高在天天线| 日韩av在线大香蕉| 国产69精品久久久久777片| 国产精品人妻久久久久久| 亚洲精品日韩在线中文字幕 | 久久精品久久久久久噜噜老黄 | 国产成人一区二区在线| 国产高潮美女av| 国产高清视频在线观看网站| 久久午夜福利片| 又爽又黄无遮挡网站| 91在线观看av| 又粗又爽又猛毛片免费看| 色综合亚洲欧美另类图片| 日韩欧美精品免费久久| 99久国产av精品| 一个人看视频在线观看www免费| 人人妻,人人澡人人爽秒播| 成年女人永久免费观看视频| av女优亚洲男人天堂| 大又大粗又爽又黄少妇毛片口| 国产精品久久久久久久电影| 久久国产乱子免费精品| 舔av片在线| 国产成人aa在线观看| 美女被艹到高潮喷水动态| 国产午夜福利久久久久久| 日日撸夜夜添| 日日摸夜夜添夜夜爱| 国产精品久久久久久久久免| 国产av麻豆久久久久久久| 12—13女人毛片做爰片一| 国产私拍福利视频在线观看| 岛国在线免费视频观看| 99久国产av精品| 欧美不卡视频在线免费观看| 欧美bdsm另类| 国产一区二区三区在线臀色熟女| 国产人妻一区二区三区在| 国产精品一区二区免费欧美| 少妇高潮的动态图| 热99在线观看视频| 欧美成人a在线观看| 精品久久久久久成人av| 久久久久国产网址| eeuss影院久久| 少妇的逼好多水| a级毛色黄片| 赤兔流量卡办理| 午夜影院日韩av| 有码 亚洲区| 亚洲欧美日韩高清专用| 亚洲欧美日韩高清在线视频| av天堂中文字幕网| 男女做爰动态图高潮gif福利片| а√天堂www在线а√下载| 97超级碰碰碰精品色视频在线观看| 中文字幕av在线有码专区| 在线观看av片永久免费下载| 波多野结衣高清作品| 男人的好看免费观看在线视频| 国内久久婷婷六月综合欲色啪| 免费av观看视频| 在线观看av片永久免费下载| 国产一区二区亚洲精品在线观看| 国产人妻一区二区三区在| 黄色欧美视频在线观看| 99视频精品全部免费 在线| 国产精品av视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 成年av动漫网址| 亚洲成人av在线免费| avwww免费| 毛片一级片免费看久久久久| 亚洲精品乱码久久久v下载方式| 国语自产精品视频在线第100页| 亚洲欧美成人综合另类久久久 | 在线国产一区二区在线| 大型黄色视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 麻豆成人午夜福利视频| 十八禁国产超污无遮挡网站| 欧美三级亚洲精品| 亚洲美女黄片视频| 亚洲精品粉嫩美女一区| 能在线免费观看的黄片| 亚洲在线观看片| 大型黄色视频在线免费观看| 国产精品爽爽va在线观看网站| 久久久久九九精品影院| 最近最新中文字幕大全电影3| 小说图片视频综合网站| 欧美日韩精品成人综合77777| 69人妻影院| 国产久久久一区二区三区| 国产精品亚洲一级av第二区| 国产精品一及| 在线观看午夜福利视频| 大又大粗又爽又黄少妇毛片口| 久久欧美精品欧美久久欧美| 亚洲精品国产成人久久av| 一卡2卡三卡四卡精品乱码亚洲| 国产探花极品一区二区| 少妇被粗大猛烈的视频| 卡戴珊不雅视频在线播放| 如何舔出高潮| 国产毛片a区久久久久| 99久久中文字幕三级久久日本| 日本欧美国产在线视频| 国产精品电影一区二区三区| 神马国产精品三级电影在线观看| a级毛片a级免费在线| 欧美3d第一页| 免费观看精品视频网站| 蜜臀久久99精品久久宅男| 欧美一区二区亚洲| 欧美另类亚洲清纯唯美| 国产老妇女一区| 国产高清视频在线观看网站| 亚洲成a人片在线一区二区| 波野结衣二区三区在线| 国产成年人精品一区二区| 国产精品嫩草影院av在线观看| 国产精品人妻久久久久久| 国产男靠女视频免费网站| 国产成人a∨麻豆精品| 国产色婷婷99| 午夜福利18| 色播亚洲综合网| 人妻丰满熟妇av一区二区三区| av在线观看视频网站免费| 日韩av在线大香蕉| 亚洲国产精品久久男人天堂| 欧美日韩一区二区视频在线观看视频在线 | 免费无遮挡裸体视频| 久久久久久久久久黄片| 非洲黑人性xxxx精品又粗又长| 少妇高潮的动态图| 日日撸夜夜添| 夜夜夜夜夜久久久久| av福利片在线观看| 亚洲内射少妇av| 国产69精品久久久久777片| 国产黄片美女视频| 床上黄色一级片| 国产v大片淫在线免费观看| 成人特级黄色片久久久久久久| a级毛片a级免费在线| 免费黄网站久久成人精品| 国产高清视频在线播放一区| 欧美日韩乱码在线| 国产大屁股一区二区在线视频| 大型黄色视频在线免费观看| 搡老妇女老女人老熟妇| 精品久久久久久久久久久久久| 国产成人精品久久久久久| 国产精品免费一区二区三区在线| 美女 人体艺术 gogo| 国产黄片美女视频| 亚洲国产欧美人成| 99热这里只有精品一区| 国产免费一级a男人的天堂| 两个人视频免费观看高清| 亚洲在线自拍视频| 日韩欧美精品v在线| 99久久精品国产国产毛片| 国产伦在线观看视频一区| 国产精品久久久久久精品电影| 国产人妻一区二区三区在| 亚洲欧美日韩东京热| 有码 亚洲区| 精品久久久久久久久久久久久| 午夜福利高清视频| 91久久精品国产一区二区三区| 特大巨黑吊av在线直播| 身体一侧抽搐| 男女视频在线观看网站免费| 可以在线观看毛片的网站| 久久久久久九九精品二区国产| 两个人视频免费观看高清| 亚洲国产精品国产精品| 成年女人看的毛片在线观看| 97碰自拍视频| 国产一区二区在线观看日韩| www.色视频.com| eeuss影院久久| 久久欧美精品欧美久久欧美| 国产三级中文精品| 久久久午夜欧美精品| 全区人妻精品视频| 看免费成人av毛片| 国产av不卡久久| 国产91av在线免费观看| 22中文网久久字幕| av在线播放精品| av专区在线播放| 菩萨蛮人人尽说江南好唐韦庄 | 日韩大尺度精品在线看网址| 精品久久久久久久久亚洲| av女优亚洲男人天堂| 婷婷亚洲欧美| 九九热线精品视视频播放| 国产乱人偷精品视频| 亚洲乱码一区二区免费版| 18禁裸乳无遮挡免费网站照片| 国产成人一区二区在线| 亚洲熟妇中文字幕五十中出| 午夜精品国产一区二区电影 | 99久久精品国产国产毛片| 亚洲,欧美,日韩| 长腿黑丝高跟| 免费黄网站久久成人精品| 婷婷精品国产亚洲av在线| 亚洲四区av| 网址你懂的国产日韩在线| 久久久精品欧美日韩精品| 色噜噜av男人的天堂激情| 一本一本综合久久| 精品国内亚洲2022精品成人| www.色视频.com| 欧美成人精品欧美一级黄| 日韩精品青青久久久久久| 亚洲欧美日韩高清专用| 在线a可以看的网站| 舔av片在线| 国产探花在线观看一区二区| 九九久久精品国产亚洲av麻豆| 夜夜夜夜夜久久久久| 成人毛片a级毛片在线播放| 在现免费观看毛片| 哪里可以看免费的av片| 能在线免费观看的黄片| 国产精品亚洲一级av第二区| 精品免费久久久久久久清纯| a级一级毛片免费在线观看| 91麻豆精品激情在线观看国产| 久久草成人影院| 午夜福利在线在线| 亚洲四区av| 又黄又爽又免费观看的视频| 国产精品爽爽va在线观看网站| 中文字幕久久专区| 久久99热这里只有精品18| 免费看日本二区| 欧美又色又爽又黄视频| 亚洲国产色片| 亚洲第一区二区三区不卡| 一边摸一边抽搐一进一小说| 天天一区二区日本电影三级| 噜噜噜噜噜久久久久久91| 欧美最黄视频在线播放免费| 亚洲av.av天堂| 色av中文字幕| 亚洲成人精品中文字幕电影| 三级毛片av免费| 国产精品一区二区性色av| 最后的刺客免费高清国语| 免费无遮挡裸体视频| 精品免费久久久久久久清纯| 麻豆一二三区av精品| 黑人高潮一二区| 国产色爽女视频免费观看| 男女做爰动态图高潮gif福利片| 极品教师在线视频| 久久精品91蜜桃| 亚洲婷婷狠狠爱综合网| 国产成人freesex在线 | 十八禁国产超污无遮挡网站| 久久精品夜夜夜夜夜久久蜜豆| 晚上一个人看的免费电影| 国产精品福利在线免费观看| 国产人妻一区二区三区在| 亚洲一区二区三区色噜噜| 一边摸一边抽搐一进一小说| 欧美不卡视频在线免费观看| 91在线精品国自产拍蜜月| 欧美激情国产日韩精品一区| 国产一区二区三区在线臀色熟女| 美女被艹到高潮喷水动态| 亚洲自偷自拍三级| 最近最新中文字幕大全电影3| 男女做爰动态图高潮gif福利片| 亚洲av电影不卡..在线观看| 日本色播在线视频| 久久久国产成人免费| 亚洲国产高清在线一区二区三| av卡一久久| 三级经典国产精品| 久久6这里有精品| 国产精品国产高清国产av| 日本a在线网址| 国产精品不卡视频一区二区| 99在线人妻在线中文字幕| 国产精品一区二区免费欧美| 在线免费观看的www视频| 国产综合懂色| 亚洲美女视频黄频| 国产91av在线免费观看| 黄片wwwwww| 国产精品伦人一区二区| 97碰自拍视频| 日韩 亚洲 欧美在线| 亚洲精品日韩在线中文字幕 | 99久久无色码亚洲精品果冻| 亚洲国产日韩欧美精品在线观看| 九九热线精品视视频播放| 国产在线男女| 亚洲精品国产成人久久av| 可以在线观看毛片的网站| АⅤ资源中文在线天堂| 国产成人a区在线观看| 亚洲真实伦在线观看| 久久九九热精品免费| 免费av不卡在线播放| 少妇猛男粗大的猛烈进出视频 | 免费av毛片视频| 亚洲成人久久爱视频| 三级男女做爰猛烈吃奶摸视频| 六月丁香七月| 黄色日韩在线| 性色avwww在线观看| 欧美三级亚洲精品| 午夜精品在线福利| 12—13女人毛片做爰片一| 国产人妻一区二区三区在| 午夜福利高清视频| 最近的中文字幕免费完整| 午夜福利在线在线| 夜夜夜夜夜久久久久| 国产三级在线视频| 99热6这里只有精品| 老熟妇仑乱视频hdxx| 乱系列少妇在线播放| 日韩精品青青久久久久久| aaaaa片日本免费| 村上凉子中文字幕在线| 中出人妻视频一区二区| 日韩中字成人| 欧美高清性xxxxhd video| 日本 av在线| 尤物成人国产欧美一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 久久精品91蜜桃| 久久久久国内视频| 精品久久久噜噜| 亚洲精品色激情综合| 淫妇啪啪啪对白视频| 变态另类丝袜制服| 中文字幕av在线有码专区| a级毛片a级免费在线| 在线观看www视频免费| 九草在线视频观看| 69精品国产乱码久久久| 精品人妻熟女av久视频| 少妇高潮的动态图| 国产高清国产精品国产三级| 日本欧美视频一区| 久久久久精品久久久久真实原创| 久久久久精品性色| 九草在线视频观看| 搡老乐熟女国产| 一级爰片在线观看| 国产精品秋霞免费鲁丝片| 狂野欧美白嫩少妇大欣赏| 久久久久久久久久人人人人人人| 噜噜噜噜噜久久久久久91| 乱码一卡2卡4卡精品| 哪个播放器可以免费观看大片| 黄色视频在线播放观看不卡| 成人亚洲欧美一区二区av| 日本av手机在线免费观看| 80岁老熟妇乱子伦牲交| 五月玫瑰六月丁香| 91成人精品电影| 日本wwww免费看| 晚上一个人看的免费电影| 成人黄色视频免费在线看| 99热全是精品| 亚洲精品乱码久久久久久按摩| 高清欧美精品videossex|