硅藻检验是判断水中尸体入水时间和入水地点的基础，在寻找尸源和判断案件性质方面发挥着重要作用，但是目前各地硅藻数据库建设还处在初始阶段，数据库建设受到多种因素的影响，如没有标准的建库方法、缺少硅藻识别专家、缺少硅藻检验仪器设备、硅藻变化复杂等。在上述因素的影响下，数据库建设覆盖面少且质量不一，直接影响尸体入水时间、地点判断的参考价值和准确性。
漳卫新河为海河支流漳卫南运河山东境内的分洪河道，前身为古黄河河道，1971年改称漳卫新河。该河从山东省德州市四女寺村四女寺枢纽起，上接卫运河，沿冀、鲁边界，途径山东省武城县、德州市、宁津县、乐陵市、庆云县、无棣县，河北省吴桥县、东光县、南皮县、盐山县、海兴县，在无棣县大口河入渤海。河道全长257 km，流域面积3 144 km2。漳卫新河穿过德州市主城区，河流沿岸风景点多，为发生溺水高风险水域[1]。
本研究对山东省德州市漳卫新河德城区段4处采样点进行逐月采样，分析该水域1年内硅藻的动态变化。同时，与之前数据库建设时本河流远点采样点数据进行对比分析，研究硅藻变化的规律。
1 材料与方法
1.1 主要仪器与试剂
Leica EM ACE200低真空镀膜机（德国Leica公司），HL-6多联真空抽滤仪（珠海黑马医学仪器有限公司），Quanta650扫描电镜（美国FEI公司），X-MaxN能谱仪、INCA v5.04软件（英国Oxford公司），0.45 μm一次性有机膜抽滤杯（珠海市祥臻生物科技有限公司）。
无水乙醇（天津市科密欧化学试剂有限公司），浓硝酸（浓度65%~68%，天津福晨化学试剂厂），甲醛（浓度37%~40%，莱阳市康德化工有限公司）。
1.2 水样采集
2020年1月—12月，于每月15日和30日（2月28日）对漳卫新河德城区段天衢路桥、三八路桥（距离天衢路桥0.5 km）、东方红路桥（距离天衢路桥1.0 km）和东风路桥（距离天衢路桥1.7 km）4处采样点进行采样（图1），每处采样点共取样24份，合计96份样品。
图1
采样方法：每个采样点在距离岸边1 m、水面下1 m处共取500 mL水样，加2 mL甲醛水溶液后常温避光保存。盛装水样的容器为一次性塑料瓶。取样前用取样处河水对采样器涮洗3次，取样后用纯净水涮洗3次。
1.3 样品处理
采用强酸消解-真空抽滤-扫描电子显微镜法[2]检测水样中的硅藻，具体步骤如下。
强酸消解：将水样摇匀，取60 mL水样，置于烧杯内，加入12 mL浓硝酸，90 ℃加热30 min。
真空抽滤：按消解水样与纯水按体积比1∶5加纯水至烧杯，将稀释后的消解液倒入醋酸纤维膜滤杯中抽滤，完成后加50 mL纯水继续抽滤，使滤膜表面接近中性，然后加10 mL无水乙醇抽滤，去除膜内水分，在电热板上烘干后粘在扫描电子显微镜样品台上。
样品镀膜：将样品置于Leica EM ACE200低真空镀膜机中镀金膜，样品距离喷金口4 cm，喷金时间180 s。
1.4 统计分析
含量和种属数据采用均值±标准差（
2 结果
2.1 硅藻检出情况
表1 硅藻种属组成
Tab. 1
| 纲 | 属 |
|---|---|
中心纲 （Centricae） | 小环藻属（Cyclotella）、冠盘藻属（Stephanodiscus）、圆筛藻属（Coscinodiscus）、直链藻属（Melosira）、海链藻属（Thalassiosira） |
羽纹纲 （Pennatae） | 针杆藻属（Synedra）、舟形藻属（Navicula）、菱形藻属（Nitzschia）、异极藻属（Gomphonema）、双壁藻属（Diploneis）、脆杆藻属（Fragilaria）、茧形藻属（Amphiprora）、桥弯藻属（Cymbella）、双菱藻属（Surirella）、波缘藻属（Cymatopleura）、尾管藻属（Urosolenia）、海氏藻属（Haslea）、布纹藻属（Gyrosigma）、菱板藻属（Hantzschia）、卵形藻属（Cocconeis）、曲壳藻属（Achnanthes）、等片藻属（Diatoma）、角毛藻属（Chaetoceros）、棒杆藻属（Rhopalodia） |
表2 2020年2月28日各采样点硅藻种属百分比 (%)
Tab. 2
| 序号 | 硅藻种属 | 天衢路桥 | 三八路桥 | 东方红路桥 | 东风路桥 |
|---|---|---|---|---|---|
| 1 | 针杆藻属 | 26.1 | 24.3 | 33.2 | 45.7 |
| 2 | 舟形藻属 | 36.5 | 47.7 | 35.5 | 29.8 |
| 3 | 小环藻属 | 24.5 | 21.5 | 23.2 | 16.5 |
| 4 | 菱形藻属 | 4.2 | 1.8 | 2.1 | 1.9 |
| 5 | 异极藻属 | 0.2 | 0.4 | 0 | 0.4 |
| 6 | 双壁藻属 | 0.3 | 0.5 | 0.2 | 0.2 |
| 7 | 脆杆藻属 | 3.2 | 0.7 | 1.3 | 1.8 |
| 8 | 茧形藻属 | 0.2 | 0.2 | 0.2 | 0 |
| 9 | 桥弯藻属 | 1.3 | 1 | 0.8 | 1.3 |
| 10 | 双菱藻属 | 0 | 0 | 0.5 | 0.1 |
| 11 | 冠盘藻属 | 0 | 0 | 0.8 | 0 |
| 12 | 圆筛藻属 | 1 | 1.4 | 1.1 | 1 |
| 13 | 波缘藻属 | 0.1 | 0 | 0 | 0.1 |
| 14 | 尾管藻属 | 0.6 | 0 | 0.2 | 0.2 |
| 15 | 海氏藻属 | 0 | 0 | 0.2 | 0 |
| 16 | 直链藻属 | 0.1 | 0.3 | 0 | 0.1 |
| 17 | 海链藻属 | 1.3 | 0.2 | 0.3 | 0.4 |
| 18 | 布纹藻属 | 0.3 | 0.1 | 0.2 | 0.2 |
| 19 | 菱板藻属 | 0 | 0.1 | 0.2 | 0 |
| 20 | 卵形藻属 | 0.1 | 0 | 0 | 0 |
| 21 | 曲壳藻属 | 0 | 0 | 0 | 0.1 |
| 22 | 等片藻属 | 0.1 | 0 | 0 | 0 |
| 23 | 角毛藻属 | 0 | 0 | 0 | 0 |
| 24 | 棒杆藻属 | 0 | 0 | 0 | 0 |
2.2 硅藻含量的变化
漳卫新河每个月4处采样点两次采样的平均含量和种属变化（表3）显示，不同月份硅藻的平均含量在1 054~13 041个/10 mL，种属构成为8~16个，全年硅藻含量为（3 879±3 551）个/10 mL，硅藻种属为（11.9±2.1）种。2~12月的硅藻含量和种属与1月之间差异有统计学意义（P<0.05）。对平均含量进行描述性分析，最大值为8月硅藻含量，最小值为12月含量。表4为2月各采样点的具体含量，2月15日和2月28日东风路桥采样点硅藻含量均高于其他采样点。
表4 2月份各采样点硅藻含量
Tab. 4
| 月份 | 含量/（个·10 mL-1） | 种属/种 |
|---|---|---|
| 1月 | 1 444±1 094 | 13.6±3.2 |
| 2月 | 1 755±2 2481） | 15.8±1.61） |
| 3月 | 3 881±27 6081） | 14.1±2.51） |
| 4月 | 1 371±1 0901） | 13.8±1.51） |
| 5月 | 2 113±2 2441） | 12.6±2.71） |
| 6月 | 1 199±7511） | 10.5±2.91） |
| 7月 | 3 057±2 3391） | 10.1±2.11） |
| 8月 | 13 041±3 1431） | 11.1±1.91） |
| 9月 | 4 842±6 6531） | 10.6±2.61） |
| 10月 | 8 051±9 2971） | 10.9±2.61） |
| 11月 | 4 743±3 8751） | 10.5±2.31） |
| 12月 | 1 054±1 2751） | 8.6±1.91） |
（个）
硅藻的种属个数也因月份的变化而不同，硅藻平均种属个数在2月达15种以上，在12月不到9种。
2.3 硅藻种属检出个数与硅藻含量的关系
将硅藻含量作为自变量（x）、种属个数作为因变量（y）进行线性回归分析，模型公式为y=11.206+0.000 2x，模型R2为0.078。模型通过F检验（F=7.917，P=0.006<0.05），硅藻含量的回归系数值为0.000 2（t=2.814，P=0.006<0.05）。
2.4 漳卫新河不同地点硅藻群落对比
计算漳卫新河天衢路采样点与三八路采样点、东方红路采样点、东风路采样点、乐陵采样点、宁津采样点5处采样点的Sprensen相似系数（表5）显示，距离2 km以内的3处采样点相似系数很高，不能进行有效区分；天衢路采样点与距离40 km和80 km的2处采样点有明显的差别，能对距离较远处水样进行有效区分。
表5 各采样点Sprensen相似系数
Tab. 5
| 采样点 | Sprensen相似系数 |
|---|---|
| 三八路 | 0.956 52 |
| 东方红路 | 0.978 72 |
| 东风东路 | 0.978 72 |
| 乐陵 | 0.736 84 |
| 宁津 | 0.588 24 |
2.5 漳卫新河德城区段硅藻种群组成相似性
图2
图2
NMDS分析漳卫新河硅藻种群组成的相似性结果
A：各采样点；B：各季节。
Fig. 2
NMDS results for diatom community similarity in Zhangweixin River
表3 每个月硅藻含量和种属的变化
Tab. 3 Monthly change of average diatom contentand species
（n=8，
表4 2月份各采样点硅藻含量 (个)
Tab. 4
| 采样点 | 2月15日 | 2月28日 |
|---|---|---|
| 天衢路桥 | 904 | 1 077 |
| 三八路桥 | 605 | 1 215 |
| 东方红路桥 | 619 | 933 |
| 东风路桥 | 7 276 | 1 408 |
3 讨论
河流浮尸案件的尸体溯源一直是公安工作中的重点，随着DNA技术的发展，溯源工作有了很大的提高。但水中尸体的入水位置推测仍然是阻碍案件侦破和定性的难点。硅藻检验技术的发展，为解决水中尸体入水位置提供了新的思路。通过硅藻检验比对分析具有理论基础和可行性，李峥等[14]已将此方法成功应用于实际案例。
3.1 漳卫新河德城区段硅藻种群组成特点
本研究发现，漳卫新河生态系统较为稳定，各采样点的距离较近，生态环境相似，各种影响硅藻变化的因素相近，因此硅藻群落结构的相似度较大。不同季节的环境条件（如温度）不同，影响硅藻变化的因素不同，水域中硅藻的群落组成因季节变化而改变。漳卫新河的硅藻种群变化主要为硅藻种属相对丰度的变化。各硅藻种属含量随时间的变化发生改变，除主要优势种属外的其他硅藻种属，在部分月份未检出，如棒杆藻属在2月未检出。4个采样点有一定的距离因而有环境差异，部分采样点会有变化趋势不同步的情况。部分时间采样点受外部环境干扰较大，如2月的平均硅藻含量普遍较低（1 000个/10 mL左右），但2月15日东风东路采样点硅藻含量达7 276个/10 mL。硅藻受环境因素影响，其分布具有明显的区域性特征，不同流域硅藻种属构成是根据硅藻检验推测尸体入水地点的前提[14-16]。Sprensen相似系数是用于比较两样本种属构成相似程度的常用指标[14]，相似系数越高两样本构成相似度越高。在确定溺死案例中，可以计算肺组织与尸体中及上下游硅藻种属构成之间的Sprensen相似系数。相似系数越高的水域是入水地点的可能性更大，通过分析比较Sprensen相似系数，可为推测落水溺亡位置提供方向。
3.2 硅藻数据库建设的价值及注意事项
在水中尸体检验实践中，取水样时水域的时空环境可能已经发生变化，若去取上下游流域的水样比较耗时耗力，且取到当时的水样无法代表水域的全部硅藻组成，因此一般只对尸体器官中检出的硅藻与发现尸体处当时的水样中硅藻进行比对。硅藻数据库建立后，器官硅藻检验结果可以与库中该河流硅藻数据进行比对。与同一河流不同流域硅藻数据进行比对，可以为判断死者的落水位置提供参考；与同一河流不同季节的硅藻数据进行比对，可以对水中尸体的死亡时间判断提供依据[5]。在鉴定实践中，肺组织中检出硅藻与水样硅藻优势种属组成是判断入水地点的重要依据，案件中会遇到硅藻组成相似性很高但肺中检出的个别硅藻在水样中未检出的情况。
在建库中应注意：（1）选择采样水域，调研建库地区水域概况和水中尸体多发水域，覆盖当地主要河流、湖泊、水库和溺水高发水域。（2）科学设置河流两采样点距离，可以通过硅藻相似系数区分合适距离。（3）建库样品量大，购置可以批量处理样品的水样消解、抽滤、镀膜等前处理设备。（4）利用硅藻识别统计自动化设备[15]，可提高数据统计效率。
科学成熟的硅藻数据库是支撑破案的基础，建成完善的数据库可以为案件节省人力、物力、财力，提高办案效率，提供侦查方向。本研究获得德州市漳卫新河硅藻的变化规律，可用于当地硅藻数据库的建设，利用好硅藻数据库，可以在水中尸体鉴定中发挥巨大作用。
参考文献
—2005）编委会
德州市水域中硅藻分布的研究
[J].
Survey into distribution of diatoms in local waters of Dezhou city
[J].
硅藻检验方法综述
[J].
Advances in the research of diatom test
[J].
硅藻检验方法及其在法医学溺死诊断中的应用现状
[J].
Diatom test： Development and application status in forensic diagnosis of drowning
[J].
溺死诊断的研究进展
[J].
Research progress on diagnosis of drowning
[J].
南四湖浮游硅藻季节性演替及其与环境因子的相关性
[J].
Seasonal succession of planktonic diatom community and its relationship with environmental factors in Nansi Lake
[J].
京杭大运河济宁段秋冬季硅藻分布及法医学应用
[J].
Distribution of diatoms in Jining section of Beijing-Hangzhou canal in autumn and winter and its forensic application
[J].
广东珠江水域硅藻种类调查
[J].
The investigation of diatom species composition and database of diatom scanning electron microscope in the Pearl River of Guangdong province
[J].
宁波市三江流域夏季硅藻分布
[J].
Diatoms distribution in Ningbo three-river watershed during summer
[J].
成都市主城区水中尸体多发河流区段硅藻分布
[J].
Distribution of diatoms in main sections of urban district rivers with drowning-prone in Chengdu
[J].
上海市浦东新区川杨河水域硅藻分布及其法医学应用
[J].
Distribution of diatoms in Chuanyang River of Pudong New Area of Shanghai and its forensic application
[J].
硅藻检验推测溺亡位置1例
[J].
The location of drowning inferred by diatom testing： A case report
[J].
Contribution to the determination of the place of death by drowning–A study of diatoms' biodiversity in Douro river estuary
[J].
肺组织硅藻定量分析在实验动物溺死诊断中的应用
[J].
Application of quantitative analysis of diatoms in lung tissue for the diagnosis of drowning of experimental animals
[J].
植物区系地理研究中的重要参数——相似性系数
[J].
Coefficient of similarity -- An important parameter in floristic geography
[J].
基于人工智能硅藻自动化识别系统的实际案例应用
[J].
Application of artificial intelligence automatic diatom identification system in practical cases
[J].
