吴 丹,吴凤平

(河海大学商学院,江苏 南京 210098)

流域初始水权配置主要包括 2个层次[1-2]:第1层次,即流域初始水权在各区域(涉及流域内省、市和县3级行政区域)之间的配置;第2层次,即区域初始水权在各用水部门(涉及区域内用水行业或产业和最终用水户)之间的配置。目前,针对流域初始水权第1层次的配置,国内外学者已进行了系统性的探索研究,国外学者主要根据具体国家和地区的社会制度、水资源情况、文化传统,沿袭了具有各自历史合理性的配置方法[3-9],通过立法和过程模拟进行分配研究;国内学者主要通过系统地构建一套水权配置指标体系,建立综合指标评价模型[1,10-11]对其进行分配研究。

针对流域初始水权第 2层次的配置,国内外学者研究的焦点主要集中于水权配置应遵循的原则与优先序位规则,通过构建目标规划模型[12-16]进行分配研究。其中,学者们均认为水权配置遵循的原则包括人的基本生活用水优先原则、保障社会稳定和粮食安全原则、重视生态及保护环境原则、经济效益原则、产业结构合理原则。针对不同用水行业初始水权配置结果是否合理,不同行业能否接受配置结果,鲜见学者进行系统探讨。吴凤平等[17]提出了一种对水权配置方案进行和谐性判别的方法,并在和谐性判别的基础上研究了一种能有效调整不同地区、不同行业初始水权配置结果的分配方法。但构建的各行业判别指标过于单一,未充分体现行业的水权配置结果与其发展目标之间是否存在一定的匹配关系。

鉴于国内外相关文献,为保障不同用水行业初始水权配置结果的合理性,满足不同行业水权配置结果与其发展目标之间的匹配关系,将流域内各区域的生活、生态环境以及生产等用水行业作为初始水权配置对象,从系统角度构建面向行业的初始水权配置系统,在剖析初始水权配置系统内涵与特征的基础上,提出初始水权配置系统优化思想,结合系统中各用水行业的发展目标构建初始水权配置系统诊断指标体系,建立改进的理想解模型,对流域内各系统中不同用水行业的初始水权配置结果进行诊断,并根据诊断结果结合不同用水行业的发展目标,基于交互式决策理论,将目标规划模型与逐步折衷法相结合,建立初始水权配置系统优化模型,对流域内各系统中不同用水行业的初始水权进行同步调整,以充分体现流域内各系统中不同行业的水权配置结果与其发展目标之间保持一定的匹配关系,保障流域内各系统中不同行业的水权配置结果合理。

1 面向行业的初始水权配置系统优化思想

1.1 初始水权配置系统的内涵与特征

面向行业的初始水权配置系统是由流域内各区域的生活、生态环境以及生产等相互作用、相互影响、相互制约的不同用水行业构成的具有一定结构、目标与功能的系统,是在水资源管理制度和法律法规体系的前提保障下,充分考虑流域内各区域的水资源与社会、经济以及生态环境之间相互协调与制约关系,结合生活、生态环境以及生产等用水行业的发展目标,为防止不同用水行业之间的恶性竞争用水问题,以区域产业结构布局的调整与优化为核心,对各用水行业初始水权进行配置的系统。其中,生活、生态环境、第一产业、第二产业以及第三产业等用水行业既是系统不可分割的有机组成部分,又是相对独立的子系统。

面向行业的初始水权配置系统既要统筹兼顾生活、生态环境、生产的用水需求,又要限制高耗水产业的发展,在保障社会稳定和粮食安全的基础上,将水资源配给低耗水、高效益的生产用水行业,最终实现区域产业结构布局的调整与优化,提高区域的社会经济综合效益,以水资源的可持续高效利用支撑区域经济社会的可持续发展。

1.2 初始水权配置系统优化的总体思想

