DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 3. 遗传算法运算 遗传算法主要通过五个步骤对复杂问题进行运算，包括编码、适应度函数、选择、交叉、变异从而得到最接近最优解的近似解，能够对非线性目标问题进行优化，同时对多个个体进行计算，有效避免陷入局部最优的情况出现。 (一) 编码 通过将目标个体转化为染色体的形式并以字符串表示，将该染色体放置在约束条件所形成的解空间中进行运算。编码直接决定了后续步骤的实现方法，常用的编码方法包括Grey编码、实数编码、多级参数编码等等。本论文采用的是实数编码，不需要进行编码的转换。 (二) 适应度函数 适应度函数是评价群体中个体好坏的标准，是进行自然选择的唯一依据，一般是由目标函数加以变化得到，根据求解最大值还是最小值问题的不同，适应度的选择也不同。 如6，3，11点需求量分别为4，1，2吨，车辆负重为6吨。 当配送路线计算到11时，超过负重上限，此时得到第一条配送路线，对比第一个需求点与哪个配送中心的距离最接近后，以此为起点，同理寻找返回的配送中心点，并计算其里程，从点11开始计算第二条配送路线，重复上述过程。 (三) 选择操作 选择操作是从就群体中选择优良个体组成新的种群，得到下一代种群，个体的选择概率与适应度相关，个体适应度值越高，被选中的概率越大。遗传算法选择操作有轮赌盘法、锦标赛发等多种方法。 (四) 交叉操作 交叉操作指从种群个体中最近选择两个染色体，通过两个染色体的交换组合，把优秀的基因遗传给子代，从而产生新的优秀个体。因为个体采用实数编码，所以交叉操作采用实数交叉法，第K个染色体 a k 和第L个染色体 a l 在j位的交叉操作方法，b为[0，1]区间的随机数。 随机生成2个染色体如下： 4 12 6 7 2 5 10 14 13 19 3 8 11 3 12 13 9 5 4 14 6 11 10 78 1 2 随机生成[0，1]之间的数，当这个数字小于交叉概率Pc时，与临近的染色体进行交叉，例如为0.3，然后在[1，5]之间随机选取一个整数为1，就从染色体中位于第一位的基因进行交叉，将两者找到相同的基因进行交换得到： |3|126 7 2 5 10 14 13 1 9 |4| 8 11 |4|12 13 9 5 |3|14 6 11 10 7 8 1 2 然后从第二位的基因继续进行交换，循环10次，从而得到交叉之后的染色体，并将产生的个体存储到子代的种群中，然后对产生的子代进行适应度计算，当子代的适应度比父代要差时，保留父代，保证不会出现退化解。 (五) 变异操作 变异操作的主要目的是保持种群多样性。变异操作从种群中随机选取一个个体，选择个体中的一点进行变异从而产生更加优秀的个体，统称下一代的群体为子代，上一代群体为父代。 将交叉后的子代进行变异操作，随机生成[0，1]之间的数，当其取值大于变异概率Pm时，对该染色体进行变异操作，随机生成一个数组，将该数组的第一位和第二位的数字所对应的交叉后产生的染色体的基因的位置进行调换。 (六) 终止条件 一般情况下，有以下几种情况发生时计算会停止：达到目标时，达到设置的时间，达到设置的迭代次数，在相同多代内没有任何变化时，系统自动停止。本文使用的方法是设置规定的迭代次数，当未达到最大次数时会重复上面的运算，然后选择适应度最好的个体进行输出。
system: [SYSTEM]

user: 4. 用遗传算法研究配送中心路径问题 其配送中心的坐标分别为(40，50)、(60，40)，其他需求点的具体情况如下图2、图3所示。其中需求点的时间约束条件如表1所示。 图2. 配送点需求量及配送时间要求 图3. 配送点距需求点距离分布 Table 1 Demand point time constrain 需求点序号 最早到达时间 最晚到达时间 1 4:00 5:00 2 3:30 6:00 3 4:30 5:50 4 3:20 5:00 5 3:00 4:30 6 4:30 6:00 7 3:20 4:30 8 3:30 5:50 9 4:30 5:48 10 3:36 5:30 11 4:06 5:30 12 3:45 5:00 13 3:06 5:15 14 3:10 5:00 15 3:30 5:00 16 3:12 5:48 17 3:12 4:48 18 4:24 6:48 表1. 需求点时间约束 各个需求点距离配送中心的距离如图3所示，如果与配送中心1的距离小于配送中心2的距离第三列为1，否则为2。 在求出各个配送点之间的距离之后，使用matlab进行计算，设定种群大小为150，迭代次数为200次，交叉概率为0.8，变异概率为0.3。 通过运算，迭代过程如图4所示。 图4. 迭代过程 得到优化后路线如表2，表3所示。 Table 2 Before optimizatio 起点 需求点 终点 车辆1 2 15-4-17 2 车辆2 1 5-2 2 车辆3 1 9-12 2 车辆4 2 13-11 2 车辆5 1 8-7-14 2 车辆6 2 18-10 2 车辆7 2 16-1 2 车辆8 1 6-3 1 表2. 优化前路线 Table 3 After optimizatio 起点 需求点 终点 车辆1 1 9-11 2 车辆2 1 7-8 2 车辆3 1 5-6 2 车辆4 1 2-10 2 车辆5 2 13-12 2 车辆6 2 17-16 2 车辆7 2 14-15-18 2 车辆8 1 3-4-1 1 表3. 优化后路线表 优化前、后路线如图5、图6所示。 图5. 优化前路线图 图6. 优化后示意图 优化前该市配送中心一共需要8辆车，但很大程度上无法满足需求点的时间窗要求，加入惩罚成本后，总成本为64，124元，在通过是用遗传算法进行优化之后，不仅没有增加车辆的使用，而且能在很大程度上满足客户对于时间的要求，加上时间窗惩罚成本后总成本为24，541元，在配送成本上有显著的降低，达到了对配送中心配送路线优化的目的。
system: [SYSTEM]

user: 5. 结论 本文利用了遗传算法的全局性和自适应性对NP问题的最优解进行寻找，以石家庄市某物流企业所提供数据作为基础数据，着重考虑用户的满意度，在相同成本的基础上，尽可能的满足客户的需求，通过建立模型并进行求解，科学地规划配送路线，将优化前后的路线方案进行对比，结果表明获得了较好的优化方案。 本文目前只针对企业较突出的问题进行研究，未来会通过加入更多约束条件和优化算法完善内容，比如应用神经网络算法、数据分析模型进行需求预测，帮助提前安排用车，利用改进的遗传算法提高路径优化效率等。
system: [SYSTEM]

