水環(huán)境治理PPP項目的多主體行為演化博弈

摘要：基于水環(huán)境治理PPP項目多主體不同策略組合下的收益函數(shù)，構(gòu)建政府主體、經(jīng)濟主體和社會主體多主體演化博弈模型，并通過MATLAB軟件進(jìn)行數(shù)值模擬和參數(shù)敏感性分析.結(jié)果表明：在一定條件下復(fù)制動態(tài)系統(tǒng)可以達(dá)到積極規(guī)制，自主推進(jìn)，主動參與的理想穩(wěn)定狀態(tài)；激勵措施可以加快自主推進(jìn)、主動參與策略收斂速度，但激勵強度并非越強越好；演化穩(wěn)定策略受懲罰力度和損失系數(shù)γ影響，對損失系數(shù)β不敏感.

關(guān)鍵詞：水環(huán)境治理； PPP模式； 演化博弈； 利益相關(guān)者； 政府規(guī)制

中圖分類號： X 522文獻(xiàn)標(biāo)志碼： A 文章編號： 1000－5013（2023）02－0233－10

Multi－Agent Behavior Evolutionary Game of Water Environment Treatment PPP Projects

SONG Wenjie， ZHAN Zhaoxi

（College of Civil Engineering， Huaqiao University， Xiamen 361021， China）

Abstract： Based on the income function of water environment treatment PPP projects under different strategy combinations of multi－agent， a multi－agent evolutionary game model was established among government agent， economic agent and social agent， and numerical simulation and parameter sensitivity analysis were carried out by MATLAB software. The results show that the replication dynamic system can achieve the ideal stable state of active regulation， independent promotion and active participation under certain conditions. The incentive measures can accelerate convergence speed of autonomous promotion and active participation strategies， but the incentive intensity is not as better as it gets stronger. The evolution stability strategy is affected by penalty intensity and loss coefficient γ， but it isn′t insensitive to loss coefficient β.

Keywords： water environment treatment； PPP model； evolutionary game； stakeholders； government regulation

在構(gòu)建人與自然生命共同體背景下，生態(tài)修復(fù)、河道整治、污水治理等水環(huán)境綜合治理是當(dāng)前乃至未來我國環(huán)境保護(hù)的重要工作之一，但是巨額投資和復(fù)雜治理技術(shù)使項目治理陷入困境［1］.政府和社會資本合作（PPP）模式的應(yīng)用為水環(huán)境治理提供了新的融資途徑［2］，在應(yīng)用前沿科技［3］和先進(jìn)管理模式［4］的同時，破解了項目建設(shè)資金的瓶頸問題.然而，水環(huán)境治理PPP項目在實施過程中，仍存在項目參與主體不同利益訴求和利益沖突的共性，影響項目治理績效［5］.

現(xiàn)有研究指出，識別治理項目的主要利益相關(guān)者及其層次關(guān)系和互相影響程度［6－7］，并將利益相關(guān)者之間的利益訴求納入決策是提高治理成效和實現(xiàn)可持續(xù)性目標(biāo)的關(guān)鍵［8］.在水環(huán)境治理過程中，通過采用網(wǎng)絡(luò)化治理手段，強調(diào)了地方政府在水環(huán)境治理中的主導(dǎo)作用，但政府過多干預(yù)并不利于環(huán)境治理效益的提升［9］.為了加強共同管理、提升治理成效，應(yīng)權(quán)衡不同利益相關(guān)者之間的利益訴求［10］.有學(xué)者應(yīng)用演化博弈理論分析不同群體之間的互動機制［11］，探索利益主體利益博弈及策略選擇的動態(tài)過程.通過對政府間、企業(yè)間、政府與企業(yè)間等主體行為分析，認(rèn)為政府的監(jiān)管成本、力度以及懲罰值是影響系統(tǒng)穩(wěn)定策略的關(guān)鍵［12］，當(dāng)政府和社會資本積極合作獲得較高收益時，可以實現(xiàn)主體間合作共贏［13］.現(xiàn)有研究肯定了政府的主導(dǎo)作用和項目主體之間的博弈行為對治理績效的影響，但是忽視了多元利益主體之間的行為博弈和影響因素分析，尤其是對公眾、營利性設(shè)施使用者等社會主體的利益訴求缺乏深入分析.社會主體作為特殊的利益相關(guān)者群體，是水環(huán)境治理效果的直接感受者和受益者［14］，在動態(tài)考量水環(huán)境治理PPP項目多主體間的利益關(guān)系與行為博弈時，理應(yīng)納入分析框架.當(dāng)前有少數(shù)學(xué)者對水環(huán)境治理PPP項目多主體行為博弈展開了研究，但也僅是從價值共創(chuàng)［15］角度去探索博弈主體對環(huán)境治理公私合作項目的影響，而如何將政府規(guī)制與外部機制相結(jié)合研究水環(huán)境治理PPP項目的演化博弈問題尚未涉及.基于此，本文對水環(huán)境治理PPP項目多主體在政府規(guī)制下策略選擇的作用機理和演化路徑進(jìn)行研究.

