我国是渔业生产大国，水产品总产量超过世界的三分之一。水产品市场的充足供给，为“粮食安全”保驾护航。以《谁来养活中国》一书震动了中国和世界。美国生态经济学家莱斯特·布朗在接受《环球时报》采访时提出：中国对世界的贡献是计划生育和淡水渔业。

在渔业生产取得重大成就的同时，渔获物的保鲜、贮藏和运输显得尤为重要。渔业用冰的需求逐年扩大。渔业用冰按用途分为商业用冰和工业用冰。商业主要用于连锁超市、酒店、餐饮、休闲渔业山庄等；工业主要用于水产品加工、海洋捕捞、远洋渔业等相关行业。

随着中西方文化的融合，西方生活方式逐渐被中国消费者接受，对越来越多的深海捕捞渔获物，人们习惯直接食用。用冰块作为直接食用渔获物保鲜的手段，已经普遍应用在酒店、餐饮和娱乐场所，甚至进入家庭。每逢夏日，为了防止水产品腐烂变质，一家中型超市每天至少需要2 500 kg冰。为保证渔获物的质量，海洋捕捞业必须要用冰块降温。船上保鲜方式一般可分为冰藏、冷海水保鲜、盐藏，以及微冻、冻结等。沿岸、近海拖网渔船主要用0℃左右的冰藏保鲜；中小型围网渔船主要采用-1 ℃左右的冷海水保鲜；离渔港较远的作业地拖网渔船或渔业基地船采用平板冻结；部分渔船采用轻盐低温保鲜或-3 ℃左右的微冻保鲜。我国海域面积广阔，总海域面积为472.7 km²，沿海渔船约30万艘，海洋捕捞的用冰市场前景十分可观。

随着我国渔业的发展，面临着庞大的用冰市场需求。迄今为止，我国渔业用冰需求量尚未得到充分统计、论证，指导制冰行业生产的市场需求数据欠缺。制冰企业无法根据市场需求调节生产规模；涉渔第一、二、三产业的用冰量结构和需求发展趋势不明确，制冰企业对用冰消费市场结构和发展趋势不掌握，直接造成制冰企业生产存在盲目性。本文基于数理统计推理模型，对我国渔业现在和未来用冰需求量进行估算分析，旨在促进渔业用冰市场的规范化、合理化发展。

1 变量与检验 1.1 变量选择和数据处理

本文选取了1995—2016年这20年间的样本数据。我们将制冰总量(Y)作为因变量，将海水产品产量(FI)、水产品加工量(SE)、水产品消费量(TH)作为自变量。为了消除数据存在异方差的问题，我们对所有的变量取对数。下面简单介绍数据来源以及定义。

制冰总量即为渔业用冰总量，分为涉渔第一产业用冰海水产品用冰量、涉渔第二产业用冰水产品加工用冰量、涉渔第三产业用冰水产品消费用冰量。

海水产品用冰量属于涉渔第一产业用冰量的范围。海水产品产量包括国内海洋捕捞产量、远洋渔业产量和海水养殖产量，数据来源于农业部《中国渔业统计年鉴》(历年)。其中国内海洋捕捞产量指利用捕捞渔具在国内海域直接从事捕捞生产所获得的渔获物重量。按捕捞海域分，可分为渤海、黄海、东海、南海捕捞产量。按捕捞渔具分，可分为拖网、围网、刺网、张网、钓具和其他渔具捕捞产量。远洋渔业产量是指由各远洋渔业企业和各生产单位按我国远洋渔业项目管理办法组织的远洋渔船(队)在非我国管辖水域(外国专属经济区水域或公海)捕捞的水产品产量。海水养殖产量是指利用滩涂、浅海、港湾及陆上海水水体，通过人工投放苗种或天然纳苗，并经人工饲养管理所获得的水产品总量。

水产品加工用冰量属于涉渔第二产业用冰量的范围。水产品加工量是指以水产品为原料，采用各种食品贮藏加工、水产综合利用技术和工艺所生产的产品重量。如冷冻冷藏品、腌制品、干制品、熏制品、罐头食品、各种生熟小包装食品以及鱼油、鱼肝油、多烯脂肪酸制剂、饲料鱼粉、藻胶、碘和贝壳工艺品等。其中：用冰保鲜的加工品主要是水产冷冻品和冷冻加工品，数据来源于农业部《中国渔业统计年鉴》(历年)。水产冷冻品是指为了保鲜，将水产品进行冷冻加工处理后得到的产品，包括冷冻品和冷冻加工品，但不包括商业冷藏品。冷冻品泛指未改变其原始性状的粗加工产品，如冷冻全鱼、全虾等。冷冻加工品是指采用各种生产技术和工艺，改变其原始性状、改善其风味后制成的产品，如冻鱼片、冻虾仁、冷冻烤鳗和冻鱼籽等。

