應(yīng)用于地震災(zāi)害人員損失評(píng)估的人員分布統(tǒng)計(jì)方法試驗(yàn)

賴 敏，范柱子，劉學(xué)華，劉 全

(1 四川省地震局，四川 成都 610041；2 雅安市石棉縣防震減災(zāi)局，四川 石棉 625400；3 綿陽(yáng)市平武縣防震減災(zāi)局，四川 平武 622500；4 雅安市蘆山縣防震減災(zāi)局，四川 蘆山 625600)

應(yīng)用于地震災(zāi)害人員損失評(píng)估的人員分布統(tǒng)計(jì)方法試驗(yàn)

賴 敏1，范柱子2，劉學(xué)華3，劉 全4

(1 四川省地震局，四川 成都 610041；2 雅安市石棉縣防震減災(zāi)局，四川 石棉 625400；3 綿陽(yáng)市平武縣防震減災(zāi)局，四川 平武 622500；4 雅安市蘆山縣防震減災(zāi)局，四川 蘆山 625600)

鑒于基于人口數(shù)據(jù)庫(kù)統(tǒng)計(jì)資料的地震災(zāi)害人員損失模型的不足，進(jìn)行了更小空間尺度和時(shí)間窗口的人員分布統(tǒng)計(jì)試驗(yàn)。利用已有研究成果，對(duì)該數(shù)據(jù)進(jìn)行了模擬演算。表明了該統(tǒng)計(jì)工作對(duì)于提高基于人口數(shù)據(jù)庫(kù)統(tǒng)計(jì)資料的地震災(zāi)害人員損失模型的應(yīng)用精度是有益的。

人員分布；震害人員損失快速評(píng)估；模擬計(jì)算；人口數(shù)據(jù)庫(kù)

地震是突發(fā)的自然災(zāi)害，一次破壞性地震往往在極短時(shí)間內(nèi)造成十分嚴(yán)重的人員傷亡。在最短時(shí)間內(nèi)進(jìn)行地震造成的接近實(shí)際狀況的人員傷亡評(píng)估，是地震應(yīng)急和政府決策的重要依據(jù)，其重要意義不言而喻。地震災(zāi)害快速評(píng)估系統(tǒng)實(shí)際上是地震災(zāi)情的盲評(píng)估，其評(píng)估結(jié)果精度取決于收集到的基礎(chǔ)數(shù)據(jù)以及所采用的模型公式精度。作為地震災(zāi)害人員傷亡快速評(píng)估的人口分布研究相對(duì)較少，諸多地震災(zāi)害人員傷亡快速評(píng)估模型應(yīng)用的人員分布(密度)數(shù)據(jù)來(lái)源均應(yīng)用全國(guó)人口普查的數(shù)據(jù)[1-8]。從中華人民共和國(guó)統(tǒng)計(jì)局公布的第六次地方人口普查公報(bào)以及2010年5月23日頒布的中華人民共和國(guó)國(guó)務(wù)院令(第576號(hào))第八條規(guī)定：“人口普查每10年進(jìn)行一次，尾數(shù)逢0的年份為普查年度。兩次普查之間，進(jìn)行一次比例為全國(guó)1%的人口抽樣調(diào)查”可知，全國(guó)人口普查數(shù)據(jù)在實(shí)際應(yīng)用中存在兩方面問(wèn)題：一是數(shù)據(jù)空間分辨率低(縣級(jí))；二是時(shí)間分辨率低(5～10年)。隨著城市化進(jìn)程的開(kāi)展，流動(dòng)人口的增加，進(jìn)一步制約了地震傷亡人口評(píng)估的精度[7]。所以，獲取接近真實(shí)狀況的人口分布統(tǒng)計(jì)結(jié)果，對(duì)于地震災(zāi)害人員損失快速評(píng)估是十分必要的。筆者嘗試了時(shí)間間隔為月的人口統(tǒng)計(jì)數(shù)據(jù)的試驗(yàn)方法。

1 技術(shù)思路和試驗(yàn)方法

首先獲取當(dāng)?shù)亟y(tǒng)計(jì)時(shí)間以前的戶籍人員統(tǒng)計(jì)資料，進(jìn)行人員分布統(tǒng)計(jì)，獲得當(dāng)月該地域的人員分布數(shù)據(jù)資料。將每月人員分布數(shù)據(jù)與公式(1)項(xiàng)的數(shù)據(jù)相比較，即得到當(dāng)月的人員分布變化數(shù)據(jù)資料。

