环境科学  2024, Vol. 45 Issue (5): 2748-2756   PDF    
市政污水中吗啡来源辨析
邵雪婷, 赵悦彤, 蒋冰, 裴伟, 李彦莹, 谭冬芹, 王德高     
大连海事大学环境科学与工程学院, 大连 116026
摘要: 采用分析市政污水方法可以辨别吗啡合法或非法使用情况. 由于污水中吗啡来源复杂, 区分不同来源的贡献是吗啡使用调查的关键问题. 选择典型城市, 在2022年10月至2023年3月期间采集61个代表性污水处理厂的262份进水样品, 分析了包括吗啡、可待因、蒂巴因、罂粟碱、那可丁和单乙酰吗啡在内的6种生物碱的浓度, 同时结合罂粟壳中生物碱的含量, 应用比值法和正定矩阵因子分解(PMF)模型对污水中的生物碱进行源解析. 污水中共检出5种生物碱, 但海洛因的代谢产物单乙酰吗啡未检出. 其中, 吗啡和可待因的浓度明显高于那可丁、罂粟碱和蒂巴因. 通过构建两种比值, 即可待因/ (吗啡+ 可待因)和那可丁/ (那可丁+ 可待因), 能够很好地定性区别罂粟壳使用来源. PMF的结果表明, 污水中吗啡主要来源于医用吗啡(44.9%)、罂粟壳药物(43.7%)和可待因(9.4%). 污水中吗啡来源具有显著月份变化趋势. 在新冠疫情和甲流暴发期间, 含有罂粟壳和可待因药物使用为主要来源, 而吗啡镇痛药物使用来源相对稳定. 清单分析与PMF源解析的结果具有一致性, 表明模型的可靠性, 同时也表明该城市吗啡主要是合法来源, 不存在滥用情况.
关键词: 污水      吗啡      比值法      正定矩阵因子分解(PMF)模型      罂粟壳     
Source Apportionment of Morphine in Wastewater
SHAO Xue-ting , ZHAO Yue-tong , JIANG Bing , PEI Wei , LI Yan-ying , TAN Dong-qin , WANG De-gao     
College of Environmental Sciences and Engineering, Dalian Maritime University, Dalian 116026, China
Abstract: It is a new approach to identify legal or illegal use of morphine through information on municipal wastewater. However, the sources of morphine in wastewater are complex, and distinguishing the contribution of different sources has become a key issue. A total of 262 influent samples from 61 representative wastewater treatment plants in a typical city were collected from October 2022 to March 2023. The concentrations of morphine, codeine, thebaine, papaverine, noscapine, and monoacetylmorphine were analyzed in wastewater and poppy straws. Combined with the proportion of alkaloids in poppy straws, the source analysis of alkaloids in wastewater was analyzed using the ratio method and positive matrix factorization model (PMF). Only five alkaloids were detected in wastewater, and monoacetylmorphine, a metabolite of heroin, was not detected. The concentrations of morphine and codeine were significantly higher than those of noscapine, papaverine, and thebaine. By constructing the ratios of codeine/(morphine + codeine) and noscapine/(noscapine + codeine), the source of poppy straw could be qualitatively distinguished. The PMF results showed that three sources of morphine for medical use, poppy straw, and codeine contributed 44.9%, 43.7%, and 9.4%, respectively. The different sources varied in these months due to the COVID-19 and influenza A outbreaks, in which the use of drugs containing poppy straws and codeine was the main source, whereas the use of morphine analgesics remained relatively stable. Inventory analysis further demonstrated the reliability of the source contributions from the PMF model, and morphine was not abused in this city.
Key words: wastewater      morphine      ratio      positive matrix factorization(PMF) model      poppy straw     

吗啡是阿片中提炼出来的一种生物碱, 被广泛用于临床疼痛的治疗. 然而, 长期大量使用吗啡可能会导致生理和精神依赖性. 吗啡在我国既是一种麻醉药品, 也被刑法定为毒品, 因此调查我国吗啡使用情况非常重要. 污水流行病学是一种监测药物使用的新方法, 已广泛应用于药物使用和毒品滥用调查[1 ~ 5]. 通过分析市政污水中吗啡信息, 可以快速、客观和真实地调查区域人群吗啡使用规律. 但是污水中的吗啡来源十分复杂, 包括海洛因、可待因和吗啡这3种阿片类物质以及罂粟壳的使用, 都可以通过人体代谢产生吗啡进入到污水中[6]. 如何准确定量吗啡的合法和非法使用是污水分析面临的难题.