水产品消费用冰量属于涉渔第三产业用冰量的范围。水产品消费量是指城乡居民食用水产品的消费量，数据来源于国家统计局《中国统计年鉴》(历年)。具体包括：城乡居民的家庭消费和社会消费。其中：城乡居民的家庭消费为其次要构成，约占水产品消费总量的44%；社会消费指餐馆、饭店消费和请客送礼等形式其他消费，为其主要构成，占水产品消费总量的56%，社会消费部分所占比例有上升之势。城乡居民食用消费的水产品主要是冷冻品、鲜活水产品和半成品、熟制品、干制品等加工产品。其中，鲜活水产品和加工水产品愈来愈受消费者欢迎，所占比例逐年提高。为方便水产消费品运输，用冰量逐年提高。

1.2 计量检验

在进行回归结果分析之前，先检验各变量间是否存在异方差，再检验各变量的平稳性，最后进行协整检验。通过一系列的检验，最终能够有效保证回归参数估计量具有良好的统计性质，借助Eviews 6.0软件来进行计量检验，结果如下。

1.2.1 异方差检验

假设总体回归函数中的随机误差项满足同方差性，即它们都有相同的方差。如果这一假定不满足，即随机误差项具有不同的方差，则称线性回归模型存在异方差性。

根据表 1的检验结果，我们可以得出，P=0.12>0.05，接受原假设，即线性回归模型不存在异方差。

1.2.2 平稳性检验

为了寻找出海水产品产量、水产品加工量、水产品消费量与制冰总量之间是否存在稳定关系，首先对其进行平稳性检验。借助Eviews 6.0分析软件，采用经典的ADF单位根检验的方法对变量的对数及其差分数列做平稳性检验，结果如表 2。

检验类型(c, t, n)中的c代表检验平稳性时估计方程中的截距项，为0代表剔除截距项；第二项t代表时间趋势项，为0代表剔除时间趋势项；括号中第三项代表自回归滞后期。原假设是被检验的变量有一个单位根，*代表在5%的显著水平上拒绝原假设，**代表在1%的显著水平上拒绝原假设。

从表 2的检验结果中可以得出，变量LNY、LNFI、LNSE、LNTH的t统计量的伴随概率均小于0.05，检验结果是平稳的。

1.2.3 协整检验

由上可知，序列是平稳的，因此对各变量进行协整检验，验证变量是否存在长期均衡关系，使各个变量间的关系趋于可靠，检验结果如表 3。

表 3检验结果表明，模型中的单整变量之间具有协整关系，各变量可以通过线性关系实现平稳。

2 建模与实证分析 2.1 建立线性回归模型

经过一系列的检验，本文将实证模型进行最小二乘法(OLS)回归，判断海水产品产量、水产品加工量和水产品消费量是否对渔业用冰量产生影响。

假设模型为LNY=A×LNFI+B×LNSE+C×LNTH+D, 其中Y代表渔业用冰量，FI代表海水产品产量、SE代表水产品加工量、TH代表水产品消费量，A、B、C分别是自变量FI、SE、TH对因变量Y的影响系数大小，D为截距项，方程的线性回归结果如表 4。

上述回归方程的表格形式可以写成一般形式：

式中：R²=0.994, =0.993，F=956.362。

回归结果显示，调整后的R²为0.993，表明该模型拟合效果很好，各变量的t统计量检验结果较好。LNFI在1%的显著水平下系数为0.155，表示海洋水产品产量每增加10%，就会推动用冰量增加0.155%；LNSE在1%的显著水平下系数为0.115，表示冷冻品产量每增加10%，用冰量会增加0.115%；LNTH在1%的显著水平下系数为0.387，表示冻品产量每增加10%，用冰量会增加0.387%。

2.2 实证分析 2.2.1 1995—2016年渔业具体产业用冰量估算