1 策略分析與模型假設(shè)

1.1 演化博弈主體與行動策略分析

水環(huán)境治理PPP項目主要涉及政府主體、經(jīng)濟主體、社會主體等3大主體，在不考慮其他約束情況下，將政府主體、經(jīng)濟主體和社會主體視為1個完整的系統(tǒng)，該系統(tǒng)中3個主體是有限理性的，不受外在因素影響，僅分析博弈主體的現(xiàn)有行為策略.

在水環(huán)境治理PPP項目實施過程中，政府發(fā)揮著主導(dǎo)作用，此主導(dǎo)作用有異于傳統(tǒng)政府的“大包大攬”，而是政府在解決水環(huán)境污染問題上擔(dān)負(fù)主要責(zé)任，是不可輕忽的主導(dǎo)力量，把控著項目的總體利益.

首先，政府是政策制度和規(guī)則的制定者.我國水環(huán)境污染防治和水環(huán)境保護(hù)自覺意識尚未形成，PPP模式引入水環(huán)境治理領(lǐng)域的時間也較短，強化規(guī)制制度建設(shè)對我國水環(huán)境治理PPP項目的推進(jìn)至關(guān)重要，而水環(huán)境治理PPP項目的復(fù)雜性也決定了此類項目需要制度和規(guī)則作為強有力的后盾.

其次，政府是監(jiān)督者和管理者.政府作為水環(huán)境治理PPP項目的主要參與主體，具有不可推卸的監(jiān)督與管理職責(zé).因此，在積極規(guī)制的情況下，政府主體通過環(huán)保宣傳、政策優(yōu)惠、專項資金或補貼等［16］激勵和扶持措施激發(fā)經(jīng)濟主體和社會主體的積極性和自主性.然而，政府主體通過利益權(quán)衡或者僅是為了政績需求，可能會處于消極規(guī)制狀態(tài)，未給予項目過多關(guān)注和經(jīng)濟補貼.

水環(huán)境治理PPP項目的經(jīng)濟主體主要為社會資本和項目公司，與政府主體構(gòu)成合作伙伴關(guān)系.作為典型的“經(jīng)濟人”，獲得合理的經(jīng)濟效益是經(jīng)濟主體參與治理項目的首要考慮因素.隨著社會環(huán)境保護(hù)意識的增強，社會認(rèn)可度和公司形象成為企業(yè)的核心競爭力之一，經(jīng)濟主體會根據(jù)項目成本、項目收益、風(fēng)險預(yù)測及政府支持力度投資水環(huán)境治理PPP項目.在政府規(guī)制下，經(jīng)濟主體為獲得更多補貼，積極落實各項制度政策，配合政府工作做好項目推進(jìn).然而，經(jīng)濟主體自身“經(jīng)濟人”的屬性可能會使其采取非合作行為或投機行為［17］，未按照合同約定或違背謹(jǐn)慎工程和運營慣例，被動地去完成水環(huán)境治理任務(wù).

水環(huán)境治理PPP項目的社會主體主要為周邊群眾和營利性設(shè)施使用者.水環(huán)境問題與社會主體聯(lián)系緊密，直接影響著周邊群眾的日常生活質(zhì)量，社會主體參與環(huán)境治理也是社會發(fā)展的必然趨勢，《環(huán)境保護(hù)公眾參與辦法》《環(huán)境影響評價法》《水污染防治法》等法律條文也明確了群眾參與環(huán)境影響評價、水污染防治等生態(tài)環(huán)境治理的權(quán)力.不同于政府主體和經(jīng)濟主體，社會主體在項目治理中既無確切的工作任務(wù)，也無實質(zhì)的資金投入.然而，社會主體作為水環(huán)境治理效果的直接感受者和受益者，為了維護(hù)賴以生存的自然環(huán)境，維護(hù)自身合法利益，社會主體會主動參與水環(huán)境治理項目，但也存在社會主體基于多一事不如少一事的原則而選擇置身事外的可能.

1.2 模型的基本假設(shè)

假設(shè)1：政府主體策略集合為{積極規(guī)制，消極規(guī)制}，經(jīng)濟主體策略集合為{自主推進(jìn)，被動履責(zé)}，社會主體策略集合為{主動參與，不參與}.分別用x，y，z表示政府主體、經(jīng)濟主體和社會主體選擇積極規(guī)制、自主推進(jìn)和主動參與策略的概率，x，y，z∈［0，1］且均為時間（t）的函數(shù)；分別用（1-x），（1-y），（1-z）表示政府主體、經(jīng)濟主體和社會主體選擇消極規(guī)制、被動履責(zé)和不參與策略的概率.

