尸溫相關因素推斷死亡時間的研究

楊宇雷，吳榮奇，費 耿，林中圣，施 群，肖 碧，王黎揚，葛延昌，孟航，馬開軍，陳新

（上海市公安局物證鑒定中心 上海市現(xiàn)場物證重點實驗室，上海200083）

尸溫相關因素推斷死亡時間的研究

楊宇雷，吳榮奇，費 耿，林中圣，施 群，肖 碧，王黎揚，葛延昌，孟航，馬開軍，陳新

（上海市公安局物證鑒定中心 上海市現(xiàn)場物證重點實驗室，上海200083）

目的 統(tǒng)計分析影響尸體溫度下降多因素和死亡時間的相關性，探索尸溫推斷死亡時間的創(chuàng)新應用。方法 建立影響尸體溫度（直腸溫度）下降相關因素的采集標準，收集實踐工作中的157例明確死亡時間真實案例的相關數(shù)據(jù)，根據(jù)實際工作經(jīng)驗結合統(tǒng)計學的基本要求對各因素進行量化評分，利用EXCEL和SPSS軟件對數(shù)據(jù)進行處理，采用多元線性回歸的方法統(tǒng)計分析各影響系數(shù)同死亡時間的相關性。結果 獲得了具有統(tǒng)計學意義的回歸方程，Y=25.993+0.04X1+0.172X2+0.88X3+0.047X4+0.373X5+0.347X6-0.766X7，決定系數(shù)R2=0.876。結論 該方法為尸溫推斷死亡時間的創(chuàng)新應用，經(jīng)測試可用于實際工作。

尸體溫度；死亡時間；影響因素；量化評分；多元線性回歸

Abstract：Objective To statistically analyze the correlation between factors related to body temperature drop and postmortem interval，and to establish an innovative method for estimating postmortem interval by body temperature.Methods The collection standard of the factors influencing body temperature drop was established.Data of 157 cases which had exact postmortem intervals were collected.The related factors were quantitatively scored according to practical experience and basic statistic requirements.The data were processed by EXCEL and SPSS software，and the correlations between the factors and the postmortem intervals were analyzed with multiple linear regression.Results A statistically significant regression equation was obtained：Y=25.993+0.04X1+0.172X2+0.88X3+0.047X4+0.373X5+0.347X6-0.766X7.The correlation coefficient R2was 0.876.Conclusion The study presented an innovation application of body temperature for estimating postmortem interval，which had been tested in the practical work.

Key words：body temperature；postmortem interval；related factor；quantitative score；multiple linear regression

死亡時間是指法醫(yī)檢查尸體時刻距死者死亡時刻的間隔時間，死亡時間推斷是法醫(yī)實際工作中需要首先解決的重要任務之一，也一直是法醫(yī)學鑒定的難點和研究的熱點[1]。由于死亡時間推斷的影響因素繁多，還尚無一種能完全克服各種影響因素或者忽略影響因素的死亡時間推斷方法[2]，而在公安法醫(yī)的實際工作中，死亡時間的推斷往往要滿足于案件偵破等需要，快速、準確是一切技術手段的首要要求。

本研究就是以“簡便實用，快速準確”為目的，根據(jù)大量的文獻資料和工作經(jīng)驗，選擇最有可能精確判斷死亡時間的尸體溫度變化規(guī)律，統(tǒng)籌考慮環(huán)境溫度、個體情況、衣著情況、載體情況、死亡原因、環(huán)境通風情況、濕度等影響因素，根據(jù)實際工作經(jīng)驗創(chuàng)新地量化重要的影響因素，進而通過一系列的統(tǒng)計分析尋求回歸方程并加以驗證，尋找出影響尸溫下降相關因素的系數(shù)，以實現(xiàn)利用尸溫推斷死亡時間這個技術手段的應用創(chuàng)新。

1 材料和方法

1.1 材料

2007～2011年，上海市公安局物證鑒定中心法醫(yī)室根據(jù)統(tǒng)一的規(guī)范收集的157例已知死亡時間的案例，并根據(jù)實際情況，列舉了實踐工作中常見的各種情況，按照統(tǒng)計分析的要求結合工作經(jīng)驗進行量化評分。