user: 摘 要 居民参与药品回收积极性低是目前制约着过期药品回收活动的最大阻碍，本文拟通过对过期药品的回收激励策略进行研究，来提高居民对过期药品回收活动的积极性。本文首先通过实证研究找出了影响居民参与回收活动的关键因素，对政府今后可采取的激励政策进行了理论分析。其次，将“系统动力学”模型引入过期药品回收领域，构建了量化规范的模型，并将模型应用于江西省瑞安市的过期药品回收的实际情形之中，对现实中的过期药品回收系统运行进行模拟和仿真，考察了不同的回收激励政策的实施效果，为政府日后制定有效科学的激励策略提供参考。 关键词 :系统动力学，政府，家庭过期药品，激励策略 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 绪论 随着人民生活水平的提高，一些常备药成为了家庭之中的必需品之一。而居民对于家中的过期药品，最普遍的处理方法就是当做普通的生活垃圾一样丢弃，这恰恰是处理过期药品最不正确的做法。对于过期药品来说，最合适的处理方法是对其进行集中回收和处理。目前我国一些大型医药企业如“广药白云山”等公司，曾自发地开展过居家过期药品回收活动，但由于过期药品的回收需要高额的运营成本来支持且短期内无法获得收益，因此企业缺乏长期回收过期药品的原动力。况且企业主导下的过期药品回收还存在着回收面窄，居民信任度不高，回收流程不规范等问题 [ 1 ]。 政府部门本就肩负着维护生态环境，保证人民生命和财产安全的责任，因此政府应自觉担负起回收居家过期药品的重任，自上而下地统筹协调好回收工作。近年来政府一直在积极开展过期药品的回收工作，将其视为保障和改善民生的举措之一，并在回收点设置、回收规范等方面取得了进展。但同时不可否认的是居民的参与积极性不高，回收率较低仍然是制约我国过期药品回收工作发展的最大障碍 [ 2 ]。 居民是回收过程的最关键主体之一，而居民的积极性不高是制约回收工作开展的最大障碍。基于这个原因，本文首先对影响居民回收意愿的因素进行了初步的筛选，随后通过调查问卷的方式进一步进行验证，得到影响居民回收意愿的关键因素。经过筛选以及问卷调查，得到了如下四个影响居民回收意愿和行为的关键因素： 1) 经济因素：居民在处理过期药品时，追求的是自身利益最大化，因此政府若是给予居民一定的补贴，会刺激他们的回收意愿；另一方面，非法的药品回收的价格若是提高，那么居民在利益的驱使下会将药品卖给不法商贩 [ 3 ]。 2) 环保意识与知识水平：以往的研究已经证实，居民的环保意识和知识水平对个体废弃物处置行为存在显著作用关系 [ 4 ]，是否拥有相关回收和回收项目的足够知识是家庭在参与回收项目上重要的决定因素 [ 5 ]。这两点在一定程度上决定了居民对于过期药品危害认识的程度。若居民充分意识到了过期药品的危害，肯定会自发地参与到回收工作中。 3) 回收便利性：回收的便利性高，则节约了居民的时间成本和精力，那么居民自然会乐意参与到回收行动中来。回收便利性可视为由居民参与药品回收所付出的时间成本、精力成本及心理成本所组成的多重构面。当回收便利性较高，节省了居民参与过期药品回收的时间和精力，其回收意愿也会提升；相反地，当回收便利性较差，其回收意愿则会降低 [ 6 ]。 4) 回收点服务水平：回收点工作人员是与居民直接接触的群体，工作人员的工作态度和服务水平直接影响了居民的回收体验和效果，居民的体验不好在一定程度上也会影响回收积极性 [ 7 ]。因此，回收渠道对消费者的回收体验与回收成效有着重要的影响。 针对上述四个关键因素，设计出如下具有针对性的策略： 1) 增设药品回收点；2) 加强回收点的监管；3) 对居民实施回收补贴；4) 加大回收宣传力度。 由于在回收系统中，存在着居民、政府、正规药品回收点、私人回收点等多个主体，各个主体之间相互影响制约，形成了动态的、具有多重反馈的复杂系统，单凭主观的臆想和经验判断是无法正确理清主体之间的关系并判断回收策略实施的效果的。因此本文将引入系统动力学理论来对过期药品回收系统进行动态仿真，观察在不同的策略下系统的动态变化，使得结果更加科学可靠。
system: [SYSTEM]

user: 2. 系统动力学模型 上一章从理论的角度探讨分析了家庭过期药品回收的策略，本章将详细介绍系统动力学的思想和方法研究家庭过期药品回收策略的可行性，为后续的模型建立打下基础。 2.1. 系统动力学的原理与特点 系统动力学(System Dynamics，简称“SD”)是一门以反馈控制理论为基础，以计算机仿真技术为手段，适用于研究复杂系统的结构、功能与行为之间的辩证对立统一关系的科学 [ 8 ]，是系统科学与管理科学的一个重要分支。系统动力学强调系统、整体的观点，认为系统行为主要是由系统内部要素间的反馈机制决定的，极力从微观结构入手把握系统内部结构、参数及总体功能，进而分析和研究系统的动态行为 [ 9 ]。系统动力学的主要特点表现在以下几个方面： 1) 系统动力学研究问题注重从因果机制出发。因果关系是客观现象之间引起与被引起的关系，是现象之间普遍联系的表现。从系统内部要素间的因果关系出发，将复杂的现象化简为一个个反馈环，才能理清现象产生和发展的脉络，从而了解系统形成的内在原因与机制。 2) 系统动力学模型的建立基于系统的反馈结构。不同于参数型模型着眼于高质量统计数据的获取，系统动力学更注重通过系统结构来把握系统的行为，即模型的建立是基于系统的反馈结构而不是统计相关性。因此系统动力学模型也被称为结构依存型模型，其突出优点在于能处理非线性、高阶次、多重反馈的复杂问题，在数据不足或某些参量难以量化时仍可以反馈环为基础可以做一些研究，适用于解决对数据精度要求不高的、复杂的社会经济问题。 3) 系统动力学可实现对真实系统的模拟。在计算机仿真技术的支持下，系统动力学模型可以成为社会、经济、生态等复杂系统的“实验室”，实现对真实系统的模拟，从而对系统的长期发展趋势进行定量研究。 4) 系统动力学能够进行政策仿真。系统动力学模型可以考察系统在不同参数或不同策略因素输入下的行为变化，从而实现对系统在不同情境下的动态仿真，使得决策者可通过观察模拟结果寻求解决问题的更优策略。 2.2. 系统动力学解决问题的主要步骤 系统动力学是一种定性分析与定量分析相结合的方法，其中，定性分析是定量分析的依据，定量分析是定性分析的具体化。从定性-定量的研究过程来看，利用系统动力学解决问题的过程大致可分为五步。第一步，对系统进行总体分析，界定系统边界，聚焦具体的研究问题并明确解决该问题需要达到怎样的目的；第二步，对系统进行结构分析，把握系统要素间的因果关系，确定系统总体与局部的反馈机制；第三步，进一步区分变量性质，建立变量间的数学关系，从而建立规范、定量的系统动力学模型，并对模型参数进行估计与确定；第四步，对系统动力学模型进行检验和评估，并在此过程中修改调试模型，确保模型的有效性；第五步，利用模型进行实际系统模拟与政策仿真分析。 2.3. 应用SD研究家庭过期药品回收策略的可行性分析 系统动力学应用于家庭过期药品回收领域研究的可行性是由被研究系统的特性和系统动力学的特殊优势决定的。本文所要研究的政府主导下的家庭过期药品回收系统主要具有以下几个特性： 1) 家庭过期药品回收系统具有复杂性。家庭过期药品回收系统包含众多要素，各要素之间相互关联、相互作用，使得系统内部形成多重反馈环，结构较为复杂，因而系统行为可能具有反直观性，即行为方式可能与多数人所预期的结果不一致； 2) 家庭过期药品回收系统是一个非线性的系统。各项回收激励策略对过期药品回收量的影响通常是高度非线性的，面对这些非线性关系，传统的数学方法没有办法很好地描述他们之间的定量关系； 3) 社会经济系统运行通常具有政策滞后性，即系统行为不会立刻对政策变化做出反应。因而家庭过期药品回收策略的实施与其效果出现之间往往也存在较长的时间延迟； 4) 家庭过期药品回收系统具有时变性。随着时间的变化，系统内、外部影响因素也在变化。时间不同，输出的结果也不同，系统动态性较为显著。 正是由于以上几个特点，使得家庭过期药品回收系统非常适合用系统动力学来进行研究。系统动力学可以依据系统内部各因素之间的因果关系，构造出能反映非线性、多重反馈关系和时间延迟的动态模型，并利用计算机仿真技术实现系统行为的再现。同时，可以改变系统参数和决策变量来考察系统行为的动态变化，从而更直观地了解不同回收激励策略的实施效果，为政府设计更优的回收方案提供参考。
system: [SYSTEM]