假設(shè)2：水環(huán)境治理PPP項目自身給政府主體帶來的基本經(jīng)濟效益為R1，基本成本為C1；政府主體選擇積極規(guī)制策略時，可獲得額外收益E1，需要額外付出成本B1（E1gt;B1）；對選擇自主推進(jìn)策略的經(jīng)濟主體提供稅收優(yōu)惠和經(jīng)濟補貼（激勵強度），記為J1；對選擇主動參與策略的社會主體給予獎勵，記為J2；對不配合政府工作的經(jīng)濟主體進(jìn)行處罰，記為K1；當(dāng)政府主體采取消極規(guī)制策略時，因經(jīng)濟主體被動履責(zé)而導(dǎo)致不可預(yù)見風(fēng)險發(fā)生而帶來的損失，記為S1.

假設(shè)3：水環(huán)境治理PPP項目自身給經(jīng)濟主體帶來的基本經(jīng)濟效益為R2，基本成本為C2；當(dāng)經(jīng)濟主體選擇自主推進(jìn)策略時，可獲得額外收益E2，需要額外付出成本B2（E2gt;B2）；當(dāng)社會主體選擇不參與策略，經(jīng)濟主體選擇自主推進(jìn)策略時給本方帶來的損失L1的概率（損失系數(shù)），記為γ；當(dāng)社會主體選擇主動參與策略，經(jīng)濟主體選擇被動履責(zé)策略時，會帶來損失S2.

假設(shè)4：水環(huán)境治理PPP項目自身給社會主體帶來的基本經(jīng)濟效益為R3；當(dāng)社會主體選擇主動參與策略時，可獲得額外收益E3，需要付出成本C3；當(dāng)經(jīng)濟主體選擇被動履責(zé)策略，社會主體選擇主動參與策略時給本方帶來的損失L2的概率，記為β；當(dāng)經(jīng)濟主體選擇自主推進(jìn)策略，社會主體選擇不參與策略時，會帶來損失S3.

考慮到一方主體因另外兩方主體采取正向行為時會獲得額外的經(jīng)濟效益，設(shè)n，m，g∈（0，1）分別為政府主體、經(jīng)濟主體、社會主體的收益系數(shù)，令政府主體因經(jīng)濟主體和社會主體都采取正向行為獲得的額外經(jīng)濟效益為n（R2+R3）；經(jīng)濟主體因政府主體和社會主體都采取正向行為獲得的額外經(jīng)濟效益為m（R1+R3）；社會主體因政府主體和經(jīng)濟主體都采取正向行為獲得的額外經(jīng)濟效益為g（R1+R2）.

2 多主體演化博弈模型的構(gòu)建

2.1 模型的構(gòu)建

基于假設(shè)1～4，水環(huán)境治理PPP項目中政府主體、經(jīng)濟主體、社會主體的策略選擇可形成8個策略組合，即G1＝{積極規(guī)制，自主推進(jìn)，主動參與}，G2＝{積極規(guī)制，自主推進(jìn)，不參與}，G3＝{積極規(guī)制，被動履責(zé)，主動參與}，G4＝{積極規(guī)制，被動履責(zé)，不參與}，G5＝{消極規(guī)制，自主推進(jìn)，主動參與}，G6＝{消極規(guī)制，自主推進(jìn)，不參與}，G7＝{消極規(guī)制，被動履責(zé)，主動參與}，G8＝{消極規(guī)制，被動履責(zé)，不參與}.結(jié)合模型的基本假設(shè)，可得多主體演化博弈收益情況，如表1所示.

在水環(huán)境治理PPP項目中，設(shè)政府主體采用積極規(guī)制策略時的期望收益為Ux1，采用消極規(guī)制策略時期望收益為U1-x1，平均期望收益為U1，則Ux1，U1-x1，U1分別為

根據(jù) Malthusian方程［18］，水環(huán)境治理PPP項目政府主體復(fù)制動態(tài)方程W（x）為

在水環(huán)境治理PPP項目中，設(shè)經(jīng)濟主體采用自主推進(jìn)策略時的期望收益為Uy2，采用被動履責(zé)策略時的期望收益為U1-y2，平均期望收益為U2，可得Uy2，U1-y2，U2分

將均衡點帶入雅可比矩陣J，可得演化博弈均衡點的特征值，如表2所示.

由于D9=（x*，y*，z*）對應(yīng)的雅可比矩陣跡為0，其特征值之和為0，不存在3個特征值都為負(fù)，D9=（x*，y*，z*）非ESS，故不在表2中討論.