面向行业的初始水权配置系统优化,是在确定流域内各系统中不同行业初始水权配置结果的基础上,通过对流域内各系统中不同行业的水权配置结果与其发展目标之间的匹配关系的诊断,以验证不同用水行业初始水权配置结果的合理性。若配置结果不合理,则根据诊断结果,结合不同用水行业的发展目标,通过流域内各系统中不同行业之间的多轮交互协商,对不同行业的初始水权进行同步调整,以保障流域内各系统中不同行业的水权配置结果与其发展目标之间保持一定的匹配关系,实现不同用水行业初始水权的优化配置。

面向行业的初始水权配置系统优化过程分为2个阶段,第1阶段为初始水权配置系统诊断过程,第2阶段即初始水权配置系统优化过程。第1阶段是在确定流域内各系统中不同行业初始水权配置结果的基础上,构建初始水权配置系统诊断指标体系,建立初始水权配置系统诊断模型,对流域内各系统中不同行业的水权配置结果与其发展目标之间的匹配关系进行诊断,以验证不同用水行业水权配置结果的合理性。第2阶段是基于不同用水行业在不同发展时期的用水需求常常处于不断调整、修正、演化的状态,根据诊断结果,结合不同用水行业的发展目标建立初始水权配置系统优化模型,通过对流域内各系统中不同用水行业的水权配置结果进行同步调整,充分体现流域内各系统中不同行业的水权配置结果与其发展目标之间的匹配关系。

2 面向行业的初始水权配置系统诊断

2.1 初始水权配置系统诊断指标体系

面向行业的初始水权配置系统诊断是对流域内各系统中不同用水行业初始水权配置结果的合理性进行的诊断,主要体现流域内各系统中不同行业的水权配置结果与其发展目标之间的匹配关系。目前,学者们主要针对水资源承载力评价指标体系、水资源配置系统协调发展评价指标体系、水资源配置系统可持续发展评价指标体系以及水资源配置系统协同度评价指标体系等水资源综合评价指标体系进行研究探讨。基于此,根据各区域不同用水行业的初始水权配置结果,结合流域内各系统中不同用水行业的发展目标,借鉴水资源综合评价指标体系,采用文献研究法、频度统计法、主成分分析法以及专家咨询法等理论分析方法,在全面深入各区域调研、与各区域利益相关者访谈、征求流域各区域管理机构和水权配置领域专家意见的基础上,确定影响流域内各系统中生活、生态环境、农业、工业以及第三产业水权配置的诊断指标,构建面向行业的初始水权配置系统诊断指标体系,见图1。

图1 面向行业的初始水权配置系统诊断指标体系

从图1可知,影响流域内各系统中不同行业水权配置的诊断指标主要是结合生活、生态环境、农业、工业以及第三产业等不同用水行业的发展目标,以体现流域内各系统中不同行业的水权配置结果与其发展目标之间的匹配关系。基于初始水权配置系统诊断指标体系,分别以生活、生态环境、农业、工业以及第三产业等用水行业水权配置结果为研究对象,通过构建诊断模型,可进一步验证流域内各系统中不同用水行业水权配置结果的合理性。

2.2 初始水权配置系统诊断模型

结合初始水权配置系统中生活、生态环境、农业、工业以及第三产业等用水行业的诊断指标,采用改进的理想解模型,对流域内各系统中不同用水行业水权配置结果的合理性进行诊断。改进的理想解模型构建的具体步骤如下:

步骤1分别以生活、生态环境、农业、工业以及第三产业等用水行业为研究对象,结合面向行业的初始水权配置系统诊断指标体系,确定流域内各系统中不同用水行业水权配置结果的诊断指标。

设系统di(即流域内第i区域水权配置系统,i=1,2,…,n)第j(j=1,2,…,5)用水行业水权配置结果的诊断指标体系可表示为Yi={Yi1,…,Yi5},其中,生活诊断指标Yi1={yi11,yi12,yi13},生态环境诊断指标Yi2={yi21,yi22,yi23},农业诊断指标Yi3={yi31,yi32,yi33,yi34},工业诊断指标Yi4={yi41,yi42,yi43},第三产业诊断指标Yi5={yi51,yi52,yi53}。针对流域内各系统中同一用水行业水权配置结果的诊断指标进行指标规范化处理,构造加权指标值。即