user: 3. 模型设计 本章将建立过期药品回收系统的动力学模型，分析、量化各要素间的关系，为之后的仿真模拟打下基础。 3.1. 模型假设 1) 假设本研究不涉及过期药品产生的因素，只研究药品的回收过程；即每户每年产生的过期药品数量是一定的； 2) 所有回收后的药品均被全部销毁，不会再次流入市场； 3) 系统内的回收点均持续经营，不存在撤出系统的情况； 4) 居民家中未被正规或非法回收的过期药品，若干年后随生活垃圾被居民丢弃。 5) 为了定量化各要素，将用各项策略的成本来表现政府对该项激励策略的投入力度，用过期药品回收量来表现激励策略的实施效果。另外，居民环保意识需要长期积累和影响，因此本文选取宣传总费用作为正规药品回收量的影响因素之一。 3.2. 过期药品回收因果反馈模型 因果回路图是一种定性描述系统结构的方法，用图示方法探究系统的动态形成内因和要素间的反馈关系 [ 10 ]。基于前文的分析和假设，绘制出家庭过期药品回收系统的因果回路图(图1)。 图1. 家庭过期药品回收系统因果回路图 3.3. 过期药品回收存量流量模型 系统存量流量图(图2)以图示的形式描述变量间的量化关系与控制过程，明确系统的反馈与动态变化规律。由于药品回收量与其影响因素之间呈非线性关系，无法使用清晰的表达式进行描述。因此，本文借鉴了以往学者研究垃圾回收问题的处理方式，引入了若干“影响因子”变量 [ 9 ]，通过建立表函数关系式来体现相关变量对于每年药品回收量的影响程度。图中速率变量“正规药品回收速率”表示每年通过正规渠道回收的家庭过期药品数量(单位：盒(瓶)/年)，是本文研究的重点考察变量。 图2. 家庭过期药品回收系统存量流量模型 为了模型仿真，必须对各变量间的数学关系进行描述，即根据系统流图写出系统动力学方程。本文建立的政府主导居家过期药品回收体系存量流量模型可划分为居家过期药品产生、回收与处置与政府回收经费分配三个子模型，各子模型的方程体系如下： 1) 居家过期药品产生子模型： ① 总人口数 = INTEG (人口增长数，总人口数初始值) ② 人口增长数 = 总人口数*人口增长率 ③ 家庭总户数 =总人口数/平均家庭户规模 ④ 居家过期药品增长速率 = 家庭总户数*平均每户家庭每年产生的过期药品数 2) 居家过期药品回收及处置子模型 ① 居民家中过期药品数量 = INTEG (居家过期药品增长速率 − 正规药品回收速率 − 随意丢弃过期药品速率 − 非法药品回收速率，居民家中过期药品数量初始值) ② 正规药品回收速率 = 居民家中过期药品数量*正规药品回收率初值*回收点数量对正规回收的影响因子*回收点监管力度对正规回收的影响因子*政府宣传力度对正规回收的影响因子*经济因素对正规回收的影响因子 ③ 回收点数量对正规回收的影响因子 = WITH LOOKUP (正规回收点数量) ④ 政府宣传力度对正规回收的影响因子 = WITH LOOKUP (政府投入的总宣传经费) ⑤ 经济因素对正规回收的影响因子 = WITH LOOKUP (非法回收单价与正规回收补贴差额) ⑥ 正规药品回收点回收过期药品的数量 = INTEG (正规药品回收速率，正规药品回收点回收过期药品的数量初始值) ⑦ 正规回收点数量 = INTEG (正规药品回收点增速，正规回收点数量初始值) ⑧ 非法药品回收速率 = 居民家中过期药品数量*非法药品回收率初值*经济因素对非法药品回收的影响因子*政府宣传力度对非法药品回收的影响因子 ⑨ 经济因素对非法药品回收的影响因子 = WITH LOOKUP (非法回收单价与正规回收补贴差额) ⑩ 政府宣传力度对非法药品回收的影响因子 = WITH LOOKUP (政府投入的总宣传经费) ⑪ 非法药品回收商贩回收过期药品的数量 = INTEG (非法药品回收速率，非法药品回收商贩回收过期药品的数量初始值) ⑫ 随意丢弃过期药品速率 = DELAY1 (居家过期药品增长速率 − 正规药品回收速率 − 非法药品回收速率，延迟时间) ⑬ 随意丢弃过期药品的数量 = INTEG (随意丢弃过期药品速率，随意丢弃过期药品的数量初始值) 3) 政府回收经费分配子模型 ① 政府每年投入的回收点设置经费 = 正规药品回收点增速*平均每个药品回收点的设置费用 ② 政府每年投入的回收点监管经费 = 正规回收点数量*政府对每个回收点每年投入的监管经费 ③ 政府每年投入的药品销毁经费 = 单位药品销毁成本*正规药品回收速率 ④ 政府每年投入的回收补贴经费 = 政府给居民的单位药品回收补贴*正规药品回收速率 ⑤ 政府每年投入的宣传经费 = 政府每年投入居家过期药品回收工作的总经费 − 政府每年投入的药品销毁经费 − 政府每年投入的回收补贴经费 − 政府每年投入的回收点设置经费 − 政府每年投入的回收点监管经费 ⑥ 政府投入的总宣传经费 = INTEG (政府每年投入的宣传经费，政府投入的总宣传经费初始值)
system: [SYSTEM]