1.1.1 材料收集規(guī)范：

（1）主要收集檢驗時測得的直腸溫度大于環(huán)境溫度的案例；

（2）詳細記錄案發(fā)當時的天氣情況；

（3）詳細記錄死者所處環(huán)境的實時溫度、載體情況、環(huán)境通風情況等。載體需要詳細描述材質(zhì)；環(huán)境需要區(qū)分室內(nèi)室外，空間大小，有無空調(diào)，有無可感風速等；

（4）尸體檢驗時詳細記錄個體情況、衣著情況、載體情況，死亡原因。個體情況包括死者年齡，性別，胖瘦等情況，其中胖瘦情況需記錄腹部脂肪最厚處的脂肪厚度；衣著情況需要詳細記錄衣服的長短袖情況，衣服的件數(shù)，布料的材質(zhì)，有無填充物及其材料等；載體情況需要記錄詳細的載體名稱，主要分成泥石類、液體類、木質(zhì)類、布棉類；死亡原因記錄大失血、顱腦損傷、機械性窒息等情況，大失血要區(qū)分內(nèi)出血或者外出血；

（5）準確測量、詳細記錄尸體溫度，本研究采用直腸溫度表示。

（6）實際死亡時間的追蹤。根據(jù)案情反映、醫(yī)院病史、監(jiān)控信息、證人證言、犯罪嫌疑人供述等綜合判斷分析確切的死亡時間。

1.1.2 影響因素常見情況列舉

（1）個體情況：著重考慮對尸體溫度下降影響最大的個體胖瘦情況，據(jù)有關研究報道，人體的胖瘦與腹壁的脂肪厚度成正相關性[3]，本研究即采用腹壁脂肪的厚度來體現(xiàn)個體情況中人體的胖瘦。單位為：cm。

（2）環(huán)境溫度：以現(xiàn)場實測的溫度表示。單位為：℃。

（3）衣著情況、載體情況、死亡原因和通風情況的列舉如下（見表1）：本研究所使用的各個案例情況均符合相應的法律法規(guī)情況和相關的保密條款。

表1 影響尸體溫度下降因素的列舉表

1.2 方法

1.2.1 測量方法的統(tǒng)一：

（1）環(huán)境溫度的測量。經(jīng)計量單位校準過的水銀溫度計和煤油溫度計各一根，水銀溫度計放置于尸體邊，水銀球距周圍物體15cm以上，煤油溫度計置于尸體所處空間的中心位置，取兩次測量值的平均數(shù)。

（2）直腸溫度的測量。經(jīng)計量單位校準過的水銀溫度計和煤油溫度計各一根，先由水銀溫度計經(jīng)肛門向上插入10cm，放置5min，然后讀取數(shù)值，緊接著用煤油溫度計再測一次，取兩次測量值的平均數(shù)。

1.2.2 影響因素的評分

（1）個體情況：腹壁脂肪最厚處厚度的實際測量值。單位為：cm。

（2）環(huán)境溫度：以現(xiàn)場實測的溫度表示。單位為：℃。

根據(jù)實際工作經(jīng)驗判斷，將具有近似保暖功能和熱傳導功能的衣著和載體歸類并同通風情況和死亡時間量化評分（見表2）。

表2 影響尸體溫度下降因素的評分表

1.2.3 數(shù)據(jù)的整理

將表2的評分數(shù)值代入157例實際案例中，X1=個體情況、X2=環(huán)境溫度、X3=衣著情況、X4=載體情況、X5=死亡原因、X6=環(huán)境通風情況、X7=直腸溫度，得出Y=死亡時間。

1.2.4 統(tǒng)計學分析

采用SPSS17.0統(tǒng)計軟件包對數(shù)據(jù)進行多元線性回歸分析。檢驗水準α=0.05。

2 結果

2.1 模型效果檢驗（復相關系數(shù)R的檢驗和自相關問題的診斷）（見表3）

復相關系數(shù)R=0.936，決定系數(shù)R2=0.876。

F=150.060，F(xiàn) 檢驗結果概率 p=0.000＜0.001，差異有統(tǒng)計學意義，回歸模擬的擬合是有效的。

