游泳池水CYA/FC比值与ORP对HPC超标的筛查价值:基于GEE的聚集效应校正分析

彭飞, 李方茹, 刘军

安徽预防医学杂志 ›› 2026, Vol. 32 ›› Issue (3) : 252-257.

PDF(948 KB)
PDF(948 KB)
安徽预防医学杂志 ›› 2026, Vol. 32 ›› Issue (3) : 252-257. DOI: 10.19837/j.cnki.ahyf.2026.03.015
卫生监督与监测

游泳池水CYA/FC比值与ORP对HPC超标的筛查价值:基于GEE的聚集效应校正分析

作者信息 +

Screening utility of the cyanuric acid/free chlorine ratio and oxidation-reduction potential for heterotrophic plate count exceedance in swimming pools: a clustered data analysis using generalized estimating equations

Author information +
文章历史 +

摘要

目的 评估亳州市游泳池水质监督监测数据的聚集效应,并在广义估计方程(generalized estimating equations, GEE)校正后,探索氰尿酸/游离性余氯(cyanuric acid/free chlorine, CYA/FC)比值及氧化还原电位(oxidation-reduction potential, ORP)对菌落总数(heterotrophic plate count, HPC)超标的筛查应用价值。方法 采集2024年7月及2025年6—7月亳州市游泳场所水样。采用混合线性模型估计组内相关系数(intraclass correlation coefficient, ICC),采用GEE开展关联分析,以场所为单位的整群Bootstrap法估计比值比(odds ratio, OR)及受试者工作特征曲线(receiver operating characteristic curve, ROC)曲线下面积(area under curve, AUC)的95%置信区间(confidence interval, CI)。结果 共在40家游泳场所采集水样231份。各指标ICC在0.191~0.780之间(P均<0.05),HPC ICC为0.297,提示存在显著聚集效应。FC合格率为21.65%,其中超过上限者占66.67%。GEE分析显示,CYA/FC>20与HPC超标呈正相关(OR=2.690,95%CI:1.348~10.283),并在高FC组(FC>1.0 mg/L)展现出最佳筛查效能(AUC=0.918,阴性预测值=0.974)。ORP在不同切点下的95%CI均包含1,未显示统计学显著性。结论 亳州市游泳池水质监测数据存在显著聚集效应,小簇样本条件下结合整群Bootstrap法可提高统计推断的稳健性。CYA/FC比值>20与HPC超标风险相关,尤其在高FC场所具有较好的排除式筛查潜力;ORP未显示出具有统计学显著性的筛查价值。上述发现需在更大样本中进一步验证。

Abstract

Objective To assess clustering effects in routine swimming pool water quality surveillance data from Bozhou, and to evaluate the screening utility of the cyanuric acid/free chlorine (CYA/FC) ratio and oxidation-reduction potential (ORP) for heterotrophic plate count (HPC) exceedance after adjustment for clustering using generalized estimating equations (GEE). Methods The water samples were collected from swimming facilities in Bozhou in July 2024 and June-July 2025. Intraclass correlation coefficients (ICCs) were estimated using a linear mixed model. Associations of the CYA/FC ratio and ORP with HPC exceedance were assessed using GEE, and cluster-level bootstrap (1 000 replications) was applied to estimate 95% confidence interval (CI) for odds ratio (OR) and area under the receiver operating characteristic curve (AUC). Results A total of 231 water samples were collected from 40 swimming facilities. ICCs ranged from 0.191 to 0.780 (all P<0.05), with an ICC of 0.297 for HPC, indicating notable clustering. FC compliance rate was only 21.65%, with 66.67% exceeding the upper limit. GEE analysis showed that a CYA/FC ratio >20 was positively associated with HPC exceedance (OR=2.690, 95%CI: 1.348-10.283) and demonstrated the best screening performance in the high-FC group (>1.0 mg/L) (AUC=0.918, NPV=0.974). Cluster-bootstrap 95%CI for ORP crossed 1.0 at both cutoffs, indicating no statistical significance. Conclusion Swimming pool surveillance data exhibited notable clustering; cluster-level bootstrap may enhance inference robustness under small-cluster conditions. A CYA/FC ratio >20 was associated with HPC exceedance risk and showed promising rule-out screening potential, particularly in the high-FC group. ORP did not demonstrate statistically significant screening utility. These findings warrant validation in studies with larger sample sizes.