模型中所有的損益變量均大于0，且E1gt;B1，E2gt;B2，由于D1=（0，0，0），D4=（0，0，1），D6=（1，0，1）存在特征值恒大于0，即E1-B1+K1+S1gt;0，E2-B2+S2gt;0，E2-B2+J1+K1+S2gt;0，不滿足李雅普諾夫穩(wěn)定性理論，故只對D2=（1，0，0），D3=（0，1，0），D5=（1，1，0），D7=（0，1，1），D8=（1，1，1）這5個均衡點的穩(wěn)定性進(jìn)行分析，存在以下5種情況.

情況1：當(dāng)E2-B2-γL1+J1+K1lt;0，E3-C3-βL2+J2lt;0時，由李雅普諾夫穩(wěn)定性理論可知，均衡點D2=（1，0，0）為ESS，此時，政府主體選擇積極規(guī)制策略，經(jīng)濟主體選擇被動履責(zé)策略，社會主體選擇不參與策略.

情況2：當(dāng)E1-B1-J1lt;0，E2-B2-γL1gt;0，E3-C3+S3lt;0時，由李雅普諾夫穩(wěn)定性理論可知，均衡點D3=（0，1，0）為ESS，此時，政府主體選擇消極規(guī)制策略，經(jīng)濟主體選擇自主推進(jìn)策略，社會主體選擇不參與策略.

情況3：當(dāng)E1-B1-J1gt;0，E2-B2-γL1+J1+K1gt;0，E3-C3+J2+S3lt;0時，由李雅普諾夫穩(wěn)定性理論可知，均衡點D5=（1，1，0）為ESS，此時，政府主體選擇積極規(guī)則策略，經(jīng)濟主體選擇自主推進(jìn)策略，社會主體選擇不參與策略.

情況4：當(dāng)E1-B1-J1-J2lt;0，E3-C3+S3gt;0時，由李雅普諾夫穩(wěn)定性理論可知，均衡點D7=（0，1，1）為ESS，此時，政府主體選擇消極規(guī)制策略，經(jīng)濟主體選擇自主推進(jìn)策略，社會主體選擇主動參與策略.

情況5：當(dāng)E1-B1-J1-J2gt;0，E3-C3+J2+S3gt;0時，由李雅普諾夫穩(wěn)定性理論可知，均衡點D8=（1，1，1）為ESS，此時，政府主體選擇積極規(guī)制策略，經(jīng)濟主體選擇自主推進(jìn)策略，社會主體選擇主動參與策略.

通過上述分析可知，演化穩(wěn)定策略（圖1）受多個參數(shù)影響，且不同約束條件下對應(yīng)的ESS不同.為了進(jìn)一步驗證博弈模型結(jié)果的合理性及有效性，通過MATLAB軟件對演化穩(wěn)定策略和參數(shù)敏感性進(jìn)行仿真模擬.

3 多主體演化博弈仿真分析

3.1 博弈模型結(jié)果驗證

1） 當(dāng)復(fù)制動態(tài)系統(tǒng)ESS為D2=（1，0，0）時，

在 E2-B2-γL1+J1+K1lt;0且 E3-C3-βL2+J2lt;0的條件下，對仿真數(shù)據(jù)進(jìn)行賦值，結(jié)果如表3所示.

將初始值帶入MATLAB軟件，得到情況1的演化軌跡，如圖2所示.由圖2可知：在 E2-B2-γL1+J1+K1lt;0且 E3-C3-βL2+J2lt;0的條件下，任意（x，y，z）的仿真結(jié)果都趨于D2=（1，0，0），這說明隨著系統(tǒng)演化，經(jīng)濟主體和社會主體逐漸趨于選擇被動履責(zé)和不參與策略，此時，政府主體為了推動水環(huán)境治理PPP項目實施，提高治理效果，最終會選擇積極規(guī)制策略.

2） 當(dāng)復(fù)制動態(tài)系統(tǒng)ESS為D3=（0，1，0）時，在E1-B1-J1lt;0，E2-B2-γL1gt;0且E3-C3+S3lt;0的條件下，對仿真數(shù)據(jù)進(jìn)行賦值，結(jié)果如表4所示.

