喬飛楊, 王丹璐, 陳 昱, 雷利民, 趙秀閣*, 陶 燕, 杜 偉, 張 婷

1.中國環(huán)境科學(xué)研究院, 環(huán)境基準(zhǔn)與風(fēng)險評估國家重點實驗室， 北京 100012 2.蘭州大學(xué)資源環(huán)境學(xué)院， 甘肅 蘭州 730000 3.中國環(huán)境科學(xué)學(xué)會， 北京 100082

空氣污染會產(chǎn)生不良的健康效應(yīng)進(jìn)而使整個社會遭受經(jīng)濟損失[1]，研究[2-4]表明，人體暴露于被污染的空氣中會增加健康風(fēng)險，嚴(yán)重威脅支氣管、肺和其他器官. 人體暴露空氣污染的健康風(fēng)險主要取決于空氣中污染物濃度水平以及人的環(huán)境暴露行為模式，即暴露參數(shù)，而呼吸量則是開展空氣污染相關(guān)的人群健康風(fēng)險評估的首要參數(shù)[5-6]，是決定健康風(fēng)險評估科學(xué)性和準(zhǔn)確性的關(guān)鍵因子. 呼吸量(Inhalation rate, IR)是指在一定溫度下，人體單位時間內(nèi)呼入空氣的體積[7]. 目前，國內(nèi)外獲得呼吸量的方法主要有直接測量法[8]、雙水標(biāo)記法[9]、心率回歸法[10]和能量代謝法[11]等，其中，直接測量法的結(jié)果準(zhǔn)確但操作繁瑣復(fù)雜；雙水標(biāo)記法準(zhǔn)確度高但測定周期長，且不能反映呼吸量實時性；心率回歸法和能量代謝法因較前兩種方法更簡便、快捷，且適用于大樣本人群呼吸量估算而被廣泛采用.

20世紀(jì)10到20年代，為深入研究人體生理機能，在肺活量的研究中出現(xiàn)了以呼吸量(當(dāng)時被稱為分鐘通氣量)為主要測量參數(shù)的研究方法[12-13]. 20世紀(jì)80年代起，人體呼吸量進(jìn)入了快速研究階段. 1981年，國際輻射防護(hù)委員會(ICRP)通過時間-活動-通氣量調(diào)查方法，收集了全球13 327名個體的身高、體重等生理參數(shù)，并估算了成人、兒童(10歲)、嬰兒(1歲)和新生兒的每分鐘通氣量，獲得了呼吸量參考值，但并未對其有效性和準(zhǔn)確性進(jìn)行說明[14]；1991年，Shamoo等[15]發(fā)現(xiàn)，心率可以作為估算呼吸量的指標(biāo)，提供一種指示呼吸量的新方法，并進(jìn)行了驗證；1993年，Layton[11]提出了一種基于耗氧量(與能量消耗相關(guān))和通氣當(dāng)量(與每分鐘耗氧量體積有關(guān))的估算呼吸量的新方法，即人體能量代謝法，美國環(huán)境保護(hù)局1997年通過此方法對其暴露參數(shù)手冊中呼吸量的推薦值進(jìn)行了更新[16]，這也標(biāo)志著呼吸量的研究進(jìn)入了成熟階段. 進(jìn)入21世紀(jì)，我國學(xué)者對比了能量代謝法發(fā)現(xiàn)，Shizgal-Rosa法更適合估算我國健康人群的基礎(chǔ)代謝率[17-18]；同時，心率回歸法也發(fā)展迅速，基于暴露情景的不同人群呼吸量估算精確度也在不斷提升[19-23]. 目前，各國暴露參數(shù)手冊中呼吸量的估算均采用了能量代謝法，但均未對其適用性和準(zhǔn)確性進(jìn)行評估[24-25]. 筆者根據(jù)2018年6—7月在甘肅省某縣開展的兒童呼吸量調(diào)查研究，基于現(xiàn)場實測獲取的身高、體重、心率等參數(shù)，應(yīng)用7種呼吸量模型估算呼吸量，基于兒童生長發(fā)育規(guī)律對估算結(jié)果進(jìn)行了分析，采用Bland-Altman方法評價了模型的一致性[26-27]，基于呼吸量實測值并結(jié)合平均偏差(MB)、標(biāo)準(zhǔn)化分?jǐn)?shù)偏差(MFB)和標(biāo)準(zhǔn)化分?jǐn)?shù)誤差(MFE)評價了模型的準(zhǔn)確性，利用篩選模型估算了調(diào)查區(qū)縣兒童的呼吸量推薦值，以期為今后因地制宜地開展兒童與空氣相關(guān)健康風(fēng)險評估及管理提供科學(xué)依據(jù).