目前, 环境中源解析技术主要包括清单分析[7 ~ 9]、扩散模型[10]和受体模型[11]. 清单分析是通过管理部门调查吗啡以及相关含有吗啡药物使用的主要数据, 结合区域人口数量, 估算区域人均使用量. 该方法可以从区域整体反映不同来源药物对吗啡数值的影响. 其次是扩散模型. 通过调查吗啡排放的主要单位, 如医院和制药企业的废水, 结合相关区域污水管道数据, 从源头逐级分析污水吗啡负荷, 辨析吗啡主要使用源对污水处理厂吗啡数据的影响, 适合源头明确的调查工作. 受体模型是来源解析中最常用的方法, 主要分为定性方法和定量方法. 定性方法有比值法[12]、特征化合物法和指纹谱图等[13]. 比值法具有简单易行的优点, 常用于定性分析[14]. 定量方法主要包括主成分分析[15]、非负约束因子模型[16]、化学质量平衡模型[17]以及正定矩阵因子分解模型(PMF)[18 ~ 20]等. PMF方法不需要测量源成分谱, 可利用约束条件解析各类源特征以及贡献率, 在源解析中被广泛应用. 与扩散模型相比, 受体模型无需考虑吗啡在污水管道中传输、扩散、吸附和降解等复杂的化学和物理过程, 一般通过源排放和受体浓度的化学成分分析推断各种源的浓度贡献率, 是源识别的主流技术方法.

污水中的吗啡主要来源于医用吗啡代谢, 可待因药物代谢, 海洛因代谢以及罂粟壳药物[21]和食品[22, 23]非法添加的使用. 根据污水中吗啡的潜在来源, 选择和分析罂粟壳的生物碱吗啡、可待因、蒂巴因、罂粟碱、那可丁[24]以及海洛因的代谢产物单乙酰吗啡, 结合比值法和PMF, 分析不同源使用对污水中吗啡的贡献, 识别并定量主要使用源, 以期为制定吗啡科学管制措施提供基础数据.

1 材料与方法 1.1 化学试剂

甲醇(HPLC ≥ 99.8%)购买自美国Merck公司;甲酸和乙酸购买自百灵威公司;Oasis HLB固相萃取柱(60 mg, 3 mL)购买自美国Waters公司. 蒂巴因、罂粟碱及那可丁购自天津阿尔塔公司. 吗啡及其内标吗啡-D3、可待因及其内标可待因-D6、6-单乙酰吗啡及其内标6-单乙酰吗啡-D3购自美国A-ChemTek公司. 可替宁及其内标可替宁-D3购买自美国ChromaDex和Cerilliant公司.

1.2 样品采集和前处理

本研究选择一个典型城市61个代表性污水处理厂. 采样时间2022年10月至2023年3月. 共采集了262份24 h混合进水样品. 采样后运回实验室立即分析.

污水样品使用固相萃取(SPE)方法进行分析. 污水样品(50 mL)经过离心后用过滤膜去除固体和悬浮颗粒物, 加入内标. 然后样品加到预先活化的固相萃取柱中, 待水样萃取完成后, 使用真空泵抽干5 min, 再使用4 mL甲醇和8 mL的5%甲醇氨来洗脱萃取物质并收集洗脱液. 用旋转蒸发器将洗脱液蒸至近干. 甲醇水溶液(10%甲醇)复溶至250 μL, 上机分析.

污水样品(40 mL)加入一定量的标准生物碱, 充分混合后静置12 h, 样品离心和过滤, 分别在水相和颗粒相中加入内标, 水相样品使用SPE方法进行分析, 颗粒相样品使用超声萃取方法进行分析.

罂粟壳样品使用研钵研磨成粉末后精确称量样品粉末0.10 g于锥形瓶中[25]. 加入20 mL含5%乙酸的20%甲醇溶液. 超声提取30 min后加入离心试管离心10 min, 吸取上清液, 上机分析.

1.3 仪器分析方法

使用超高效液相色谱-三重四级杆串联质谱法(Triple QuadTM 4500, 美国AB SCIEX公司)用于样品中生物碱的定量分析. 色谱柱温40 ºC;液相色谱流动相的水相A为0.1%甲酸水溶液, 流动相B为甲醇;流速0.4 mL·min-1, 进样体积3 µL. 色谱柱为Kinetex XB-C18 (2.1 mm × 100 mm, 2.6 µm). 质谱主要参数设置为电喷雾离子源(ESI+), 多离子反应监测(MRM)模式;离子喷射电压(IS)5 500 V, 离子源温度(TEM)500 ℃;雾化气(GS1)压力(50 × 6.9) kPa, 辅助加热气(GS2)压力(55 × 6.9) kPa, 气帘气(CUR)(35 × 6.9) kPa, 碰撞气(CAD)(9 × 6.9) kPa. 具体质谱参数如表 1所示. 生物碱的颗粒相含量、水相浓度、颗粒物-水分配系数(Kd)、方法检出限(MDL)和回收率如表 2所示. 生物碱在水相中浓度比例在82%~95%, 前处理方法满足分析的要求.

表 1 目标物质的质谱参数1) Table 1 Mass spectra parameters of target compounds

表 2 生物碱的水相浓度、颗粒物相含量、颗粒物-水分配系数(Kd)、回收率和方法检出限(MDL) Table 2 Concentrations in water and solid phase, partition coefficients(Kd), recovery, and method detection limits of all alkaloids

1.4 PMF模型原理与计算

在市政污水中吗啡等生物碱来源解析中, 应用PMF计算是通过对市政污水(受体)的生物碱数据进行分析, 提取若干因子, 利用特征化合物将因子识别为不同的源类, 再通过多元线性回归计算不同因子(源类)对污水中吗啡样品的贡献.