式中:xijk为系统di第j用水行业水权配置结果的第k项加权指标值;yijk为影响系统di第j用水行业水权配置结果的第k项诊断指标值;Zijk表示对流域内各系统中同一用水行业的诊断指标进行归一化处理后的指标值;wijk为yijk的权重,可充分利用规范化评价的先验信息以及专家咨询法等相结合的综合赋权法予以确定;k=1,…,L,其中L为系统di第j用水行业的指标总个数(j=1,2,4,5时L=3;j=3时,L=4)。

步骤2确定系统di第j用水行业水权配置的第k项诊断指标的理想值和负理想值。式中:K为效益型指标的下标组成的集合;K′为成本型指标的下标组成的集合。

步骤3计算系统di第j用水行业水权配置结果的诊断指标到指标理想值的距离。

式中为系统第j用水行业水权配置结果的诊断指标到其理想值的距离为系统di第j用水行业水权配置结果的诊断指标到其负理想值的距离。

步骤4计算系统di第j用水行业水权配置结果的诊断指标对其指标理想值的相对接近度。

式中:Cij为系统di第j用水行业水权配置结果的诊断指标对其指标理想值的相对接近度,0≤Cij≤1;α为决策者对诊断指标到其理想值距离的偏好权重;β为决策者对诊断指标到其负理想值距离的偏好权重 ,通常取 α=β=0.5。

步骤5检验系统di第j用水行业水权配置结果与其发展目标之间的匹配关系。可用公式表示为

式中:Wij为系统Si第j用水行业的配水量(Wi1,Wi2,Wi3,Wi4,Wi5分别为系统Si的生活、生态环境、农业、工业以及第三产业的配水量);WDij为系统Si第j用水行业的需水量;Eij为系统di第j用水行业的水权配置效率(Ei1,Ei2,Ei3,Ei4,Ei5分别为系统di的人均生活配水量单位绿地面积配水量、单位灌溉面积配水量、单位第二产业产值配水量、单位第三产业产值配水量);Pi,Si,TAi,TIi,TSi分别为系统di的人口总数、绿地面积、灌溉面积、工业产值、第三产业产值;ηj为流域内各系统中第j用水行业(Sij,)的配水效率差异阈值,可根据现状年国内外典型流域的配水方案,依据流域内各系统不同行业的特点,分析不同系统同一行业的配水效率差值,通过专家咨询的方法取其平均值予以确定。充分体现了流域内各系统中不同用水行业配水效率差异的合理性。

对于流域内各系统中的生活、生态环境、农业、工业以及第三产业等用水行业的水权配置结果,如果均满足式(5),则认为流域内各系统中不同用水行业的水权配置结果合理;反之,则需要对流域内各系统中不同用水行业的水权配置结果进行调整。

3 面向行业的初始水权配置系统优化

根据初始水权配置系统的诊断结果,若流域内各系统中不同用水行业的水权配置结果未通过系统诊断,则必须通过各系统中不同用水行业之间的交互协商对初始水权配置系统进行优化。基于此,按照初始水权配置系统中不同用水行业初始水权配置的优先序位,在构建不同用水行业初始水权配置的目标规划模型基础上,基于交互式决策理论,将交互式决策方法中的逐步折衷法[18]与目标规划模型相结合,建立初始水权配置系统优化的交互式目标逐步折衷耦合模型,通过多轮交互式迭代过程对流域内各系统中不同用水行业的初始水权进行适用性调整,以保障流域内各系统中不同用水行业的水权配置结果通过系统诊断。初始水权配置系统中不同用水行业初始水权配置的优先序位如下:①优先满足人的基本生活用水需求;②优先保障粮食生产用水需求;③保障河道外生态环境的用水需求;④为提高整体的国民经济效益,按照第二产业、第三产业以及第一产业的用水序位,合理保障三次产业的用水需求;⑤确定三次产业结构比例,以促进三次产业之间的协调发展。