將初始值帶入MATLAB軟件，得到情況2的演化軌跡，如圖3所示.由圖3可知：在E1-B1-J1lt;0，E2-B2-γL1gt;0且E3-C3+S3lt;0的條件下，任意（x，y，z）的仿真結(jié)果都趨于D3=（0，1，0），這是因為政府主體對自主推進(jìn)項目的經(jīng)濟主體提供稅收優(yōu)惠和經(jīng)濟補貼高于積極規(guī)制下可以獲得的額外經(jīng)濟利潤，因此，政府主體最終選擇消極規(guī)制策略.對于社會主體而言，此時主動參與項目成本較大，因此，選擇不參與策略.然而，經(jīng)濟主體此時獲得的額外經(jīng)濟利潤大于社會主體不參與項目時，經(jīng)濟主體選擇自主推進(jìn)策略給本方帶來的損失，因此，經(jīng)濟主體選擇自主推進(jìn)策略.

3） 當(dāng)復(fù)制動態(tài)系統(tǒng)的ESS為D5=（1，1，0）時，

在E1-B1-J1gt;0，E2-B2-γL1+J1+K1gt;0且E3-C3+J2+S3lt;0的條件下，對仿真數(shù)據(jù)進(jìn)行賦值，結(jié)果如表5所示.

將初始值帶入MATLAB軟件，可得情況3的演化軌跡，如圖4所示.由圖4可知：在E1-B1-J1gt;0，E2-B2-γL1+J1+K1gt;0且E3-C3+J2+S3lt;0的條件下，任意（x，y，z）的仿真結(jié)果趨于D5=（1，1，0），說明當(dāng)社會主體主動參與水環(huán)境治理PPP項目成本過大時，即使政府主體選擇積極規(guī)制策略，經(jīng)濟主體選擇自主推進(jìn)策略，社會主體選擇主動參與的概率也會逐漸降低，最終趨于0.

4） 當(dāng)復(fù)制動態(tài)系統(tǒng)ESS為D7=（0，1，1）時，

在 E1-B1-J1-J2lt;0且E3-C3+S3gt;0的條件下，對仿真數(shù)據(jù)進(jìn)行賦值，結(jié)果如表6所示.

將初始值帶入MATLAB軟件，可得情況4的演化軌跡，如圖5所示.在 E1-B1-J1-J2lt;0且E3-C3+S3gt;0的條件下，任意（x，y，z）仿真結(jié)果都趨于D7=（0，1，1），這說明當(dāng)政府主體積極規(guī)制下獲得的額外經(jīng)濟利潤低于給予經(jīng)濟主體和社會主體的補貼時，政府主體選擇積極規(guī)制的概率會逐漸趨于0；而當(dāng)采取積極行為獲得的額外收益過高時，即使政府主體選擇消極規(guī)制策略，經(jīng)濟主體和社會主體最終依舊會選擇自主推進(jìn)策略和主動參與策略.

5） 當(dāng)復(fù)制動態(tài)系統(tǒng)ESS為D8=（1，1，1）時，在E1-B1-J1-J2gt;0且E3-C3+J2+S3gt;0的條件下，對仿真數(shù)據(jù)進(jìn)行賦值，結(jié)果如表7所示.

將初始值帶入MATLAB軟件，可得情況5的演化軌跡，如圖6所示.由圖6可知：在E1-B1-J1-J2gt;0且E3-C3+J2+S3gt;0的條件下，任意（x，y，z）的仿真結(jié)果都趨于D8=（1，1，1），即在確保水環(huán)境治理PPP項目利益相關(guān)者采取積極行為可獲得足夠利益的情況下，政府主體選擇積極規(guī)則策略，經(jīng)濟主體選擇自主推進(jìn)策略，社會主體選擇主動參與策略的概率均趨近于1.

3.2 政府主體激勵強度對復(fù)制動態(tài)系統(tǒng)穩(wěn)定性的影響

政府主體為了激發(fā)經(jīng)濟主體和社會主體的內(nèi)在活力，提高水環(huán)境治理效率，通常對參與主體提供優(yōu)惠政策，如進(jìn)行經(jīng)濟補貼或稅收減免.根據(jù)上文假設(shè)，對仿真數(shù)據(jù)進(jìn)行賦值（表8），并設(shè)x=y=z=0.65.

1） 不同J1取值對策略選擇的影響.

在其他參數(shù)不變的條件下，令J1=10，15，20，35，可得J1變化時復(fù)制動態(tài)系統(tǒng)的演化軌跡圖，如圖7所示.

由圖7可知：隨著J1不斷增大，整個系統(tǒng)從D2=（1，0，0）向D5=（1，1，0）趨近，說明政府主體對經(jīng)濟主體的激勵強度提升在一定程度上可以促使其在水環(huán)境治理過程中產(chǎn)生積極行為.然而，在（20，35）之間存在一個臨界值J*1，當(dāng)J1gt;J*1時，政府主體和經(jīng)濟主體的演化軌跡呈上下波動形式，復(fù)制動態(tài)系統(tǒng)不會穩(wěn)定在中心點，政府主體和經(jīng)濟主體的策略選擇處于不斷變化中，沒有演化穩(wěn)定策略.因此，不宜盲目增大激勵強度.