PMF模型的基本原理是将浓度数据矩阵(X)分解成因子贡献矩阵(G)、因子成分矩阵(F)和残差矩阵(E), 计算公式如下:

(1)

式中, n为样本个数, m为特征化合物个数, p为因子个数(即来源个数). PMF模型对每个因子的载荷和得分做非负约束, 即G≥ 0, F≥ 0.

PMF模型基于加权最小二乘法进行限定和迭代计算, 使得目标函数Q实现最小化, 从而获得最优源解析结果, 计算公式如下:

(2)

式中, xabgacfcb分别为XGF矩阵中的因素, uab为不确定度. 分析生物碱的不确定度决定了其在PMF模型中的权重, 不确定度越低则其所占权重就越大. 不确定度具体计算公式如下:

(3)

式中, EF为生物碱的误差分数, MDL为生物碱方法检出限. 当生物碱浓度低于或等于其相应MDL时, 其不确定度取MDL的5/6, 而相应的生物碱浓度则以1/2的MDL替代.

本文利用美国EPA推荐的PMF5.0软件对采样期间的污水中生物碱的来源进行解析. 源解析的结果既要符合模型计算要求也要符合实际情况. 在PMF计算过程中, 使用2~7个因子进行多次优化计算, 最终在因子数量为6时, 能合理解释来源类别并且解析残差均在-3.0~3.0之间, 5种生物碱总浓度的模型预测值与实测值很接近, 所选择的主因子能够充分解释原始数据, 表明计算结果较为可靠. PMF解析结果还需要符合实际情况要求, 包括经过识别的源谱特征和相对应源贡献具有合理性. 本研究通过源谱特征与实测结果进行对比、源贡献变化趋势以及清单分析结果比较3种方法用于评价模型结果合理性.

1.5 吗啡使用量计算

根据污水中吗啡浓度, 污水厂流量, 代谢校正因子和服务人口数量计算人均吗啡使用量. 服务人口数量可以使用尼古丁代谢产物可替宁计算. 吗啡人均使用量的计算公式如下:

(4)

式中, mMOR为吗啡人均使用量mg·(d·1 000人)-1, ρMOR为污水中吗啡浓度(ng·L-1);ρCOT为污水中可替宁浓度(ng·L-1);mCOT为可替宁人均产生量, 取值700 mg·(d·1 000人)-1[7]fMOR为吗啡的校正因子, 取值1.29[26].

2 结果与讨论 2.1 污水中生物碱浓度

污水中吗啡、可待因、那可丁、罂粟碱和蒂巴因的检出率分别为82.4%、72.1%、30.2%、17.6%和14.9%. 海洛因的代谢产物单乙酰吗啡未检出, 说明该地区的污水中吗啡没有来自海洛因的贡献. 检出浓度由高到低依次为:吗啡 > 可待因 > 那可丁 > 蒂巴因 > 罂粟碱, 对应的浓度均值分别为34.4、19.0、0.450、0.440和0.330 ng·L-1. 污水中的吗啡和可待因浓度与我国城市污水中浓度水平相当[27 ~ 33], 但是远低于欧美国家污水中的浓度[9, 34, 35], 也表明了我国吗啡的使用量处于较低水平.

污水中生物碱的浓度随时间变化如图 1所示. 从2022年10月开始, 吗啡浓度逐渐增加, 2022年12月达到最高峰, 然后开始降低, 2023年3月又开始增加. 中国疾病预防控制中心对全国哨点医院新冠和流感阳性率统计结果发现, 2022年12月和2023年3月中旬, 我国新冠病毒和流感病毒阳性率超过50%. 受到新冠疫情和甲流的影响, 2022年12月和2023年3月含有罂粟壳的药物使用量增加, 使生物碱的检出率和浓度显著高于其他月份.

图 1 2022年10月至2023年3月污水样品中5种生物碱浓度 Fig. 1 Concentrations of five alkaloids in wastewater collected from October 2022 to March 2023

2.2 比值法定性分析

比值法是环境污染源识别中简单有效的方法. 污水中5种生物碱检出率和浓度有所差异, 其中吗啡、可待因和那可丁检出率和浓度相对较高, 因此围绕3种生物碱构建比值判据.

2.2.1 罂粟壳中生物碱浓度和比值

从污水中生物碱的时间变化趋势可以发现, 含有罂粟壳药物使用是影响污水中吗啡浓度的重要因素. 为了得到罂粟壳中相应的比值数据, 对罂粟壳中的5种生物碱进行了分析. 结果表明, 罂粟壳中5种生物碱全部检出, 其检出含量由高到低依次为:吗啡 > 可待因 > 蒂巴因 > 罂粟碱 > 那可丁, 对应的含量均值分别为3.06、0.535、0.019 5、0.019 3和0.016 7 mg·g-1. 吗啡和可待因含量最高[25], 所占质量分数分别为84%和15%, 而那可丁、罂粟碱和蒂巴因含量都小于1%. 罂粟壳中可待因/(吗啡+可待因)值为0.15, 那可丁/(那可丁+可待因)值为0.03, 这些值可以作为判断罂粟壳来源的参考数值.

2.2.2 污水中生物碱比值

