CN112633645B - 一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法 - Google Patents
一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法 Download PDFInfo
- Publication number
- CN112633645B CN112633645B CN202011451687.2A CN202011451687A CN112633645B CN 112633645 B CN112633645 B CN 112633645B CN 202011451687 A CN202011451687 A CN 202011451687A CN 112633645 B CN112633645 B CN 112633645B
- Authority
- CN
- China
- Prior art keywords
- social
- water resource
- matrix
- water
- index
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 192
- 230000008901 benefit Effects 0.000 title claims abstract description 142
- 238000000034 method Methods 0.000 title claims abstract description 65
- 230000000694 effects Effects 0.000 title claims abstract description 22
- 238000005065 mining Methods 0.000 title claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 70
- 238000011156 evaluation Methods 0.000 claims abstract description 59
- 238000004519 manufacturing process Methods 0.000 claims abstract description 45
- 238000004458 analytical method Methods 0.000 claims abstract description 29
- 238000011160 research Methods 0.000 claims abstract description 18
- 238000005516 engineering process Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000009189 diving Effects 0.000 claims description 9
- 230000008020 evaporation Effects 0.000 claims description 8
- 238000001704 evaporation Methods 0.000 claims description 8
- 238000003973 irrigation Methods 0.000 claims description 8
- 230000002262 irrigation Effects 0.000 claims description 8
- 238000000611 regression analysis Methods 0.000 claims description 8
- 230000008878 coupling Effects 0.000 claims description 7
- 238000010168 coupling process Methods 0.000 claims description 7
- 238000005859 coupling reaction Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- 238000013139 quantization Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 238000001556 precipitation Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 238000004445 quantitative analysis Methods 0.000 abstract description 3
- 238000004451 qualitative analysis Methods 0.000 abstract description 2
- 238000011158 quantitative evaluation Methods 0.000 abstract 1
- 230000008859 change Effects 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 235000013399 edible fruits Nutrition 0.000 description 3
- 239000008235 industrial water Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- JKFYKCYQEWQPTM-UHFFFAOYSA-N 2-azaniumyl-2-(4-fluorophenyl)acetate Chemical compound OC(=O)C(N)C1=CC=C(F)C=C1 JKFYKCYQEWQPTM-UHFFFAOYSA-N 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 229910021612 Silver iodide Inorganic materials 0.000 description 1
- 230000003416 augmentation Effects 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 229940045105 silver iodide Drugs 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Economics (AREA)
- Computational Mathematics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Marketing (AREA)
- General Engineering & Computer Science (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Probability & Statistics with Applications (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法,运用物元分析理论,选取指标建立判断矩阵,并以此划分经典域和节域,通过关联度大小作为定性评价社会经济效益的数据来源;同时构建适应河源干旱区特征的C‑D生产函数,利用能值理论将不同要素全部转化为太阳能,若不同要素之间存在共线性,则利用岭回归解决要素间共线性的问题,最终根据所得函数求出年际水资源效益具体值,作为定量评价结果;分析研究区水资源挖潜和高效利用情况,分析投入和产出,得到研究区社会经济效益核算结果。本发明运用AHP和熵权法计算指标权重,克服主观性影响,定性和定量相结合避免了单一评价存在的偶然性问题,使得核算结果更加准确。
Description
技术领域
本发明涉及水资源利用领域,具体是一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法。
背景技术
近年来随着生产力的不断进步,对河源干旱区水资源挖潜和高效利用的研究不断深入,一定程度上缓解了地区所面临的水资源压力。在河源干旱区开展降雨径流挖潜,主要途径是利用人工增雨技术,通过补充某些形成降雨的必要条件,提高降雨的发生概率和降雨量,从而达到增加降雨量的效果;高效利用指提高水资源利用效率。由于地区土质特性造成田间渗漏严重,在田间灌区系统中,渗漏量和蒸发量占到了田间用水总量的65%以上,采用陶瓷滴灌等节水灌溉技术减少农作物无效蒸发以及减少渗漏量,构建坡/塬面径流高效利用技术体系,运用逐级汇集、统一调节的坝系水利用技术,可以有效提高水资源利用效率。
当前对于社会经济效益评价方法主要集中在以下几个:一般均衡理论;影子价格;模糊数学法;能值理论;边际机会成本等,但都存在不足。如一般均衡理论需要的数据量过大,针对河源区及干旱区而言本就数据量缺少,使用此方法不利于得到准确的结果;影子价格是对稀缺状况下水资源与社会经济效益之间关系做出说明,使用情况单一;模糊数学法虽然使用非常方便,但其指标权重的确定方法不统一,易带来结果的不确定性;能值理论存在计算能值转化率计算时,不同地区量值不一,不具有代表性;边际机会成本在进行效益核算时没有相对性且研究较为片面。
针对上述降雨径流挖潜与水资源高效利用的现状,本发明主要运用物元分析和生产函数来进行地区水资源挖潜和高效利用效果的社会经济效益核算,从定性和定量的角度来分析,使得到的计算结果更为准确。物元分析在指标权重计算部分,运用AHP法和熵权法耦合的方式得到计算结果,在利用生产函数定量核算社会经济效益时,运用能值方法将其全部转化为太阳能值,解决指标无法统一的问题。
发明内容
本发明的目的在于结合河源干旱区水资源挖潜和高效利用社会经济效益的影响因素,通过建立相应的指标体系,运用物元法和生产函数提出一种河源干旱区水资源挖潜和高效利用社会经济效益核算方法,能够为降水径流挖潜和水资源高效利用提供合理的评价方法,可为水资源管理提供参考。
本发明提供一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法,包括
步骤一、基于物元分析理论对河源干旱区水资源的社会经济效益进行定性评价,具体步骤如下:
(1)确定社会经济效益评价指标及物元分析变量;
选取水质(X1)、水贫困指数(X2)、人均GDP(X3)、农田灌溉亩均用水量(X4)、人均生活用水量(X5)、万元GDP用水量(X6)、水利相关行业投资(X7)总计七个指标建立水资源社会经济效益评价指标体系;物元分析中所描述的对象N、特征变量C和特征量值X组成社会经济效益评价物元R=(N、C、X),对象N为地区水资源社会经济效益,特征变量为所述七个指标;
(2)确定水资源社会经济效益评价的经典域物元矩阵和节域物元矩阵作为评价矩阵;
根据水资源社会经济效益评价的可拓性,将其划分为N01,N02,N03,N04,N05个等级,即Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ级,由此建立水资源社会经济效益评价的经典域物元矩阵R01、R02、R03、R04、R05和节域物元矩阵RP作为评价矩阵;
(3)利用AHP和熵权法耦合确定指标权重;
通过步骤(1)和步骤(2)得到的评价矩阵,依据特征量值X的取值,利用AHP法和熵权法计算得到相应的指标权重;
(4)确定综合关联度及水资源社会经济效益趋势分析
运用各项指标属于各质量等级的关联度公式和待评单元属于各质量等级的综合关联度公式得到相应的评价指标关联度计算结果,进而得到相应的评价等级,得到社会经济效益的趋势;
步骤二、基于生产函数对河源干旱区水资源的社会经济效益进行定量核算,具体步骤如下:
(1)构建相应产业效益核算的C-D生产函数
将研究区水资源的社会经济效益分为农牧业、工业和服务业三大类,并以此为研究对象构建对应的效益核算C-D生产函数,农牧业方面选取农牧业增加值、固定资产投资、劳动力、用水总量四个指标;工业方面选取工业增加值、固定资产投资、劳动力、用水总量四个指标;服务业方面选取第三产业增加值、固定资产投资、劳动力、用水总量四个指标;
C-D生产函数的基本表达式为
Q=A0(1+λ)t KαLβ
式中:Q为产出,A0为常数;λ为技术进步系数;t为年份序列;K为资本投入;L为劳动力投入;
考虑了水资源的C-D生产函数表示为:
Q=A0(1+λ)t KαLβWγ
式中:W为GDP单位用水量;γ为水资源产出弹性,将上式对数线性化,即
ln Q=ln A0+t ln(1+λ)+αln K+βln L+γln W
(2)依据能值理论进行指标统一量化,利用能值具有将所有指标全部转化为太阳能的特性,将上述四个指标全部转化为统一度量下的数值;
(3)考虑投入要素间的相关性
在出现|XTX|≈0,(XTX)-1不存在时给矩阵XTX加上一个正常数矩阵CI,其中C>0,I为单位矩阵,使得改变后的矩阵XTX+CI较矩阵XTX的奇异程度大幅度减小,C为岭参数,即采用岭回归分析进行修正,对包括时间变量在内的五个指标做岭回归分析,设置岭参数C在[0,1]之间变化,每次增加步长0.02,根据其他参数的稳定情况得到最终C的值;
(4)核算单位水对应的社会经济效益
将C-D生产函数线性化后求W的偏导数确定单位产出的边际效益BW:
(5)计算准则层的指标权重,并计算单方水综合效益
将水资源效益分为农牧业、工业、服务业三大类;权重分别为D1、D2、D3,则综合效益为
BW总=BW1D1*BW2D2*BW3*D3
步骤三、水资源挖潜和高效利用效果的社会经济效益核算,具体步骤如下:
(1)挖潜
某地区进行人工增雨后的技术参数为:假设原有的年平均水资源量为X1,经过挖潜之后年平均水资源量为X2,则有增加的水资源量为△X,水资源在运输和储存过程中存在一定的蒸散发和管损,设定蒸散发系数为ζ,管损系数为μ,同时考虑进行挖潜这一措施的成本为C1,则增加的水资源效益为Y1:
(2)高效利用
某地区运用水资源高效利用技术后的技术参数为:假设原有地区年平均水资源利用量为Z1,经过高效利用之后,水资源使用效率提高了X%,高效利用技术革新的成本为C2,则增加的水资源效益为Y2:
本发明具有如下有益效果:
本发明基于物元分析理论和C-D生产函数等,提出了河源干旱区水资源挖潜和高效利用效果的社会经济效益核算方法;构建了对应的核算指标和核算模型;并以此开展了相关核算工作,充分考虑了地区特异性以及地区现有的发展现状,运用能值理论克服了过往研究中存在的指标无法统一量化的困难,定性和定量相结合的方法,使得结果更加准确可靠,对于后续开展此类研究提供科学依据,进一步完善了水资源社会经济效益核算在河源干旱区的研究。
附图说明
图1是本发明河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法其中一个实施例的技术路线图;
图2是本发明河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法的流程图;
图3是本发明实施例中农牧业方面单方水效益曲线图;
图4是本发明实施例中工业方面单方水效益曲线图;
图5是本发明实施例中服务业方面单方水效益曲线图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明针对研究区的特点选取相应的指标,利用物元分析理论对河源干旱区水资源效益做出定性评价,得到研究区某一时间段内水资源社会经济效益状况趋势分析,而后运用生产函数对效益做出定量核算,得到水资源效益的核算函数和效益随时间的变化图,同时依据挖潜和高效利用两个要素对水资源所产生的不同影响,得到对应的效益核算结果。
请参阅图1及图2,本发明河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法其中一个实施例的具体技术方案如下:
步骤一:基于物元分析理论对河源干旱区水资源的社会经济效益进行定性评价
首先介绍物元分析的理论基础:
物元分析理论以研究处理矛盾问题的思维过程并将其数学形式化为核心。“物元理论”是可拓理论的两大重要支柱之一,和“可拓数学”构成了其理论体系的“硬核”。物元,是把事物的质和量有机联系在一起的重要概念,也是物元理论的首要基本概念,“事物”、“特征”及“量值”构成了物元三要素。从数学角度,给定事物的名称N,它关于特征C的量值为X,以有序3元R=(N,C,X)组作为描述事物特征的基本元,简称物元。
河源干旱区水资源社会经济效益综合评价中质量等级的归属问题实质上是一个典型的矛盾问题,利用物元分析法可以建立事物多指标性能参数的效益评定模型,以定量的数值较完整地反映效益的综合水平。
物元分析方法步骤如下:
①建立物元矩阵。根据事物特征选取多个特征量实测值,建立相应物元矩阵R0(待评环境单元),则
式中,P0为待评单元;Ci为待评单元的第i项特征,i=1,2,…,n;Xi为关于P0的量值,即对待评单元第i项特征进行分析的原始数据。
②划定经典域和节域。根据国家标准或拟定新的准则,建立经典域和节域矩阵。经典域矩阵为:
式中,Nj为第j个质量等级;Ci为质量等级Nj的第i个特征量,i=1,2,…,n;Xji为第j个质量等级中关于第i项特征(C)所规定的量值范围—经典域,节域矩阵为:
式中,P为质量等级的全体;Xpi为P关于Ci所取的量值范围;Ci为质量评价的第i个特征量。
③计算关联函数。首先确定各项指标属于各质量等级的关联度:
P(Xi,Xji)=|Xi-(aji+bji)/2|-(bji-aji)/2 (5)
P(Xi,Xpi)=|Xi-(api+bpi)/2|-(bpi-api)/2 (6)
|Xji|=|aji-bji| (7)
Kj(Xi)为第i项特征(指标)属于第j级的关联度,然后确定待评单元属于各质量等级的综合关联度:
式中,Wi为第i项特征的权重;Kj(P0)为待评单元(P0)属于第j级的综合关联度。
④待评单元质量等级评定。若Kj=max{Kj(P0)},j=1,2,…,m,则P0质量等级为第j级.当0<Kj(P0)<1时,表明评价单元符合标准对象的要求,其值愈大,符合程度愈好.当-1<Kj(P0)<0时,表示评价单元不符合某级评价标准的要求,但具备转化为标准对象的条件,其值愈小,愈易转化。当Kj(P0)<-1时,表示评价单元不符合某级评价标准的要求,而且又不具备转化为标准对象的条件。
所述步骤一具体步骤如下:
步骤1.1:确定社会经济效益评价指标及物元分析变量;
选取水质(X1)、水贫困指数(X2)、人均GDP(X3)、农田灌溉亩均用水量(X4)、人均生活用水量(X5)、万元GDP用水量(X6)、水利相关行业投资(X7)总计七个指标建立水资源社会经济效益评价指标体系;物元分析中所描述的对象N、特征变量C和特征量值X组成社会经济效益评价物元R=(N、C、X),对象N为地区水资源社会经济效益,特征变量为所述七个指标;
步骤1.2:确定水资源社会经济效益评价的经典域物元矩阵和节域物元矩阵作为评价矩阵;
根据水资源社会经济效益评价的可拓性,将其划分为N01,N02,N03,N04,N05个等级,即Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ级,由此建立水资源社会经济效益评价的经典域物元矩阵R01、R02、R03、R04、R05和节域物元矩阵RP作为评价矩阵;
步骤1.3:利用AHP和熵权法耦合确定指标权重(1)AHP法计算指标权重
①构造判断矩阵A
根据3标度理论,对指标进行两两比较,构造判断矩阵A为:
其中:
②计算判断矩阵的最优传递矩阵
最优传递矩阵R为:
式中:dik=exp(rik),因其具有良好的一致性,故又叫一致性矩阵。
③由判断矩阵D,计算权重系数
根据判断矩阵,计算出该层各指标关于上层次某指标的优先权重,运用乘积方根法计算出特征向量的值,即:
ω=[ω1,ω2,…ωn]T (13)
式中:ω=[ω1,ω2,…ωn]T为特征向量,作为该层次n个指标的权重向量。
(2)熵权法确定指标权重
①标准化处理评价指标
在选定的评价指标中,既有产生正面影响的指标,即正向指标;又有产生负面影响的指标,即负向指标,因此需要标准化处理初始数据使其具备可比性。
对于正向指标,标准化处理公式为:
对于负向指标,归一化处理公式为:
②熵值的确定
熵值由下式计算得到:
式中:Hi为评价指标i的熵值,在此假定当fij=0时,fij ln fij=0;gij为原始评价指标的归一化化矩阵,k为玻尔兹曼常量,n为研究对象的总个数。熵值越小,说明该指标对于决策者来说越能提供有用的信息。
③确定评价指标的熵权
在得到熵值H,之后,通过下式确定评价指标的熵权:
(3)确定综合权重
将上述确定的AHP法权重值ωi和熵权法权重值ci进行权重耦合,则综合权重为Wi:
步骤1.4、确定综合关联度及水资源社会经济效益趋势分析
首先利用上述依据评价指标计算得到的综合关联度,然后根据待评单元质量等级评定规则,通过每一年计算数据对应的定性评价结果,分析相对应的社会经济效益变化趋势。
步骤二、基于生产函数对河源干旱区水资源的社会经济效益进行定量核算,具体步骤如下:
步骤2.1:构建相应产业效益核算的C-D生产函数
经济学增长理论认为,经济增长的主要因素是各生产要素的增加,经济增长取决于各生产要素的贡献。某生产要素对经济发展的贡献率一般指该项生产要素增长对经济增长的直接影响。生产函数就是描述生产过程中投入的生产要素和其他要素组合所产生的可能最大产出量之间的依存关系的数学关系表达式,一般记为
Q=f(X1,X2...Xn) (22)
式中:Q为产出;X1,X2,…,Xm为投入的各生产要素。
C-D生产函数的基本表达式为
Q=A0(1+λ)tKαLβ (23)
式中:A0为常数;λ为技术进步系数;t为年份序列;K为资本投入;L为劳动力投入。
生产函数模型中投入量和产出量之间的依存关系普遍存在于各种生产过程之中。水资源是国民经济可持续发展的重要资源基础,也是重要的生产要素,所以在生产函数中不仅要考虑劳动力投入L和资本投入K这2种生产要素,也有必要将水资源视为一种生产要素,将产品中的“水量投入”(非水的经济投入)列为生产函数的自变量之一,由此分析水资源对经济发展的贡献。
考虑了水资源的C-D生产函数可表示为
Q=A0(1+λ)tKαLβWγ (24)
式中:W为GDP单位用水量;γ为水资源产出弹性。将上式对数线性化,即
ln Q=ln A0+t ln(1+λ)+αln K+βln L+γln W (25)
这样就确立了以水资源作为一种投入要素的生产函数,可以利用水资源产出弹性求水资源的边际效益。
步骤2.2:基于能值理论的指标统一量化处理
从上述生产函数可以看出,不同的投入要素间量纲并不统一,而能值理论则做到很好的解决这一问题,其将所有要素都转化为太阳能J,让不同类指标全部转化为J这一统一度量。
EM=τ×B (26)
式中:EM为能值(sej);τ为能值转化率(sej/J或sej/g);B为能量或物质的质量(J或g)
步骤2.3:考虑投入要素间的相关性
运用多元线性回归最重要的就是不同指标之间的互不相关,若存在相关性,此时使用最小二乘估计计算无法得到有效的估计系数,模型预测没有任何意义。因此需要采用岭回归来客服多重共线的问题。岭回归就是在出现|XTX|≈0,(XTX)-1不存在时给矩阵XTX加上一个正常数矩阵CI,使得改变后的矩阵相较原来的奇异程度大幅度减小,C为岭参数。
步骤2.4:核算单位水对应的社会经济效益
弹性分析可以明确表示各种要素对国民产出的影响程度,所以可以利用弹性系数求得水资源对各产业的边际效益。将上式线性化后求W的偏导数就可以确定单位产出的边际效益BW。
弹性是因变量的相对变化与自变量相对变化的比值,其绝对值大小可显示出经济因变量相对变动对经济自变量相对变动的反映程度,其符号的正负可反映出水资源投入量追加时对经济增长影响的趋势,计算效益时取其绝对值。
步骤2.5:计算准则层的指标权重,并计算单方水综合效益
将水资源效益分为农牧业、工业、服务业三大类;权重分别为D1、D2、D3,则综合效益为
BW总=BW1D1*BW2D2*BW3*D3 (28)
步骤三、水资源挖潜和高效利用效果的社会经济效益核算,具体步骤如下:
(1)挖潜
根据上述背景资料显示,某地区进行人工增雨后的技术参数为:假设原有的年平均水资源量为X1,经过挖潜之后年平均水资源量为X2,则有增加的水资源量为△X,水资源在运输和储存过程中存在一定的蒸散发和管损,设定蒸散发系数为ζ,管损系数为μ,同时考虑进行挖潜这一措施的成本为C1,则增加的水资源效益为Y1。
(2)高效利用
根据上述背景技术资料显示,某地区运用水资源高效利用技术后的技术参数为:假设原有地区年平均水资源利用量为Z1,经过高效利用之后,水资源使用效率提高了X%,高效利用技术革新的成本为C2,则增加的水资源效益为Y2。
本发明充分考虑水资源挖潜和高效利用两个因素,以水资源的社会经济效益来表征水资源挖潜和高效利用的效果,构建一种河源干旱区域水资源挖潜与高效利用效果的社会经济效益的核算方法,可为河源干旱区水资源挖潜与高效利用评价提供支撑。
实例分析
本发明以果洛藏族自治州作为研究区为实施例,以验证本发明的效果。具体技术方案如下:
步骤一、基于物元分析理论对河源干旱区水资源的社会经济效益进行定性评价
(1)确定社会经济效益评价指标及物元分析变量
针对研究区的特点,选取水质(X1)、水贫困指数(X2)、人均GDP(X3)、农田灌溉亩均用水量(X4)、人均生活用水量(X5)、万元GDP用水量(X6)、水利相关行业投资(X7)总计七个指标建立水资源社会经济效益评价指标体系。
物元分析中所描述的对象N、特征变量C和特征量值X组成社会经济效益评价物元R=(N、C、X),在本例中对象N为地区水资源社会经济效益,特征变量为水质、水贫困指数等七个指标,特征量值X为表格中的具体数值。采用2013-2019年的数据进行区域水资源社会经济效益评价,见表1。
表1研究区2013-2019评价指标取值
(2)确定水资源社会经济效益评价的经典域和节域
根据水资源社会经济效益评价的可拓性,将其划分为N01,N02,N03,N04,N05个等级,即Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ级。由此建立水资源社会经济效益评价的经典域物元矩阵R01、R02、R03、R04、R05和节域物元矩阵RP分别为:
(3)通过熵权法和AHP法确定指标权重
通过步骤(1)和步骤(2)得到评价矩阵,依据表1的数据利用AHP法和熵权法计算得到相应的指标权重,通过熵权法和AHP法得到指标权重如下:
Wi=(0.1147,0.0655,0.1719,0.0918,0.2748,0.1912,0.0901)
(4)确定综合关联度及水资源社会经济效益趋势分析
运用公式4(各项指标属于各质量等级的关联度公式)和公式8(待评单元属于各质量等级的综合关联度公式)得到相应的评价指标关联度计算结果,结果如下表2所示:
表2水资源社会经济效益评价结果
利用上表的综合关联度具体数据,得到相应的评价等级,可以看出社会经济效益的趋势,从而定性评价研究区的水资源社会经济效益,根据所得结果显示,研究区域水资源社会经济效益呈现逐步上升趋势。
步骤二、基于生产函数对河源干旱区水资源的社会经济效益进行定量核算
(1)构建相应产业效益核算的C-D生产函数
将研究区水资源的社会经济效益分为农牧业、工业和服务业三大类,并以此为研究对象构建对应的效益核算C-D生产函数。农牧业方面选取农牧业增加值、固定资产投资、劳动力、用水总量四个指标;工业方面选取工业增加值、固定资产投资、劳动力、用水总量四个指标;服务业方面选取第三产业增加值、固定资产投资、劳动力、用水总量四个指标,具体数据见表3,表4,表5。
表3 2009-2019果洛州农牧业水资源社会经济效益核算指标
表4 2009-2019果洛州工业水资源社会经济效益核算指标
表5 2009-2019果洛州服务业水资源社会经济效益核算指标
(2)依据能值理论进行指标统一量化,利用能值具有将所有指标全部转化为太阳能的特性,将上述四个指标全部转化为统一度量下的数值,因本例中存在有一定的共线性,故采用岭回归分析进行修正,对包括时间变量在内的五个解释变量做岭回归分析,设置岭参数C在[0,1]之间变化,每次增加步长0.02,根据其他参数的稳定情况得到最终C的值。
经过上述步骤可得:
Q农牧业=3.88(1+0.06)tK-0.02L-0.16W-0.13;
Q工业=1.93(1+0.08)tK0.415L0.255W0.107;
Q服务业=2.85(1+0.11)tK0.728L0.045W0.026。
由此可知,农牧业方面每增加1%的产值中有0.13%是由水资源贡献的;同理,工业方面有0.107%;服务业方面有0.026%,根据上式27可以计算出对应产业的单方水效益。根据所得结果绘制如图3-5所示效益曲线,由图3-5可以看出各个行业水资源效益都在上升,与物元分析所得结果一致。
(3)考虑投入要素间的相关性
在出现|XTX|≈0,(XTX)-1不存在时给矩阵XTX加上一个正常数矩阵CI,使得改变后的矩阵相较原来的奇异程度大幅度减小,C为岭参数,即采用岭回归分析进行修正,对包括时间变量在内的五个指标做岭回归分析,设置岭参数C在[0,1]之间变化,每次增加步长0.02,根据其他参数的稳定情况得到最终C的值;
(4)核算单位水对应的社会经济效益
将C-D生产函数线性化后求W的偏导数确定单位产出的边际效益BW:
(5)计算准则层的指标权重,并计算单方水综合效益
将水资源效益分为农牧业、工业、服务业三大类;权重分别为D1、D2、D3,则综合效益为
BW总=BW1D1*BW2D2*BW3*D3
步骤三、水资源挖潜和高效利用效果的社会经济效益核算
根据现有研究进展得到的数据显示,果洛地区进行声波和传统碘化银人工增雨效果为增加地区可用水量的10%,同时通过滴灌和其他节水措施并行,减少地区用水量在8%左右。
首先根据上式(27)计算出地区单方水资源效益为:5.72元/m3;81.97元/m3;7.77元/m3;利用上面计算得到的权重,三者分别为0.6491;0.0719;0.2790,得到综合效益BW总为3.925元/m3。可以看出地区工业水资源价值较高,因为地区主要工业产值来自于青海威斯特铜业,占比超过了80%,随着技术的不断改进,用水量效率提高,这也就导致了地区工业水资源效益较高。再利用式(28)、(29)可以得到地区通过水资源挖潜和高效利用增加的效益为:
挖潜:
高效利用:
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何属于本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (3)
1.一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法,其特征在于:包括如下步骤:
步骤一、基于物元分析理论对河源干旱区水资源的社会经济效益进行定性评价,具体步骤如下:
(1)确定社会经济效益评价指标及物元分析变量;
选取水质(X1)、水贫困指数(X2)、人均GDP(X3)、农田灌溉亩均用水量(X4)、人均生活用水量(X5)、万元GDP用水量(X6)、水利相关行业投资(X7)总计七个指标建立水资源社会经济效益评价指标体系;物元分析中所描述的对象N、特征变量C和特征量值X组成社会经济效益评价物元R=(N、C、X),对象N为地区水资源社会经济效益,特征变量为所述七个指标;
(2)确定水资源社会经济效益评价的经典域物元矩阵和节域物元矩阵作为评价矩阵;
根据水资源社会经济效益评价的可拓性,将其划分为N01,N02,N03,N04,N05个等级,即Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ级,由此建立水资源社会经济效益评价的经典域物元矩阵R01、R02、R03、R04、R05和节域物元矩阵RP作为评价矩阵;
(3)利用AHP和熵权法耦合确定指标权重;
通过步骤(1)和步骤(2)得到的评价矩阵,依据特征量值X的取值,利用AHP法和熵权法计算得到相应的指标权重;
(4)确定综合关联度及水资源社会经济效益趋势分析
运用各项指标属于各质量等级的关联度公式和待评单元属于各质量等级的综合关联度公式得到相应的评价指标关联度计算结果,进而得到相应的评价等级,得到社会经济效益的趋势;
步骤二、基于生产函数对河源干旱区水资源的社会经济效益进行定量核算,具体步骤如下:
(1)构建相应产业效益核算的C-D生产函数
将研究区水资源的社会经济效益分为农牧业、工业和服务业三大类,并以此为研究对象构建对应的效益核算C-D生产函数,农牧业方面选取农牧业增加值、固定资产投资、劳动力、用水总量四个指标;工业方面选取工业增加值、固定资产投资、劳动力、用水总量四个指标;服务业方面选取第三产业增加值、固定资产投资、劳动力、用水总量四个指标;
C-D生产函数的基本表达式为Q=A0(1+λ)tKαLβ
式中:Q为产出,A0为常数;λ为技术进步系数;t为年份序列;K为资本投入;L为劳动力投入;
考虑了水资源的C-D生产函数表示为:
Q=A0(1+λ)tKαLβWγ
式中:W为GDP单位用水量;γ为水资源产出弹性,将上式对数线性化,即
ln Q=ln A0+tln(1+λ)+αln K+βln L+γln W
(2)依据能值理论进行指标统一量化,利用能值具有将所有指标全部转化为太阳能的特性,将上述四个指标全部转化为统一度量下的数值;
(3)考虑投入要素间的相关性
在出现|XTX|≈0,(XTX)-1不存在时给矩阵XTX加上一个正常数矩阵CI,其中C>0,I为单位矩阵,使得改变后的矩阵XTX+CI较矩阵XTX的奇异程度大幅度减小,C为岭参数,即采用岭回归分析进行修正,对包括时间变量在内的五个指标做岭回归分析,设置岭参数C在[0,1]之间变化,每次增加步长0.02,根据其他参数的稳定情况得到最终C的值;
(4)核算单位水对应的社会经济效益
将C-D生产函数线性化后求W的偏导数确定单位产出的边际效益BW:
(5)计算准则层的指标权重,并计算单方水综合效益
将水资源效益分为农牧业、工业、服务业三大类;权重分别为D1、D2、D3,则综合效益为
BW总=BW1D1*BW2D2*BW3*D3
步骤三、水资源挖潜和高效利用效果的社会经济效益核算,具体步骤如下:
(1)挖潜
某地区进行人工增雨后的技术参数为:假设原有的年平均水资源量为X1,经过挖潜之后年平均水资源量为X2,则有增加的水资源量为△X,水资源在运输和储存过程中存在一定的蒸散发和管损,设定蒸散发系数为ζ,管损系数为μ,同时考虑进行挖潜这一措施的成本为C1,则增加的水资源效益为Y1:
(2)高效利用
某地区运用水资源高效利用技术后的技术参数为:假设原有地区年平均水资源利用量为Z1,经过高效利用之后,水资源使用效率提高了X%,高效利用技术革新的成本为C2,则增加的水资源效益为Y2:
2.如权利要求1所述的河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法,其特征在于,步骤一中利用AHP和熵权法耦合确定指标权重具体为:
(1)AHP法计算指标权重
①构造判断矩阵A
根据3标度理论,对指标进行两两比较,构造判断矩阵A为:
其中:
②计算判断矩阵的最优传递矩阵
最优传递矩阵R为:
式中:dik=exp(rik),因其具有良好的一致性,故又叫一致性矩阵;
③由判断矩阵D,计算权重系数
根据判断矩阵,计算出该层各指标关于上层次某指标的优先权重,运用乘积方根法计算出特征向量的值,即:
ω=[ω1,ω2,…ωn]T
式中:ω=[ω1,ω2,…ωn]T为特征向量,作为该层次n个指标的权重向量;
(2)熵权法确定指标权重
①标准化处理评价指标
对于正向指标,标准化处理公式为:
对于负向指标,归一化处理公式为:
②熵值的确定
熵值由下式计算得到:
式中:Hi为评价指标i的熵值,在此假定当fij=0时,fij ln fij=0;gij为原始评价指标的归一化化矩阵,k为玻尔兹曼常量,n'为研究对象的总个数;
③确定评价指标的熵权
在得到熵值H,之后,通过下式确定评价指标的熵权:
(3)确定综合权重
将上述确定的AHP法权重值ωi和熵权法权重值ci进行权重耦合,则综合权重为Wi:
3.如权利要求1所述的河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法,其特征在于,步骤二中依据能值理论进行指标统一量化具体为:
EM=τ×B
式中:EM为能值,单位为sej;τ为能值转化率,单位为sej/J或sej/g;B为能量或物质的质量,单位为J或g。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011451687.2A CN112633645B (zh) | 2020-12-10 | 2020-12-10 | 一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011451687.2A CN112633645B (zh) | 2020-12-10 | 2020-12-10 | 一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112633645A CN112633645A (zh) | 2021-04-09 |
CN112633645B true CN112633645B (zh) | 2023-05-02 |
Family
ID=75309768
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011451687.2A Active CN112633645B (zh) | 2020-12-10 | 2020-12-10 | 一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112633645B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113344440B (zh) * | 2021-06-30 | 2023-05-23 | 西北农林科技大学 | 一种地下水承载力评估模型的构建方法 |
CN113723704B (zh) * | 2021-09-07 | 2023-04-18 | 东莞理工学院 | 一种基于连续和分级混合数据的水质快速预测方法 |
CN116307936B (zh) * | 2023-05-17 | 2023-08-18 | 长江水利委员会长江科学院 | 用水总量的分解方法、系统、电子设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014142907A (ja) * | 2012-12-25 | 2014-08-07 | Nippon Telegr & Teleph Corp <Ntt> | 経済被害額評価装置及び方法及びプログラム |
CN109243178A (zh) * | 2018-11-14 | 2019-01-18 | 上海应用技术大学 | 一种不良气候条件下城镇公路交通安全分析与评价方法 |
CN109377093A (zh) * | 2018-11-22 | 2019-02-22 | 国网湖南省电力有限公司 | 配电网项目投资效益评价方法 |
-
2020
- 2020-12-10 CN CN202011451687.2A patent/CN112633645B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014142907A (ja) * | 2012-12-25 | 2014-08-07 | Nippon Telegr & Teleph Corp <Ntt> | 経済被害額評価装置及び方法及びプログラム |
CN109243178A (zh) * | 2018-11-14 | 2019-01-18 | 上海应用技术大学 | 一种不良气候条件下城镇公路交通安全分析与评价方法 |
CN109377093A (zh) * | 2018-11-22 | 2019-02-22 | 国网湖南省电力有限公司 | 配电网项目投资效益评价方法 |
Non-Patent Citations (1)
Title |
---|
基于物元可拓模型的区域用水效率控制综合评价;王晓琳等;水电能源科学;第35卷(第3期);30-34 * |
Also Published As
Publication number | Publication date |
---|---|
CN112633645A (zh) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112633645B (zh) | 一种河源干旱区水资源挖潜与高效利用效果的社会经济效益核算方法 | |
CN102183621B (zh) | 水产养殖溶解氧浓度在线预测方法及系统 | |
CN108876209A (zh) | 一种考虑相对产量的红壤稻田土壤肥力评价方法 | |
CN113377966B (zh) | 一种基于知识图谱的水利工程调度规程推理方法 | |
CN111462450B (zh) | 一种考虑降雨空间异质性的山洪预警方法 | |
CN111984700B (zh) | 一种基于水热耦合平衡原理的月径流预测方法 | |
CN113592278A (zh) | 考虑非期望产出的sbm水环境承载力评价方法 | |
Song et al. | Catastrophe progression method based on MK test and correlation analysis for assessing water resources carrying capacity in Hubei province | |
Wu et al. | Assessment of water resources carrying capacity based on fuzzy comprehensive evaluation–case study of Jinan, China | |
Wang et al. | Study on the spatial distribution of water resource value in the agricultural system of the Yellow River Basin | |
CN109002946B (zh) | 一种河湖补水的“两库—两站”系统水资源优化调度方法 | |
Zhou et al. | Regional water resource security in China based on a new fuzzy method with combination weighting | |
Chen et al. | Analysis of the characteristics and driving forces of water footprint productivity in paddy rice cultivation in China | |
CN105096216A (zh) | 一种快速计算水电站发电量的方法 | |
CN115439022B (zh) | 养殖场粪污还田后环境污染的风险评估方法及系统 | |
CN115358587A (zh) | 区域多部门协同的基础设施规划方法及系统 | |
CN115496382A (zh) | 一种区域农业生态发展动态评估方法、系统及可存储介质 | |
Wu et al. | Research on water environmental capacity accounting of the Yongzhou Section of Xiangjiang River Basin based on the SWAT-EFDC coupling model | |
CN108364136B (zh) | 一种基于证据推理的水资源短缺风险分析方法及系统 | |
Yanjuan et al. | Exploitation, utilization status, and development prospects of water resources in Ningbo City | |
CN112906313A (zh) | 基于系统协同理论的生态流量计算方法 | |
CN112801339A (zh) | 一种基于电力大数据的地区经济发展预测方法 | |
Shi et al. | Evaluation and drivers of green agricultural water use efficiency in Yangtze River Economic Zone | |
CN110187072A (zh) | 一种基于贴进度模型的水环境质量评价方法 | |
CN115421217B (zh) | 一种基于降雨预报信息的洪水过程及其类型预报方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |