中国综合社会调查抽样设计方案.doc

《中国综合社会调查抽样设计方案.doc》由会员分享，可在线阅读，更多相关《中国综合社会调查抽样设计方案.doc（19页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、中国综合社会调查中国综合社会调查（CGSSCGSS）第二期（第二期（2010-20192010-2019）抽样方案抽样方案2目目 录录一、调查背景.2二、调查目标总体.2三、抽样设计原则.2四、抽样设计中的几个问题.3（一）关于分层.3（二）各阶段抽样单元.3（三）样本量的界定及分配.3五、具体设计.3（一）必选层的样本抽取.3（二）抽选层的样本抽取.3六、最终样本的构成.3七、样本权数的确定.3（一）设计权数.3（二）实际执行情况调整.3（三）最终权数的确定.3八、主要估计量.3（一）各小层内部目标量的估计.3（二）全国目标量的估计.3九、估计量的方差估计.3十、附录：抽选层 100 个样本

2、初级单元.33中国综合社会调查抽样设计方案中国综合社会调查抽样设计方案（2010）一、调查背景一、调查背景我国是世界上人口最多的国家，随着我国市场经济的发展，社会正在发生巨大的变革，经济发展使得对全国社会多方面信息的需求日益突出。及时、全面、客观地了解我国收入、医疗卫生、教育、失业保障等方面的信息，将为国家宏观调控政策的制定和企业的更好发展提供保障，对于促进城乡社会的协调发展、保持国民经济可持续发展都有着重大的现实意义。中国综合社会调查是一项全国性的大型调查项目，调查总体是全国城镇居民和农村居民。本次调查的总体要求：（1）能够全面了解我国城乡社会发展情况；（2）能够对城市群体和农村群体进行对比

3、分析；（3）能够在地理概念或者区域发展水平方面体现我国社会发展的地域差异性。本方案着重介绍此次调查项目的抽样方案设计以及样本数据的权数调整方法。设计方案充分考虑了全国及不同地域估计的需要，对调查总体进行了科学、细致的分层，在一定程度上提高了估计精度；特别是在城乡样本配比、设立自我代表层、样本数据加权等问题的处理上体现了本方案的科学性、高效性、可操作性等特点。二、调查目标总体二、调查目标总体此次调查的目标总体为全国 31 个省、自治区、直辖市（不含港澳台）的所有城市、农村家庭户。三、抽样设计原则三、抽样设计原则首先，作为全国性的抽样调查，整体方案必须是严格的概率抽样，要求样本对全国及某些指定的城

4、市或地区具有代表性。其次，抽样方案必须保证具有4较高的效率，即在相同的样本量条件下，方案设计应使抽样误差尽可能小，调查精度尽可能高。最后，方案必须具有较强的可操作性，不仅便于具体抽样的实施，也要求便于后期的数据处理。四、抽样设计中的几个问题四、抽样设计中的几个问题（一）关于分层（一）关于分层根据本次调查的研究需要，将调查总体分为两大类：一是必选层，该层总体为入选大城市的市辖区家庭户；二是抽选层，该层总体为除去必选层市辖区以外全国所有家庭户。1、必选层根据调查需要及以往调查经验，本次调查对那些发展处于国内领先水平的大城市将特殊对待，将该类城市市辖区家庭户作为单独一层进行设计，作为必选层。对于大城

5、市的界定，着眼于直辖市、省会城市和副省级城市共 36 座城市，从经济水平、教育水平及城市开放性程度等角度出发，选取 GDP、拥有教师总数、外国直接投资（FDI）实际使用外资金额这三个总量指标进行考察，采用因子分析方法确定排名前五的城市进入必选层。需要特别说明的是：该层的调查对象为这些城市的市辖区居民，该层最终调查单元均划为城市家庭户。2、抽选层抽选层的调查总体由必选层以外的城市、农村家庭户组成。为了便于在后期数据分析中采用平衡半样本进行方差估计，对初级抽样单元的分层划分打破省级地域限制，进一步增加分层的层数。将抽选层划分为区层和县层（包含县级市和县） ，采用人口密度、非农业人口比重和人均地区生

6、产总值三个指标，在区层和县层中分别进行因子分析，得到区层和县层内各个区县的综合因子得分；在对综合因子得分进行排名的基础上将区层进一步分为 19 层，县层划进一步分为 31 层，抽选层共计细分划分为 50 个层。（二）各阶段抽样单元（二）各阶段抽样单元本次调查采用分层三阶段概率抽样，视所在层情况，各阶段抽样单元略有5不同，详见表 1。表表 1 各阶段抽样单元各阶段抽样单元第一阶段抽样单元第二阶段抽样单元第三阶段抽样单元必选层必选层街道居委会家庭户抽选层抽选层区、县级市、县居委会、村委会家庭户这样设计的原因在于：对于必选层，选择街道作为初级抽样单元可以细化抽样框，使得样本点相对分散，有利于总体信息

7、的采集，避免由于抽样框过粗而导致样本有偏。对于抽选层，全国区、县级市、县的数量较多，以其作为初级抽样单元比较合适。（三）样本量的界定及分配（三）样本量的界定及分配1、目标样本量目标样本量是指在一定精度要求下，实现调查目标所必需的样本量。根据以往调查经验，本次调查设定目标样本量为 12000 户，其中必选层 2000 户，抽选层 10000 户。后续部分涉及的样本量分配均立足于目标样本量。2、各阶段样本量分配由于必选层与抽选层设计的各阶段抽样单元略有不同，因此，分别介绍两者各阶段样本量的分配。对于必选层，该层总样本量为 2000 户，计划抽取 40 个初级抽样单元（街道） ，每个初级抽样单元（P

8、SU）抽取 2 个二级抽样单元（居委会） ，每个二级抽样单元（SSU）中抽取 25 个家庭户。对于抽选层，该层总样本量为 10000 户，计划抽取 100 个 PSU（区、县级市、县） ，每个 PSU 中抽取 4 个 SSU（居委会、村委会） ，每个居委会（村委会）中抽取 25 个家庭户。最终，本方案共需抽取 140 个 PSU，480 个 SSU。3、样本城乡分配为了能切实反映全国范围城乡社会的真实情况，并在后续研究中能够对城乡进行对比研究，需要保证样本在城乡分配上与实际情况一致。首先明确本方案中的城乡样本概念，本方案设计基于这样的假设：居委会的家庭户为城市居6民，村委会的家庭户为农村居民。

9、统计资料显示，目前我国城市常住人口数与农村常住人口数基本持平，由于城市居民主体的各方面差异相对明显，方差较大，因此将样本量的城乡分配比例确定为 64。根据第一阶段样本量的分配结果，必选层中共抽取 80 个居委会，共计2000 户，因此抽选层内城乡家庭户数需分别为 5200 和 4800 户才能满足 6：4的要求。由于抽选层每个 PSU 下抽取 4 个 SSU（居委会或村委会） ，每个 SSU内最终抽样单元的目标样本量均为 25，因此对城乡样本比例的控制，主要是使得抽选层居委会与村委会下的城乡样本比例达到 5200：4800，也即抽选层中的居委会与村委会样本个数比约为 208：192。为了实现样

10、本二级单元 208：192 的目的，需要根据样本初级单元的城市化水平（非农业人口比重）分配样本居委会和村委会的数量。这里采用分别在区层和县层样本初级单元内，根据各个区县的城市化水平（用非农人口比重表示）进行分段，形成若干个区间，然后根据不同的区间对样本区县中的居委会、村委会个数进行分配的方法来实现样本居委会和村委会的比例要求。4、接触样本量如果回答率达到 100%，则调查时需要的接触样本量即为有效样本量，但现实中无回答现象不可避免。要克服这个困难，可以采取替换样本的方法，但样本替换在实际操作中存在缺陷，因此本方案采取利用膨胀系数扩大样本量的方法，对第三阶段样本量进行放大。根据往年调查经验，发达

11、城市的市辖区居民由于种种原因，回答率在 50% 左右，即膨胀系数在 2 左右，因此在必选层每个二级单元抽取 50 户家庭，该层接触样本量扩大至 4000；抽选层的居民群体回答率高于必选层，但其内部还有差异，大体上城市居民的回答率在 65% 左右，农村居民的回答率高于城市居民，大致在 85%左右，因此，对于抽选层，在每个居委会抽取 38 户，在每个村委会内抽取 30 户。五、具体设计五、具体设计（一）必选层的样本抽取（一）必选层的样本抽取71、必选层入样城市的确定首先对全国 36 座城市（包括直辖市、省会城市、副省级城市）市辖区的GDP、拥有教师总数、外国直接投资（FDI）实际使用外资金额这三项

12、指标进行因子分析，最终确定 5 个城市进入必选层。这里给出进行因子分析之后，综合得分排名前五的城市名单（见表 2） 。表表 2 36 座城市中前座城市中前 5 座城市座城市排名排名城市名称城市名称1上海2北京3广州4深圳5天津2、初级单元的抽取根据最新的全国行政规划，这 5 座城市市辖区总数为 67 个，如果以 67 个市辖区为初级阶段抽样框，初级单元个数有限，可能会导致样本在分布上过于集中，影响样本对该层总体的代表性，为了使样本点分布较为分散，将初级抽样单元细化为街道，由此，必选层的抽样框为入选城市市辖区的街道，且该层人口规模均只采用城市人口。以必选层入选城市市辖区的街道作为抽样框，以各街道

13、的城市人口规模为辅助信息，采取与人口规模成比例的 PPS 抽样抽取 40 个街道作为该层的PSU，这里：最终抽样单元为家庭户，理论上应该以街道的户数作为辅助信息，但目前这方面信息缺乏，因而用街道人口数作近似处理。具体步骤说明如下：设共有个初级单元（街道） ，第 个初级单元的人口Ni规模为（） ，首先按地理区域将街道进行排序并产生累计人口数iM1,2,iN列，即排序第一位二级抽样单元为，对应的累积人口为；排序第二位的1M1M单元人口数为，对应的累积人口为，依次进行，则排序第 的单元2M12MMi对应的人口数为，对应的累积人口为（） ，需要注意的是：iM1ii jM1,2,iN8这里的 是指排序后

14、单元对应序号。接下来，对第一个单元赋以共个i11M1M代码；对第二个单元赋以共个代码，对第 个单元赋1121MMM 2Mi以共个代码，人口总数为。1111iijj jjMM iM0 1Ni iMM若已知初级单元样本量为，则首先可以确定抽样间隔，理论上，nk0Mkn实际中一般取最接近于的整数。接下来从范围内随机地产生一个整数0M n1k作为抽样起点，则代码为所对应的单元即为被抽中的初r,(1)r rkrnk级单元，且第 个抽样单元被抽中的概率为。i0i iMZM3、二级单元的抽取在确定初级抽样单元（街道）后，二级单元抽样框为样本街道内所有居委会名单，同时收集各个居委会人口数。二级单元的抽取方法：

15、在入选街道内，采用与各居委会人口规模成比例的系统 PPS 抽样抽取 2 个居委会。具体抽取过程同初级单元部分。4、最终单元的抽取最终单元的抽取在每个入选的居委会中进行，在每个居委会内，按家庭户的门牌号进行排序，采取等概率系统抽样抽出 50 户家庭作为最终调查单元。这里，实际接触样本量在调查时视回答率状况分批投放，以达到既满足目标样本量，同时各二级样本单元（居委会）下最终样本（家庭户）数量差异不至于悬殊。必选层样本抽取情况见表 3。表表 3 必选层样本抽取情况必选层样本抽取情况区域初级单元数居委会数目标样本量实际接触样本量必选层40402=8040225=200040250=4000（二）抽选层

16、的样本抽取（二）抽选层的样本抽取抽选层的初级抽样单元总体为：除去必选层城市市辖区以外的区、县级市、县，根据中国统计年鉴 2009的行政区划设置信息，除去必选层 67 个区后，9抽选层共有 792 个区和 2003 个县级市和县。1、抽样框的构建为了实现样本分布与总体分布的一致，根据全国分县市人口统计资料2008 年的户籍人口统计，区层和县层的人口规模比约为 7：18，差异较大，为了区县两层样本初级单元分配数量差异，采用与人口规模平方根成比例的分配方法，将抽选层中的 50 个层划分为为 19 个区层和 31 个县层。在区层和县层中，以人口密度、非农业人口比例、人均地区生产总值三个重要指标对区县进

17、行因子分析，首先利用综合因子得分法对区层和县层的区县进行排序，然后根据区县的个数基本等分为 19 层和 31 层。抽选层中区层共有792 个初级单元不能被 19 整除，将余数 13 归入区层的最后一层，结果为前 18个区层每层含 41 个初级单元，第 19 层含 54 个初级单元。县层做法类似，2003个初级单元不能被 31 整除，将余数 19 归入县层的最后一层中，结果为前 30 个县层每层含 64 个初级单元，第 31 层含 83 个初级单元。具体分配结果见附录excel 表。2、初级单元的抽取在上述 50 个小层中，以各个初级单元的综合因子得分排序，各个初级单元的人口数为辅助信息，按照与

18、各初级单元人口数成比例的系统 PPS 抽样方式分别在各小层中抽取 2 个样本区或县。3、二级单元的抽取在每个样本区或县中，抽取 4 个村委会或居委会，二级单元的具体抽取要根据初级单元的类别来确定，这是本方案实现城乡样本配比的关键环节，根据入样初级单元的类别确定应在该初级单元内抽取的居委会、村委会数目。一般说来，城市化水平较高的区（县） ，居委会个数较多、村委会个数较少，因此，城市化水平不同的区县，居委会与村委会抽取个数之间的配比亦不同，经测算，具体的居委会、村委会分配标准见表 4。表表 4 初级单元内二级单元的分配标准初级单元内二级单元的分配标准分配标准分配标准居委会居委会村委会村委会95%及

19、以上4050%95%3115%50%221015%以下13在每个入样初级单元内部分别构建居委会、村委会抽样名单，同时收集该居、村委会人口数。根据上表的分配标准，以人口数作为规模辅助信息，按照与人口规模成比例的 PPS 抽取相应个数的居委会、村委会，确保每个初级单元下有 4 个二级单元。4、最终户单元的抽取最终户单元的抽取在每个入选的居委会或村委会中进行，要求每个二级单元内部达到 25 户的目标样本，最终抽选层目标样本量为 10000。但在实际调查中要涉及接触样本量的膨胀问题，正如在样本量界定中提到的，最终接触样本量在每个居委会（村委会）内必须进行扩大，由于城乡居民的回答率不同（大约分别在 65

20、%、85% 左右） ，因此将居委会内的接触样本量扩大至 38 户，每个村委会内的接触样本量扩大至 30 户。具体抽取方法仍采用等概率系统抽样。5、户内调查对象的抽取在入选的户内，列出所有 18 岁及以上人口，随机抽取一人，作为最终的调查对象。六、最终样本的构成六、最终样本的构成本次调查目标样本量为 12000，必选层样本量为 2000，抽选层样本量为10000。在必选层中，抽取 40 个初级单元，每个初级单元内抽取 2 个二级单元，在每个二级单元内目标样本量达到 25 户；在抽选层中，共抽取 100 个初级单元，400 个二级单元，其中居委会 208 个，村委会 192 个。全部目标样本城乡比