user: 4. 回收激励策略的仿真分析 浙江省瑞安市是我国较早开始实施过期药品回收工作的城市之一，于2008年成立居家过期药品回收服务工作小组，按照街道统筹、社区自愿相结合的原则进行过期药品回收网络布局，并于2011年基本实现了全市15个镇(街道)区域全覆盖。按照系统动力学的建模方法，将模型运用到瑞安市的实际药品回收系统中，通过政策仿真探究政府如何采用更好的激励策略来开展药品回收工作。 4.1. 参数估计及模型检验 对存量流量模型中的参数进行估计和确定是建模的重要部分。由于系统动力学模型行为的模式与结果主要取决于模型的结构，而不是参数值的大小，因此在进行参数估计时，只要保证参数具有恰当的准确性即可，无需追求参数的高精确度 [ 8 ]。本文选取的模拟时间边界为2008年~2028年，步长为1年，采用的参数确定方法主要有回归分析、趋势分析、表函数法等。数据的主要来源是瑞安市统计年鉴及瑞安市食品药品监管局对于居家过期药品回收工作阶段性成果的总结报告，同时也通过社会调查及文献研究的方法进行了相关资料的获取和整理。 模型的有效性直接影响了仿真结果的价值，因此在对模型进行仿真之前首先应该对模型进行检验，验证模型是否能够很好地描述现实系统的行为规律，保证仿真结果的可靠性。分别对模型进行了一致性检验，历史性检验以及极端条件检验后，验证了该模型能够有效如实地反映现实系统的行为规律。 4.2. 政策仿真分析 政策仿真是系统动力学的一个重要功能，即通过改变模型中某些关键变量，分析该变量所代表的政策变化对于系统行为的影响。本节将利用Vensim Dss软件进行政策仿真模拟，从“正规药品回收速率”这一变量的动态变化情况来考察不同激励措施的实施效果。 4.2.1 . 增设正规过期药品回收点 假设从2017年开始，瑞安市政府增设药品回收点，观察回收便利性提高后药品回收速率的变化情况。 图3. 增设回收点激励策略的仿真结果 由上图3可知，增设药品回收点的激励政策对于提升药品回收率的效果并不明显。这是由于瑞安市早在2011年，就已基本实现了全市15个镇(街道)区域的回收网络全覆盖，因而再增设回收点对于居民回收行为的激励效果并不明显。另外，从图中可观察到，该项激励策略的回收曲线在后期低于了策略实施前的初始情形，说明很多回收点将处于闲置状态，造成资源的浪费。 4.2.2 . 加大对正规药品回收点的监管力度 假设自2017年起，瑞安市政府增加对每个回收点年监管经费，观察政府加大对正规药品回收点的监管力度后，正规药品回收速率的变化情况。 图4. 加大对回收点监管力度的仿真结果 由上图4可知，加大回收点的监管力度对于居民参与药品回收的激励效果颇为显著，并且对着监管力度的逐渐加大，激励效果更加明显。这是因为目前瑞安市的药品回收点普遍存在着自身管理不善的问题，这在一定程度上打击了居民参与药品回收的积极性。在这种情况下，加大对回收点的监管力度，逐步提高回收点服务水平，从而有效地提升过期药品回收量。 4.2.3. 对居民实行回收补贴 此前瑞安市开展过期药品回收活动的主要形式是“无偿回收”，但是过期药品的回收奖励不宜采用直接的金钱补贴，政府可采用发放生活用品等的间接奖励。为了方便核算，采取0.5元、1元作为政府的单位药品补贴成本，观察政府对居民实行回收补贴后正规药品回收速率的变化情况。 图5. 对居民实行回收补贴激励策略的仿真结果 由上图5可知，当政府对于居民实行回收补贴的单位成本为0.5元时，对回收量的提升作用并不明显，而当单位成本增加到1元时，回收量将得到迅猛增长。居民作为理性的经济人，在参与药品回收活动时自然要考虑收益与成本。当补贴较少时，对公众没有很大的吸引力，而当单位药品的补贴提高，居民参与的积极性得以激发。从长期来看对居民实行回收补贴的策略成本过高，在不追加总回收经费的情形下，会大大削弱政府对于原有宣传工作的投入力度，从而使得总体激励效果处于不断减缓的趋势。从图上可以看出在策略实施的后期回收量甚至低于实施策略前。 4.2.4 . 加大过期药品回收宣传力度 由于实际工作中政府投入的总经费有限，上述激励策略的实施，都会在一定程度上制约原有宣传经费的支配和使用，因而这些激励措施的实施效果也综合了宣传力度变化对于过期药品回收量的影响。下面假定2017年起瑞安市政府在原有回收经费的基础上追加对于宣传工作的经费投入，单独考察加大过期药品回收宣传力度这一激励策略的实施效果(图6)。 图6. 加大回收宣传力度激励策略的仿真结果
system: [SYSTEM]

user: 5. 结论 本文对政府主导下的过期药品回收系统展开了较为全面的研究，针对现阶段居民对回收工作的积极性不高的现实问题，制定了相应的激励策略，并运用模型将激励策略带入到瑞安市的实际情景之中，观察激励策略的实施效果。通过改变参数输入进行仿真，观察不同的激励措施的实施效果，以此为依据对未来瑞安市政府制定居家过期药品的激励政策提供参考意见。 依据上述仿真结果，可对瑞安市政府制定居家过期药品回收激励策略提出相关建议如下： 当前情形下，不宜盲目增设过多回收点，以免造成回收点闲置、资源浪费的现象；应将关注焦点集中在对现有回收点的监管上，缓解由于其管理不善造成的居民回收积极性遭受打击的情况；在实行回收补贴方面，现阶段居民对于经济因素较为敏感，短期内补贴政策可有效促进过期药品回收量的增长，但从长期来看，经济手段的激励效果不断减弱，因而不能过度依赖。在回收宣传方面，提高居民环保健康意识是一个长期战略任务，随着时间的推移，加大回收宣传力度对于促进居家过期药品回收的激励效果会越来越明显。因此，若瑞安市政府注重居家过期药品回收工作长期发展后劲的增强，则不能因为短期内的回报率较低而忽视对于宣传工作的投入。
system: [SYSTEM]

user: 摘 要 焦石坝地区上奥陶统五峰组-下志留统龙马溪组下段发育一套优质页岩层，页岩品质横向分布较为稳定。通过近几年涪陵页岩气田的开发实践，发现水平井单井产量与埋深、微构造形态、岩层倾角以及岩层曲率等构造因素关系密切，但不同区块对页岩气单井产量起控制作用的主控因素存在差异。从构造因素出发对JD区块的单井产量差异进行分析，明确了该区块控制产量的主要构造因素，认为埋深以及岩层倾角是制约页岩气产能的关键因素，而微构造以及岩层曲率仅在局部区块对页岩气产能起到控制作用。 关键词 :单井产量，构造变形，埋深，地层倾角，微构造形态，焦石坝构造 Copyright © 2018 by authors， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 受深水陆棚相制约，大焦石坝地区晚奥陶世五峰期-早志留世龙马溪期发育一套厚度较为稳定的富碳、高硅质优质页岩，是页岩气“成烃控储”的基础 [ 1 ] [ 2 ] 。同时，成岩、成气以来的良好保存条件使得焦石坝主体箱状构造页岩气富集且高产，成为中国最大的页岩气田 [ 3 ] [ 4 ] [ 5 ] 。在焦石坝主体箱状构造获取页岩气商业开发后，如何实现页岩气的稳产及扩大勘探领域是亟需解决的问题。 JD区块位于焦石坝箱状背斜的北西翼北段，目前，钻井揭示该区具有较好的页岩气勘探潜力 [ 6 ] [ 7 ] [ 8 ] [ 9 ] 。构造变形是制约中国海相页岩气富集的核心因素之一，立足于构造变形的差异性，构建页岩气产能与构造变形的响应关系，查明构造对页岩气富集的关键控制作用，以有效推进大焦石坝地区二期页岩气的产建发展。
system: [SYSTEM]

user: 2. 区域地质背景 大焦石坝位于中上扬子过渡带、齐岳山断层中段西缘，为湘鄂西隔槽式褶皱向川东隔挡式褶皱转换处 [ 10 ] [ 11 ] [ 12 ] 。大焦石坝地区古生代为克拉通型沉积，中新生代经历了雪峰内陆拆离造山改造和内陆构造叠加。在江南雪峰陆内造山作用下，受湘鄂西-川东递进式扩展变形，焦石坝箱状构造发育 [ 10 ] 。垂向上，受寒武系、志留系、雷口坡组-嘉陵江组3套滑脱层影响，发生多层系滑脱作用，形成上、中、下3套构造形变层，并发生构造解体。 焦石坝箱状背斜形态完整，东西两侧受北东向和近南北向两组逆断层夹持围限而成，地表主要出露三叠系。JD区块位于焦石坝箱状背斜向川东构造过渡处，构造上为焦石坝背斜西翼转折端及西倾斜坡，东侧以吊水岩断层为界，西侧和北侧为珍溪场宽缓向斜，南部被近南北向的乌江断层改造，内部主要发育西倾的天台场1号断层和天台场2号断层。
system: [SYSTEM]

user: 3. 构造变形特征 受江南-雪峰造山带北西-南东方向递进构造变形作用，大焦石坝地区自东向西变形逐渐变弱(图1)。受滑脱层逐级升高，大焦石坝地区由前寒武系区域性滑脱层转变为寒武系滑脱层 [ 13 ] ，因此在焦石坝箱状背斜西缘形成西倾斜坡构造，并发生反冲作用，西倾天台场断层发育(图2)。 图1. 大焦石坝地区北西-南东方向测线地震解释剖面图 图2. JD区块奥陶系五峰组底界反射层构造与单井产量叠合图 JD区块构造变形总体较弱，可以划分为6个局部构造区(图2)：以北西方向倾没的天台场1号鼻状构造(④)为主，另外还包括焦石坝背斜西翼(①)、天台场2号断鼻(②)、吊水岩北斜坡(③)、沿江鞍部(⑤)以及乌江1号断背斜(⑥)。各分区局部构造特征差异明显，以断鼻、斜坡局部构造类型为主；主控断层天台场1、2号断层断面较缓，为层间断层，伴生断层不发育，具有活动强度较弱的特点，平面上断距南大北小，鼻状构造附近断距基本小于100 m。
system: [SYSTEM]

user: 4. 页岩气产能与构造变形响应关系 JD区块断层不发育，构造保存条件好，整体具有较好的产能。但受差异构造作用，不同区页岩气产能亦具一定差异。 4.1. 页岩气产能特征 目前，①~④号构造分区的高部位和翼部均有完钻试气井，水平井轨迹方向以平行构造走向为主，也有少量基本垂直于构造走向的。总体来看，①号分区内单井平均产量最高，②号分区次之，③号和⑥号分区内单井平均产量较低。同一构造分区内的单井产量差异也较大，如②号构造分区内单井最高测试产量是最低测试产量的4倍多，并且相邻单井的产量也出现较大差异(图2)。 4.2. 页岩气产能与构造变形响应关系 勘探显示，JD区块页岩气产能与埋深、微构造形态、岩层倾角以及岩层曲率密切相关。 1) 埋深。JD区块总体埋深变化较大，已试气井的埋深平均在3000~4000 m。通过对JD区块单井埋深与破裂压力、单井产量的交会分析表明，单井的埋深越大，破裂压力越大，单井产量越低(图3)。JD区块①号分区内构造埋深最小，单井的平均产量最高，④号分区内埋深较大，其单井的平均产量较低。YY1井位于②号分区鼻状构造的高部位，构造埋深为最小，穿行地层构造平缓，该井的测试产量为该分区的最高值。 图3. JD区块单井埋深与破裂压力(a)、单井产量(b)的交会图 2) 微构造形态。在埋深基本相似的情况下，构造形态的差异也会导致单井产量的差异。正向构造主要表现为张应力，地应力相对较小，负向构造以压应力为主，叠加上构造应力后，地应力往往较大，压裂施工难度会明显高于正向构造，造成单井产量偏低。YY4井与相邻水平井产量相比出现异常低值，该井水平段在相邻井中最长，①~③号分区的穿行率与其他井差异不大，分析原因认为该现象与井水平段穿行的局部构造形态的变化有关。该井水平段A靶点穿行于斜坡构造，地层具正向变形特征，到了B靶点附近穿行于两条逆冲断层夹持的洼槽构造部位，地层具负向变形特征(图2)，局部构造形态的变化代表着应力已经发生改变，说明井水平段从A靶点到B靶点应力的非均质性较强，影响了压裂改造，造成产量偏低。 3) 岩层倾角。统计分析表明，在地质条件、埋深基本相当的情况下，单井的产量与地层平均倾角为负相关关系，即地层倾角越小，垂直埋深差越小，产量越高。在鼻状构造北东方向低部位钻探的YY3井其埋深相对较大，穿行地层产状较陡，地层倾角超过25˚，水平段的方向基本垂直构造方向，靶点的垂直埋深差异大，超过了800 m，测试产量为该区块最低值，YY2井与YY3井构造部位相同，但垂直埋深较小，故产能较高(图4)。 图4. JD区块奥陶系五峰组底埋深与单井产量叠合平面图 4) 岩层曲率。构造大尺度裂缝的发育不利于页岩气保存，还会造成钻井的泥浆漏失及压裂改造效果差。实际中采用叠后岩层曲率属性值的大小及平面展布特征来表征构造裂缝的发育程度。统计分析表明，单井产量与岩层曲率呈负相关关系，但不是主要控制因素，仅在局部区域控制单井产能。
system: [SYSTEM]

user: 5. 结论 1) JD区块位于焦石坝箱状背斜西翼向川东构造过渡处，构造上表现为向西倾斜的单斜构造，仅受西倾的天台场断层局部改造，局部发生低幅挠曲，整体构造改造弱，为较为稳定的西倾斜坡。 2) 埋深与岩层倾角是制约JD区块单井产能的主要控制因素，埋深大、岩层倾角小，有利于页岩气的保存，因此产量高。微构造形态以及岩层曲率是该区单井产量差异的非主控因素，在埋深相当的情况下，宽缓的正向构造形态有利于压裂改造，该部位井的测试产量高，同时岩层曲率在局部区域控制单元产能。
system: [SYSTEM]

user: 摘 要 隔夹层构成的渗流屏障对地下油水运移特别是残余油的形成与散布具有较大的控制作用，隔夹层研究是提高油气采收率、最高效开发油气资源的关键因素。研究表明，西江24-1区隔夹层对剩余油分布控制作用强，其遮挡导致油田开发中后期有大量残余油分布，利用测井解释成果对影响低渗透层的主要因素进行分析，能够清晰地识别钙质、泥质和物性3类隔夹层。在一定沉积微相内，夹层具有稳定对比、随机分布的特征。 关键词 :隔夹层，测井响应，渗透率，分布模式，演变规律 Copyright © 2018 by author， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 油区地质特征研究 西江24-1区地层包括古近系，新近系文昌组、珠海组、珠江组、韩江组、粤海组、万山组，第四系，含油层段主要分布在新近系珠江组，厚778 m，共划分出23个油层，各油层之间均被较厚的泥岩分开 [ 1 ] [ 2 ] [ 3 ] 。研究区是一个位于主断层控制下的逆牵引背斜构造，形态简单，圈闭面积小，构造幅度低，走向北西-南东向，闭合面积0.87~4.32 km 2 ，闭合高度8~47.5 m [ 4 ] [ 5 ] 。区内发育退积充填低位体系域和海进体系域沉积，包括滨岸体系的滩砂、滩砂水道、沿岸坝、滨砂和辫状三角洲体系的分流河道、水下分流河道、河口坝、风暴席状砂等微相类型，可作为良好的储集层 [ 6 ] [ 7 ] 。
system: [SYSTEM]