自相關問題的診斷 ：DW值一般要求1.5～2.5h，殘差與自變量互為獨立。回歸模型的DW值為1.888，說明該模型無自相關的問題，此模型可以被使用。

2.2 回歸系數(shù)的檢驗及結果（見表4）

2.3 建立回歸方程

根據(jù)以上回歸系數(shù)的檢驗及結果，建立回歸方程如下：

復相關系數(shù)R=0.936，決定系數(shù)R2=0.876。

Y為死亡時間，X1為個體胖瘦情況 （脂肪厚度表示），X2為環(huán)境溫度，X3為衣著情況，X4為載體情況，X5為環(huán)境通風情況，X6為死亡原因，X7為直腸溫度。

2.4 組內(nèi)考核和組外驗證

組內(nèi)考核，誤差在1h以內(nèi)的為60.6%，誤差在2h以內(nèi)的為87.9%，誤差在3h以內(nèi)的為95.5%；組外驗證3個案例，顯示2個案例誤差在1h以內(nèi)，1個案例誤差在2h以內(nèi)。

3 討論

根據(jù)尸溫的下降速度來推測死亡時間是最早且較準確的方法之一，國內(nèi)外法醫(yī)學家對其的研究一直沒有停止過，也找出了一定的規(guī)律性，但在實際檢案中往往做不到簡便而準確，要么是沒有充分考慮多種因素的影響而得出簡單的線性公式，要么是各種條件下多個多元回歸方程式[4-5]，不利于實戰(zhàn)中快速準確地得到結果。尸體溫度下降產(chǎn)生影響的這六個因素里，各個因素在尸體溫度下降的過程中占有的比重究竟有多大，一直是困擾著我們的一個問題。

事實證明，尸體溫度下降的規(guī)律和影響因素之間的關系研究，僅僅通過實驗室研究的手段是完全不夠的。規(guī)律的掌握或者模擬不僅僅只能依靠實驗室工作，實踐證明有些規(guī)律雖然人們不知道他確切的聯(lián)系性，但是可以通過大量數(shù)據(jù)的統(tǒng)計研究，盡量接近或者模擬這種規(guī)律性。公安法醫(yī)工作的性質(zhì)決定了我們有得天獨厚的收集各類數(shù)據(jù)及獲得確切死亡時間的優(yōu)勢，于是就有了該課題的設想和實施。

表3 模型效果檢驗匯總2）

表4 回歸系數(shù)1）的檢驗及結果

根據(jù)統(tǒng)計分析的基本要求，影響尸體溫度下降的六個指標，即個體條件、環(huán)境溫度、衣著與履蓋物、載體、環(huán)境通風情況及死亡原因，除去個體中的胖瘦情況、環(huán)境溫度本身即以數(shù)值表示，其他因素需要進行數(shù)字化評分。每一類影響因素本身包含有多種情況，也許不勝枚舉，但在實際工作中我們發(fā)現(xiàn)每類影響因素中均有常見的情況，同時根據(jù)經(jīng)驗判斷將有些對尸體溫度下降有較類似作用的常見情況進行歸類，然后再加以列舉，最后根據(jù)各種常見情況的影響力大小，依據(jù)工作中的經(jīng)驗判斷對于每項進行賦值量化。

經(jīng)采用SPSS軟件，多元逐步回歸死亡時間Y和個體情況 X1、環(huán)境溫度X2、衣著及覆蓋物情況X3、載體情況X4、環(huán)境通風情況X5、死亡原因X6、直腸溫度X7之間的相關性，得到回歸方程：

復相關系數(shù)R=0.936，決定系數(shù)R2=0.876。

從回歸方程的復相關系數(shù)和方差分析結果表明，本課題采用的死亡時間同影響因素之間的多元回歸方程是有統(tǒng)計學意義的。從組內(nèi)考核和組外驗證情況來看，該回歸方程計算出的死亡時間能夠較準確反映出真實的死亡時間，且相應的誤差基本能夠在公安工作需要的允許范圍之內(nèi)，具有一定的實際應用價值。