21、为64，基本与我国实际情况相吻合。本次调查的最终接触样本量为 17664，其中必选层 4000，抽选层 13664。在必选层中，在初级单元、二级单元个数不变的情况下，将每个二级单元内部的接触样本量增加至 50 户；在抽选层中，初级单元数 100、每个初级单元内二级单元数为 4，在每个入样居委会内将接触样本量扩大至 38 户，在每个入样村委会内部将接触样本量扩大至 30 户，因此，抽选层接触样本量为 13664，其中城市居民样本为 7904，农村居民样本为 5760。11七、样本权数的确定七、样本权数的确定（一）设计权数（一）设计权数设计权数是指每个样本单元在调查总体中所代表的单元数，它等于该样

22、本单元入样概率的倒数。由前面确定的抽样方法，必选层与抽选层样本的设计权数形式相似，但表达的内容有所区别。统一的设计权数公式可表达为：（）125hijhhi hijk hihijDMMwMM1,2,3,51h 其中，下标 1 代表设计权数，代表必选层，分别对应于1h 2,3,51h 抽选层中的 50 个小层。代表人口规模总数，代表每个二级单元中包含的MD三级单元数（户数） 。其他各符号及下标的含义分别解释如下：在必选层中， 代表第 个初级单ii元（街道） ，代表第个二级单元（居委会） ，代表第个三级单元（家庭户）jjkk；在抽选层中， 代表第 个初级单元（区、县） ，代表第个二级单元（居委iij

23、j会、村委会） ，代表第个三级单元（家庭户） 。kk（二）实际执行情况调整（二）实际执行情况调整本调查为了保证目标样本量，对接触样本量进行扩大，调查的执行力度将影响调查样本的权数，因此需要再次进行权数调整。设某个居委会（村委会）最后样本调查的执行量为，则该家庭户的实际执行调整权数为：n225hijkwn（三）最终权数的确定（三）最终权数的确定每个最终调查单元的最终权数等于每个家庭户的初始设计权数与执行情况调整权数的乘积，即：12hijkhijkhijkwww其中，（ ）代表第小层，1 代表设计权数，2 代表执行h1,2,3,51h h调整权数。八、主要估计量八、主要估计量12利用加权调整后得到

24、的最终权数，可以得到各种类型目标量的估计公hijkw式。在这里，我们记第小层第 个样本区、县（第 1 小层为样本街道）第个hij样本居委会、村委会中第个家庭户中被调查者的目标变量为。khijky（一）各小层内部目标量的估计（一）各小层内部目标量的估计1、总体总值的估计Yrhhhijkhijk SYwy其中，是第小层内回答单元所组成的集合。rhSh2、总体均值的估计Yrhrhhijkhijk S h hijk Swy Yw3、总体比例的估计Prhhrhhijk SC h hijk Sw Pw其中，是第 小层内具有某种特征的单元集合。hCh4、总体比率估计量的估计rhrhhijkhijk S h

25、hijkhijk Syw Rxw其中，、为该比率估计量所涉及的两指标的样本回答值。hijkxhijky5、域估计针对特定群体的估计定义示性函数10hijkI 不属于该群体属于该群体该特定群体均值估计量为：13rhrhhijkhijkhijk S Th hijkhijk SIyw YIw（二）全国目标量的估计（二）全国目标量的估计1、总体总值的估计Y0rhzhijkhijk hShMYwyM1,2,3,51h 其中，表示全国总人数，表示第层的人口总数。0MhMh2、总体均值的估计Y0rhzh hShMYYM3、总体比例的估计P0rhhrhhijk hSChhijk hSMwMPw 其中，是第 层

26、内具有某种特征的单元集合。hCh4、总体比率估计量的估计00rhrhhijkhijk hShhijkhijk hShMywMRMxwM其中，、为该比率估计量涉及的两指标的样本回答值。hijkxhijky5、域估计针对特定群体的估计定义示性函数10hijkI 不属于该群体属于该群体该特定群体均值估计量为：1400rhrhhijkhijkhijk hSh Thijkhijk hShMIywMYMIwM九、估计量的方差估计九、估计量的方差估计由于本方案采取的是分层三阶段复杂抽样，采用常规的方差计算方法，过程十分繁杂，且需要用到大量的辅助信息。结合本次调查将初级单元划分为 50个层，每层抽取 2 个初