根据回归后的方程LNY=0.155×LNFI+0.115×LNSE+0.387×LNTH+10.14，估算出1997—2016这20年间的渔业用冰量，具体见表 5。

由表 5可知，从渔业具体产业来看，1995—2016这20年间，涉渔第一、二、三产业的用冰量总体上均呈递增趋势。其中：涉渔第一产业用冰量预估增长87.52%，年均增长3.04%；涉渔第二产业用冰量预估增长402.69%，年均增长7.99%；涉渔第三产业用冰量预估增长163.87%，年均增长4.73%。

就渔业各产业用冰结构而言，以2016年为例，涉渔第三产业用冰占比最高为47.09%，其次是涉渔第一产业用冰占比为40.74%，涉渔第二产业用冰占比最低为12.17%，如表 6。

从总量来看，1995—2016年间的渔业用冰量预估增长138.13%，年均增长4.22%，而全国水产冷库制冰量实际增长70.44%，年均增长2.57%。我们将全国水产冷库实际制冰量(Y)与估算的用冰量(Y1)进行拟合，结果见图 1。

从图 1看出，全国水产冷库实际制冰量趋势线(Y)呈上升趋势，与估算的用冰量趋势线(Y1)方向一致，表明该估算与实际相契合。图 1中预估趋势线Y1紧紧围绕实际制冰线Y上下波动，且Y1略高于Y，表明实际全国水产冷库制冰量的波动会影响渔业用冰量，从估算结果来看，渔业用冰量总是略高于实际制冰量，表明当前渔业制冰量供给尚有不足，而渔业用冰市场需求较大。

2.2.2 2020—2030年渔业用冰量预测

根据1997—2016年的用冰量估算数据，我们对2020—2030年的用冰量进行预测，以确定未来10年的用冰需求。假设Y_t+k=A1+A2*K，其中下标t表示时间，这里以2016年为基础年。K为滞后期，表示从基年到预测当年的期数，预测结果如表 7。

从预测方程的形式得知，用冰量总体呈逐年递增趋势，利用二次指数平滑法对原序列进行拟合，得到用冰量平滑序列和原序列折线图，如图 2所示。

最后，对2020—2030年的渔业用冰量进行预测，具体结果如表 8。

3 结论

通过建立经典的线性回归模型，估算出1995—2016年渔业用冰量，并与全国水产冷库制冰量趋势线相拟合，得出两者趋势一致，表明估算得到的结果相对客观准确，符合我国渔业用冰的实际情况。其次，通过异方差检验、平稳性检验、协整检验等一系列计量检验，表明回归模型不存在异方差，变量之间存在稳定及长期均衡关系，因而估算结果具有严格的统计意义。从总体趋势上看，我国的渔业用冰量不断增长，说明我国渔业生产规模正呈现不断扩大的趋势。

进一步研究表明，我国涉渔第一产业预估用冰量占比较大为40.74%，但年均增长率为3.04%，低于总用冰量的年均增长率4.22%，表明涉渔第一产业后期的用冰需求增长会放缓。涉渔第二产业预估用冰量在三个产业中占比最小为12.17%，但年均涨幅最高，达到7.99%，远高于总用冰量的年均增长率4.22%，说明涉渔第二产业用冰需求呈现快速增长的趋势。涉渔第三产业预估用冰量占比较大为47.09%，年均增长率为4.73%，略高于总用冰量的年均增长率4.22%，表明涉渔第三产业后期的用冰需求增长呈稳定增长趋势。预计到2030年，渔业总用冰量达到1927万吨。

综上，随着我国渔业生产规模的扩大，渔业用冰需求量将不断增长。渔业生产(涉渔第一产业)和渔业流通与服务业(涉渔第三产业)用冰需求量增长速度低于渔业工业和建筑业(涉渔第二产业)。我国渔业十三五规划明确提出“要推进水产加工业转型升级，重点发展水产品精深加工”。2016年，我国水产加工品总量同比增长3.50%，高于海水产品产量同比增长率(2.36%)1.14个百分点。考虑海洋捕捞生产规模逐年缩减，预计两者增长率的差距将逐年加大。水产品加工业规模的壮大，直接拉动了涉渔第二产业的用冰需求。在涉渔第二产业中，片冰使用量较大，建议制冰企业加大片冰生产规模。同时，要保障制冰用水清洁，以防止用冰过程中细菌滋生，影响水产加工品的品质。