交互式目标逐步折衷耦合模型是逐步折衷法与目标规划模型相结合的耦合模型,其将流域内各系统中不同用水行业的发展目标细分为3个子集,即:①必须加以改进现有水平的目标集;②至少应保持现有水平的目标集;③可适当减少现有水平的目标集。通过多轮交互迭代过程,不断调整流域内各系统中不同用水行业的目标值与约束条件,然后建立新的目标规划函数,以修正各系统不同用水行业的配水量,从而获得较满意的水权配置结果。

交互式目标逐步折衷耦合模型有2个核心思想:第一,充分考虑到决策者的偏好,用 Δqij修正约束条件;第二,为避免新模型与旧模型的解重复,用新的约束条件替代旧模型的约束条件,同时,不改变用水行业水权配置的优先级别。根据交互式目标逐步折衷耦合模型的思想,通过流域内各系统中不同用水行业之间的交互协商,假定用水行业的目标调整规则为:①系统Si的生活用水行业、生态环境用水行业的发展规划目标至少应保持不变;②系统Si的农业用水行业在保障粮食总产量的基础上,其配水量Wi3可适度降低,降低的量为 Δqi3(流域内其他系统的行业配水量若需调整,则可类似考虑),建立系统优化的交互式目标逐步折衷耦合模型,即

式中为系统Si第j用水行业经水权调整后得到的解;Zi1,Zi2,Zi3,Zi5,Zi7,Zi8分别为系统Si生活用水的目标值、粮食总产量的目标值、河道外生态环境用水的目标值、第三产业生产总值的目标值、第一产业与第二产业产业结构比的目标值、第二产业与第三产业产业结构比的目标值;ai,bi,ci,ei分别为系统Si每立方米水的粮食产量、第二产业增加值、第三产业增加值、第一产业增加值,分别为系统Si第m个目标超过Zim目标值的正偏差量,以及未达到Zim目标值的负偏差量。

采用GLPS目标规划软件,对式(7)进行求解,对于流域内各系统的不同用水行业可获得系统Si不同用水行业的新的水权配置结果。通过多轮交互迭代,最终保障流域内各系统中不同用水行业的水权配置结果通过系统诊断。

4 实例分析

以大凌河流域为例,在确定规划年2030年该流域各区域的初始水权配置方案的基础上,对各区域不同行业的配水量进行诊断与优化。大凌河是一条独流入海的河流,大凌河流域是内蒙古自治区、辽宁省、河北省3个省级行政区域的跨省流域,包括辽宁省的锦州市、阜新市、朝阳市、盘锦市、葫芦岛市5个地级市13个县(市、区)及内蒙古自治区通辽市、赤峰市和河北省承德市等部分地区。规划年2030年大凌河流域各地级市的社会经济发展指标值可通过大凌河流域各地市的水资源公报及《水权制度建设试点经验总结——大凌河流域初始水权制度建设资料汇编》[19]以及调研等方式予以得到。大凌河流域多年平均可分配水资源总量为12.494亿m3,流域内各区域的水权初始配置方案见表1。

表1 大凌河流域各区域的水权初始配置方案 万m3

根据表1,以生活、生态环境、农业、工业以及第三产业等用水行业为研究对象,结合改进理想解模型,对大凌河流域各区域不同用水行业的水权配置结果进行诊断,依据区域用水行业的特点,确定生活、生态环境、农业、工业以及第三产业配水效率的差值阈值η1=20m3/人,η2=0.15万m3/hm2,η3=0.45万m3/hm2,η4=0.004m3/元,η5=0.0002m3/元。大凌河流域各区域不同用水行业的诊断指标特征值见文献[20]。

根据式(1)～(4)得到流域各区域中同一用水行业发展目标的相对贴近度,见表2。

表2 大凌河流域各区域同一用水行业发展目标的相对贴近度

