黃 毅

(貴州省水利水電勘測設(shè)計(jì)研究院有限公司，貴陽 550002)

0 引 言

近幾十年來，我國的水利水電事業(yè)獲得高速發(fā)展，建設(shè)很多大規(guī)模的水電水利工程，在供水、灌溉、發(fā)電、防洪等方面發(fā)揮了極大的環(huán)境效益、社會(huì)效益、經(jīng)濟(jì)效益，也為我國國民經(jīng)濟(jì)的快速、穩(wěn)定、持續(xù)發(fā)展提供了重要保障與支撐[1]。與此同時(shí)，這些工程也帶來了大批非自愿移民，帶來了很大的經(jīng)濟(jì)與人文負(fù)擔(dān)。這些水庫移民數(shù)量眾多，其中很大一部分是農(nóng)村移民。而且隨著不斷建設(shè)的水電水利工程，水庫移民的數(shù)量還將持續(xù)增加[2]。對于水庫移民來說，他們祖祖輩輩生活的地方被淹沒，失去其賴以生存的生計(jì)與生產(chǎn)資源，為國家的發(fā)展做出了十分巨大的貢獻(xiàn)，因此必須重視水庫移民的安置與補(bǔ)償問題。水庫移民通常數(shù)量龐大，且涉及到社會(huì)重建、生活生計(jì)、生態(tài)環(huán)境、區(qū)域經(jīng)濟(jì)發(fā)展、工程建設(shè)等多個(gè)方面，一旦處理不好水庫移民的安置與補(bǔ)償問題，可能會(huì)導(dǎo)致生態(tài)環(huán)境惡化、社會(huì)不安定、人民生活困難、區(qū)域經(jīng)濟(jì)結(jié)構(gòu)被遭到破壞等嚴(yán)重后果[3]。在水庫移民的安置與補(bǔ)償中，最重要的就是制訂其長期補(bǔ)償方式，使水庫移民能夠獲得穩(wěn)定、長期的經(jīng)濟(jì)收入，保障他們能夠獲得穩(wěn)定的生活，實(shí)現(xiàn)水庫移民安置后逐步致富、長期發(fā)展的目標(biāo)，避免移民安置區(qū)與庫區(qū)的不穩(wěn)定因素，同時(shí)避免產(chǎn)生移民遺留問題，真正實(shí)現(xiàn)移民群眾與工程項(xiàng)目的雙贏，維護(hù)社會(huì)穩(wěn)定、構(gòu)建和諧社會(huì)。在該背景下，對水庫移民長期補(bǔ)償方式進(jìn)行經(jīng)濟(jì)效益方面的分析，具有十分重要的應(yīng)用價(jià)值與理論意義。對于水庫移民長期補(bǔ)償方式經(jīng)濟(jì)效益分析問題，很多研究人員與學(xué)者從不同角度展開了研究，獲得了多樣化的研究成果。綜合現(xiàn)有研究成果，設(shè)計(jì)一種基于Logistic模型的水庫移民長期補(bǔ)償方式經(jīng)濟(jì)效益分析方法，實(shí)現(xiàn)水庫移民長期補(bǔ)償方式高效、準(zhǔn)確的經(jīng)濟(jì)效益分析。

1 基于Logistic模型的水庫移民長期補(bǔ)償方式經(jīng)濟(jì)效益分析方法設(shè)計(jì)

1.1 分析數(shù)據(jù)獲取與數(shù)據(jù)統(tǒng)計(jì)性描述

針對典型區(qū)域的水庫移民，按照移民比例分布，通過分層整群隨機(jī)抽樣法獲取分析數(shù)據(jù)。具體來說，需要獲取兩個(gè)階段的分析數(shù)據(jù)，分別為遷徙前、遷徙并執(zhí)行長期補(bǔ)償方案后[4]。

在遷徙前，通過隨機(jī)抽樣方式對樣本水庫移民進(jìn)行確定；在遷徙并執(zhí)行長期補(bǔ)償方案后，通過追蹤調(diào)查方法對樣本水庫移民進(jìn)行各方面的追蹤。

針對樣本水庫移民相關(guān)指標(biāo)的基本數(shù)據(jù)信息實(shí)施統(tǒng)計(jì)性描述，具體包括就業(yè)狀況、教育、人均純收入、人均居住面積、恩格爾系數(shù)、居住環(huán)境等指標(biāo)。

將就業(yè)狀況的基本數(shù)據(jù)信息分為以下幾個(gè)類別：兼業(yè)安置就業(yè)、非農(nóng)安置就業(yè)、農(nóng)業(yè)安置就業(yè)[5]。

將教育的基本數(shù)據(jù)信息分為以下幾個(gè)類別：大專及以上教育、高中或中專教育、初中教育、小學(xué)教育、小學(xué)以下。