1 調(diào)查地區(qū)與研究方法

1.1 調(diào)查地區(qū)

調(diào)查縣隸屬于甘肅省慶陽市，介于107°56′20″E～108°38′08″E、35°14′40″N～35°36′18″N之間，屬于黃土高原溝壑區(qū)，黃土層厚90～150 m；地勢自東北向西南傾斜，平均海拔1 460 m. 該縣氣候受季風(fēng)影響明顯，為溫帶大陸性季風(fēng)半濕潤氣候. 截至2017年末，全縣共有各類學(xué)校143所，其中小學(xué)86所、初中8所；小學(xué)在校兒童 15 872 人，初中 4 849 人[28].

1.2 調(diào)查方法

1.2.1抽樣方法

根據(jù)最小樣本量計算公式〔見式(1)〕[7]，每層需要最小樣本量118人.

(1)

式中：n為每層最小樣本量，人；Uα/2為顯著性水平為95%相應(yīng)的標(biāo)準(zhǔn)正態(tài)差，取1.96；σ為體重的標(biāo)準(zhǔn)差，根據(jù)預(yù)調(diào)查結(jié)果，調(diào)查兒童的體重為(39.42±12.58)kg，σ取12.58 kg；δ為允許誤差，問卷調(diào)查法中δ一般取10%；μ為體重的算數(shù)平均值，取39.42 kg；deff為設(shè)計效應(yīng)值，取3.0.

調(diào)查以城鄉(xiāng)(城市、農(nóng)村)、性別(男、女)和年齡(6～<9歲、9～<12歲、12～<15歲)為主要分層因素進(jìn)行分層抽樣，根據(jù)總樣本量計算公式〔見式(2)〕，需樣本 1 574 個，此次實際調(diào)查人數(shù)為 3 303 人，滿足調(diào)查樣本量的要求.

N=(n×q)/(1-p)

(2)

式中：N為總樣本量，人；p為失訪率，取10%；q為分層因素的乘積，為12層(2×2×3).

第一階段，抽取調(diào)查學(xué)校. 抽取學(xué)校的教學(xué)質(zhì)量及規(guī)模等均為當(dāng)?shù)刂械人?，此次共隨機抽取學(xué)校10所，其中，獨立小學(xué)6所，城市和農(nóng)村分別為2所和4所；獨立中學(xué)3所，城市和農(nóng)村分別為1所和2所，1所九年制學(xué)校(即包含小學(xué)和初中)(見圖1).

圖1 調(diào)查縣調(diào)查學(xué)校分布Fig.1 The geographical location of the distribution of schools

第二階段，抽取調(diào)查班級. 抽取班級的教學(xué)質(zhì)量盡量為學(xué)校的中等水平，在所選擇學(xué)校的每個年級隨機抽取2個班(排除重點班或最差班).

第三階段,抽取調(diào)查對象. 抽取所選班級的所有兒童.

1.2.2現(xiàn)場測量

采用身高體重計(江蘇蘇宏醫(yī)療器械有限公司，RGT-120/160)根據(jù)《國民體質(zhì)測定標(biāo)準(zhǔn)手冊》[29]中身高和體重測量方法對兒童的身高和體重進(jìn)行測量，保留小數(shù)點后2位.

采用夾式脈搏血氧儀(Beurer GmbH，德國，PO30型)測定受試者心率. 由于部分受試兒童年齡較小，手指發(fā)育尚未達(dá)到儀器的光感檢測條件，從而無法顯示讀數(shù)，故對此類受試者采用人工計數(shù)方式. 根據(jù)崔焱[30]提供的測量方法，測量前待受試者安靜休息片刻后，以坐姿等候測量. 測量者用右手食指、中指和無名指的指端，按在被測者的橈動脈表面，清晰觸到動脈波動時開始計時，計錄60 s (以秒表計時)內(nèi)脈搏波動次數(shù). 由于此次受試者均為健康人群，脈搏與心率一致，故采用測量脈搏的方式測定心率.

采用肺通氣量儀(青島精誠儀器儀表有限公司，F(xiàn)T-1型)對受試者的坐姿呼吸量進(jìn)行測量. 測量前將傳感器插入面具吸氣口，再將面具戴在口鼻處，面具的邊緣貼緊顏面. 打開儀器，受試者正常呼吸5 min后，關(guān)閉儀器并記錄平均呼吸量(單位為L/min).

1.3 呼吸量模型估算

采用7種呼吸量估算模型，并結(jié)合不同活動強度水平下的能量消耗量與基礎(chǔ)代謝率比值對甘肅省某縣兒童的呼吸量進(jìn)行估算(見表1). 其中，能量代謝法利用呼吸量與能量消耗率和空氣攝取率成正比的關(guān)系計算受試者的呼吸量，Shizgal-Rosa法的原理與能量代謝法一致，但每日能量消耗的計算方法不同；心率回歸法(包括Greenwald法、Zuurbier法、Cozza法、Ramos法和Vale法)可以通過直接測量心率并結(jié)合回歸模型來估算呼吸量，由于缺少兒童心率回歸估算模型，因此均采用已有的成人心率回歸模型，其中Greenwald法采用肺活量與心率估算呼吸量，肺活量的計算采用文獻(xiàn)[32-33]報道方法.

表1 用于估算兒童呼吸量的數(shù)學(xué)模型

1.4 統(tǒng)計分析

1.4.1數(shù)據(jù)錄入與清洗

采用EpiData軟件進(jìn)行雙錄入. 若錄入結(jié)果一致，則直接入庫；錄入結(jié)果不一致則將該樣本數(shù)據(jù)轉(zhuǎn)到質(zhì)控人員進(jìn)行二次錄入. 采用SPSS 14.0軟件，以性別、年齡和城鄉(xiāng)等為關(guān)鍵變量，對每個變量的缺失值進(jìn)行標(biāo)記和分析，剔除邏輯錯誤和非法值.

1.4.2模型的一致性檢驗

采用SPSS 14.0軟件對樣本數(shù)據(jù)進(jìn)行正態(tài)性檢驗，用MedCalc 15.1軟件對7種模型進(jìn)行Bland-Altman分析，評價其一致性.

1.4.3呼吸量的評估與估算

將篩選的最適模型作為呼吸量估算方法，采用Crystal Ball軟件進(jìn)行蒙特卡洛模擬估算調(diào)查區(qū)縣兒童呼吸量. 基于兒童呼吸量實測值，利用平均偏差(MB)、標(biāo)準(zhǔn)化分?jǐn)?shù)偏差(MFB)和標(biāo)準(zhǔn)化分?jǐn)?shù)誤差(MFE)評價各模型，這3種參數(shù)可通過量化估算值與實測值的差異來反映模型評估的可靠性，參數(shù)值越接近0，說明模型越可靠，計算公式：

(3)

(4)

(5)

式中，Ci為估算值，C0為測量值，K為參與對比的有效樣本數(shù).

1.5 質(zhì)量控制

調(diào)查質(zhì)量控制貫穿于方案修訂、技術(shù)培訓(xùn)、物資

準(zhǔn)備和預(yù)調(diào)查以及現(xiàn)場調(diào)查和數(shù)據(jù)處理分析等全過程. 現(xiàn)場調(diào)查過程中，采用調(diào)查員自查、質(zhì)量控制員復(fù)查等措施，保證漏項率、邏輯錯誤率和測量不清率低于15%，其中身高、體重由兩名測量員完成，質(zhì)控員抽取一定比例的調(diào)查對象進(jìn)行復(fù)核，復(fù)核一致率高于95%. 采用SPSS 14.0軟件對缺失和異常數(shù)據(jù)進(jìn)行清洗.

1.6 質(zhì)量評價

1.6.1問卷有效性

此次共調(diào)查在校兒童 3 303 人，全部應(yīng)答，調(diào)查兒童應(yīng)答率100%；以性別、年齡和城鄉(xiāng)等為關(guān)鍵變量進(jìn)行數(shù)據(jù)清理，最終用于分析的有效樣本數(shù) 2 866 份，問卷有效率86.77%；以性別、年齡、身高、體重、心率等關(guān)鍵變量進(jìn)行復(fù)核，關(guān)鍵變量應(yīng)答率均為100%.

1.6.2問卷可靠性

抽取1 045人對問卷中的20個問題進(jìn)行復(fù)核詢問，結(jié)果一致率為99.9%.

2 結(jié)果與討論

2.1 調(diào)查兒童基本特征

2.1.1樣本分布