將上述公式(2)項(xiàng)的結(jié)果應(yīng)用于地震災(zāi)害人員傷亡快速評(píng)估的模擬計(jì)算。我們選取了四川省三個(gè)縣的六個(gè)村/場(chǎng)鎮(zhèn)(見(jiàn)圖1及表1)，作為統(tǒng)計(jì)試驗(yàn)點(diǎn)，進(jìn)行了關(guān)于春節(jié)期間人口數(shù)量變化的統(tǒng)計(jì)試驗(yàn)。鑒于2017年春節(jié)時(shí)間為1月27日(除夕)，所以確定每月25日為統(tǒng)計(jì)時(shí)間，作為參照也統(tǒng)計(jì)當(dāng)時(shí)、當(dāng)?shù)氐膽艏丝跀?shù)量。統(tǒng)計(jì)時(shí)間是2016年12月～2017年2月。六個(gè)村/場(chǎng)鎮(zhèn)參與統(tǒng)計(jì)試驗(yàn)的基本信息見(jiàn)表1。平武縣和蘆山縣統(tǒng)計(jì)點(diǎn)的村民經(jīng)濟(jì)收入以傳統(tǒng)務(wù)農(nóng)、外出務(wù)工為主。石棉縣統(tǒng)計(jì)點(diǎn)村民的經(jīng)濟(jì)收入以傳統(tǒng)務(wù)農(nóng)、當(dāng)?shù)仄髽I(yè)打工或外出務(wù)工為主。統(tǒng)計(jì)點(diǎn)地形條件是蘆山縣為丘陵，平武縣和石棉縣為山區(qū)。

2 人口統(tǒng)計(jì)試驗(yàn)數(shù)據(jù)結(jié)果

2016年12月-2017年2月的人口統(tǒng)計(jì)數(shù)據(jù)結(jié)果分別見(jiàn)表2～表7。將表2～表7的人口變化率作圖，形成圖2。

表1 參與統(tǒng)計(jì)試驗(yàn)的村/鎮(zhèn)信息表

表2 蘆山縣大川鎮(zhèn)三江村人數(shù)統(tǒng)計(jì)數(shù)據(jù)

注：1學(xué)生—在校學(xué)生及教職員工；2 機(jī)關(guān)—村(鎮(zhèn))級(jí)機(jī)關(guān)；3 工商戶—從事商業(yè)活動(dòng)的所有商戶；4 務(wù)工—在當(dāng)?shù)仄髽I(yè)工作的所有人員；5 務(wù)農(nóng)—排除前4項(xiàng)參與者的所有人員；6 流動(dòng)—未在當(dāng)?shù)厣罨蚬ぷ鞯乃腥藛T以及1、2、4項(xiàng)由于放假滯留在當(dāng)?shù)氐娜藛T。

表3 蘆山縣雙石鎮(zhèn)雙河村人數(shù)統(tǒng)計(jì)數(shù)據(jù)

表4 平武縣高村鄉(xiāng)場(chǎng)鎮(zhèn)人數(shù)統(tǒng)計(jì)數(shù)據(jù)

表5 平武縣虎牙鄉(xiāng)場(chǎng)鎮(zhèn)人數(shù)統(tǒng)計(jì)數(shù)據(jù)

表6 石棉縣栗子坪鄉(xiāng)紫馬村人數(shù)統(tǒng)計(jì)數(shù)據(jù)

圖1 統(tǒng)計(jì)點(diǎn)位置

表7 石棉縣永和鄉(xiāng)白馬村人數(shù)統(tǒng)計(jì)數(shù)據(jù)

表8 各統(tǒng)計(jì)點(diǎn)人口變化比率的相對(duì)誤差(%)

圖2 2016年12月-2017年2月人口變化趨勢(shì)

從圖2及表8可見(jiàn)，若以2016年12月的人口變化率作為基準(zhǔn)的話，3個(gè)月的各統(tǒng)計(jì)點(diǎn)的人口數(shù)量變化率的相對(duì)誤差均遵循：基準(zhǔn)-高-基準(zhǔn)的形態(tài)。這表明，人口變化較大的時(shí)間出現(xiàn)在春節(jié)期間。

3 將統(tǒng)計(jì)數(shù)據(jù)應(yīng)用于實(shí)際震例的模擬試驗(yàn)

文獻(xiàn)[7]提出了計(jì)算地震人員死亡的評(píng)估公式(見(jiàn)原文的公式8)：

(1)

關(guān)于地震期間人員分布狀況，本文引用文獻(xiàn)[7]表2，見(jiàn)表9。

表9 汶川地震前6度以上地區(qū)人員分布狀況表