將人均純收入的基本數(shù)據(jù)信息分為以下幾個(gè)類別：2 000元以下、2 000～3 000元、3 000～4 000元、4 000～5 000元、5 000元以上。

將人均居住面積的基本數(shù)據(jù)信息分為以下幾個(gè)類別：10～20 m2、20～30 m2、30～40 m2、40～50 m2、50～60 m2、60 m2以上。

將恩格爾系數(shù)的基本數(shù)據(jù)信息分為以下幾個(gè)類別：20%～30%、30%～40%、40%～50%、50%～60%。

將居住環(huán)境分為以下幾個(gè)類別：惡劣、良好、優(yōu)越，調(diào)查各類別比例并實(shí)施統(tǒng)計(jì)性描述。

分別對各類別的比例進(jìn)行調(diào)查，對其進(jìn)行統(tǒng)計(jì)性描述[6]。

同時(shí)，對遷徙前與遷徙并執(zhí)行長期補(bǔ)償方案后水利水電工程的補(bǔ)償投入與補(bǔ)償收益數(shù)據(jù)進(jìn)行收集，并對數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì)。

1.2 指標(biāo)數(shù)據(jù)預(yù)處理

對于獲取的指標(biāo)原始數(shù)據(jù)，對其實(shí)施以下預(yù)處理：數(shù)據(jù)清洗、標(biāo)準(zhǔn)化處理。

其中，數(shù)據(jù)清洗主要使用SVD算法填充空值，也就是對原始數(shù)據(jù)矩陣實(shí)施降維處理，補(bǔ)全原始數(shù)據(jù)矩陣[7]。

降維處理的具體步驟如下：

將原始數(shù)據(jù)矩陣T設(shè)為一個(gè)m×n矩陣，對該矩陣中保留的信息量進(jìn)行確定。用能量閾值表示該矩陣中保留的信息量，其計(jì)算公式具體如下：

(1)

式中：Pω為前ω個(gè)大奇異值的對應(yīng)平方和函數(shù)；W為全部奇異值的對應(yīng)平方和函數(shù)[8]。

通過能量閾值對ω進(jìn)行求取，具體公式如下：

(2)

式中：ω為滿足上式的最小整數(shù)值。

根據(jù)ω值，在原始數(shù)據(jù)矩陣T中對前ω個(gè)大奇異值組進(jìn)行去除，構(gòu)建新的數(shù)據(jù)矩陣∑ω。

對于數(shù)據(jù)矩陣∑ω，根據(jù)填充后矩陣的奇異值分解結(jié)果，在A和B中分別選取前ω個(gè)奇異向量，構(gòu)建新的矩陣Aω和Bω，獲取降維后即填充后的數(shù)據(jù)矩陣。

其中，在進(jìn)行填充后矩陣的奇異值分解之前，需要填充空白元素，主要通過相關(guān)列項(xiàng)目評分實(shí)際平均值對原始數(shù)據(jù)矩陣中的空白元素進(jìn)行代替，獲取一個(gè)稠密矩陣。將填充后的矩陣用T*來表示，其奇異值分解公式具體如下：

(3)

降維后即填充后的數(shù)據(jù)矩陣具體如下：

T*=Aω∑ωBω

(4)

標(biāo)準(zhǔn)化處理能夠?yàn)橛兄娜徊煌烤V與內(nèi)涵的各指標(biāo)帶來可合成性。對于福利水平與變動(dòng)正相關(guān)的指標(biāo)，通過下式進(jìn)行標(biāo)準(zhǔn)化處理：

(5)

式中：x為指標(biāo)賦值；xmin為指標(biāo)最小值；xmax為指標(biāo)最大值[10]。

對于福利水平與變動(dòng)負(fù)相關(guān)的指標(biāo)，通過下式進(jìn)行標(biāo)準(zhǔn)化處理：

(6)

就此完成指標(biāo)原始數(shù)據(jù)的預(yù)處理。

1.3 水庫移民長期補(bǔ)償方式經(jīng)濟(jì)效益分析模型構(gòu)建

基于Logistic模型構(gòu)建水庫移民長期補(bǔ)償方式經(jīng)濟(jì)效益分析模型，該模型具體可以分為兩個(gè)子模型：①遷徙前的經(jīng)濟(jì)效益分析模型；②遷徙并執(zhí)行長期補(bǔ)償方案后的經(jīng)濟(jì)效益分析模型[11]。