2022年10月至2023年3月污水中两种比值如图 2所示. 新冠疫情暴发期开始前2022年10月和11月点分布比较相似, 既有来自镇痛药物单独使用, 也有来自含有罂粟壳止咳药的使用. 2022年12月新冠疫情和2023年3月甲流暴发期主要围绕在罂粟壳指纹比值附近, 表明此期间污水中3种物质主要来源于含有罂粟壳药物使用. 2023年1月样品可以分为两部分, 一部分围绕罂粟壳指纹区域, 说明新冠疫情影响仍在继续, 另一部分样品则远离罂粟壳指纹区域, 说明部分人群已经康复. 2023年2月的大部分样品都远离罂粟壳指纹区域, 表明新冠疫情放缓, 止咳药物使用减少. 通过污水中指纹比值的时间变化趋势可以发现, 含有罂粟壳的药物使用对污水中吗啡的影响较大.

横竖线交叉点区域为罂粟壳指纹区 图 2 2022年10月至2023年3月污水样品中可待因/(吗啡+可待因)和那可丁/(那可丁+可待因)值分布 Fig. 2 Values of codeine/(morphine + codeine) and noscapine/(noscapine + codeine) in poppy straw and wastewater collected from October 2022 to March 2023

2.3 PMF定量分析 2.3.1 因子载荷和浓度贡献

图 3表示不同生物碱对6类因子浓度贡献和贡献率. 根据每个因子的特征组成, 对6类因子进行来源识别. 因子1中全部生物碱都有贡献, 这5种生物碱是罂粟壳的特征成分[21, 36], 可以归为罂粟壳使用源. 因子2中以可待因为主, 吗啡也有一定贡献. 这两种物质都是可待因的主要代谢产物, 可以归为可待因药物使用源. 可待因药物是强效中枢性镇咳药[37]. 因子3以吗啡为主, 可以归为吗啡药物使用源. 吗啡药物主要用于中重度疼痛镇痛, 如手术、严重创伤、烧伤和晚期癌症. 因子4以罂粟碱为主. 罂粟碱作为一种非特异性的平滑肌松弛药, 主要用于治疗脑血栓、肺栓塞、肢端动脉痉挛症及动脉栓塞性疼痛. 因此, 因子4可以归为罂粟碱药物使用源. 因子5以那可丁为主, 可以归为那可丁药物使用源. 那可丁主要作为镇咳和平喘药在临床中广泛使用[38]. 因子6以蒂巴因为主, 可以归为蒂巴因源. 与那可丁和罂粟碱不同, 由于蒂巴因具有强烈的惊厥毒性而未被应用于临床, 没有单独的药物使用来源, 主要作为合成多种阿片药物的中间体使用, 因此市政污水中蒂巴因可能来源制药企业生产排放或者含有罂粟壳药物使用[39].

图 3 5种生物碱对6个因子贡献浓度和贡献率 Fig. 3 Percent of five alkaloids load on six factors and contribution to the concentrations in wastewater

2.3.2 源贡献时间变化趋势

PMF得到的6个源对每种生物碱的贡献率如图 4所示. 污水中吗啡主要有3个来源, 最主要来源是吗啡药物(44.9%)和罂粟壳(43.7%)的使用, 而可待因的贡献率较低(9.4%). 可待因的主要来源是本身可待因的使用(71.5%), 罂粟壳贡献率较低(27.5%). 另外浓度较低的3种生物碱主要来源于罂粟壳和本身药物的使用或者排放, 罂粟壳贡献率在27.5%~53.3%波动. 总体上, 污水中的5种生物碱主要来源于罂粟壳的使用(38.6%), 其次是吗啡(32.1%)和可待因(22.1%)使用, 其余3种物质贡献率低于4%.

图 4 污水中不同来源对5种生物碱贡献率 Fig. 4 Source contribution on five alkaloids in wastewater

污水中吗啡来源贡献具有显著月份变化趋势, 如图 5所示. 2022年12月新冠疫情和2023年3月甲型流感暴发时间段内, 罂粟壳和可待因镇咳药物为主要来源, 吗啡药物使用贡献率最低为25%和27%. 而未集中发生疫情的时间, 医用吗啡贡献率相对较高在31%~45%之间波动. 特别是2022年10月, 疫情暴发前贡献率最高. 定性和定量的源解析结构都表明, 罂粟壳使用是市政污水中吗啡的重要来源.

图 5 污水中6种来源贡献的时间变化趋势 Fig. 5 Trend of relative contribution from six sources

2.3.3 污水中不同来源吗啡浓度变化趋势

污水中3个主要来源的吗啡浓度变化时间趋势如图 6所示. 罂粟壳来源的吗啡在污水中浓度呈现出先增加后减少然后增加的趋势. 受到新冠疫情和甲流的影响, 2022年12月和2023年3月罂粟壳来源吗啡浓度均值分别为43.2 ng·L-1和29.4 ng·L-1, 显著高于其他月份. 可待因来源吗啡也表现出类似的规律, 主要原因是在此期间罂粟壳和可待因作为治疗咳嗽的药物[40]被大量使用. 而来自镇痛药物吗啡的使用表现出较稳定的趋势, 12月和3月浓度与其他月份之间不存在显著差异性(P > 0.05), 说明吗啡作为镇痛药物使用不受新冠和甲流影响. 污水中不同来源吗啡的浓度时间变化趋势进一步证明模型计算的合理性.