2） 不同J2取值對策略選擇的影響.

在其他參數(shù)不變的條件下，令J2=8，12，16，28，可得J2變化時復(fù)制動態(tài)系統(tǒng)演化軌跡圖，如圖8所示.

由圖8可知：（8，12）之間存在一個臨界值J*2，（16，28）之間存在一個臨界值J**2，J*2lt;J**2；當(dāng)J2lt;J*2時，系統(tǒng)向D2=（1，0，0）趨近；當(dāng)J*2lt;J2lt;J**2時，系統(tǒng)向D8=（1，1，1）趨近，隨著J2的增大，系統(tǒng)收斂速度逐漸加快；當(dāng)J2gt;J**2時，政府主體和社會主體的演化軌跡上下震動，復(fù)制動態(tài)系統(tǒng)不會穩(wěn)定在中心點，政府主體和社會主體的策略選擇處于不斷變化中，沒有演化穩(wěn)定策略.由此可知，增大J2不僅可以激發(fā)社會主體參與水環(huán)境治理項目的熱情，還可以促進(jìn)經(jīng)濟主體積極主動投身于水環(huán)境治理PPP項目，實現(xiàn)項目合作共贏.然而，與J1取值類似，應(yīng)當(dāng)把J2控制在合理取值范圍之內(nèi)，以防因J2過高導(dǎo)致復(fù)制動態(tài)系統(tǒng)不穩(wěn)定.

3.3 政府主體懲罰力度對復(fù)制動態(tài)系統(tǒng)穩(wěn)定性的影響

在其他參數(shù)不變的條件下，令K1=20，25，30，35，可得K1變化時復(fù)制動態(tài)系統(tǒng)的演化軌跡圖，如圖9所示.由圖9可知：（25，30）之間存在一個臨界值K*1，當(dāng)K1lt;K*1時，復(fù)制動態(tài)系統(tǒng)向D2=（1，0，0）趨近；當(dāng)K*1lt;K1時，系統(tǒng)向D5=（1，1，0）趨近，且隨著K1的增大，經(jīng)濟主體策略概率趨近于1的收斂速度逐漸加快.因此，增大K1可以避免經(jīng)濟主體選擇消極行為，但此時經(jīng)濟主體的行為轉(zhuǎn)變難以帶動社會主體參與水環(huán)境治理PPP項目的積極性.

3.4 損失系數(shù)對復(fù)制動態(tài)系統(tǒng)穩(wěn)定性的影響

1） 損失系數(shù)γ對策略選擇的影響.

在其他參數(shù)不變的條件下，令γ=0.14，0.16，0.18，0.20，可得不同損失系數(shù)γ下復(fù)制動態(tài)系統(tǒng)的演化軌跡圖，如圖10所示.由圖10可知：（0.16，0.18）之間存在一個臨界值γ*，當(dāng)γ*lt;γ時，系統(tǒng)向D2=（1，0，0）趨近；當(dāng)γlt;γ*時，系統(tǒng)向D5=（1，1，0）趨近，且隨著γ的減小，經(jīng)濟主體策略概率趨近于1的收斂速度逐漸加快，說明γ減小可推動經(jīng)濟主體選擇積極策略.

2） 損失系數(shù)β對策略選擇的影響.

在其他參數(shù)不變的條件下，令β=0.24，0.26，0.28，0.30，可得不同損失系數(shù)β下復(fù)制動態(tài)系統(tǒng)的演化軌跡圖，如圖11所示.由圖11可知：系統(tǒng)對損失系數(shù)β不敏感，此時，系統(tǒng)演化穩(wěn)定策略仍為D2=（1，0，0），這說明當(dāng)經(jīng)濟主體選擇被動履責(zé)策略時，即使減小社會主體主動參與項目給本方帶來的風(fēng)險，也不會改變他們的意愿.

4 結(jié)論與建議

采用演化博弈分析法構(gòu)建多主體行為演化博弈模型，分析政府主體、經(jīng)濟主體和社會主體在博弈過程中的演化均衡策略，并通過MATLAB軟件進(jìn)行數(shù)值仿真模擬，可得以下3個結(jié)論.

1） 從演化博弈角度看，水環(huán)境治理PPP項目政府主體、經(jīng)濟主體和社會主體構(gòu)成的復(fù)制動態(tài)系統(tǒng)穩(wěn)定性受多個參數(shù)影響，但在基本條件E1gt;B1，E2gt;B2的情況下，可通過調(diào)節(jié)系統(tǒng)參數(shù)達(dá)到{積極規(guī)制，自主推進(jìn)，主動參與}的理想穩(wěn)定狀態(tài).