遷徙前經(jīng)濟(jì)效益分析模型的變量選取具體如下：遷徙前水庫移民的就業(yè)狀況、教育、人均純收入、人均居住面積、恩格爾系數(shù)、居住環(huán)境、工程項(xiàng)目補(bǔ)償投入的描述性統(tǒng)計(jì)數(shù)據(jù)，將以上數(shù)據(jù)導(dǎo)入到SPSS Statistics 軟件內(nèi)，實(shí)施Logistic回歸分析，構(gòu)建遷徙前經(jīng)濟(jì)效益分析模型[12]。

其中，各項(xiàng)參數(shù)的擬合系數(shù)分別取各種數(shù)據(jù)的解釋變量，構(gòu)建的遷徙前經(jīng)濟(jì)效益分析模型具體如下：

(7)

式中：PLW為遷徙前水庫移民生活狀態(tài)穩(wěn)定、工程項(xiàng)目補(bǔ)償投入較低的概率；PLO為遷徙前水庫移民生活狀態(tài)不穩(wěn)定、工程項(xiàng)目補(bǔ)償投入較高的概率；x1-x6分別為遷徙前水庫移民的就業(yè)狀況、教育、人均純收入、人均居住面積、恩格爾系數(shù)、居住環(huán)境的描述性統(tǒng)計(jì)數(shù)據(jù)；λ1-λ7為7種描述性統(tǒng)計(jì)數(shù)據(jù)的對應(yīng)解釋變量；w為工程項(xiàng)目補(bǔ)償投入的決定系數(shù)[13]。

模型中各維數(shù)數(shù)據(jù)的對應(yīng)VIF、容差以及條件指數(shù)具體見表1。

表1 模型中各維數(shù)的對應(yīng)VIF、容差以及條件指數(shù)(遷徙前)

根據(jù)表1的VIF、容差以及條件指數(shù)數(shù)據(jù)，遷徙前經(jīng)濟(jì)效益分析模型中的數(shù)據(jù)不存在多重共線性，說明模型參數(shù)偏差較小，可以投入使用[14]。

遷徙并執(zhí)行長期補(bǔ)償方案后的經(jīng)濟(jì)效益分析模型的變量選取具體如下：遷徙并執(zhí)行長期補(bǔ)償方案后水庫移民的就業(yè)狀況、教育、人均純收入、人均居住面積、恩格爾系數(shù)、居住環(huán)境、補(bǔ)償收益數(shù)據(jù)的描述性統(tǒng)計(jì)數(shù)據(jù)，同樣將以上數(shù)據(jù)導(dǎo)入到SPSS Statistics 軟件內(nèi)，實(shí)施Logistic回歸分析，構(gòu)建遷徙并執(zhí)行長期補(bǔ)償方案后的經(jīng)濟(jì)效益分析模型[15]。

其中，模型的擬合系數(shù)分別取各種數(shù)據(jù)的解釋變量，構(gòu)建的遷徙并執(zhí)行長期補(bǔ)償方案后的經(jīng)濟(jì)效益分析模型具體如下：

(8)

式中：PER為遷徙并執(zhí)行長期補(bǔ)償方案后水庫移民經(jīng)濟(jì)情況好、補(bǔ)償收益較高的概率；PED為遷徙并執(zhí)行長期補(bǔ)償方案后水庫移民經(jīng)濟(jì)情況差、補(bǔ)償收益較低的概率；x1′-x6′分別為遷徙并執(zhí)行長期補(bǔ)償方案后水庫移民的就業(yè)狀況、教育、人均純收入、人均居住面積、恩格爾系數(shù)、居住環(huán)境的描述性統(tǒng)計(jì)數(shù)據(jù)；ε1-ε7為7種描述性統(tǒng)計(jì)數(shù)據(jù)的對應(yīng)解釋變量；ζ為遷徙并執(zhí)行長期補(bǔ)償方案后補(bǔ)償收益的決定系數(shù)[16]。

模型中每個(gè)維數(shù)數(shù)據(jù)的對應(yīng)VIF、容差以及條件指數(shù)具體見表2。

表2 模型中每個(gè)維數(shù)的VIF、容差和條件指數(shù)

根據(jù)表2的VIF、容差以及條件指數(shù)數(shù)據(jù)，遷徙并執(zhí)行長期補(bǔ)償方案后的經(jīng)濟(jì)效益分析模型中的數(shù)據(jù)也不存在多重共線性，說明模型可以投入使用。

通過比較兩種模型，即可獲得水庫移民長期補(bǔ)償方式的經(jīng)濟(jì)效益情況。

2 經(jīng)濟(jì)效益分析實(shí)證分析

2.1 水電站長期補(bǔ)償安置計(jì)劃數(shù)據(jù)

選擇某水電站水庫移民實(shí)施設(shè)計(jì)的基于Logistic模型的水庫移民長期補(bǔ)償方式經(jīng)濟(jì)效益分析方法的實(shí)證分析。該水電站的總庫容為19 560×104m3，于2009年投入使用。該水電站于2010年開始實(shí)施水庫移民的長期補(bǔ)償安置計(jì)劃。在其長期補(bǔ)償安置計(jì)劃的實(shí)施中，主要采用3種安置方式：兼業(yè)安置、非農(nóng)安置以及農(nóng)業(yè)安置。

該水庫的長期補(bǔ)償安置計(jì)劃數(shù)據(jù)見表3。

表3 實(shí)驗(yàn)水庫的長期補(bǔ)償安置計(jì)劃數(shù)據(jù)

收集以上移民遷徙前、遷徙并執(zhí)行長期補(bǔ)償方案后的分析數(shù)據(jù)，利用設(shè)計(jì)方法針對該水庫的移民展開長期補(bǔ)償方式經(jīng)濟(jì)效益分析實(shí)驗(yàn)。

2.2 兼業(yè)安置方式經(jīng)濟(jì)效益分析的實(shí)證分析

在兼業(yè)安置方式中，水庫移民的經(jīng)濟(jì)效益分析結(jié)果見圖1。

圖1 水庫移民的經(jīng)濟(jì)效益分析結(jié)果(兼業(yè)安置方式)

根據(jù)圖1的水庫移民經(jīng)濟(jì)效益分析結(jié)果，兩種擬合曲線的數(shù)值差異較大，表明遷徙并執(zhí)行長期補(bǔ)償方案后兼業(yè)安置方式的水庫移民獲得了較大的經(jīng)濟(jì)效益，說明設(shè)計(jì)方法在兼業(yè)安置方式上有著良好的經(jīng)濟(jì)效益分析效果。

項(xiàng)目工程方的經(jīng)濟(jì)效益分析結(jié)果見圖2。

圖2 項(xiàng)目工程方的經(jīng)濟(jì)效益分析結(jié)果(兼業(yè)安置方式)

根據(jù)圖2的項(xiàng)目工程方經(jīng)濟(jì)效益分析結(jié)果可以發(fā)現(xiàn)，兩種擬合曲線的數(shù)值差異很大，表明遷徙并執(zhí)行長期補(bǔ)償方案后,項(xiàng)目工程方在水庫淹沒面積處與耕地農(nóng)用地占用面積處獲得了巨大的經(jīng)濟(jì)效益。實(shí)驗(yàn)結(jié)果同樣說明，設(shè)計(jì)方法在兼業(yè)安置方式上有著良好的經(jīng)濟(jì)效益分析效果。

2.3 農(nóng)業(yè)安置方式與非農(nóng)安置方式經(jīng)濟(jì)效益分析的實(shí)證分析

對農(nóng)業(yè)安置方式與非農(nóng)安置方式進(jìn)行綜合經(jīng)濟(jì)效益分析，其中水庫移民的經(jīng)濟(jì)效益分析結(jié)果具體見圖3。

根據(jù)圖3的水庫移民經(jīng)濟(jì)效益分析結(jié)果，兩種擬合曲線的數(shù)值整體差異較大，表明遷徙并執(zhí)行長期補(bǔ)償方案后農(nóng)業(yè)安置方式與非農(nóng)安置方式的水庫移民獲得了較大的經(jīng)濟(jì)效益。說明設(shè)計(jì)方法有著良好的長期補(bǔ)償方式經(jīng)濟(jì)效益分析效果。

圖3 水庫移民的經(jīng)濟(jì)效益分析結(jié)果(農(nóng)業(yè)/非農(nóng)安置方式)

在林地占用面積處、建設(shè)用地占用面積處以及未利用地占用面積處，項(xiàng)目工程方的經(jīng)濟(jì)效益分析結(jié)果見圖4。

圖4 項(xiàng)目工程方的經(jīng)濟(jì)效益分析結(jié)果(農(nóng)業(yè)/非農(nóng)安置方式)

根據(jù)圖4的分析結(jié)果可知，兩種擬合曲線的數(shù)值差異較大，表明遷徙并執(zhí)行長期補(bǔ)償方案后項(xiàng)目工程方在這幾處獲得了很大的經(jīng)濟(jì)效益。實(shí)驗(yàn)結(jié)果表明，設(shè)計(jì)方法有著良好的長期補(bǔ)償方式經(jīng)濟(jì)效益分析效果。

3 結(jié) 語

在對水庫移民長期補(bǔ)償方式經(jīng)濟(jì)效益分析進(jìn)行研究的過程中，應(yīng)用Logistic模型設(shè)計(jì)水庫移民長期補(bǔ)償方式經(jīng)濟(jì)效益分析模型，實(shí)現(xiàn)了各種長期補(bǔ)償安置方式的經(jīng)濟(jì)效益分析，取得了一定研究成果。