調(diào)查兒童的樣本分布如表2所示. 由表2可見：共獲得 2 866 個有效樣本，其中城市 1 327 人，農(nóng)村 1 539 人，分別占46.30%和53.70%；男性和女性分別有 1 501 和 1 365 人，各占52.37%和47.63%. 從年齡分布看，6～<9歲、9～<12歲和12～<15歲分別有449、1 106 和 1 311 人，占比分別為15.67%、38.59%和45.74%.

表2 調(diào)查兒童的樣本分布

2.1.2基本生理特征

調(diào)查兒童生理參數(shù)的描述性統(tǒng)計如表3所示. 由表3可見：從身高來看，隨著年齡的增長，兒童身高逐漸增長，同年齡段城市兒童身高均高于農(nóng)村. 城市6～<9歲、9～<12歲和12～<15歲男童的平均身高分別為141.86、145.31和158.31 cm，女童分別為142.57、144.64和154.77 cm；農(nóng)村6～<9歲、9～<12歲和12～<15歲男童的平均身高分別為129.38、142.86和153.26 cm，女童分別為128.21、142.64和151.04 cm. 從體重來看，隨著年齡的增長，體重逐漸增加，同年齡段男童體重均大于或與女童持平. 城市6～<9歲、9～<12歲和12～<15歲男童的平均體重分別為35.64、38.51和46.77 kg，女童分別為35.90、37.10與45.18 kg；農(nóng)村6～<9歲、9～<12歲和12～<15歲男童的平均體重分別為28.79、35.33和43.80 kg，女童分別為28.67、35.01和42.53 kg. 從心率來看，城市6～<9歲、9～<12歲和12～<15歲男童平均心率分別為87、86和84次/min，女童分別為87、85和86次/min；農(nóng)村6～<9歲、9～<12歲和12～<15歲男童平均心率分別為87、82和81次/min,女童分別為88、84和81次/min. 與衣風(fēng)等[34-35]研究結(jié)果一致，兒童的心率基本隨年齡的增長而逐漸下降，這是因為心臟收縮壓隨著年齡增長不斷提高，而每搏輸出量絕對值隨之增加，因此不斷發(fā)育的兒童心臟收縮力逐漸增強，心率也隨之逐漸下降.

表3 調(diào)查地區(qū)兒童基本生理參數(shù)描述性統(tǒng)計

2.2 呼吸量估算模型兒童適用性評估

注：縱坐標(biāo)表示兩種方法呼吸量計算值的比值，橫坐標(biāo)表示兩種方法呼吸量計算值的平均值.圖2 7種估算模型的Bland-Altman分析散點圖Fig.2 Bland-Altman plots of 7 models for estimating inhalation rates

2.2.1一致性分析

經(jīng)正態(tài)和方差齊性檢驗，7種估算模型的計算結(jié)果數(shù)據(jù)行為良好，均滿足Bland-Altman分析條件. 分別以Shizgal-Rosa法的呼吸量計算值與其他6種模型呼吸量計算值的平均值為橫坐標(biāo)，二者的比值為縱坐標(biāo)，繪制Bland-Altman分析散點圖，并統(tǒng)計各模型95%一致性界限范圍(見圖2). 各模型落在95%一致性界限外的樣本數(shù)均小于5%，表明Shizgal-Rosa法的計算值與其他模型計算值的偏差均較小. 在Bland-Altman分析中，95%一致性界限范圍越小，一致性越好. 由圖2 可見：Shizgal-Rosa法與能量代謝法的95%一致性界限范圍(0.57～0.75)最小；在心率回歸法中，Ramos法與Shizgal-Rosa法的95%一致性界限范圍(0.18～0.73)最小，Vale法的范圍(0.25～1.17)最大. 相較于能量代謝法，Shizgal-Rosa法與心率回歸法對兒童呼吸量的估算一致性較差.

2.2.2基于實測值的估算模型評估