user: 2. 隔夹层识别研究 西江24-1区的隔夹层实际上是各油层内部分布的夹层，尽管部分夹层厚度较大、分布面积较广，也构成隔层，但考虑其处在油层之内，所以统称为隔夹层。区内储层物性以中高孔、中高渗为特征，难以借用其他油田隔夹层分析方法和标准来识别该区隔夹层。 储层岩心分析结果表明，当渗透率小于120 mD、孔隙度小于16%时，含油率低于5%，该特征虽不能反映工区储层隔夹层划分标准，但表明在中高孔、中高渗储层中，相对较低的渗透率仍能影响油气的分布。以此为参考，将测井渗透率小于120 mD的层位暂作为低渗透层，将区内的隔夹层分为3类：钙质隔夹层、泥质隔夹层、物性隔夹层。其中，钙质隔夹层的渗透率小于10 mD，钙质质量分数大于10%；泥质隔夹层的渗透率小于10 mD，泥质质量分数大于30%；物性隔夹层的泥质质量分数为19%~31.3%之间，钙质质量分数小于10%。
system: [SYSTEM]

user: 3. 隔夹层分布模式 3.1. 单层砂体 1) 厚层分流河道砂体顶部常见钙质隔夹层，底部粗砂岩中少见，仅发育一些物性隔夹层，厚度一般为1~1.5 m。 2) 河口坝砂体主体部位隔夹层不发育(图1)，上部局部发育钙质隔夹层，但厚度小，分布少；侧翼发育钙质或物性隔夹层。 图1. 河口坝砂体上部发育的钙质隔夹层(1X井H12层) 3.2. 复合砂体 复合砂体常由垂向上多个期次的单砂体叠置而成，可划分为不同成因单元砂体复合和同成因单砂体复合，不同类型复合砂体隔夹层表现形式不同。研究区常见复合砂体有多期河道砂体复合、多期河口坝砂体复合、河口坝与河道砂体复合。多期河道叠加形成泥质隔夹层、物性隔夹层和钙质隔夹层，主河道附近隔夹层发育较薄，河道边缘隔夹层多且厚度大(图2)；河口坝砂体之间形成泥质隔夹层(图3)；河口坝与河道砂体之间形成物性隔夹层，厚度大(图4)。 图2. 多期河道叠加形成泥质隔夹层 图3. 河口坝砂体之间形成隔夹层 图4. 河口坝与河道砂体之间形成物性隔夹层
system: [SYSTEM]

user: 4. 隔夹层演化规律 平面展布上河道分叉区域砂体厚度相对较大，水动力条件和沉积构造组合类型复杂，泥质沉积量大，隔夹层发育。水下分流河道边缘部位的水动力条件弱，使得泥质成分增加，而河道中央水动力运动相对较强，构成泥质沉积缩减。研究区主要在水下分流河道分叉处以及河道边缘发育隔夹层，在单一水下分流河道及河道主流线部位的隔夹层发育较少或不发育。 河口坝沉积早期或河口坝边缘由于沉积能量较弱，往往发育被泥晶碳酸盐胶结的钙质夹层。多期河口坝砂体形成时，每期旋回早期水体较深，三角洲向湖盆推进，形成一套三角洲前缘砂体，同时期可容空间减小，水体变浅；次一级旋回早期水体加深，则可容空间陡然增大，河口位置倒退，在前一期三角洲前缘砂体之上沉积了前三角洲泥岩，形成泥质隔夹层，若被后期改造则易形成物性夹层。河道砂体与河口坝砂体叠加时，易形成物性隔夹层。
system: [SYSTEM]

user: 5. 结论 1) 西江24-1区储层可识别出钙质、泥质和物性3类隔夹层，并建立了分类识别标准。 2) 单期河道砂体顶部易发育钙质隔夹层，多期河道叠加时，主河道不发育隔夹层，河道边部和前缘易发育泥质、钙质隔夹层。单期河口坝主体无夹层或在顶部发育薄层钙质隔夹层，侧翼易形成钙质隔夹层；多期河口坝叠加易形成泥质或物性隔夹层。河道与河口坝之间易形成物性隔夹层。
system: [SYSTEM]

user: 摘 要 基于东濮凹陷桥口构造带油气地质特征，为获得油气地质储量，应用成藏动力学理论和技术方法，对区域内油气成藏动力学条件、油气富集规律进行了系统研究，分析了油气成藏模式和分布特点。 关键词 :成藏动力学，成藏模式，富集规律，桥口构造带 Copyright © 2018 by authors， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 区域油气地质概况 桥口构造带位于东濮凹陷黄河南探区中央隆起带北端，是受兰聊和黄河断层控制的继承性背斜构造，包括桥口和新霍两个构造，东西宽约10 km，南北长20 km，面积约200 km 2 。东临葛岗集洼陷，北东倾没于前梨园洼陷，西与海通集洼陷相邻，是多个洼陷中间的隆起带，形成多源供油、供气的态势，为沙河街组一段，二段下亚段，三段上、中、下亚段均含油气的油气富集区 [ 1 ] 。
system: [SYSTEM]

user: 2. 油气成藏动力学特征 2.1. 地温场特征分析 根据实测地温资料 [ 1 ] 可以发现，桥口地区的地温梯度在3.0℃/100 m以上，地表温度约为14.6℃。桥口构造的地温梯度稍高，约3.42℃/100 m (图1)，新霍构造的地温梯度比桥口构造稍低，约为3.08℃/100 m。同时，由于局部断层活动和局部热的封隔体的存在，地温梯度呈现复杂的分布格局，从隆起 → 洼陷呈现出逐渐降低的趋势。 2.2. 压力场特征分析 桥口地区的压力场垂向上分为3套：沙河街组一段~三段上亚段，为常压系统，压力系数多数分布在0.96~1.06之间；沙三段中亚段为过渡带，压力系数为1.06~1.3左右；沙三段下亚段为高压异常带，压力系数普遍大于1.3，局部压力系数达到了1.5。平面上，由于断层活动的复杂性与多期性，形成有规律的压力平面分布特征：东部，沙三段下亚段的异常高流体压力场中形成高压异常带；中部，沙三段中亚段压力过渡带系统中的流体压力介于高压异常带和正常流体压力系统之间；西部，常压流体系统中形形成正常压力系统(图2)。 图1. 桥口构造温度–深度关系图 图2. 桥口地区压力系统垂向变化 2.3. 流体势对油气分布的影响 沙二段下亚段、沙三段上亚段流体来自于前梨园和葛岗集次洼，油气从洼陷向构造顶部运移，油源充足 [ 2 ] [ 3 ] 。在靠近洼陷的地方，原油密度自东向西变小，进入断层发育带，原油密度自东向西变大；沙三段中亚段主流体来自于葛岗集次洼，由于李屯断裂的控制，大量油气不能穿过断层向构造高部位运移。因此，油气主要来自于次流体方向，即从孟岗集次洼向构造高部位运移，由于次流体的势能小，流体运移速度缓慢，因而沙三段中亚段含油面积和丰度均不及沙二段下亚段和沙三段上亚段。沙三段下亚段的油气大部分来自于东北方向的前梨园次洼(图3)。 图3. 桥口地区流体分布图 2.4. 油气运移输导体系分析 通过对桥口地区流体输导网络的分析认为，区内断层、砂体及层序界面构成三维输导体系网络，存在沿黄河断层展布的断层–岩性输导脊、沿李屯断层和玉皇庙断层展布的断层-岩性输导脊，对油气由洼陷向构造高部位的运移提供了良好的通道。在以岩性与层序界面输导体为主的地区，油气以面状运移为主，在构造-岩性输导脊部位，油气以脊状运移为主，油气在桥口构造高部位聚集，并在构造东部形成天然气气藏(图4)。 图4. 桥口地区输导体系平面配置示意图
system: [SYSTEM]

