海洋旅游资源非使用价值支付意愿与家庭收入水平的关系研究——以舟山、三亚为例

肖建红 李小嵩

1.青岛大学旅游与地理科学学院 青岛266071

2.青岛大学经略海洋研究中心 青岛266071

3.青岛大学商学院 青岛266071

评估海洋旅游资源的非使用价值（non-use value）可以为海洋旅游资源的保护和管理提供参考。条件价值评估法（contingent valuation method，CVM）是最常被采用的非使用价值评估方法，其通过假想情景引出受访者对非市场物品的支付意愿（willingness to pay，WTP）[1]。为了提出更合理、精准的海洋旅游资源保护政策，通常需要通过各种回归模型对受访者WTP 的影响因素进行研究。在对受访者WTP 影响因素的相关研究中，学者们常采用最小二乘回归、Tobit、Logit 等模型研究WTP 均值的影响因素，但是这些模型都假定对于不同WTP 群体，解释变量对被解释变量的影响是相同的，不能清晰地反映WTP 的条件分布全貌，而分位数回归模型（quantile regression，QR）能更加全面地了解解释变量对WTP 条件分布的影响[2]。在研究收入对WTP 的影响时，关注WTP 整体分布的研究较少，而通过分位数回归模型可以了解受访者家庭收入水平对不同群体WTP 的影响程度的差异，从而为不同案例地不同海洋旅游资源的保护政策提供参考[3]。

分位数回归模型最早由Koenker 和Bassett 提出[2]，在绝大多数回归模型都关注被解释变量的条件均值的情况下，分位数回归提供了对解释变量X 和被解释变量Y 的分位数之间线性关系的估计方法，与普通均值回归相比，分位数回归使用残差绝对值的加权平均作为最小化的目标函数，能避免受极端值的影响，适合具有异方差的模型，使回归结果更加稳健。分位数回归对被解释变量的条件分布刻画得更加细致，能给出被解释变量条件分布的大体特征。每个分位点上的回归都被赋予条件分布上某个特殊点（中央或尾部）的一些特征；把不同的分位点回归集中起来就能提供一个关于被解释变量条件分布得更完整的统计特征描述。分位数回归并不要求很强的分布假设，在扰动项非正态分布的情形下，分位数估计量可能比最小二乘估计量更为有效。目前分位数回归模型已被学者用于环境物品影响因素的相关研究中，代表性研究成果主要有：Krishnamurthy 等运用分位数回归模型研究了11个国家家庭绿色能源WTP 的影响因素，发现收入对WTP 分布的中心处有显著的影响[4]；O'Garra 等则运用分位数回归模型研究了受访者对于在伦敦引入氢公交车以减少空气和噪声污染WTP 的影响因素，结果表明分位数回归对分析CVM 数据有效，增强了对支付意愿决定因素的理解[5]；Notaro等运用分位数回归模型研究了游客对保护加尔达湖特伦蒂诺湖岸柏树WTP 的影响因素，表明分位数回归可以分析条件分布下WTP 的影响因素[6]。国内也有学者运用分位数回归模型研究WTP 的影响因素，代表性成果主要有：郑海霞等运用分位数回归模型，研究了居民对改善密云水库水质WTP 的影响因素，揭示了不同支付层次、支付能力和收入水平的群体对环境政策的需求和偏好[3]；李佳佳等采用CVM 评估了居民对两种假想医疗保险的WTP，并在此基础上运用分位数回归模型研究了WTP 的影响因素，进一步证实了低值WTP 和高值WTP 的影响因素不同[7]；黄颖等运用分位数回归模型研究了在各个分位点上农户对政策性森林保险WTP 的影响因素，发现农户的WTP 在受教育程度上具有明显分化[8]；孙聪等运用分位数回归模型研究了住房支出的影响因素，识别了群体差异[9]；朱宁等运用分位数回归模型，实证探讨城镇居民家庭鸡蛋价格承受能力以及在蛋价高位、品牌与普通鸡蛋价差小时的消费倾向[10]。

目前对海洋旅游资源的研究多以单一资源为评估对象[11-14]，而鲜有评估多个案例地、多种海洋旅游资源非使用价值的研究成果。本文以舟山和三亚两个案例地的海滩、海洋文化和珊瑚礁3 种海洋旅游资源作为研究对象，采用最小二乘回归模型和分位数回归模型，研究居民和游客对保护海洋旅游资源WTP 的影响因素，尤其是公众家庭收入对受访者WTP 的影响，以期为分位数回归模型在海洋旅游资源领域的应用以及对海洋旅游资源政策的制定提供参考。