為了全面評估估算模型的可靠性，按照城鄉(xiāng)、性別和年齡等分層因素，從此次調(diào)查兒童中隨機抽取799人(其中，男性411人，女性388人)進(jìn)行呼吸量實測，依據(jù)兒童呼吸實際測量結(jié)果對模型估算結(jié)果進(jìn)行評估(見表4). 由表4可見：Shizgal-Rosa法的MB、MFB和MFE均為最小，且與能量代謝法的MB、MFB和MFE的差值分別為3.21 L/min、0.30和0.19；心率回歸法中，Ramos法、Cozza法和Zuurbier法估算的呼吸量在全年齡段均遠(yuǎn)高于實測值，Vale法與Greenwald法的估算結(jié)果在心率回歸法中與實測值較為接近，說明這兩種方法穩(wěn)定性相對較好，但Vale法僅用心率作為唯一自變量來建立模型，可能并不完全適用于兒童呼吸量的估算. 一方面是由于單個變量具有偶然性，若發(fā)生較大波動會嚴(yán)重影響模型結(jié)果；另一方面，由于兒童各方面機能發(fā)育并不完善，其心率不夠穩(wěn)定[35]，亦會對模型的穩(wěn)定性及結(jié)果的準(zhǔn)確性產(chǎn)生影響. 整體來看，以能量代謝為基礎(chǔ)的模型估算表現(xiàn)要好于心率回歸法. Shizgal-Rosa法估算的呼吸量與實測值最為相近，其MB、MFB和MFE分別為0.08 L/min、0.03和0.14，在估算結(jié)果中均為最優(yōu).

表4 呼吸量估算模型計算結(jié)果評估情況

2.3 不同呼吸量估算模型計算結(jié)果分析

基于不同呼吸量估算模型的兒童呼吸量估算結(jié)果如表5所示. 由表5可見，各模型的呼吸量估算結(jié)果差異明顯，采用Shizgal-Rosa法估算男童和女童的呼吸量范圍分別為4.79～6.35和5.08～6.09 L/min，能量代謝法為7.85～10.07和7.24～8.79 L/min，Vale法為10.09～12.00和10.20～12.14 L/min，Ramos法為16.78～19.36和14.17～16.33 L/min，Cozza法為14.22～17.07和14.45～17.29 L/min，Zuurbier法為15.88～18.46和12.04～14.19 L/min，Greenwald法為8.76～15.34和8.22～13.44 L/min. 從城鄉(xiāng)分層看，不同估算模型對于城市和農(nóng)村兒童呼吸量估算結(jié)果相差不大. 從性別方面來看，各模型估算的同年齡段男性的呼吸量高于女性. Shizgal-Rosa法的估算結(jié)果顯示，無論是城鄉(xiāng)還是男女，呼吸量均隨著年齡的增長而增加，能量代謝法估算結(jié)果規(guī)律與其基本一致，不同分層之間呼吸量差異較為明顯. 有研究[36]表明，呼吸量基本隨體重的升高而增加，即正常發(fā)育的兒童隨著年齡的增長，其呼吸量不斷增加. 而Zuurbier法、Cozza法、Ramos法和Vale法估算結(jié)果均呈6～9歲年齡段兒童的呼吸量大于其他年齡段的特征，違背了兒童的生長發(fā)育規(guī)律. 進(jìn)一步對不同估算模型的估算值分布進(jìn)行對比(見圖3)，發(fā)現(xiàn)這些方法較其他方法估算結(jié)果范圍大且離群值多. 由于缺少針對兒童的心率呼吸量估算模型，故選用成人心率回歸法進(jìn)行計算，且僅有一個自變量，以上均給估算結(jié)果帶來了較大的不確定性. Greenwald法估算結(jié)果呈現(xiàn)的規(guī)律與能量代謝法基本一致，但其估算結(jié)果總體高于能量代謝法，即高于實測值，雖然該方法在心率作為自變量的基礎(chǔ)上又加入了身高、體重、年齡3個自變量，但仍有待進(jìn)一步驗證. 綜上，Shizgal-Rosa法更適合我國兒童，尤其是筆者調(diào)查縣兒童呼吸量的估算.

圖3 基于不同估算模型的不同分層兒童呼吸量箱形圖Fig.3 Box Figure of inhalation rates of children with different stratification based on different estimation models

表5 基于調(diào)查兒童生理參數(shù)不同估算模型的計算結(jié)果

2.4 研究地區(qū)兒童呼吸量推薦值估算

基于調(diào)查縣兒童呼吸量估算模型適用性、一致性和可靠性評價結(jié)果，利用分析模型中表現(xiàn)最好的Shizgal-Rosa法，采用Crystal Ball進(jìn)行蒙特卡洛模擬(循環(huán) 10 000 次)調(diào)查縣兒童休息狀態(tài)下的短期呼吸量推薦值，依據(jù)HJ 877—2017《暴露參數(shù)調(diào)查技術(shù)規(guī)范》[7]獲得了調(diào)查縣不同年齡段兒童的長期呼吸量及其他狀態(tài)下的短期呼吸量推薦值(見表6).