经诊断可知,行业水权配置结果存在不合理性。结合优化模型式(7),优化大凌河流域各区域不同用水行业的配水量,使流域各区域的用水行业水权配置结果与其发展目标之间保持一定的匹配关系。大凌河流域各区域的水权优化配置方案见表3。

表3 大凌河流域各区域的水权优化配置方案 万m3

表3中,赤峰和通辽第三产业水权配置结果虽未满足式(5),但其已达到第三产业需水量,因此可认为其第三产业的水权配置结果合理,即流域各区域的用水行业水权配置结果与其发展目标之间保持一定的匹配关系,流域各区域不同用水行业的水权配置结果合理。

5 结 语

在构建面向行业的初始水权配置系统、剖析初始水权配置系统内涵及其特征的基础上,提出了面向行业的初始水权配置系统优化的思想,建立了初始水权配置系统诊断指标体系与模型,通过系统优化的交互式目标逐步折衷耦合模型,实现了流域各区域的不同用水行业水权合理配置。基于大凌河流域各区域不同行业的水权配置方案,进一步验证了模型的合理性。

:

[1]吴凤平,葛敏.水权第一层次初始分配模型[J].河海大学学报:自然科学版,2005,33(2):216-219.

[2]葛敏,吴凤平.水权第二层次初始分配模型[J].河海大学学报:自然科学版,2005,33(5):592-594.

[3]GOPALAKRISHNAN C.The doctrine of prior appropriation andits impact on water development:a critical survey[J].American Journal of Economics and Sociology,1973,32(1):61-72.

[4]HOWE C W,SCHURMEIER D R,SHAW W D.Innovative approaches to water allocation:the potential for water markets[J].Water Resources Research,1986,22(4):439-445.

[5]MATHER J.Water resources[M].New York:John Wiley&Sons,Inc.,1984.

[6]WILLIAM G.Waterlaw[M].2nded.Michigan:Lewis Publishers,Inc.,1988.

[7]KEL MAN J,KELMAN R.Waterallocationfor economic production in a semi-arid region[J].WaterResources Development,2002,18(3):391-407.

[8]TEERINK J R,MASAHIRO N.美国日本水权水价水分配[M].刘斌,高建恩,王仰仁,译.天津:天津科学技术出版社,2000.

[9]ROBERT B,EDWYNA H.Efficiency gainsfrom water markets:empirical analysis of water move in Australia[J].Agricultural Water Management,2008,95:391-399.

[10]陈艳萍,吴凤平,吴丹.基于模糊优选和TOPSIS法的流域初始水权分配模型[J].河海大学学报:自然科学版,2009,37(4):467-471.

[11]佟金萍,王慧敏,牛文娟.流域水权初始分配系统模型[J].系统工程,2007,25(3):105-110.

[12]王治.关于建立水权与水市场制度的思考[N].中国水利报,2001-12-25(4).

[13]曾勇,杨志峰,刘静玲.一种区域水资源初始产权配置模型[J].资源科学,2005,27(2):28-32.

[14]周念来,纪昌明.区域水权初始配置模型研究[J].水电能源科学,2007,25(3):6-8.

[15]董增川.水资源规划与管理[M].北京:中国水利水电出版社,2008.

[16]吴丹,吴凤平,陈艳萍.面向行业的初始水权配置系统模型构建[J].水利水电科技进展,2010,30(5):29-32.

[17]吴凤平,葛敏.基于和谐性判断的交互式水权初始分配方法[J].河海大学学报:自然科学版,2006,34(1):104-107.

[18]汪敏,王丽铮,朱美琪.一种交互式船型优化方法的探讨及应用[J].武汉理工大学学报:交通科学与工程版,2002,26(1):99-102.

[19]高而坤,党连文.水权制度建设试点经验总结:大凌河流域初始水权制度建设资料汇编[G].北京:中国水利水电出版社,2008.

[20]吴风平,吴丹,陈艳萍.流域初始水权配置系统方案诊断模型[J].系统工程,2010,28(4):24-29.