user: 3. 成藏模式及油气分布特点 3.1. 成藏模式 桥口构造为洼中低隆起，聚集了来自周围多个洼陷地油气资源。桥口构造带复杂的断裂系统、砂体及层序界面为油气的运移提供了良好的输导通道，油气沿输导通道(断层、砂体、层序界面等)由洼陷向构造位置较高的桥口构造带运移，遇到有利圈闭则形成油气藏，或继续向上倾方向运移，纵向上沿断层向浅层运移，直至封堵条件好的圈闭中形成油气藏。 东营运动是整个东濮凹陷的一次主要的油气运移时间，油气主要是沿断层向上部层系运移 [ 4 ] 。东营运动后是油气二次生成和运移期，生成的油气主要沿生油岩在层系内的低部位向构造高部位运移，在环绕生油洼陷较近的圈闭中形成背斜油气藏和断鼻油气藏。 3.2. 油气藏特点及油气分布规律 该区的玉皇庙、桥东、新东断层大多数时间具有开启性，为东部洼陷生成的油气向构造顶部、翼部运移的通道 [ 4 ] [ 5 ] 。黄河二级、李屯断层系列西倾配套断层具有封堵性，能使大量油气在断鼻圈闭内聚集起来，形成有工业价值油气藏。受断层多、断块小、含油气层系多等因素影响，桥口构造形成复杂断块群油气藏，其特征表现为：垂向上含油层系多，纵向上油层分布受断块控制，主要为沿断层上升盘呈“刷子型”分布的成藏模式，即具有逢油必断，高点富集，油层分布范围小，垂向上多套油水系统，平面上叠加连片的特点。 从圈闭评价情况来看，二级断层控制着油气的分布，不同断块含油层系不一样，富集程度也不一样。油气差异聚集明显，构造高部位油气聚集最好，为油藏，腰部次之，为油气藏，低部位差，为气藏。
system: [SYSTEM]

user: 4. 结论 1) 桥口构造区油气资源丰富，运移通道良好，在构造高部位及浅层圈闭有利区易形成油气藏。 2) 区内的油气分布具有逢油必断，高点富集，油层分布范围小，垂向上多套油水系统，平面上叠加连片的特点。依据构造部位的由高至低，依次发育油藏、油气藏、气藏。
system: [SYSTEM]

user: 充填体压缩性能的试验性研究 肖雄飞，孟浩，赵明，苗健 中国煤炭地质总局勘查研究院采空所，北京 作者简介：肖雄飞(1988-)，男，工程师，主要从事采空塌陷地质灾害防治。 收稿日期：2018年2月7日；录用时间：2018年5月14日；发布日期：2018年12月15日
system: [SYSTEM]

user: 摘 要 为保障充填开采对抑制地表沉陷影响的高效性，以冀中能源邢台煤矿为工程背景，对充填体变形多个影响因素展开分析，就充填矸石粒径大小的影响进行了试验性研究。研究表明，影响地表沉陷的主要因素之一是充填体压缩率；矸石粒径大小与充填体压缩率有一定的相关性，可通过优化选择矸石粒径来调节矸石的压缩率，以此来控制地表沉陷。 关键词 :地表沉陷，矸石粒径，充填体压缩率，充填开采 Copyright © 2018 by authors， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 综采高产高效充填工作面充填体的压缩性能是矿井控制地表沉陷和预防地表移动变形的一个重要性质 [ 1 ] [ 2 ] [ 3 ] [ 4 ] [ 5 ] 。它不仅受到充填材料本身的影响，而且受到充填配料方式和工作面作业方式的影响。因此，井下综采中的高产、高效矸石充填对地表沉陷的控制应该从充填体压缩性能上入手，减小压缩率，加强对顶板下沉的控制。一般来说，矸石综采充填方法可以减少环境污染和最大潜能地发挥工作面综采设备的能力，对提高工作面的产量和效率意义重大。但是充填后较大的矸石压缩率将会导致工作面充填体对覆岩控制程度不够，给工作面的顶板管理以及地表沉陷的控制带来较大的挑战，不利于工作面的稳产，经济效益的提高和建下采煤、采动的要求。因此，应根据地质、经济、技术、矿压、管理等因素综合考虑，通过完善充填体颗粒大小、提高充填体强度等一系列措施减少充填体压缩率。充填体压缩率的变化影响着上覆岩层的移动变化，充填体压缩率越大，覆岩的移动范围就越大，对地表沉陷与移动的控制效果越差 [ 5 ] - [ 15 ] 。 以冀中能源集团金牛能源股份有限公司邢台煤矿为工程背景，通过实验室试验研究矸石粉煤灰混合充填料中矸石粒径的大小对充填体压缩率的影响，以便在将来对煤矿充填采煤地表沉陷控制进行理论指导。
system: [SYSTEM]

user: 2. 工程概况 充填采煤即通过综合机械化设备将充填材料置于采空区内，抑制上覆岩层的移动、变形，进而有效地控制地表移动变形 [ 16 ] 。充填采煤控制地表沉陷的因素有充填前顶底板移近量、充填体欠接顶量、充填体压缩率、底板浮煤压缩率等。研究表明，前3个因素是影响地表沉陷的主要因素。笔者主要通过研究充填体的压缩率入手，强调充填体的压缩率大小对地表的影响，加强对地表沉陷与移动变形的控制。其中，充填体压缩率对地表沉陷影响示意图如图1所示。 图1. 充填体压缩率对地表沉陷影响示意图
system: [SYSTEM]

user: 3. 充填体压缩影响性分析 结合邢台煤矿工作面生产情况，在有利于顶板稳定和控制地表下沉量的前提下合理提出影响充填体压缩率的因素。充填采煤法类似于以岩性较软的充填体和岩性较硬的煤层进行替换。在顶板覆岩的压力下充填体就会发生相对于顶底部的压缩变形，对地表沉陷和变形造成影响 [ 17 ] 。其中，充填材料的含水量或泌水性、充填体颗粒大小以及充填体的强度会影响充填体压缩率的变化。 3.1. 充填材料粒径大小对充填体压缩变形的影响 对于矸石充当单一充填材料，国内已研究过不同矸石粒径应力、应变的关系。通过试验可知，对矸石充填体施加同样的轴向载荷，矸石的粒径越大，也将产生较大的压缩变形，即粒径越大，充填体压缩率越高。产生该现象的主要原因是大粒径的矸石材料，其颗粒间具有较大的孔隙和较高的孔隙度，在其受轴向力产生压实过程中，这些孔隙将被压缩和填塞，以致产生较大的应变。 3.2. 充填材料含水量对充填体压缩变形影响 通常情况下矸石充填材料含水量会对充填体的压缩变形造成一定程度的影响 [ 18 ] 。但一般的矸石都是就近取材，改变其含水量的难度较大，故对其含水量不作重点研究。通过对7608工作面的现场测试，充填材料的平均含水量为9.18%。 3.3. 压实方式对充填体压缩变形的影响 充填工作面倾角越大，充填的矸石在自重力作用下越容易密实，其充填效果亦较好 [ 19 ] 。提高矸石的初始压实率能有效减小等价采高，更好地控制地表沉陷与变形。具体表现为捣实机加强对矸石的作用，使其尽量充满充填区域。 3.4. 充填采煤提高充填体压缩性能的措施 充填体压缩率还受充填体强度和围岩压力的影响，充填体强度越大，压缩率越小 [ 6 ] 。矸石粉煤灰充填与传统单一矸石充填不同，矸石作为充填骨料，粉煤灰则用来当作辅料减少充填料中的空隙，以便提升充填体压缩性能。提高充填体的压缩性能可以通过减小压缩率的方式来实现，例如可以调节充填体材料颗粒大小和压实方式，还可以进行注浆补强，同时也要控制好充填的成本，合理进行井下充填。采煤过程中，需要不断改善工作条件，采用更多的完善方案，在满足安全要求的前提下最大限度地提高经济效益，选择相应合适的措施。
system: [SYSTEM]