由表6可見：調(diào)查縣6～<9歲、9～<12歲和12～<15歲兒童的長期呼吸量估算值分別為9.52、12.45和13.52 m3/d，其中，男童分別為9.52、12.35和14.54 m3/d，女童為10.04、12.46和11.99 m3/d；調(diào)查縣6～<9歲、9～<12歲和12～<15歲兒童休息狀態(tài)下的短期呼吸量估算值分別為3.89、4.55和5.52 L/min，靜坐狀態(tài)下分別為4.67、5.46和6.62 L/min，輕度運動狀態(tài)下分別為7.78、9.10和11.04 L/min，中度運動狀態(tài)下分別為15.56、18.20和22.08 L/min，劇烈運動狀態(tài)下分別為38.90、45.50和55.20 L/min. 與同年齡段我國西北地區(qū)兒童休息狀態(tài)下短期呼吸量推薦值[24]比較發(fā)現(xiàn)：該縣男童休息狀態(tài)下的短期呼吸量(6～<9歲、9～<12歲和12～<15歲分別為3.88、4.52和5.62 L/min)略低于推薦值(4.30、4.90和5.70 L/min)，而女童則略高. 其他狀態(tài)下的短期呼吸量和長期呼吸量的估算值也與推薦值存在類似差異. 造成這種差異的原因：一方面，此次調(diào)查是在我國兒童暴露參數(shù)發(fā)布后約5 a開展的，有研究[37-38]表明兒童的生長發(fā)育水平與我國社會經(jīng)濟發(fā)展和兒童營養(yǎng)狀況等密切關(guān)系；另一方面，此次調(diào)查只選取了一個示范縣，不能代表整個西北地區(qū)兒童，估算值可能與推薦值存在一定差異.

年齡段∕歲性別短期呼吸量∕(L∕min)休息靜坐輕度運動中度運動劇烈運動估算值推薦值1)估算值推薦值1)估算值推薦值1)估算值推薦值1)估算值推薦值1)長期呼吸量∕(m3∕d)估算值推薦值1)6～<99～<1212～<15男3.88 4.30 4.66 5.20 7.77 8.60 15.53 17.20 38.83 43.00 9.52 10.50女4.10 3.80 4.92 4.60 8.20 7.60 16.40 15.20 41.00 38.00 10.04 9.30 小計3.89 4.10 4.67 4.90 7.78 8.20 15.56 16.30 38.90 40.80 9.52 10.00 男4.52 4.90 5.42 5.90 9.03 9.90 18.07 19.70 45.17 49.40 12.35 13.50 女4.55 4.50 5.46 5.40 9.10 9.00 18.20 17.90 45.50 44.80 12.46 12.20 小計4.55 4.70 5.46 5.70 9.10 9.50 18.20 18.90 45.50 47.30 12.45 12.90男5.62 5.70 6.74 6.80 11.23 11.30 22.47 22.60 56.17 56.50 14.54 14.60 女5.21 4.90 6.25 5.80 10.42 9.70 20.83 19.50 52.08 48.70 11.99 11.20 小計5.52 5.30 6.62 6.40 11.04 10.60 22.08 21.20 55.20 53.10 13.52 13.10

3 結(jié)論與建議

3.1 結(jié)論

a) 各估算模型呼吸量估算值差異顯著，能量代謝法與Shizgal-Rosa法估算結(jié)果均符合兒童生長發(fā)育規(guī)律；心率回歸法6～<9歲兒童的呼吸量估算結(jié)果大于其他年齡段，不適合調(diào)查縣兒童呼吸量估算.

b) Bland-Altman分析結(jié)果表明，Shizgal-Rosa法與心率回歸法的一致性較差；基于實測值的評估模型可靠性評估顯示，Shizgal-Rosa法估算結(jié)果較為準(zhǔn)確，為所選模型中最優(yōu)，最適合該研究調(diào)查縣兒童的呼吸量估算.

c) 基于Shizgal-Rosa法，采用蒙特卡洛模擬調(diào)查縣各年齡段兒童呼吸量，結(jié)果顯示該縣男童長期呼吸量與不同狀態(tài)下的短期呼吸量略低于西北地區(qū)推薦值，女童略高于西北地區(qū)推薦值.

3.2 建議

a) 我國兒童呼吸量估算模型建立迫在眉睫. 當(dāng)前我國兒童呼吸量暴露參數(shù)的推薦值采用美國的能量代謝法估算，但結(jié)果表明該方法并不完全適合我國兒童，估算結(jié)果與實測值存在一定差異. 因此，需要開展兒童呼吸量研究，基于兒童呼吸量實測數(shù)據(jù)構(gòu)建適合我國兒童的呼吸量與日常體力活動及各項生理參數(shù)的估算模型，提高兒童空氣相關(guān)健康風(fēng)險評估結(jié)果的科學(xué)性和準(zhǔn)確性.