2） 從激勵角度看，政府主體提高獎勵力度可促使經(jīng)濟主體和社會主體采取正向積極行為.其中，增加對社會主體的激勵強度可以同時加快自主推進(jìn)、主動參與策略的收斂速度.然而，系統(tǒng)穩(wěn)定性對激勵強度較為敏感，政府主體需要控制獎勵力度，當(dāng)超過臨界值時，激勵效果反而會降低.

3） 從懲罰與損失角度看，政府主體實施懲罰規(guī)制在一定程度上可以約束經(jīng)濟主體策略選擇，干預(yù)消極行為；損失系數(shù)γ的減小對經(jīng)濟主體策略選擇具有正向推動作用，而損失系數(shù)β對演化博弈影響較弱.因此，弱化社會主體消極行為對經(jīng)濟主體的負(fù)面影響，促使政府主體選擇積極規(guī)制策略是提高項目治理績效的關(guān)鍵點.

為了優(yōu)化政府規(guī)制，提高水環(huán)境治理PPP項目的治理成效，提出以下3個建議.

1） 激發(fā)企業(yè)內(nèi)在創(chuàng)新力.融合生態(tài)環(huán)境導(dǎo)向的開發(fā)（EOD）模式理念，政府可以鼓勵相關(guān)治污主體在完成水環(huán)境治理的同時，結(jié)合當(dāng)?shù)靥厣_發(fā)特色產(chǎn)業(yè)鏈，如生態(tài)農(nóng)業(yè)、生態(tài)旅游、康養(yǎng)院等，從而提高項目經(jīng)濟效益，確保基本條件E1gt;B1，E2gt;B2在治理項目實施過程中處于恒成立狀態(tài).

2） 完善政府獎懲規(guī)制.政府規(guī)制部門應(yīng)制定詳細(xì)的水環(huán)境治理PPP項目特征指標(biāo)、運營維護(hù)和績效考核指標(biāo)等，以此作為獎懲標(biāo)準(zhǔn)，對于履責(zé)或主動參與的，可給予經(jīng)濟獎勵、經(jīng)營性項目使用優(yōu)惠、其他項目政策補償?shù)燃畋Ｕ?；對于違背合同協(xié)議或采取投機行為的，可進(jìn)行高力度經(jīng)濟懲罰，嚴(yán)重者在一定時間范圍內(nèi)禁止其參與此類項目，從而約束消極行為.同時，要在項目建設(shè)前期測算J1，J2的臨界值，采取動態(tài)監(jiān)管手段，將激勵強度控制在一定閾值內(nèi)，防止出現(xiàn)復(fù)制動態(tài)系統(tǒng)失衡現(xiàn)象.

3） 優(yōu)化多主體協(xié)同治理機制.為了促進(jìn)經(jīng)濟主體策略概率趨近于1，在制定規(guī)制策略時可以把弱化社會主體消極行為對經(jīng)濟主體的負(fù)面影響作為提高項目績效的切入點，通過環(huán)保宣傳提高社會主體的環(huán)保意識；建立多渠道參與機制，通過交互性新媒體網(wǎng)絡(luò)傳播形成網(wǎng)絡(luò)傳遞效應(yīng)，對經(jīng)濟主體消極行為形成嚇阻.

參考文獻(xiàn)：

［1］宋民雪，劉德海，尹偉巍.經(jīng)濟新常態(tài)、污染防治與政府規(guī)制： 環(huán)境突發(fā)事件演化博弈模型［J］.系統(tǒng)工程理論與實踐，2021，41（6）：1454－1464.DOI：10.12011/SETP2019－2338.

［2］AN Xiaowei，LI Huimin，WANG Lunyan，et al.Compensation mechanism for urban water environment treatment PPP project in China［J］.Journal of Cleaner Production，2018，201：246－253.DOI：10.1016/j.jclepro.2018.08.003.

［3］張亞瓊，何楠，楊絲雯，等.基于演化博弈的水生態(tài)PPP項目第三方監(jiān)管策略［J］.水資源保護(hù)，2021，37（4）：166－172.DOI：10.3880/j.issn.1004－6933.2021.04.024.

［4］LI Huimin，XIA Qing，WEN Shiping，et al.Identifying factors affecting the sustainability of water environment treatment public－private partnership projects［J］.Advances in Civil Engineering，2019，2019（11）：1－15.DOI：10.1155/2019/7907234.

［5］呂天宇，趙微，李娜，等.基于農(nóng)民滿意度的農(nóng)村水環(huán)境治理績效研究： 以湖北省部分縣市的調(diào)查為例［J］.水土保持通報，2020，40（6）：137－145，152.DOI：10.13961/j.cnki.stbctb.2020.06.020.