按上式引用表9的數(shù)據(jù)試算結(jié)果為，地震災(zāi)害的盲評(píng)估地震死亡人員數(shù)量為83 489人。汶川地震中，汶川地震中實(shí)際死亡人數(shù)為69 227人，失蹤18 194人，總數(shù)為87 421人。評(píng)估死亡人數(shù)與實(shí)際死亡人數(shù)的相對(duì)誤差為4.709%。將表3的雅安市蘆山縣雙石鎮(zhèn)雙河村三個(gè)月的統(tǒng)計(jì)人口變化比例結(jié)果，乘以表9的不同烈度地區(qū)的人口分布，按上面介紹的方法計(jì)算，結(jié)果見(jiàn)表10。

表10 不同人員分布比例模擬汶川地震人員損失評(píng)估結(jié)果

從表10可見(jiàn)，不同統(tǒng)計(jì)時(shí)間的人員分布將導(dǎo)致地震災(zāi)害的盲評(píng)估地震死亡人員數(shù)量最大誤差為14.39%。這表明，不同時(shí)間的人口數(shù)據(jù)統(tǒng)計(jì)結(jié)果將導(dǎo)致附加的地震災(zāi)害的盲評(píng)估地震死亡人員數(shù)量誤差。

4 結(jié)語(yǔ)

通過(guò)上述的試驗(yàn)可見(jiàn)，僅僅引用人口數(shù)據(jù)庫(kù)的數(shù)據(jù)進(jìn)行地震災(zāi)害人員損失計(jì)算，是不可靠的。若要提高基于人口分布(密度)的地震災(zāi)情人員損失評(píng)估模型的精度，必須切實(shí)進(jìn)行更低統(tǒng)計(jì)地域空間密度和更小時(shí)間窗口的人員分布統(tǒng)計(jì)數(shù)據(jù)的基礎(chǔ)統(tǒng)計(jì)工作。否則，各種基于人口庫(kù)數(shù)據(jù)資料進(jìn)行的地震災(zāi)害人員損失的模型研究工作均是不適宜的。

[1] 蔡宗文，等.基于居民地分布的震害快速評(píng)估方法[J].中國(guó)地震，2007,23(4)：410-415.

[2] 危福泉，等.基于人口統(tǒng)計(jì)數(shù)據(jù)的區(qū)域快速震害評(píng)估方法[J].地震學(xué)報(bào)，2008,30(5)：518-524.

[3] 吳海燕，等.1999年以來(lái)地震生命損失研究的新進(jìn)展[J].震害防御技術(shù)，2014,34(1)：90-102.

[4] 李曉杰，等.強(qiáng)震人員損失回歸預(yù)測(cè)方法[J].西北地震學(xué)報(bào)，2012,34(1)：44-48.

[5] 李媛媛，等.地震人員傷亡評(píng)估方法研究[J].災(zāi)害學(xué)，2014,29(2)：223-227.

[6] 肖光先.震后災(zāi)害損失快速評(píng)估[J].災(zāi)害學(xué)，1991,6(4)：12-17.

[7] 韓貞輝，等.人口、房屋數(shù)據(jù)空間化及其震害快速評(píng)估的應(yīng)用[J].地震地質(zhì)，2013,35(4)：894-906.

[8] 徐國(guó)棟，等.汶川地震損失的快速評(píng)估[J].地震工程與工程振動(dòng)，2008,28(6)：74-83.

MonthDistributionStatisticalMethodExperimentUsedtothePersonnelLossAssessmentduetoEarthquakes

LAI Min1，F(xiàn)AN Zhuzi2，LIU Xuehua3，LIU Quan4

(1.Sichuan Earthquake Agency, Sichuan Chengdu 610041;2.Earthquake Disaster Reduction Bureau of Shimian, Sichuan Shimian 625400;3.Earthquake Disaster Reduction Bureau of Pingwu, Sichuan Pingwu 622500; 4.Earthquake Disaster Reduction Bureau of Lushan, Sichuan Lushan 625600, China)

In view of the deficiency of earthquake disaster personnel loss model based on population library statistics, the personnel distribution statistics of smaller space scale and time window are conducted. Using the existing research results, the data is simulated. It is indicated that the statistical work is necessary to improve the assessment accuracy of earthquake disaster personnel loss model based on population statistics.

personnel distribution；damage rapid assessment; simulation calculation；population database

2017-06-16；

2017-10-11

四川省地震局地震科技專項(xiàng)LY1704。

賴敏(1957-)，男，四川省成都市人.

P315.941

B

1001-8115(2017)04-0021-03

10.13716/j.cnki.1001-8115.2017.04.006