b) 我國兒童呼吸量推薦值亟待更新. 基于我國兒童呼吸量估算模型對暴露參數(shù)手冊中兒童呼吸量推薦值進(jìn)行更新，能夠更為準(zhǔn)確地反映我國兒童在不同運動強度下活動模式的變化特征，并為今后我國開展兒童環(huán)境暴露及健康風(fēng)險研究打下基礎(chǔ).

參考文獻(xiàn)(References)：

[1] 李惠娟,周德群,魏永杰.空氣污染的健康經(jīng)濟損失評價研究進(jìn)展[J].環(huán)境科學(xué)研究,2020,33(10):2421-2429.

LI Huijuan,ZHOU Dequn,WEI Yongjie.Review on health economic loss assessment of air pollution[J].Research of Environmental Sciences,2020,33(10):2421-2429.

[2] GALLAGHER P,LAZARUS W,SHAPOURI H,etal.Cardiovascular disease-risk benefits of clean fuel technology and policy:a statistical analysis[J].Energy Policy,2010,38(2):1210-1222.

[3] SHI Liuhua,ZANOBETTI A,KLOOG I,etal.Low-concentration PM2.5and mortality:estimating acute and chronic effects in a population-based study[J].Environmental Health Perspectives,2016,124(1):46-52.

[4] ROVIRA J,DOMINGO J L,SCHUHMACHER M.Air quality,health impacts and burden of disease due to air pollution (PM10,PM2.5,NO2and O3):application of AirQ+ model to the Camp de Tarragona County (Catalonia,Spain)[J].Science of the Total Environment,2020,703:135538.

[5] ZHOU Maigeng,LIU Yunning,WANG Lijun,etal.Particulate air pollution and mortality in a cohort of Chinese men[J].Environmental Pollution,2014,186:1-6.

[6] SHANG Yu,SUN Zhiwei,CAO Junji,etal.Systematic review of Chinese studies of short-term exposure to air pollution and daily mortality[J].Environment International,2013,54:100-111.

[7] 環(huán)境保護(hù)部. HJ 877—2017暴露參數(shù)調(diào)查技術(shù)規(guī)范[S].北京:中國環(huán)境科學(xué)出版社,2017.

[8] SHEPHARD R J.A critical examination of the Douglas bag technique[J].Journal of Physiology,1955,127(3):515-524.

[9] LIFSON N,GORDON G B,MCCLINTOCK R,etal.Measurement of total carbon dioxide production by means of D218O[J].Obesity Research,1997,5(1):78-84.

[10] SHAMOO D A,TRIM S C,LITTLE D E,etal.Improved quantitation of air pollution dose rates by improved estimation of ventilation rate[M].Total exposure assessment methodology:a new horizon,Pittsburgh PA:Air and Waste Management Association,1990:553-564.

[11] LAYTON D W.Metabolically consistent breathing rates for use in dose assessments[J].Health Physics,1993,64(1):23-36.

[12] KROGH A,LINDHARD J.The regulation of respiration and circulation during the initial stages of muscular work[J].Journal of Physiology,1913,47(1/2):112-136.

[13] STURGIS C C,PEABODY F W,HALL F C,etal.Clinical studies on the respiration:Ⅷ.the relation of syspnea to the maximum minute-volume of pulmonary ventilation[J].Archives of Internal Medicine,1922(2):236-244.

[14] ICRP (International Commission on Radiological Protection).Report of the task group on reference man[M].New York:Pergammon Press,1981.

[15] SHAMOO D A,JOHNSON T R,TRIM S C,etal.Activity patterns in a panel of outdoor workers exposed to oxidant pollution[J].Journal of Exposure Analysis & Environmental Epidemiology,1991,1(4):423-438.

[16] US EPA.Exposure factors handbook[R].Washington DC:Office of Research and Development,1997.

[17] 梁潔,蔣卓勤,何玉敏,等.中國健康成人基礎(chǔ)代謝率估算公式的探討[J].中國校醫(yī),2008,22(4):372-374.

LIANG Jie,JIANG Zhuoqin,HE Yumin,etal.Estimated equations of basal metabolic rate in Chinese healthy adults[J].Chinese Journal of School Doctor,2008,22(4):372-374.