1 数据来源与计量经济模型

1.1 数据来源

通过CVM 调查了居民和游客对保护舟山和三亚两个案例地的海滩、海洋文化和珊瑚礁3 种海洋旅游资源的WTP，两个案例地3 种海洋旅游资源采用的问卷结构、调查设计和调查方法一致。问卷包括对待海洋旅游资源的详细描述、核心估值问题和受访者人口社会特征及不愿意支付的原因3 部分。在调查过程中，为了提高调查结果的有效性和可靠性，对假想偏差、信息偏差、抗议性反应偏差和引导方式偏差进行了控制。本文采用明确提醒受访者的方法控制假想偏差，提醒内容包括想象假想表决为真实表决，家庭预算约束，因这项计划而额外增加的成本，以及海滩、海洋文化和珊瑚礁其他替代品的存在[15-17]；采用对调查人员进行严格培训，使调查人员掌握详细调研信息，再通过调查人员面对面向受访者讲解相关信息来控制信息偏差[15]；采用设置一系列不愿意支付的原因，用于识别真零支付和抗议性零支付，被识别的抗议性零支付在进一步分析中被剔除[18-20]；目前，海洋旅游资源非市场经济价值评估中，有关支付卡式、二分式和开放式3 种投标值引导方式的研究成果均较多，研究采用了支付卡式引导方式，一方面可以有效规避二分式引导方式存在的肯定性回答偏差和起点偏差，另一方面可以有效减轻受访者的认知负担，降低不反应率和提高研究的高效性和精准性[21-22]。

调查对象为舟山和三亚的年龄为18～65 周岁的游客和当地居民。正式调查采用面对面调查方式，于2017 年7—8 月进行，调查共获得有效问卷4 763 份。舟山共获得有效问卷2 059 份，其中海滩1 195 份，海洋文化864份；三亚共获得有效问卷2 704份，其中海滩1 364份，珊瑚礁1 340份。

1.2 计量经济学模型

分位数回归模型与绝大多数关注被解释变量条件均值的回归模型相比，能避免极端值的影响，对条件分布刻画得更加细致，给出被解释变量条件分布更完整的统计特征描述，使回归结果更加稳健。分位数回归模型不要求很强的分布假设，在扰动项非正态分布的情形下，分位数估计量可能比最小二乘估计量更为有效。本文运用最小二乘回归模型和分位数回归模型研究WTP 与样本人口社会特征之间的关系，各变量的详细定义参见表1。

表1 回归模型各变量定义

建立OLS回归模型：

在式（1）中，α为常数项，βi（i=1，2，3，4，5，6）为各个解释变量的回归系数，μ为误差项，其他变量解释见表1。

建立分位数回归模型：

在式（2）和式（3）中，Xi表示解释变量的向量，βθ表示各个解释变量的回归系数，μi,θ为误差项，quantθ(WTPi|Xi)为给定Xi条件下的第θ分位数回归的WTPi，0 ＜θ＜1。

为了求解样本中第θ分位数回归，使用线性规划方法求解最小化残差绝对值的加权平均：

在式（4）中，bθ为βθ估计量。

2 结果分析

2.1 受访者人口社会特征

受访者人口社会特征如表2 所示，舟山海滩和海洋文化受访者的人口社会特征类似，在两个样本中，男女比例相当，平均年龄约为37岁，约47%的受访者接受过高等教育，家庭年收入约为14万元，超过54%的受访者是游客。三亚海滩受访者的男女比例相当，平均年龄约为34 岁，约60%的受访者接受过高等教育，家庭年收入约为13万元，约60%的受访者为游客。三亚珊瑚礁的受访者中男性居多，平均年龄约为33 岁，约66%的受访者接受过高等教育，家庭年收入约为14万元，居民和短途、中途以及长途游客的比例相当。

表2 受访者人口社会特征均值

2.2 受访者WTP分布

受访者WTP 的分布比例参见图1。舟山的样本中分别有96.7%和97.1%的受访者对保护海滩和海洋文化具有正WTP，分别有84.4%和79.1%的受访者选择200 元及以下的投标值；三亚的样本中分别有98.4%和98.7%的受访者对保护海滩和珊瑚礁具有正WTP，分别有85.2%和84.9%的受访者选择200元及以下的投标值。