图 6 污水中罂粟壳、镇痛吗啡和可待因来源的吗啡在污水中浓度变化趋势 Fig. 6 Trend of morphine concentrations from pappy straw, morphine analgesics, and codeine sources

2.3.4 吗啡主要来源使用量变化趋势

根据PMF解析的3个主要来源计算人均使用量结果如图 7所示. 罂粟壳使用源吗啡整体呈现先增高后降低再增高的趋势. 2022年10月和11月均值为3.05 mg·(d·1 000人)-1和4.50 mg·(d·1 000人)-1, 新冠疫情发生的2022年12月和2023年1月显著增加到27.7 mg·(d·1 000人)-1和9.08 mg·(d·10 00人)-1, 新冠疫情逐渐消退, 2023年2月均值降为0.240 mg·(d·1 000人)-1. 甲流暴发的3月又增加为11.4 mg·(d·1 000人)-1.

图 7 污水中罂粟壳、镇痛吗啡和可待因来源的吗啡使用量变化趋势 Fig. 7 Trend of morphine consumption from pappy straw, morphine analgesics, and codeine sources

可待因使用来源的吗啡也呈现出罂粟壳来源类似的波动趋势. 2022年10月和11月均值为0.86 mg·(d·1 000人)-1和1.26 mg·(d·1 000人)-1, 新冠疫情发生的2022年12月和2023年1月增加到2.21 mg·(d·1 000人)-1和2.26 mg·(d·1 000人)-1, 2023年2月均值降为1.71 mg·(d·1 000人)-1, 3月又升高为1.74 mg·(d·1 000人)-1. 两种来源的吗啡波动规律主要原因为罂粟壳和可待因主要作为镇咳药物在疫情发生期间被大量使用.

医学镇痛吗啡的人均使用量在采样期间保持稳定, 没有统计学意义上的显著差异(P > 0.05), 人均使用量中位数范围为1.65 ~4.35 mg·(d·1 000人)-1, 总体上中位数为3.15 mg·(d·1 000人)-1.

总体来看, 3种来源的吗啡的使用量还是低于世界平均使用水平[384.9 mg·(d·1 000人)-1][41]. 国外污水中吗啡主要来自医用吗啡的使用和海洛因的滥用[42, 43], 可待因来源贡献率较低[44]. 因为添加罂粟壳的中成药的普遍使用[45], 我国吗啡来源与国外吗啡来源存在差异.

2.3.5 吗啡使用清单分析

为了进一步验证PMF识别不同来源吗啡的可靠性, 对3种主要来源的吗啡进行清单分析. 2021年国际麻醉药管制委员会(INCB)报告[46]中国医用吗啡年使用总量为2 003 kg. 根据14亿人口估算, 医用吗啡人均使用量为3.92 mg·(d·1 000人)-1. 该计算结果与污水中医用吗啡使用量3.15 mg·(d·1 000人)-1基本一致.

2021年全国可待因使用量4 971 kg, 按照可待因6%代谢为吗啡计算[26], 可待因来源吗啡使用量为0.584 mg·(d·1 000人)-1, 也与疫情前两个月可待因来源吗啡使用量0.86 mg·(d·1 000人)-1和1.26 mg·(d·1 000人)-1相当.

2021年我国合法种植罂粟壳折算吗啡当量为8 402 kg, 用于生产可待因等生物碱的量为5 525 kg[46], 剩余吗啡2 877 kg用于生产镇咳药物, 由此计算来自罂粟壳源的吗啡人均使用量约为5.63 mg·(d·1 000人)-1, 与疫情前两个月罂粟壳来源吗啡使用量3.05 mg·(d·1 000人)-1和4.50 mg·(d·1 000人)-1接近. 因此, 根据清单分析可以发现, 吗啡3个来源使用量与污水解析结果基本一致, 再次证明了源解析结果的可靠性, 同时也表明该城市吗啡来源清晰, 不存在滥用现象.

3 结论

(1) 通过分析污水和罂粟壳中的吗啡类生物碱, 构建可待因/(吗啡+可待因)和那可丁/(那可丁+可待因)值, 能够定性区分污水中吗啡来源趋势变化规律.

(2) PMF定量分析结果表明, 城市污水中吗啡主要来源为吗啡处方药和罂粟壳, 可待因代谢贡献较低. 污水中吗啡浓度有明显时间变化趋势, 新冠疫情和甲流暴发期, 含有罂粟壳和可待因镇咳药物使用导致污水中吗啡浓度增加, 而作为镇痛药物使用的吗啡则不受疫情影响, 在污水中保持相对稳定.

(3) 清单分析结果与PMF分析基本保持一致, 进一步证明了源解析结果的可靠性, 也表明该地区吗啡主要为合法使用, 不存在非法滥用现象.

感谢山东省济宁市公安局刑事技术科学研究所杨崇俊警官对本研究的支持.