27、级单元的设计特点，我们将采用复杂样本方差估计方法中的平衡半样本方法（BRR）对目标量的方差进行估计，从而达到简化计算复杂样本方差估计的目的。假设对总体，采用分层随机抽样，每层有放回地简单随机抽取 LhhNN12 个单元，设和是第层的样本观测值（=1，） 。现从每层中抽1hy2hyhhL取一个单元，组成一个半样本，总共有个这样的半样本。L2的基于第（=1，2，）个半样本的估计为：YL2= ,sty1122 1Lhhhhh hWyy其中，=。11,0, hh 如果第层中的第一个单元被选入第个半样本其它2h11h所有个这样的估计量的平均恰好是，即L2sty2, 11 2LststLyy 假设这个较小

28、的半样本子集包含了个半样本，有：k=+ stkyv kststyyk12 ,1 styvhLhLhhhhhkhhddWWk 11)()( 21因此,如果所选择的个半样本对于所有的都满足条件khhL, 1=0 khh 1)()(15则就正好等于。从而，基于这个半样本便可实现对的估计， stkyv styvk styV而不必基于全部的个半样本。满足上式的组半样本称为平衡半样本。L2kFay 对平衡半样本进行了修正，即 Fays BRR，引入 Fay 因子，即01估计量变为：= ,sty1122 1Lhhhhh hWyy11, hh 如果第层中的第一个单元被选入第个半样本其它方差估计量变为： stk

29、yv kststyyk12 ,2)1 (1当层数 L 很大时，即使利用正交矩阵得到完全正交的平衡半样本个数已经大大减少，所需要的计算量仍然十分庞大，此时可用部分平衡半样本的方法，即 k 部分平衡半样本方法。 方差估计量是： stkyv kststyyk12 ,2)1 (1在本次调查中，抽选层在 50 个区县层中共抽取 100 个样本区县，在每个样本区县中平均抽取 4 个居委会或村委会。根据抽样设计，有两种平衡半样本方法：1.1.初级样本单元（区县）层次上的平衡半样本方法初级样本单元（区县）层次上的平衡半样本方法按照抽样设计，将抽选层的区县分成 50 层，在每一层中抽取一定数目的样本区县。若是在

30、此基础上应用平衡半方法，则需要将每一层的区县随机分成 2组。此处我们选取 60 或 64 阶的正交半样本随机组，利用分组平衡半样本方法，得到目标变量的估计量和方差估计量分别是：目标变量的估计量是= ,sty1122 1Lhhhhh hWyy其中160.51 0.51(),1(), hhh h hhhnmhmnm m 如果第层中的第一个单元被选入第个半样本其它方差估计量是： stkyv kststyyk12 ,2)1 (1此处选取 Fay 因子.0.50或相关系数和比例估计量和方差估计量计算方法同上。2.2.二级样本单元（居委会或村委会）层次的平衡半样本方法二级样本单元（居委会或村委会）层次的平

31、衡半样本方法若以 100 个区县为层，在每个层中将 4 个二级样本单元（居委会或村委会）随机分成 2 组，利用分组平衡半样本方法得到目标变量的估计量及其方差估计量。我们利用阶正交 Hadamard 矩阵得到完全正交的平衡半样本。由此，104104得到目标变量的估计量及其方差估计量分别是估计量是= ,sty1122 1Lhhhhh hWyy0.51 0.51(),1(), hhh h hhhnmhmnm m 如果第层中的第一个单元被选入第个半样本其它方差估计量是： stkyv kststyyk12 ,2)1 (1上述计算以目标变量的均值估计量为例，相关系数和比例估计量和方差估计量计算方法同上。十

32、、附录：抽选层十、附录：抽选层 100 个样本初级单元个样本初级单元省份省份市市区县区县小层小层大层大层居委居委 会会村委村委 会会浙江省杭州市上城区1区层40辽宁省沈阳市和平区1区层4017江苏省常州市钟楼区2区层40浙江省温州市鹿城区2区层31山东省济宁市市中区3区层31浙江省杭州市江干区3区层31江苏省南京市六合区4区层31浙江省杭州市余杭区4区层22吉林省长春市朝阳区5区层31山西省太原市迎泽区5区层40湖北省武汉市江夏区6区层40贵州省贵阳市云岩区6区层31河北省唐山市丰南区7区层22黑龙江省哈尔滨市道外区7区层31江苏省连云港市新浦区8区层31青海省西宁市城北区8区层31四川省成都

33、市温江区9区层31浙江省湖州市吴兴区9区层31吉林省吉林市昌邑区10区层31山东省泰安市泰山区10区层31山东省日照市岚山区11区层22西藏自治区拉萨市城关区11区层31辽宁省铁岭市银州区12区层31辽宁省葫芦岛市连山区12区层22广西省柳州市柳北区13区层31辽宁省朝阳市双塔区13区层31广东省汕头市潮南区14区层40贵州省遵义市红花岗区14区层31江西省上饶市信州区15区层31海南省海口市秀英区15区层22四川省乐山市五通桥区16区层22重庆市重庆市涪陵区16区层22四川省泸州市龙马潭区17区层22广西省玉林市玉州区17区层22河南省信阳市浉河区18区层22河南省商丘市睢阳区18区层22湖

34、北省孝感市孝南区19区层22贵州省安顺市西秀区19区层22山东省济宁市汶上县1县层31福建省漳州市东山县1县层31江苏省扬州市江都市2县层22湖北省襄樊市老河口市2县层31山东省烟台市莱州市3县层22福建省三明市永安市3县层3118山西省运城市河津市4县层22江西省宜春市丰城市4县层22河南省漯河市临颍县5县层22湖北省襄樊市宜城市5县层22吉林省吉林市桦甸市6县层22黑龙江省黑河市逊克县6县层22河北省沧州市肃宁县7县层22河南省周口市郸城县7县层13江苏省徐州市新沂市8县层22黑龙江省七台河市勃利县8县层22吉林省白城市洮南市9县层22安徽省阜阳市太和县9县层13山西省长治市襄垣县10县层

35、22湖南省邵阳市邵东县10县层13河北省邯郸市磁县11县层13浙江省宁波市奉化市11县层22四川省南充市营山县12县层22江西省鹰潭市余江县12县层22河南省周口市扶沟县13县层13陕西省延安市志丹县13县层22四川省眉山市丹棱县14县层22湖南省怀化市洪江市14县层22黑龙江省哈尔滨市双城市15县层22黑龙江省哈尔滨市延寿县15县层22安徽省蚌埠市固镇县16县层13河南省洛阳市宜阳县16县层13山东省临沂市沂水县17县层13黑龙江省齐齐哈尔市克山县17县层22宁夏自治区吴忠市同心县18县层22安徽省宣城市绩溪县18县层22湖南省邵阳市隆回县19县层13陕西省汉中市城固县19县层22江西省赣州

36、市兴国县20县层22湖北省襄樊市南漳县20县层22四川省宜宾市宜宾县21县层13福建省福州市永泰县21县层22湖北省恩施土家族苗族自治州来凤县22县层13重庆市重庆市奉节县22县层13安徽省黄山市歙县23县层13江西省赣州市会昌县23县层22内蒙古自治区呼和浩特市土默特左旗24县层13吉林省松原市扶余县24县层1319广西省崇左市大新县25县层13陕西省延安市宜川县25县层22新疆自治区喀什地区英吉沙县26县层13陕西省铜川市宜君县26县层22云南省玉溪市新平彝族傣族自治县27县层13湖南省郴州市汝城县27县层13甘肃省兰州市榆中县28县层13重庆市重庆市城口县28县层13云南省大理白族自治州剑川县29县层13广西省柳州市三江侗族自治县29县层13云南省昆明市寻甸回族彝族自治县30县层13湖南省怀化市中方县30县层13云南省昭通市巧家县31县层13甘肃省庆阳市镇原县31县层13 合计208192