［6］王麗瓊，張云峰.鄉(xiāng)村振興視閾下泉州市農(nóng)村環(huán)境多元共治有效路徑研究［J］.中國農(nóng)業(yè)資源與區(qū)劃，2019，40（8）：219－225.DOI：10.7621/cjarrp.1005－9121.20190827.

［7］

王雪瑩，張新月，雷陽，等.于橋水庫浮游動物群落結(jié)構(gòu)與水環(huán)境因子關(guān)系多元分析［J］.福建師范大學(xué)學(xué)報（自然科學(xué)版），2021，37（4）：32－40.DOI：10.12046/j.issn.1000－5277.2021.04.005.

［8］OKUMAH M，YEBOAH A S.Exploring stakeholders′ perceptions of the quality and governance of water resources in the Wenchi municipality［J］.Journal of Environmental Planning and Management，2020，63（8）：1375－1403.DOI：10.1080/09640568.2019.1663724.

［9］黃德春，宋佳，賀正齊，等.瀾滄江－湄公河環(huán)境利益合作網(wǎng)絡(luò)主體治理效益評價［J］.亞太經(jīng)濟，2019（4）：21－29，149－150.DOI：10.16407/j.cnki.1000－6052.20190819.010.

［10］COTLER H，CUEVAS M L，LANDA R，et al.Environmental governance in urban watersheds： The role of civil society organizations in Mexico［J］.Sustainability，2022，14（2）：1－26.DOI：10.3390/SU14020988.

［11］ZHAO Xin，BAI Yu，DING Lili.Incentives for personal carbon account： An evolutionary game analysis on public－private－partnership reconstruction［J］.Journal of Cleaner Production，2021，282（9）：125358.DOI：10.1016/j.jclepro.2020.125358.

［12］潘鶴思，李英，柳洪志.央地兩級政府生態(tài)治理行動的演化博弈分析： 基于財政分權(quán)視角［J］.生態(tài)學(xué)報，2019，39（5）：1772－1783.DOI：10.5846/stxb201803230578.

［13］LIU Jicai，GAO Ruolan，CHEAH C Y J，et al.Evolutionary game of investors′ opportunistic behaviour during the operational period in PPP projects［J］.Construction Management and Economics，2017，35（1/2/3）：137－153.DOI：10.1080/01446193.2016.1237033.

［14］申津羽，王煜琪，趙正.公眾參與城市水環(huán)境治理行為的影響因素： 基于北上廣三地的調(diào)查數(shù)據(jù)［J］.資源科學(xué)，2021，43（11）：2289－2302.DOI：10.18402/resci.2021.11.12.

［15］任志濤，郭亞，孫彥武.環(huán)境治理公私合作項目多方主體價值共創(chuàng)演化博弈研究［J］.軟科學(xué)，2022，36（4）：132－137.DOI：10.13956/j.ss.1001－8409.2022.04.18.

［16］高旭闊，席子云.組合措施下政府與企業(yè)排污行為演化博弈［J］.中國環(huán)境科學(xué)，2020，40（12）：5484－5492.DOI：10.19674/j.cnki.issn1000－6923.2020.0607.

［17］張飛漣，梁秀峰，張志良，等.PPP績效付費項目參與主體機會主義行為演化研究［J］.鐵道科學(xué)與工程學(xué)報，2020，17（4）：1053－1060.DOI：10.19713/j.cnki.43－1423/u.T20191012.

［18］CHEN Hongmei，WANG Jianxue，MIAO Yujun.Evolutionary game analysis on the selection of green and low carbon innovation between manufacturing enterprises［J］.Alexandria Engineering Journal，2021，60（2）：2139－2147.DOI：10.1016/J.AEJ.2020.12.015.

［19］崔萌.協(xié)同治理背景下環(huán)保信用監(jiān)管的三方演化博弈分析［J］.系統(tǒng)工程理論與實踐，2021，41（3）：713－726.DOI：10.12011/SETP2020－0480.

［20］XU Lang，DI Zhongjie，CHEN Jihong，et al.Evolutionary game analysis on behavior strategies of multiple stakeholders in maritime shore power system［J］.Ocean and Coastal Management，2021，202：105508.DOI：10.1016/J.OCECOAMAN.2020.105508.

（責(zé)任編輯：" 錢筠 英文審校： 劉源崗）

收稿日期： 2022－12－29

通信作者： 詹朝曦（1970－），男，副教授，博士，主要從事建筑經(jīng)濟、工程管理的研究.E－mail：75963701@qq.com.

基金項目： 國家社科基金一般資助項目（21BZZ069）； 福建省軟科學(xué)創(chuàng)新戰(zhàn)略研究項目（2021R0056）http：∥www.hdxb.hqu.edu.cn