图1 受访者WTP分布比例

2.3 回归结果分析

选取0.1—0.9 共9 个分位点，分位数回归结果参见表3—表6。从总体来看，在舟山和三亚的3 种海洋旅游资源的分位数回归模型中，家庭收入变量在9 个分位点上都呈现出统计显著性，且回归系数都为正值，说明受访者保护两个案例地的3 种海洋旅游资源的WTP 都与收入成正比，受访者的收入越高其WTP越大。根据收入回归系数图（图2—图5）可以得出，收入变量的回归系数随着分位点的升高而出现递增趋势，说明随着WTP 的提高，收入变化对受访者的影响变大，即收入对高WTP 群体的影响更大。通过收入变量回归系数的递增趋势，可以预测随着公众收入水平的上升，高WTP 群体WTP 的增加量更大，高WTP 群体与低WTP群体之间WTP的差距也会随之变大。

图2 舟山海滩旅游资源回归系数图

图3 舟山海洋文化旅游资源回归系数图

图4 三亚海滩旅游资源回归系数图

图5 三亚珊瑚礁旅游资源回归系数图

表3 舟山海滩旅游资源回归结果

在对同一案例地不同海洋旅游资源（海滩与海洋文化和海滩与珊瑚礁）以及不同案例地同一资源（海滩）的分位数回归结果的比较中（图6），可以观察到受访者家庭收入对WTP 的影响存在一定差距。在舟山海滩与海洋文化的比较中发现，在70%及以下分位点上受访者家庭收入的回归系数不存在差异，均较为接近，在70%以上分位点上舟山海滩样本收入变量的回归系数大于舟山海洋文化样本中的系数，这说明对于舟山海滩和海洋文化而言，收入对中低WTP 群体的影响力没有显著差异，而对于高WTP 群体，收入对舟山海滩的影响力大于对海洋文化的影响力。在三亚海滩与珊瑚礁的比较中发现，收入变量的回归系数在所有9 个分位点上都较为接近，其中在80%和90%分位点上，三亚珊瑚礁中收入变量的回归系数略大于三亚海滩中的系数，说明对于中低WTP 群体而言，收入对三亚海滩和珊瑚礁WTP 的影响力没有显著差异，而对于高WTP 群体，收入对珊瑚礁WTP 的影响力略大于其对海滩WTP 的影响力。在舟山海滩与三亚海滩的比较中，收入变量的回归系数在50%及以下分位点不存在显著差异，而随着分位点的提高，收入系数的差距呈现递增的趋势，同时，舟山海滩中的收入系数始终大于三亚海滩的系数，说明对于低WTP 群体，收入对舟山海滩和三亚海滩WTP 的影响力没有显著差异，而对于高WTP 群体，收入对舟山海滩WTP 的影响力大于对三亚海滩WTP的影响力。

图6 收入变量分位数回归系数图

最小二乘回归模型的回归结果（表3—表6）与分位数回归模型的结果一致，两个模型都表明，在两个案例地、3种海洋旅游资源的4组样本中，家庭收入均显著正向影响受访者WTP，即家庭收入水平越高，受访者的WTP 越高。但是在最小二乘回归模型和分位数回归模型中家庭收入变量的系数大小存在显著差异，分位数回归模型中家庭收入系数呈递增状态，最小二乘回归模型中的系数仅与某些分位点上的系数较为接近。如在舟山海滩的两种模型中，最小二乘回归模型中家庭收入的系数与分位数回归模型中60%分位点上的系数最为接近，而与其他分位点上的系数差距很大。这说明最小二乘回归模型仅能反映家庭收入对WTP 均值的影响，不能充分地反映收入对低分位点群体和高分位点群体WTP 的影响，因此将假设家庭收入变量对WTP 影响具有同质性的最小二乘回归模型应用到CVM 估计结果的分析中得到的信息量较少，不能反映WTP 的条件分布全貌，需要采用分位数回归模型来更精确地刻画家庭收入对不同群体WTP 的不同影响程度。