[18] 董靜梅,夏麗,俞益,等.青少年體力活動中環(huán)境健康風(fēng)險評估的呼吸暴露參數(shù)研究進(jìn)展[J].中國運動醫(yī)學(xué)雜志,2016,35(7):664-668.

[19] GREENWALD R,HAYAT M J,BARTON J,etal.A novel method for quantifying the inhaled dose of air pollutants based on heart rate,breathing rate and forced vital capacity[J].Plos One,2016,11(1):e0147578.

[20] ZUURBIER M,HOEK G,HAZEL P V D,etal.Minute ventilation of cyclists,car and bus passengers:an experimental study[J].Environmental Health,2009,8:1-10.

[21] COZZA I C,ZANETTA D M T,FERNANDES F L A,etal.An approach to using heart rate monitoring to estimate the ventilation and load of air pollution exposure[J].Science of the Total Environment,2015,520(2):160-167.

[22] RAMOS C A,REIS J F,ALMEIDA T,etal.Estimating the inhaled dose of pollutants during indoor physical activity[J].Science of the Total Environment,2015,527/528:111-118.

[23] VALE I D D,VASCONCELOS A S,DUARTE G O.Inhalation of particulate matter in three different routes for the same OD pair:a case study with pedestrians in the city of Lisbon[J].Journal of Transport & Health,2015,2(4):474-482.

[24] 環(huán)境保護(hù)部.中國人群暴露參數(shù)手冊(兒童卷:6～17歲)[M].北京:中國環(huán)境出版社,2016.

[25] US EPA.Child-specific exposure factors handbook[R].Washington DC:US EPA,2008.

[26] 陳卉.Bland-Altman分析在臨床測量方法一致性評價中的應(yīng)用[J].中國衛(wèi)生統(tǒng)計,2007,24(3):308-309.

[27] BLAND J M,ALTMAN D G.Statistical methods for assessing agreement between two methods of clinical measurement[J].Lancet,1986,327(8476):307-310.

[28] 正寧縣統(tǒng)計局.2018年正寧縣國民經(jīng)濟和社會發(fā)展統(tǒng)計公報[R].正寧:正寧縣統(tǒng)計局,2019:11-12.

[29] 國家體育總局.國民體質(zhì)測定標(biāo)準(zhǔn)手冊[M].北京:人民體育出版社,2003:8-9.

[30] 崔焱.護(hù)理學(xué)基礎(chǔ)[M].北京:人民衛(wèi)生出版社,2002:209.

[31] 任建安,李寧,黎介壽.能量代謝監(jiān)測與營養(yǎng)物質(zhì)需要量[J].中國實用外科雜志,2001,21(10):631-637.

[32] 洪建國,王乃禮,李臻,等.上海地區(qū)學(xué)齡兒童肺功能正常值及其預(yù)計值公式[J].臨床兒科雜志,1995,13(1):54-56.

[33] 馮逵,陳莉,韓少梅,等.黑龍江省漢族兒童青少年肺通氣功能正常值和預(yù)計方程式[J].基礎(chǔ)醫(yī)學(xué)與臨床,2010,30(5):480-486.

FENG Kui,CHEN Li,HAN Shaomei,etal.Lung function reference values and prediction equations in children and adolescents of the Han nationality in Heilongjiang Province[J].Basic & Clinical Medicine,2010,30(5):480-486.

[34] 衣風(fēng),任廣峰.健康兒童動態(tài)心電圖竇性心率統(tǒng)計分析[J].中國實用醫(yī)藥,2011,6(28):94-95.

[35] 李敏,王健怡,劉萬紅,等.0～14歲正常兒童心率變異性分析[J].臨床心血管病雜志,2011,27(5):374-376.

LI Min,WANG Jianyi,LIU Wanhong,etal.Analysis of heart rate variability of 0-14 years old healthy children[J].Journal of Clinical Cardiology,2011,27(5):374-376.

[36] LAYTON D W.Simulation of average energy expenditures and related inhalation rates for the U.S. population[J].Air Quality,Atmosphere & Health,2013,6(2):527-539.

[37] SHEN Tiefu,HABICHT J P,CHANG Ying.Effect of economic reforms on child growth in urban and rural areas of China[J].The New England Journal of Medicine,1996,335(6):400-406.

[38] ZONG Xinnan,LI Hui.Physical growth of children and adolescents in China over the past 35 years[J].Bulletin of the World Health Organization,2014,92(8):555-564.