user: 4. 矸石粒径大小对充填体压缩率影响的试验 通常情况下充填体强度的研究比较多一些，而关于充填体粒径的研究比较少，通过大量试验研究，得出了不同矸石粒径大小充填体对地表沉陷的影响规律。鉴于该矿是矸石和粉煤灰混合充填，与单一矸石充填有所区别，该试验研究不同矸石粒径大小对混合充填体压缩变形的影响。由7606工作面现场实际生产可知充填开采时，采用的是矸石和粉煤灰固体混合物密实充填，通过实际过程中科学的配比研究，选定矸石与粉煤灰的质量配比为10:3。可统一试验矸石和粉煤灰的质量配比，对充填材料进行单轴抗压试验，对不同矸石粒径充填体压缩率进行分析，研究相同载荷作用下充填体的压缩变形规律。 4.1. 矸石粉煤灰压缩试验设计 试验采用自制圆柱形钢筒，将制备好的矸石粉煤灰混合体放入圆柱形钢筒内，在试件上放上传递载荷的钢柱，采用压力试验机对其进行加载。加载不同粒径的矸石粉煤灰充填体，记录施加的载荷及其相应的压缩率，计算各组试件在试验中的压缩率。 用实验室粒径筛选机筛选矸石，得到0~0.6、0.6~2.5、2.5~5、5~10、10~15、15~20 mm共6组单粒径矸石。对不同粒径的矸石进行轴向压缩试验，研究各组试样的变形特征，轴向压力加载速率为1 N/s，采集的数据为轴向应力和轴向压缩变形量。试验中采用的自制模具为圆柱形钢筒，内径为8.0 cm，外径为9.0 cm，壁厚为0.5 cm，筒高16.0 cm。统计研究结果表明，试验所用钢筒内径(D)与试样最大粒径(d max )之间合理关系为：d max ≤ D/5。试验所用钢筒D为80 mm，因此矸石的最大单粒径区间为15~20 mm，基本满足要求。试验时，圆柱形钢筒侧压大致处于相同条件，可采用控制变量法，只研究轴向应力对充填体压缩率的影响。设计单轴压力分别为10、25、50、100 kN，根据圆筒内充填体受压面积，计算可得对应的轴向应力依次为2、5、10、20 MPa。试验所用筛选矸石粒径的器材、矸石粉煤灰放置在圆柱形钢筒后进行加载的示意图如图2所示。 图2. 试件加载设备示意图 4.2. 矸石粒径与充填体压缩变形相关性分析 试验的研究目的是掌握矸石粉煤灰充填体压缩特性，找出充填体矸石粒径与压缩变形的对应关系，可在实际可操作范围内选出优化的粒径方案，进一步减小充填体的压缩率，进而更好地控制地表沉陷。充填体压缩率可通过压缩率和装样高度的比值进行计算。其中，压缩率的计算公式为： V = ( H 0 − H 1 ) / H 0 × 100% 式中：V为压缩率，%；H 0 为充填体试件未压缩前的高度，mm；H 1 为充填体受轴向压力作用稳定后的高度，mm。 将试验前后充填体进行对比可知，未压缩前的充填料比较松散，颗粒较大，压缩后矸石发生破碎，颗粒变小，并且矸石和粉煤灰密实在一起，充填体随着传载钢柱往下移动了一段距离，挤压破碎后的矸石体积变小将空隙进一步减小，使充填体进一步压实，如图3所示。 图3. 钢筒内矸石粉煤灰压缩前(a)、后(b)示意图 由图4可以看出，轴向应力为20 MPa时，最大的压缩率为33%左右，最小的压缩率也有23%左右，数值有较大差异。不同矸石粒径大小的充填体充填效果不同，具体表现为矸石粒径越小充填效果越好，但并不是简单的线性关系。当粒径小于2.5 mm时，充填体压缩率的变化非常缓慢，可能由于材料之间距离较近，分子之间产生斥力，导致压缩变得困难。当不同粒径矸石充填体承受相同轴压时，压缩率随着矸石粒径的减小而减小，矸石粒径最小的充填体压缩率最小，最大的充填体压缩变形也最明显。 图4. 轴向应力20 MPa时充填体压缩率 对不同矸石粒径大小的混合料充填体进行相同的轴向加载试验可得其压缩曲线如图5所示。 图5. 不同矸石粒径充填体受力压缩变化图 对比曲线可以看出，充填体中矸石粒径大小和充填体压缩变形量呈正相关关系。施加同样的轴向载荷，矸石粒径较大的充填体，对应的压缩率也比较大，矸石粒径较小的充填体，压缩变形也会较小一些。矸石粉煤灰混合充填体在刚开始施加轴向应力载荷时，压缩率随着载荷增加而迅速增加，轴向加载从0到5 MPa之间曲线切线斜率最大，压缩率增长最快。轴向加载进行了一段时间后，充填体中的矸石粉煤灰进行了初步的压实，颗粒之间相互运动困难，活动的空间变小，主要在10 MPa以后压缩变形增速变得特别小。施加轴向载荷时间继续增加，应力也不断变大，充填体的压缩率逐渐趋向稳定，压缩率保持一个稳定值。
system: [SYSTEM]

user: 5. 结语 通过对充填料矸石粒径大小的研究表明，矸石粒径大小对充填体压缩率有一定的规律可循，可通过优化选择矸石粒径调节矸石的压缩率因素来控制地表沉陷。充填体受压缩之后压缩率越大，则控制地表及覆岩移动的效果越差，反之则控制效果越好。充填体中矸石粒径越大，对应充填体压缩率越大，抵抗压缩变形的能力越弱，要对输送、加工以及施工成本等进行综合考虑，采用合理粒径范围的矸石进行充填。当矸石的压缩率较小时，充填完成后的充填区域顶底板移近量较小，其控制地表沉陷的效果也越好，充填效果显著。
system: [SYSTEM]

user: 摘 要 研究区块位于塔里木盆地塔中隆起带桑塔木地区，地表条件复杂，地貌变化频繁，地形总体平坦，但局部差异较大，有效砂岩储层厚度较薄，平面分布差异大，非均质性强，导致静校正问题突出。层析静校正技术和分层地表一致性静校正技术能够改善速度谱质量，提高弱反射信号强度和叠加剖面同向轴连续性。该技术在研究区的应用结果表明，成果数据体沿层切片展现了沉积体的整体连续性，且井震结合标定符合度较高，改善了地震弱反射叠加剖面的品质，提高了成像质量，为该区的地震地质解释提供了正确依据。 关键词 :地震弱反射，静校正，层析静校正，地表一致性剩余静校正 Copyright © 2018 by authors， Yangtze University and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/