目前根據(jù)157例案例的統(tǒng)計分析，得出各影響因素對應的影響系數(shù)分別為：個體胖瘦情況0.040，環(huán)境溫度0.172，衣著情況0.88，載體情況0.047，環(huán)境通風情況0.373，死亡原因0.347，直腸溫度-0.766。其中個體胖瘦情況和載體情況的影響系數(shù)相對較小，僅僅從回歸方程本身來講，其他因素保持不變的情況下，個體胖瘦和載體情況的變動對結果的影響較小，該兩項影響系數(shù)似乎與實際情況和經(jīng)驗判斷有一定誤差。分析其中的原因，我們考慮如下：（1）由于樣本量的限制，目前僅有157例，各影響因素的系數(shù)準確度還有待進一步提高；（2）根據(jù)實際情況，觀察我們收集的數(shù)據(jù)本身，在個體胖瘦情況和載體情況中這兩項評分中，數(shù)據(jù)均較集中，也即大部分個體胖瘦均較一致，且上海地區(qū)主要以泥石類為主要現(xiàn)場；（3）在評分環(huán)節(jié)，特別是載體類因素中，由于缺乏對熱傳導相關數(shù)據(jù)的證實，評分中未體現(xiàn)出不同類的具體差異。綜上，受制于樣本量等因素，本課題所得出的多元回歸方程，不可否認，還是存在著一定的誤差，在實際應用中需要借助其他推斷死亡時間的方法綜合判斷，能做出更接近真實死亡時間的推斷。

本課題的創(chuàng)新應用研究，不僅僅是為了獲得一個具有統(tǒng)計學意義且可為實踐所用的多元回歸方程，更重要的是結合公安法醫(yī)在收集數(shù)據(jù)方面得天獨厚的優(yōu)勢，從課題設計思路上進行創(chuàng)新，而方程的獲得和檢驗結果的分析，更加證明本課題的設計是可以滿足尸體溫度下降因素同死亡時間之間的相關性研究，且該統(tǒng)計方法是有效的，能發(fā)揮作用的，是利用尸體溫度下降規(guī)律來推算死亡時間的創(chuàng)新應用。

為了使本課題所研究的推斷死亡時間的新方法能夠得到普遍應用，在繼續(xù)收集各種場所中各種已知明確死亡時間的案例，用于不斷補充統(tǒng)計研究的樣本量的同時，還需要對研究思路、自變量的選擇、自變量的評估方面不斷地完善。

[1]黃光照，麻永昌.中國刑事科學技術大全法醫(yī)病理學[M].北京：中國人民公安大學出版社，2002：116-133.

[2]C.Henssge，B，Madea.Estimation of the time since death.Forensic Sci Int[J].2007，(165)：182-184.

[3]劉力劍，趙文海，榮梅.成人腹壁脂肪厚度與體重指數(shù)、腰圍的相關性分析[J].中國誤診學雜志，2009，(31)：7577-7578.

[4]八省市（區(qū)）公安廳（局）尸溫研究小組.581例非正常死亡尸體直腸溫度與時間關系的研究[J].刑事技術，1984，（4）：1-10.

[5]徐長苗，陳希生.推定死亡時間的新方法[J].刑事技術，1994，（3）：23-25.

（本文編輯：秦志強）

Study on Estimating Postmortem Interval with Factors Related to Body Temperature

YANG Yu-lei，WU Rong-qi，F(xiàn)EI Geng，LIN Zhong-sheng，SHI Qun，XIAO Bi，WANG Li-yang，GE Yan-chang，MENG Hang，MA Kai-jun，CHEN Xin
（Shanghai Key Laboratory of Crime Scene Evidence，Institute of Forensic Science，Shanghai Public Security Bureau，Shanghai 200083，China）

DF795.4

A

10.3969/j.issn.1671-2072.2012.05.009

1671-2072-（2012）05-0075-04

2012-05-27

公安部應用創(chuàng)新計劃項目（2008YYCXSHSJ032）

楊宇雷（1978—）男，碩士研究生，主要從事法醫(yī)學研究工作。 E-mail：wangnianjiao2007@yahoo.com.cn。

馬開軍（1975—），男，副主任法醫(yī)師，主要從事法醫(yī)病理學和法醫(yī)臨床學研究工作。E-mail：makaijun@sina.cn。