舟山和三亚两个案例地3 种海洋旅游资源的最小二乘回归结果和分位数回归结果表明（表3—表6），性别、年龄、学历和受访者类型的结果也存在差异，一些变量在最小二乘回归模型中不显著，在分位数回归模型中却在部分分位点上显著，一些变量在最小二乘回归模型中显著，在分位数回归模型中却仅在部分分位点上显著。舟山海滩的回归结果（表3）表明：在最小二乘回归模型中，年龄和受访者类型对受访者WTP 没有显著影响，但是在分位数回归模型中，年龄在70%分位点上对受访者WTP有显著负向影响，受访者类型在60%分位点上对WTP有显著负向影响。舟山海洋文化的回归结果（表4）表明：在最小二乘回归模型中，性别和受访者类型对受访者WTP 没有显著影响，学历对受访者WTP 有显著正向影响，而在分位数回归模型中，性别在80%分位点上显著正向影响受访者WTP，受访者类型在90%分位点上显著负向影响受访者WTP，学历仅在40%分位点上显著正向影响受访者WTP。三亚海滩的回归结果（表5）表明：在最小二乘回归模型中，性别对受访者WTP 有显著正向影响，年龄、学历和受访者类型对受访者WTP 没有显著影响，而在分位数回归模型中，性别在9 个分位点上都对WTP没有显著影响，年龄在60%和70%分位点上显著负向影响WTP，学历在40%分位点上显著正向影响WTP。三亚珊瑚礁的回归结果（表6）表明：在最小二乘回归模型中，年龄与受访者WTP 有显著负相关关系，学历和受访者类型对受访者WTP 没有显著影响，而在分位数回归模型中，年龄在60%、70%和90%分位点上显著负向影响受访者WTP，学历在30%、40%和50%分位点上显著正向影响受访者WTP，受访者类型在40%、50%、60%和70%分位点上都显著影响受访者WTP。由此可见，采用最小二乘回归模型得到的结果可能只能反映部分受访者的偏好，而不能反映不同群体WTP 影响因素的全貌，得到的政策建议的针对性不强，可能达不到很好的实施效果。而分位数回归模型可以得到更加细化的不同群体WTP 的影响因素，从而得到更加具体、实施性更强的政策建议。

表4 舟山海洋文化旅游资源回归结果

表5 三亚海滩旅游资源回归结果

表6 三亚珊瑚礁旅游资源回归结果

通过最小二乘回归结果和分位数回归结果可以观察到，在两个案例地、3 种海洋旅游资源的4 组样本中，仅有家庭收入在两种模型中（分位数回归的所有分位点上）始终对受访者WTP 存在显著正向影响，而性别、学历、年龄和受访者类型则没有呈现普遍的显著影响，这说明受访者的家庭收入是受访者保护海洋旅游资源WTP的主要驱动因素。

3 结论与建议

本文以舟山海滩和海洋文化、三亚海滩和珊瑚礁为研究对象，采用最小二乘回归模型和分位数回归模型研究了居民和游客保护海洋旅游资源WTP 的影响因素，研究结果如下文。

（1）对于舟山和三亚两个案例地的海滩、海洋文化和珊瑚礁3 种海洋旅游资源，仅有家庭收入在最小二乘回归模型和分位数回归模型（所有分位点上）中都显著正向影响受访者WTP，而性别、年龄、学历和受访者类型则没有这种普遍的显著影响，这表明家庭收入是受访者保护海洋旅游资源WTP 的主要影响因素，随着家庭收入水平的升高，受访者对保护海洋旅游资源的WTP也会随之增大。

（2）家庭收入的回归系数在分位数回归模型中呈现递增趋势，说明家庭收入对高WTP 群体保护海洋旅游资源WTP 的影响更大，随着未来家庭收入水平的提高，高WTP 群体与低WTP 群体保护海洋旅游资源WTP 的差距将会增加。在未来公众收入水平迅速提高的背景下，在海洋旅游资源保护政策制定时应充分考虑家庭收入对WTP 的影响，将高收入群体作为政策制定的主要目标群体，以便能更好地对海洋旅游资源进行保护同时满足公众对其日益增长的需求。

（3）在最小二乘回归模型和分位数回归模型中，性别、年龄、学历和受访者类型的显著性存在差异；尽管家庭收入的显著性一致，但是最小二乘回归模型和分位数回归模型不同分位点上的系数大小有很大差距。

本文通过对最小二乘回归模型和分位数回归模型的对比可以得出，将假设解释变量对WTP 的影响具有同质性的最小二乘回归模型应用到海洋旅游资源非使用价值估计结果的分析中得到的信息量较少，不能反映WTP 的条件分布全貌，而分位数回归模型得到了家庭收入与不同群体WTP 之间更精确的相关关系，因此建议未来对WTP 影响因素的相关研究应优先选用分位数回归模型。