关键词

游泳池 / 水质监测 / 聚集效应 / 氰尿酸 / 游离性余氯 / 氧化还原电位 / 广义估计方程

Key words

Swimming pool / Water quality monitoring / Clustering effect / Cyanuric acid / Free chlorine / Oxidation-reduction potential / Generalized estimating equations

引用本文

导出引用
彭飞, 李方茹, 刘军. 游泳池水CYA/FC比值与ORP对HPC超标的筛查价值:基于GEE的聚集效应校正分析[J]. 安徽预防医学杂志. 2026, 32(3): 252-257 https://doi.org/10.19837/j.cnki.ahyf.2026.03.015
PENG Fei, LI Fangru, LIU Jun. Screening utility of the cyanuric acid/free chlorine ratio and oxidation-reduction potential for heterotrophic plate count exceedance in swimming pools: a clustered data analysis using generalized estimating equations[J]. Anhui Journal of Preventive Medicine. 2026, 32(3): 252-257 https://doi.org/10.19837/j.cnki.ahyf.2026.03.015
中图分类号: R123.1   

参考文献

[1]
陈凤灵, 刘慧, 潘雪梅, 等. 广东省一起因游泳引起的咽结膜热暴发调查[J]. 疾病监测, 2018, 33(8):703-706.
[2]
韩旭, 李莉, 宋风景, 等. 基于风险矩阵和排序赋值法的公共场所环境健康风险评估[J]. 环境卫生学杂志, 2025, 15(10):885-893.
[3]
孙波, 姚孝元, 潘力军, 等. 《公共场所卫生指标及限值要求》(GB 37488—2019) 标准解读[J]. 环境卫生学杂志, 2021, 11(3):268-274.
[4]
于洋, 农惠婷, 梁林涵, 等. 2021年广西六城市公共场所微生物监测结果[J]. 环境卫生学杂志, 2022, 12(12):897-902.
[5]
李艳艳, 魏昊, 周燕, 等. 2016—2021年马鞍山市八类公共场所卫生学监测结果分析[J]. 环境卫生学杂志, 2023, 13(7):516-523.
[6]
梁晓军, 张建新, 杨海兵. 游泳池水中氰尿酸的研究进展[J]. 环境卫生学杂志, 2017, 7(4):315-320.
[7]
梁晓军, 王媛, 胡斌, 等. 三氯异氰尿酸消毒游泳池水中消毒副产物监测[J]. 中国消毒学杂志, 2022, 39(8):577-580.
目的 观察使用三氯异氰尿酸消毒游泳场所时主要消毒副产物水平。方法 通过现场采样和化学分析方法,对昆山市区1家使用三氯异氰尿酸消毒的游泳场所水中消毒副产物进行监测。结果 经三氯异氰尿酸消毒的游泳池水pH值、游离性余氯、化合性余氯和尿素的合格率分别为49.1%、97.2%、54.7%和97.2%。消毒后的游泳池水浑浊度、游离性余氯、化合性余氯、总余氯指标值均高于自来水。消毒后室内外泳池水中11种消毒副产物均可检出,其中三氯甲烷、三氯乙酸和二氯乙酸浓度最高。结论 三氯异氰尿酸在满足消毒效果的前提下,三氯甲烷浓度处于较低水平,应当注意高浓度卤乙酸可能带来的风险。
[8]
O'Brien JE, Morris JC, Butler JN. Equilibria in aqueous solutions of chlorinated isocyanurate[M]// RUBINA J, Distribution.ed. Chemistry of Water Supply, Treatment, and Ann Arbor: Ann Arbor Science Publishers,1974:333-358.
[9]
蒋玉佳, 余燕, 杨晓蓉, 等. 2013—2014年昆明市游泳池水的致病菌污染研究[J]. 环境与健康杂志, 2016, 33(3):241-242.
[10]
鲁雪立, 任秀娟, 李方茹, 等. 亳州市2022—2024年游泳池水质监测结果分析[J]. 安徽预防医学杂志, 2025, 31(3):213-215,227.
[11]
Falk RA, Blatchley ER, Kuechler TC, et al. Assessing the impact of cyanuric acid on bather's risk of gastrointestinal illness at swimming pools[J]. Water, 2019, 11(6):1314.
Current regulatory codes for swimming pool disinfection separately regulate free chlorine (FC) and cyanuric acid (CYA). It is well-known that CYA affects disinfection rates by reversibly binding to FC in aqueous solutions. However, limits for these regulated parameters have neither systematically accounted for this chemistry nor been based on the risk of gastrointestinal illness. This study was intended to determine the minimum concentration of FC relative to CYA based on the risk of gastrointestinal illness from normal fecal sloughing of selected pathogens and to find a simple regulatory rule for jointly managing FC and CYA for consistent disinfection. Literature data on CYA’s effect on microbial inactivation rates were reanalyzed based on the equilibria governing hypochlorous acid (HOCl) concentration. A model was developed that considers the rates of pathogen introduction into pool water, disinfection, turbulent diffusive transport, and pathogen uptake by swimmers to calculate the associated risk of illness. Model results were compared to U.S. Environmental Protection Agency (EPA) untreated recreational water acceptable gastrointestinal illness risk. For Cryptosporidium, correlation between log inactivation and Chick–Watson Ct was far better when C refers to HOCl concentration than to FC (r = −0.96 vs. −0.06). The HOCl concentration had a small variation (± 1.8%) at a constant CYA/FC ratio for typical FC and CYA ranges in swimming pools. In 27 U.S. states, the allowed FC and CYA results in HOCl concentrations spanning more than a factor of 500. Using conservative values for a high bather load pool with 2 mg/L FC and 90 mg/L CYA, the model predicted a 0.071 annual probability of infection for Giardia, exceeding the EPA regulatory 0.036 limit for untreated recreational waters. FC and CYA concentrations in swimming pools should be jointly regulated as a ratio. We recommend a maximum CYA/FC ratio of 20.
[12]
张璐, 彭雯洁. 氧化还原电位对游泳池水微生物指标判定的价值[J]. 预防医学论坛, 2024, 30(2):102-106.
[13]
Dragneva M. Analysis of the need to examine cyanuric acid levels in swimming pools in Bulgaria[J]. J ChemTechnolMetall, 2025, 60(6):1023-1028.
[14]
王心琪, 马茜, 曹赟, 等. 泳池水氧化还原电位卫生标准制定必要性[J]. 环境卫生学杂志, 2022, 12(10):701-704,746.
[15]
Campbell MK, Thomson S, Ramsay CR, et al. Sample size calculator for cluster randomized trials[J]. Comput Biol Med, 2004, 34(2):113-125.
Cluster randomized trials, where individuals are randomized in groups are increasingly being used in healthcare evaluation. The adoption of a clustered design has implications for design, conduct and analysis of studies. In particular, standard sample sizes have to be inflated for cluster designs, as outcomes for individuals within clusters may be correlated; inflation can be achieved either by increasing the cluster size or by increasing the number of clusters in the study. A sample size calculator is presented for calculating appropriate sample sizes for cluster trials, whilst allowing the implications of both methods of inflation to be considered.
[16]
Shrout PE, Fleiss JL. Intraclass correlations: Uses in assessing rater reliability[J]. Psychol Bull, 1979, 86(2):420-428.
Reliability coefficients often take the form of intraclass correlation coefficients. In this article, guidelines are given for choosing among six different forms of the intraclass correlation for reliability studies in which n target are rated by k judges. Relevant to the choice of the coefficient are the appropriate statistical model for the reliability and the application to be made of the reliability results. Confidence intervals for each of the forms are reviewed.
[17]
Liang KY, Zeger SL. Longitudinal data analysis using generalized linear models[J]. Biometrika, 1986, 73(1):13-22.
[18]
Laird NM, Ware JH. Random-effects models for longitudinal data[J]. Biometrics, 1982, 38(4):963-974.
Models for the analysis of longitudinal data must recognize the relationship between serial observations on the same unit. Multivariate models with general covariance structure are often difficult to apply to highly unbalanced data, whereas two-stage random-effects models can be used easily. In two-stage models, the probability distributions for the response vectors of different individuals belong to a single family, but some random-effects parameters vary across individuals, with a distribution specified at the second stage. A general family of models is discussed, which includes both growth models and repeated-measures models as special cases. A unified approach to fitting these models, based on a combination of empirical Bayes and maximum likelihood estimation of model parameters and using the EM algorithm, is discussed. Two examples are taken from a current epidemiological study of the health effects of air pollution.
[19]
Cameron AC, Gelbach JB, Miller DL. Bootstrap-based improvements for inference with clustered errors[J]. Rev Econ Stat, 2008, 90(3):414-427.
[20]
中华人民共和国卫生部. 公共场所卫生检验方法第6部分:卫生监测技术规范:GB/T 18204.6—2013[S]. 北京: 中国标准出版社, 2013.
[21]
中华人民共和国国家卫生健康委员会. 公共场所卫生指标及限值要求:GB 37488—2019[S]. 北京: 中国标准出版社, 2019.
[22]
中华人民共和国卫生部. 公共场所卫生检验方法第1部分:物理因素:GB/T 18204.1—2013[S]. 北京: 中国标准出版社, 2013.
[23]
中华人民共和国国家卫生和计划生育委员会, 中国国家标准化管理委员会. 公共场所卫生检验方法第2部分:化学污染物:GB/T 18204.2—2014[S]. 北京: 中国标准出版社, 2014.
[24]
中华人民共和国国家卫生健康委员会. 生活饮用水标准检验方法第4部分:感官性状和物理指标:GB/T 5750.4—2023[S]. 北京: 中国标准出版社, 2023.
[25]
中华人民共和国国家卫生健康委员会. 生活饮用水标准检验方法第11部分:消毒剂指标:GB/T 5750.11—2023[S]. 北京: 中国标准出版社, 2023.
[26]
中华人民共和国国家卫生健康委员会. 生活饮用水标准检验方法第12部分:微生物指标:GB/T 5750.12—2023[S]. 北京: 中国标准出版社, 2023.
[27]
中华人民共和国水利部. 氧化还原电位的测定(电位测定法):SL 94—1994[S]. 北京: 中国水利水电出版社, 1995.
[28]
中华人民共和国住房和城乡建设部. 游泳池水质标准:CJ/T 244—2016[S]. 北京: 中国标准出版社, 2016.
[29]
周婷, 兰蓝, 邱建青, 等. GEE、GLMM和MLM分析卫生重复测量资料的效果比较[J]. 现代预防医学, 2017, 44(16):2881-2885,2899.
[30]
张梦楠, 倪洋, 曾强. 大气污染人群健康流行病学研究中统计学分析方法应用进展[J]. 环境卫生学杂志, 2022, 12(5):388-395.
[31]
McNeish D, Stapleton LM. Modeling clustered data with very few clusters[J]. Multivar Behav Res, 2016, 51(4):495-518.
[32]
Peduzzi P, Concato J, Kemper E, et al. A simulation study of the number of events per variable in logistic regression analysis[J]. J Clin Epidemiol, 1996, 49(12):1373-1379.
We performed a Monte Carlo study to evaluate the effect of the number of events per variable (EPV) analyzed in logistic regression analysis. The simulations were based on data from a cardiac trial of 673 patients in which 252 deaths occurred and seven variables were cogent predictors of mortality; the number of events per predictive variable was (252/7 =) 36 for the full sample. For the simulations, at values of EPV = 2, 5, 10, 15, 20, and 25, we randomly generated 500 samples of the 673 patients, chosen with replacement, according to a logistic model derived from the full sample. Simulation results for the regression coefficients for each variable in each group of 500 samples were compared for bias, precision, and significance testing against the results of the model fitted to the original sample. For EPV values of 10 or greater, no major problems occurred. For EPV values less than 10, however, the regression coefficients were biased in both positive and negative directions; the large sample variance estimates from the logistic model both overestimated and underestimated the sample variance of the regression coefficients; the 90% confidence limits about the estimated values did not have proper coverage; the Wald statistic was conservative under the null hypothesis; and paradoxical associations (significance in the wrong direction) were increased. Although other factors (such as the total number of events, or sample size) may influence the validity of the logistic model, our findings indicate that low EPV can lead to major problems.

脚注

利益冲突声明 全部作者声明无利益冲突


PDF(948 KB)

Accesses

Citation

Detail

段落导航
相关文章

/