参考文献
[1] Wang D G, Dong Q Q, Du J, et al. Using Monte Carlo simulation to assess variability and uncertainty of tobacco consumption in a city by sewage epidemiology[J]. BMJ Open, 2016, 6(2). DOI:10.1136/bmjopen-2015-010583
[2] Zheng Q D, Lin J G, Pei W, et al. Estimating nicotine consumption in eight cities using sewage epidemiology based on ammonia nitrogen equivalent population[J]. Science of the Total Environment, 2017, 590-591: 226-232. DOI:10.1016/j.scitotenv.2017.02.214
[3] Xiao Y, Shao X T, Tan D Q, et al. Assessing the trend of diabetes mellitus by analyzing metformin as a biomarker in wastewater[J]. Science of the Total Environment, 2019, 688: 281-287. DOI:10.1016/j.scitotenv.2019.06.117
[4] Shao X T, Liu Y S, Tan D Q, et al. Methamphetamine use in typical Chinese cities evaluated by wastewater-based epidemiology[J]. Environmental Science and Pollution Research, 2020, 27(8): 8157-8165. DOI:10.1007/s11356-019-07504-w
[5] Shao X T, Zhao Y T, Jiang B, et al. Evaluation of three chronic diseases by selected biomarkers in wastewater[J]. ACS ES & T Water, 2023, 3(4): 943-953.
[6] Liu S Y, Yu W J, Wang Y R, et al. Tracing consumption patterns of stimulants, opioids, and ketamine in China by wastewater-based epidemiology[J]. Environmental Science and Pollution Research, 2021, 28(13): 16754-16766. DOI:10.1007/s11356-020-12035-w
[7] 刘晓, 胡京南, 王红梅, 等. 基于多源数据融合的河南省建材行业排放清单[J]. 环境科学, 2023, 44(4): 1924-1932.
Liu X, Hu J N, Wang H M, et al. Emission inventory of building material industry in Henan Province based on multi-source data integration[J]. Environmental Science, 2023, 44(4): 1924-1932.
[8] 徐小凡, 王宝庆, 吴俊成, 等. 港口自有移动源大气污染物排放清单[J]. 环境科学, 2023, 44(3): 1357-1367.
Xu X F, Wang B Q, Wu J C, et al. Port-owned mobile source air pollutant emission inventory[J]. Environmental Science, 2023, 44(3): 1357-1367.
[9] Rice J, Kannan A M, Castrignanò E, et al. Wastewater-based epidemiology combined with local prescription analysis as a tool for temporalmonitoring of drugs trends-A UK perspective[J]. Science of the Total Environment, 2020, 735. DOI:10.1016/j.scitotenv.2020.139433
[10] 任秀龙, 胡伟, 吴春苗, 等. 华北南部重污染城市周边区域二次气溶胶的化学特征及来源解析[J]. 环境科学, 2022, 43(3): 1159-1169.
Ren X L, Hu W, Wu C M, et al. Chemical characteristics and sources of atmospheric aerosols in the surrounding district of a heavily polluted city in the southern part of north China[J]. Environmental Science, 2022, 43(3): 1159-1169.
[11] Chen H Y, Teng Y G, Li J, et al. Source apportionment of trace metals in river sediments: a comparison of three methods[J]. Environmental Pollution, 2016, 211: 28-37. DOI:10.1016/j.envpol.2015.12.037
[12] 王鸣, 陈文泰, 陆思华, 等. 我国典型城市环境大气挥发性有机物特征比值[J]. 环境科学, 2018, 39(10): 4393-4399.
Wang M, Chen W T, Lu S H, et al. Ratios of volatile organic compounds in ambient air of various cities of China[J]. Environmental Science, 2018, 39(10): 4393-4399.
[13] 田福林, 陈景文, 敖江婷. 受体模型应用于典型持久性有毒物质的来源解析研究进展[J]. 环境化学, 2009, 28(3): 319-327.
Tian F L, Chen J W, Ao J T. Research progress on source apportionment of typical persistent toxic substances based on receptor models[J]. Environmental Science, 2009, 28(3): 319-327. DOI:10.3321/j.issn:0254-6108.2009.03.001
[14] 王德高, 杨萌, 贾宏亮, 等. 原油及油制品中多环芳烃化学指纹的分布规律研究[J]. 环境污染与防治, 2018, 30(11): 62-65.
Wang D G, Yang M, Jia H L, et al. Chemical fingerprinting of polycyclic aromatic hydrocarbons in crude oil and petroleum product samples[J]. Environmental Pollution and Control, 2018, 30(11): 62-65.
[15] 曹渺, 郭昌胜, 张恒, 等. 黄河流域入海口典型区域有机磷酸酯分布特征和风险评估[J]. 环境科学, 2023, 44(3): 1378-1386.
Cao M, Guo C S, Zhang H, et al. Occurrence distribution and risk assessment of organophosphate esters in a typical area of the estuary in the Yellow River basin[J]. Environmental Science, 2023, 44(3): 1378-1386.
[16] Tian F L, Chen J W, Qiao X L, et al. Sources and seasonal variation of atmospheric polycyclic aromatic hydrocarbons in Dalian, China: factor analysis with non-negative constraints combined with local source fingerprints[J]. Atmospheric Environment, 2009, 43(17): 2747-2753. DOI:10.1016/j.atmosenv.2009.02.037
[17] 周敏. 上海大气PM2.5来源解析对比: 基于在线数据运用3种受体模型[J]. 环境科学, 2020, 41(5): 1997-2005.
Zhou M. Comparison of three receptor models for source apportionment of PM2.5 in Shanghai: using hourly resolved PM2.5 chemical composition data[J]. Environmental Science, 2020, 41(5): 1997-2005.
[18] 陈明, 王琳玲, 曹柳, 等. 基于PMF模型的某铅锌冶炼城市降尘重金属污染评价及来源解析[J]. 环境科学, 2023, 44(6): 3450-3462.
Chen M, Wang L L, Cao L, et al. pollution assessment and source analysis of heavy metals in atmospheric deposition in a lead-zinc smelting city based on PMF model[J]. Environmental Science, 2023, 44(6): 3450-3462.
[19] 高越, 吕童, 张蕴凯, 等. PMF和RF模型联用的土壤重金属污染来源解析与污染评价: 以西北某典型工业园区为例[J]. 环境科学, 2023, 44(6): 3488-3499.
Gao Y, Lü T, Zhang Y K, et al. Source apportionment and pollution assessment of soil heavy metal pollution using PMF and RF model: a case study of a typical industrial park in northwest China[J]. Environmental Science, 2023, 44(6): 3488-3499.
[20] Wang D G, Tian F L, Yang M, et al. Application of positive matrix factorization to identify potential sources of PAHs in soil of Dalian, China[J]. Environmental Pollution, 2009, 157(5): 1559-1564. DOI:10.1016/j.envpol.2009.01.003
[21] 庞晓星, 麻风华, 王清华, 等. UPLC法测定强力枇杷露中吗啡、可待因和罂粟碱的含量[J]. 中医药信息, 2011, 28(1): 35-38.
[22] 刘敏敏, 谢小霞, 陈俏, 等. 液相色谱-质谱联用法(HPLC-MS/MS)快速测定凉茶中的5种罂粟壳类生物碱[J]. 食品工业科技, 2021, 42(23): 266-272.
Liu M M, Xie X X, Chen Q, et al. Rapid determination of five alkaloids in poppy shell in Herbal tea by HPLC-MS/MS[J]. Science and Technology of Food Industry, 2021, 42(23): 266-272.
[23] 李兴根, 乔勇升, 陈伟, 等. 超高效液相色谱-串联质谱法快速测定火锅底料中的5种罂粟壳生物碱残留[J]. 分析测试学报, 2018, 37(4): 446-451.
Li X G, Qiao Y S, Chen W, et al. Rapid determination of five alkaloids of poppy husk in hot pot soup by ultra performance liquid chromatography-tandem mass spectrometry[J]. Journal of Instrumental Analysis, 2018, 37(4): 446-451. DOI:10.3969/j.issn.1004-4957.2018.04.010
[24] Gümüşçü A, Arslan N, Sarihan E O. Evaluation of selected poppy (Papaver somniferum L.) lines by their morphine and other alkaloids contents[J]. European Food Research and Technology, 2008, 226(5): 1213-1220. DOI:10.1007/s00217-007-0739-0
[25] 张晓萍, 石晓峰, 张虹艳, 等. 罂粟壳的HPLC指纹图谱建立及其5个成分的含量测定[J]. 中国药房, 2021, 32(22): 2755-2760.
Zhang X P, Shi X F, Zhang H Y, et al. Establishment of HPLC fingerprint and content determination of 5 components in papaveris pericarpium[J]. China Pharmacy, 2021, 32(22): 2755-2760. DOI:10.6039/j.issn.1001-0408.2021.22.12
[26] Khan U, Nicell J A. Refined sewer epidemiology mass balances and their application to heroin, cocaine and ecstasy[J]. Environment International, 2011, 37(7): 1236-1252. DOI:10.1016/j.envint.2011.05.009
[27] Du P, Zhou Z L, Bai Y, et al. Estimating heroin abuse in major Chinese cities through wastewater-based epidemiology[J]. Science of the Total Environment, 2017, 605-606: 158-165. DOI:10.1016/j.scitotenv.2017.05.262
[28] Zhao J, Lu J J, Zhao H J, et al. Illicit drugs and their metabolites in urban wastewater: Analysis, occurrence and consumption in Xinjiang, China[J]. Science of the Total Environment, 2022, 852. DOI:10.1016/j.scitotenv.2022.158457
[29] Zhao J, Lu J J, Zhao H J, et al. Analysis, occurrence, and consumption of substances with abuse potential in Xinjiang, China, from 2021 to 2022[J]. Science of the Total Environment, 2023, 889. DOI:10.1016/j.scitotenv.2023.164310
[30] Wang J Y, Qi L K, Hou C Z, et al. Automatic analytical approach for the determination of 12 illicit drugs and nicotine metabolites in wastewater using on-line SPE-UHPLC-MS/MS[J]. Journal of Pharmaceutical Analysis, 2021, 11(6): 739-745. DOI:10.1016/j.jpha.2021.01.002
[31] Du P, Thai P K, Bai Y, et al. Monitoring consumption of methadone and heroin in major Chinese cities by wastewater-based epidemiology[J]. Drug and Alcohol Dependence, 2019, 205. DOI:10.1016/j.drugalcdep.2019.06.034
[32] Jin H B, Yang D, Hao Y B, et al. Estimation of the psychoactive substances consumption within 12 wastewater treatment plants service areas in a certain city of Guangxi, China applying wastewater-based epidemiology[J]. Science of the Total Environment, 2021, 778. DOI:10.1016/j.scitotenv.2021.146370
[33] Wang H B, Xu B Y, Yang L, et al. Consumption of common illicit drugs in twenty-one cities in southwest China through wastewater analysis[J]. Science of the Total Environment, 2022, 851. DOI:10.1016/j.scitotenv.2022.158105
[34] Bishop N, Jones-Lepp T, Margetts M, et al. Wastewater-based epidemiology pilot study to examine drug use in the Western United States[J]. Science of the Total Environment, 2020, 745. DOI:10.1016/j.scitotenv.2020.140697
[35] Gushgari A J, Venkatesan A K, Chen J, et al. Long-term tracking of opioid consumption in two United States cities using wastewater-based epidemiology approach[J]. Water Research, 2019, 161: 171-180. DOI:10.1016/j.watres.2019.06.003
[36] 张明童, 刘东升, 李冬华, 等. 罂粟壳与蜜罂粟壳水煎液中5个生物碱类成分的对比研究[J]. 中药材, 2021, 44(9): 2079-2084.
Zhang M T, Liu D S, Li D H, et al. Comparative study of five alkaloids contents between the water decoction of papaveris pericarpium and honey-stir-baked papaveris pericarpium[J]. Journal of Chinese Medicinal Materials, 2021, 44(9): 2079-2084.
[37] 刘文炜, 高玉琼, 刘建华, 等. RP-HPLC法同时测定哮喘片中吗啡、盐酸麻黄碱、磷酸可待因的含量[J]. 药物分析杂志, 2009, 29(5): 731-734.
Liu W W, Gao Y Q, Liu J H, et al. RP- HPLC simultaneous determination of morphine, ephedrine hydrochloride and codeine phosphate in Xiaochuan tablets[J]. Chinese Journal of Pharmaceutical Analysis, 2009, 29(5): 731-734.
[38] 孙素杰. 那可丁片含量均匀度测定方法的研究[J]. 药物分析杂志, 1997, 17(2): 96-97.
Sun S J. A study on the determination content uniformity of noscapine tablets[J]. Chinese Journal of Pharmaceutical Analysis, 1997, 17(2): 96-97.
[39] 左文松, 杨凯歌, 刘继华. 液相色谱-质谱联用法测定复方甘草片中可待因和蒂巴因的含量[J]. 中国处方药, 2022, 20(3): 38-39.
Zuo W S, Yang K G, Liu J H. Determination of codeine and thebaine in compound licoricetablet by LC-MS/MS[J]. Journal of China Prescription Drug, 2022, 20(3): 38-39.
[40] 卢森华, 黎强, 梁爽, 等. 超高效液相色谱-串联质谱法同时测定强力枇杷露中5种罂粟壳生物碱含量[J]. 中国中医药信息杂志, 2019, 26(12): 62-66.
Lu S H, Li Q, Liang S, et al. Simultaneous determination of five alkaloids from opium poppy shell in Qiangli Pipa Syrup by UPLC-MS/MS[J]. Chinese Journal of Information on Traditional Chinese Medicine, 2019, 26(12): 62-66.
[41] Zarei S, Salimi Y, Repo E, et al. A global systematic review and meta-analysis on illicit drug consumption rate through wastewater-based epidemiology[J]. Environmental Science and Pollution Research, 2020, 27(29): 36037-36051. DOI:10.1007/s11356-020-09818-6
[42] Gushgari A J, Driver E M, Steele J C, et al. Tracking narcotics consumption at a Southwestern U.S. university campus by wastewater-based epidemiology[J]. Journal of Hazardous Materials, 2018, 359: 437-444. DOI:10.1016/j.jhazmat.2018.07.073
[43] Bijlsma L, Botero-Coy A M, Rincón R J, et al. Estimation of illicit drug use in the main cities of Colombia by means of urban wastewater analysis[J]. Science of the Total Environment, 2016, 565: 984-993. DOI:10.1016/j.scitotenv.2016.05.078
[44] Asimakopoulos A G, Kannan K. Neuropsychiatric pharmaceuticals and illicit drugs in wastewater treatment plants: a review[J]. Environmental Chemistry, 2016, 13(4): 541-576.
[45] 陈楠. 服用含罂粟壳药物尿液中阿片类物质的检测[J]. 中国药物依赖性杂志, 2003, 12(4): 313-320.
Chen N. To test morphine component in urine from the opium poppy species of medicine[J]. Chinese Journal of Drug Dependence, 2003, 12(4): 313-320.
[46] INCB, Report of the international narcotics control board for 2022[R]. Vienna: United Nations, 2022.