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

user: 4. 仿真实验结果及分析 仿真实验通过人为控制迭代次数的神经网络模型RMA算法程序的运行，实现了对木材干燥实验数据相对精确、理想的系统状态数学解析式，其具体表达式如下： [ d x 1 d x 2 ] = − A × [ x 1 x 2 ] + [ w 11 w 22 w 33 w 44 ] × [ tanh ( t 1 ) tanh ( t 2 ) ] 其中 d x 1 ， d x 2 分别表示温度和湿度关于时间的微分， − A 是形如 [ α 0 0 α ] 的对角阵，其中 α 为小于1的常数， ( w 11 w 22 w 33 w 44 ) 是RMA神经网络模型的连接权矩阵，其中元素具体表达式如下： w 11 = − 0.0008365 x 2 + 0.0135 x + 0.4354 w 22 = − 0.0001135 x 2 + 0.01734 x + 0.7403 w 33 = 0.0001015 x 2 − 0.0003172 x − 0.01734 w 44 = − 0.0013 x 2 + 0.0197 x + 0.8454 根据以上数学形式的描述，可以清晰地看到木材干燥试验中含水率影响因素——温度和湿度关于时间的系统解析表达式，将木材干燥过程中温度、湿度与时间的变化关系组成了类似动力系统常微分方程的形式，根据方程表达式的数学性质可以进行系统解等稳定性分析，由实验结果进行更理论化的木材含水率分析处理。并且由图3和图4已经充分证明RMA神经网络模型算法对给定的离散数据逼近性质良好，算法收敛性强，迭代次数具有可控性。仿真模拟的方法还具备因地制宜的可变性，即根据实际数据情况应用于实验室中系统函数模拟和分析，具备可检验性和推广性，对进一步研究木材干燥含水率系统函数具有十分重要意义。
system: [SYSTEM]

user: 5. 结语 本文已经证明了关于木材干燥过程中的含水率问题可以通过建立神经网络模型的方法实现数据的误差处理以及利用RMA算法计算出数据的逼近函数，在实际应用中，此算法可以帮助实验室工作人员精确且有效的处理好实验数据，极大提高木材干燥的效率。接下来，继续研究木材干燥过程中的除温度、湿度外其他影响因素并加入多层反馈神经网络模型中进行研究也是十分有意义的。
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 ] [ 3 ] 。然而，近年来随着全球气候的变暖，火灾与气候的关系问题日益突出，火灾发生多是对极端气候 (高温、干旱)的即时响应。因此，长期以来气候–火灾之间的相互作用有助于我们探寻火灾发生的气候背景及基本规律 [ 4 ] [ 5 ] [ 6 ] 。 全球气候变化已经对社会经济发展和自然生态系统造成了重大的影响，是人类面临的重大挑战之一。目前全球变化引起的全球变暖、持续干旱、夏季少雨等现象导致全球各地火灾频发，高频率、高强度的火灾导致森林退化及片段化，栖息地破坏，生物多样性下降 [ 7 ] 。近年来的厄尔尼诺现象也被认为是导致一些地区火灾发生的重要因素之一 [ 8 ] 。随着气候的周期性变化，火灾的多发年和少发年也有准周期的波动。极端高温或干旱是天然火发生及蔓延的主要气候背景，根据气候模拟，在CO 2 增加一倍的情况下，未来100年增温幅度能够达到1.4˚C~5.8˚C [ 9 ] ，因此，天然火发生与全球气候变化的关系引起越来越多研究者的关注，目前却还没有具体的针对气候变化下的火灾预防措施和应对措施。末次冰消期以来的全新世是距离人类最近的一个时期，过去气候变化诱导下火灾发生的频率与强度的研究，对我们全面了解亚轨道千年尺度上火灾的时空序列演变与分布及过去火事件如何响应气候变化等科学问题提供了重要依据，为以后火灾发生的预测提供了基本信息，对于高温或干旱期间森林火灾的预防和扑救将有重要意义。本文主要回顾古火灾研究的发展史，介绍了目前古火灾研究的替代性指标，总结末次冰消期以来火灾–气候相关性研究的突破性进展，分析全球变化背景下火灾活动并对未来火灾研究的发展进行预测。
system: [SYSTEM]

user: 2. 古火灾研究简史 古火灾研究最早始于20世纪中期。1941年，欧洲学者Iversen首次利用显微镜对陆地、湖泊沉积中的微观炭屑进行了观察分析。尽管该研究的主要目的是针对孢粉探讨人类活动对植被覆盖的影响，然而，炭屑的发现，仍然是一项开创性的工作，提出了一个新的研究领域，跨越时间空间尺度，为地质历史时期古火灾数据收集提供了可能，为古火灾研究提供了有力的证据 [ 10 ] 。20世纪60年代后，光学显微镜的广泛应用，使炭屑的鉴定分类更加明确 [ 11 ] ，为古火灾重建奠定基础。1973年Smith等第一次在海洋沉积物中发现黑碳并对其进行了测量，为古火灾研究提供了又一新的指标，该成果引起科研工作者的极大兴趣并开始致力于扩展和完善黑碳对于火灾的指示作用，试着用黑碳重建古火灾历史并反演过去的气候环境变化，古火灾研究开始萌芽 [ 12 ] 。20世纪80年代，首次“国际火历史”会议在美国召开 [ 13 ] ，标志着古火灾研究进入了新时期。1996年，“加拿大全国森林火灾研讨会”指出，保存在各个地质历史时期地层中的木炭是古火灾历史数据获取的主要途径，木炭作为古火灾代用指标得到认可。1997年，在英国召开的“古火灾环境学研讨会”对于木炭的埋葬学、古火灾和现代火灾进行了深入交流和探讨，进一步明确了古火灾指标的价值和古火灾重建的重要意义。Pyne等(1996)和Scott (2000)根据宏观木炭所提供的地质历史时期火灾数据，发现最古老的炭化植物化石出现于晚志留世–早泥盆世 [ 8 ] [ 14 ] [ 15 ] ，由此推断地球上最早的火灾发生于距今约 3.5~4.0亿年。而二叠纪木炭已广泛分布于美洲、欧洲、非洲、亚洲及澳大利亚，由此推断，二叠纪火灾已经广布全球，成为陆地生态系统的普遍现象。我国古火灾研究起步晚，主要借鉴国外已有研究成果和经验发展起来，目前仍处于探索阶段。
system: [SYSTEM]

user: 3. 火灾的替代性指标 火灾发生后通常会在周围环境中留下燃烧痕迹及产物，因此，可以通过岩石圈(磁学参数、地球化学)、生物圈(孢粉学、古生态学、炭屑及树木年代学)甚至人类圈(考古学、历史图鉴)的历史档案中获取信息重建古火灾事件 [ 14 ] 。从天然档案的古火灾替代指标中获取的信息通常会出现不连续、难以解释的情况，并且会受到沉积环境等多种因素的影响，导致古火灾历史重建与现实或多或少会有一些差异。因此，获取非常准确的古火灾序列对古火灾研究工作者来说是一个挑战 [ 16 ] 。目前对于末次冰消期火灾事件重建的研究多是以泥炭、黄土、湖泊、海洋沉积物为研究对象，以燃烧产物(炭屑、黑碳)为古火灾代用指标，结合沉积环境的氧化还原特征及微体植物化石，探讨火灾–气候–植被之间的耦合性。当然，不同指标所反映的沉积环境信息、时间跨度、火灾的时空分布及分辨率等也存在差异。 3.1. 树木年轮火灾学 树木年轮火灾学重建古火灾事件主要采用树木年轮火疤记录和历史档案资料进行分析。不同方向和高度的火疤，记录了不同时期的林火。树木年轮火疤记录可以追溯近代森林火灾历史，因而被广泛应用 [ 17 ] 。利用树轮开展林火研究不但可以记录火灾发生的不同年代，如果获取样本量的区域尺度足够大，还可以分析火灾的空间分布特征。虽然树轮年代学准确描述了过去的火灾历史，但因为其记录年代有限，对较长时间尺度上的深入分析有局限性 [ 18 ] 。 3.2. 沉积中的炭屑分析 炭屑是植物不完全燃烧产生的无机碳化合物，能在湖泊、泥炭、海洋等沉积物中长期保存下来，是反映古火灾事件的直接证据。炭屑的产量与分布取决于地质历史时期火灾发生的频率、强度、范围和变化。通过炭屑的定量统计和形态分析，可以重建连续几千甚至几万年天然火事件的历史记录，因此，炭屑作为古火灾替代指标在沉积物中被广泛应用。相对于树木火疤，炭屑不仅应用时间尺度长，能反映火灾事件的时间跨度 [ 19 ] ，且不同粒径炭屑的分布可以反映火灾源的距离，粒径较大(>100~200 μm)的炭屑通常不会迁移很远，因此可以反映区域火灾事件 [ 20 ] [ 21 ] [ 22 ] 。 3.3. 燃烧后的分子标志物 化学标志物是由不同分子组成的火灾事件的化学示踪剂，是新兴的古火灾代理指标 [ 23 ] ，这些化合物在燃烧时产生、挥发，并快速凝结在炭屑、灰烬和其他较重的颗粒甚至类似于焦油等的半流体物质当中。分子标志物也可能来源于火灾产物的分解、成岩作用及其他理化性质转换等，这种情况下通常被认为是火灾的次生产物。近几年，利用分子标志物重建古火灾事件才广泛开展起来，但仍存在一定的局限性，例如一些化学标志物主要存在于工业革命之后的沉积物中，工业革命之前的沉积物中化学标志物含量很低，难以检测 [ 24 ] ；对于不同沉积成岩环境中分子标志物的稳定性、持久性及降解机制的研究仍不成熟。因此，利用分子标志物重建古火灾事件多集中于表层沉积而很少应用于长时间尺度的研究。 3.4. 磁性矿物 火灾的发生会直接或者间接的影响周围环境的磁学参数。火烧过程中，土壤中的顺磁性铁矿物或者是燃料中的痕量铁被转化成二次强铁磁性氧化物，而这些强铁磁性氧化物是磁铁矿和磁赤铁矿的重要组分 [ 25 ] [ 26 ] 。火烧之后，集水区域磁性风化层的侵蚀可能导致亚铁磁性粒子的流入进而在湖底沉积物中形成磁层。当然，其他因素也可能导致局部区域土壤中磁性矿物含量增加，如黄土–古土壤中的磁学信号反映了富含铁矿物的土壤暴露于高温多雨条件下的成壤及风化过程，而火灾也可能引发相同的磁学信号 [ 27 ] [ 28 ] 。随处可见的趋磁细菌也可能引起土壤及沉积环境中磁学参数的改变 [ 28 ] [ 29 ] ，空气中的磁污染以及除火灾之外其他因素引起侵蚀作用的增强(尤其是基岩中含有大量的磁性矿物)都可能导致沉积物中磁学信号的峰值 [ 29 ] 。沉积物中磁学参数是重建火灾历史和古生态的一个有前景的替代指标 [ 30 ] [ 31 ] 。此外，磁学参数的深入研究有可能反映过去火灾事件发生的强度、类型及频率，为区分天然火和人为火提供新的可能性。不足之处是磁学参数容易受到外界干扰因素的影响，因此在火灾历史重建研究中应用还不是很多。 3.5. 多指标结合重建长期古火灾序列 随着研究的不断深入，越来越多的学者试图通过多指标结合理解长时间尺度上气候–火灾–植被的相互作用及火灾过后生态系统的响应模式。炭屑和孢粉结合定量分析是目前为止探讨气候，植被、火灾及人类活动之间的相互关系的重要手段之一 [ 21 ] 。然而，孢粉易受到风、水流等外力的搬运无法完全反映原地植被丰度。最近研究表明，火灾时期，沉积物中含有丰富的燃烧植硅体，燃烧植硅体峰值对应炭屑的峰值 [ 7 ] ，而植硅体受外力影响较少，可以反映原地植被面貌，因此，植硅体在恢复古气候、古植被重建古火灾方面有一定优势。中国开展恢复古火灾–古气候的研究较晚，目前运用古气候记录(炭屑、孢粉等)对火与气候变化及植被类型之间关系的综合研究还相对较少。对火灾–气候耦合性进行长期、深入、全面和系统的研究是国内古环境研究的重要发展方向之一。
system: [SYSTEM]

user: 5. 火灾及全球变化意义 全球正在经历以变暖为特征的气候变化过程。火作为全球生态系统重要的干扰因子，气候变化直接或间接影响到燃烧环境，一般情况下，气温高、降水少易引发火灾。气候变暖使很多区域呈现暖干化趋势，导致火灾频率、火灾强度和过火面积增加，还延长整个火险期；可燃物达到着火点需要一定的热量，这种热量主要来源于周围环境的温度，因此，高温还导致可燃物更加易燃 [ 75 ] 。气候变暖引发的极端气候事件(强降水、飓风、干旱和冰冻灾害等)会导致大量植被受损和死亡，地表易燃可燃物猛增，火灾发生的危险性增大。如2008年初发生在我国长江流域至江南地区百年一遇的低温雨雪冰冻灾害，导致林木大批折断，地表可燃物猛增2~10倍，平均地表可燃物载量超过50 t/hm 2 ，部分严重地段达到100 t/hm 2 以上，已超过可发生高强度林火和大火的标准(30 t/hm 2 ) [ 76 ] 。预期气候变暖情景下，未来全球大部分区域的林火天气状况会更加严峻，林火数量将比当前状况更多，林火强度更大。林火的发生亦与干旱严重程度和持续时间有关。海洋沉积物黑碳记录显示，冰期黑碳浓度远远高于间冰期，原因可能是冰期气候干冷，火灾的发生频次比较高 [ 77 ] [ 78 ] [ 79 ] 。大兴安岭林区春季特大森林火灾与火灾发生前气象因子相互关系的研究中发现，43.9%的特大森林火是由前期湿或者特湿转为后期干或者特干导致发生的，长期持续干旱容易导致特大火灾的发生 [ 80 ] 。 美国部分地区研究表明，在十年尺度上，火灾的发生受到ENSO的影响，多数火灾年倾向于在El Niño事件向La Niña事件转换时发生，火灾减少年与ENSO减弱年较同步 [ 81 ] [ 82 ] 。据统计，世界历史上78%的森林大火发生在El Niño年 [ 82 ] 。我国气候条件相对复杂，不同类型不同发展阶段的厄尔尼诺现象对不同区域的影响不同，因此，ENSO事件对我国气候的影响需要进一步研究。一般来说，厄尔尼诺引起的高温和干旱天气会导致东北地区的森林火灾变得严重。由于华北地区没有大面积片林分布，而且交通方便，人口密度大，虽然火灾天气条件严峻，但不易发生火灾 [ 83 ] 。 火不仅即时响应气候的变化，同时对气候变化具有反馈作用，火灾的发生可能导致森林碳循环的改变，使森林成为大气中CO 2 和CH 4 的源，造成大气中CO 2 和CH 4 浓度升高，除此之外，燃烧产生的固体颗粒物也会引起空气污染，植被燃烧产生的卤代烷会破坏臭氧层，从而进一步加强全球变暖趋势，形成恶性循环，干扰全球气候变化 [ 84 ] 。
system: [SYSTEM]

user: 6. 结论与研究展望 全球古火灾序列综合对比研究表明，在全新世早期，火灾发生的频率低，周期长，而在全新世中期，火灾频率增加，这与全新世早期至中期，全球升温具有一致性。然而许多研究区域在早全新世火灾频率很高，可能受局部环境及其他因素影响。晚全新世全球火灾频率普遍增加，显然受人类活动影响。目前，随着更多的火灾代用指标的发现与应用，关于火灾历史重建已经取得了突破性进展，并得到了一些有意义的结论，预计未来可能的发展主要体现在以下几个方面： 1) 综合不同地区，不同指标发现，不同时间尺度及空间尺度的古火灾序列研究得到的结论也不大一致。为更好的理解火灾历史，保证重建古火灾序列的可靠性和完整性，在半球甚至全球尺度上不同地质记录进行对比，多指标相结合弥补单一记录存在的缺陷，高分辨率地恢复过去的火灾事件，可能会成为火灾历史重建未来发展的一个重要方向。 2) 火灾代用指标分析方法、分类鉴定、埋藏学及现代过程研究可能是今后研究中一个挑战，代用指标的沉积、迁移和保存是古火灾研究的基础和依据，目前关于火灾机制的研究更多的是结合同一地质历史时期的环境状况和人类活动推断火灾发生的原因，因此，熟悉掌握代用指标的特征及属种鉴定，对于区分天然火和人工火，在火灾发生机制、强度等研究中发挥作用应该会是未来研究的趋势。 3) 近年来，世界各地火灾频繁，受人为活动干扰和气候变化的影响更加频繁。火灾的时空分布也变得更加难以预测。重建火灾历史，特别是近几千年来的火事件(火灾发生的频率、强度、大小、范围等)，有助于理解气候变化、人类活动和火事件之间的相互作用机制，揭示火活动的变化特征。 4) 近期火灾的发生与全球变化及人类活动密不可分，为最大化降低火灾引起的资源损失及环境破坏，火灾预防工作变得尤为重要。比如定期清除枯枝落叶以减少易燃可燃物；人为种植难燃树种以降低森林燃烧性；建设防火林带、防火线及雷电检测以阻止火灾蔓延并预防雷击火等；高温干旱气候时，适当投入人力物力做好相应的防范工作。
system: [SYSTEM]

user: 摘 要 本文给出n-cell模糊数近似大于等于关系的定义，并对近似大于等于关系的性质进行了详细的研究，重点讨论了这一关系在加、减、乘、除和数乘下的变化。 关键词 :近似大于等于，N-Cell模糊数，模糊数 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. 引言 在决策中，影响决策对象的因素常具有模糊性。为了对模糊决策对象进行排序，迫切需要考虑模糊数上的近似大于、近似等于以及近似大于等于等基本关系。 Buckley [ 1 ] [ 2 ] 率先提出了一维模糊数近似大于、近似等于以及近似大于等于关系的概念，并将其应用于具有模糊性决策对象的排序中。吴望名在 [ 3 ] 中提出了一维模糊数相容的概念。王绪柱在 [ 4 ] 中针对Buckley在文 [ 1 ] 中提出的近似相等关系给出了易于判定的充要条件，并讨论了一维模糊数近似相等的性质。但是，在实际应用中，一维模糊数往往无法满足需求。因此，王桂祥、吴从炘在 [ 5 ] 中提出了n-cell模糊数的概念，n-cell模糊数是一种特殊的n维模糊数，同时也是对一维模糊数的推广。王桂祥及其团队 [ 5 ] [ 6 ] 将n-cell模糊数应用于分类、模式识别、信息融合等领域中，受到了广泛的关注，因此，讨论n-cell模糊数上的近似大于等于关系也是非常重要的。 本文引入了n-cell模糊数上的近似大于等于关系，并详细讨论了这一关系的性质，我们的结论对n-cell模糊数排序问题的相关理论和应用有着重要的意义。
system: [SYSTEM]

user: 2. 预备知识 本节介绍了n维模糊数和n-cell模糊数的基本概念和性质，以及n-cell模糊数的代数运算，详细内容请读者参阅文献 [ 7 ] [ 8 ] [ 9 ] [ 10 ] 。 设u是 R n 上的模糊集，u可看作 R n → [ 0 ， 1 ] 的函数。任取 α ∈ [ 0 ， 1 ] ，我们称 [ u ] α 是u的a-截集，当 α > 0 ， [ u ] α = { x ∈ R n | u ( x ) ≥ α } ，当 α = 0 ， [ u ] 0 = { x ∈ R n | u ( x ) > 0 } ¯ 。 [ u ] 0 也称为u的支集，记作suppu。 若 u : R n → [ 0 ， 1 ] 满足以下性质(1)-(4)： 1) u是正规的模糊集，即存在 x 0 ∈ R n 使得 u ( x 0 ) = 1 ； 2) u是凸模糊集，即对任意 x ， y ∈ R n ， λ ∈ [ 0 ， 1 ] 有 u ( λ x + ( 1 − λ ) y ) ≥ min { u ( x ) ， u ( y ) } ； 3) u是上半连续函数，即 lim sup x → x 0 u ( x ) ≤ u ( x 0 ) ； 4) [ u ] 0 有界。 则称u为n维模糊数，全体n维模糊数记作 E n 。 用 a 表示 R n 中的点 ( a ， a ， ⋯ ， a ) ，用 表示模糊集 a ^ : R n → [ 0 ， 1 ] ， a ^ ( x ) = { 0 ， x ≠ a ， 1 ， x = a . 在不引起混淆的情况下，我们也用 a 表示 a ^ ，此时对任意的 α ∈ [ 0 ， 1 ] ， [ a ] α = { a } 。 王桂祥在 [ 5 ] 中引入了n-cell模糊数的概念：设 u ∈ E n ，若对任意 α ∈ [ 0 ， 1 ] ， [ u ] α = ∏ i = 1 n [ u i − ( α ) ， u i + ( α ) ] ，其中 u i − ( α ) ， u i + ( α ) ∈ R 且 u i − ( α ) ≤ u i + ( α ) ， i = 1 ， 2 ， ⋯ ， n ，则称u是n-cell模糊数。记 L ( E n ) 为 R n 上全体n-cell模糊数，显然 L ( E n ) ⊂ E n ，且 L ( E n ) ≠ E n 。 设 u ∈ L ( E n ) ，若 u i − ( 0 ) > 0 ， i = 1 ， 2 ， ⋯ ， n ，则称u为正n-cell模糊数；若 u i + ( 0 ) < 0 ， i = 1 ， 2 ， ⋯ ， n ，则称u为负n-cell模糊数。 N-cell模糊数的代数运算如下：对任意 u ， v ∈ L ( E n ) ， k ∈ R ， α ∈ [ 0 ， 1 ] ，有： [ u + v ] α = [ u ] α + [ v ] α = ∏ i = 1 n [ u i − ( α ) + v i − ( α ) ， u i + ( α ) + v i + ( α ) ] ， [ u − v ] α = [ u ] α − [ v ] α = ∏ i = 1 n [ u i − ( α ) − v i + ( α ) ， u i + ( α ) − v i − ( α ) ] ， [ k u ] α = k [ u ] α = { ∏ i = 1 n [ k u i − ( α ) ， k u i + ( α ) ] ， k ≥ 0 ， ∏ i = 1 n [ k u i + ( α ) ， k u i − ( α ) ] ， k ≤ 0 ， [ u v ] α = ∏ i = 1 n [ min { u i − ( α ) v i − ( α ) ， u i − ( α ) v i + ( α ) ， u i + ( α ) v i − ( α ) ， u i + ( α ) v i + ( α ) } ，                                 max { u i − ( α ) v i − ( α ) ， u i − ( α ) v i + ( α ) ， u i + ( α ) v i − ( α ) ， u i + ( α ) v i + ( α ) } ] ， [ u v ] α = [ min { u i − ( α ) v i − ( α ) ， u i − ( α ) v i + ( α ) ， u i + ( α ) v i − ( α ) ， u i + ( α ) v i + ( α ) } ， max { u i − ( α ) v i − ( α ) ， u i − ( α ) v i + ( α ) ， u i + ( α ) v i − ( α ) ， u i + ( α ) v i + ( α ) } ] (其中 v 为正n-cell模糊数或负n-cell模糊数)。 容易看出 u − v = u + ( − v ) 。 Buckley在 [ 1 ] 中给出了一维模糊数近似大于等于关系的概念： 对任意两个模糊数 u ， v ∈ E 1 ，任取 α ∈ ( 0 ， 1 ] ，定义 ν ( u ≥ v ) = sup x ≥ y min ( u ( x ) ， v ( y ) ) 。 1) 若 ν ( u ≥ v ) = 1 且 ν ( v ≥ u ) < α ，则称 u 在 α 水平下近似大于 v ，记为 u ≻ α v 或 v ≺ α u ； 2) 若 min ( ν ( u ≥ v ) ， ν ( v ≥ u ) ) ≥ α ，则称 u 与 v 在 α 水平下近似相等，记为 u ≈ α v ； 3) 若 u ≻ α v 或 u ≈ α v ，则称 u 在 α 水平下近似大于等于 v ，记为 u ≥ α v 或 v ≤ α u 。
system: [SYSTEM]

user: 摘 要 为了解决装配线平衡这类NP-hard问题，要求将工序合理地分配到各工作站，使各工作站工时尽可能相等以有效地提高生产效率。以某公司遥控车装配线的现状为案例，依据装配线实际情况设置约束条件，建立了数学模型，提出运用遗传算法对装配线进行优化改进，设计了算法流程以及编码、变异、交叉、解码和选择操作，并运用MATLAB编程实现了该算法。并且提出了工作站负荷率折线评价指标，更好的反映出工作站装配能力的发挥空间和装配线平衡理想程度。装配线平衡率由原来的68.6%上升到93.3%，装配线损失率由原来的31.4%下降到6.7%，各工作站负荷率更加平稳，优化效果显著。研究结果表明遗传算法结合MATLAB编程解决装配线平衡问题是可行的、有效的，并且提出的工作站负荷率折线评价指标为装配线平衡问题的研究提供了参考依据。 关键词 :装配线，平衡，遗传算法，工业工程，负荷率 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 ] 。装配线平衡问题是典型的NP-hard问题，问题的复杂度随任务数的增加呈几何级增长，难以在合理的时间完全解决 [ 3 ] 。依据陈星宇 [ 4 ] 的总结，关于研究装配生产线平衡问题大致分为如下5类：1) 工业工程方法，如程序分析、操作分析以及动作分析等；2) 精确算法，如分支定界法、割平面法、整数规划以及动态规划等；3) 启发式算法，如Johnson启发式算法、CDS启发式算法以及NEH启发式算法等；4) 基于仿真的方法，如通常用的仿真软件有EM-power，Witness以及Flexin等；5) 人工智能算法，如神经网络、模拟退火以及遗传算法等智能求解方法。本研究运用遗传算法对装配线平衡问题进行规划求解。
system: [SYSTEM]

user: 2. 装配线平衡问题模型建立 本研究在某公司遥控车装配线已知生产节拍的基础上，求最优工作站分配。要求在满足工序间优先关系(图1)的情况下把生产工序合理地、有效地分配到每个工作站，使得整条装配线工作站的数量以及装配线平衡率达到最优。依据相关研究，装配线模型目标函数 [ 5 ] 可以表示为： max f ( x ) = ∑ i = 1 m T i m × C T × 100 % 其中： T i 表示第i个工序的时间；CT表示生产节拍；m表示工作站的数量。 依据装配线平衡率公式，装配线的瓶颈时间对整个装配线效率有较大影响。装配线在一定的作业周 图1. 遥控车装配工序关系图 期内完成一个产品，即生产周期，是由装配线上作业最长的工作站决定的，无论其他工作站作业速度多快也只能在瓶颈作业内传送产品。如果瓶颈时间过长会导致装配线上其它工作站停工待料，造成一定的工时损失与浪费成本。因此降低瓶颈工序标准作业时间，使装配线各工作站生产均衡是装配线平衡率提高的关键。 模型约束条件为： 1) S i ∩ S j = ∅ ，( i ≠ j ； i ， j = 1 ， 2 ， ⋯ ， m )即同一个工序不能被同时分配在不同的工作站中。其中 S i 表示第i号工作站， S j 表示第i号工作站， ∅ 表示空集。 2) ∪ k S k = E ， k = 1 ， 2 ， ⋯ ， m ，即一道工序只能被分配在一个工作站中， S k 表示第k号工作站， E 表示全集。 3) T ( S k ) ≤ C T ，( k = 1 ， 2 ， ⋯ ， m )即各工作站作业时间不大于装配线的节拍。其中 S k 表示在符合优先关系矩阵A的前提下，各工序分配到第k个工作站的工序集合。 4) 作业元素优关系先矩阵： A = ( a i j ) n × n 。 其中， a i j = { 1               作 业 元 素 i 是 j 的 直 接 先 行 元 素 0               作 业 元 素 i 不 是 j 的 直 接 先 行 元 素 A = ( 0 1 1 1 1 1 1 1 0 0 ⋯ 1 0 0 1 1 1 1 1 1 0 0 ⋯ 1 0 0 0 1 1 1 1 1 0 0 ⋯ 1 0 0 0 0 1 1 1 1 0 0 ⋯ 1 0 0 0 0 0 1 1 1 0 0 ⋯ 1 0 0 0 0 0 0 1 1 0 0 ⋯ 1 0 0 0 0 0 0 0 1 0 0 ⋯ 1 0 0 0 0 0 0 0 0 0 0 ⋯ 1 0 0 0 0 0 0 0 0 0 0 ⋯ 1 0 0 0 0 0 0 0 0 0 0 ⋯ 1 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋱ ⋮ 0 0 0 0 0 0 0 0 0 0 ⋯ 0 )
system: [SYSTEM]

user: 3. 装配线平衡遗传算法设计 遗传算法(Genetic Algorithm，简称GA)是以生物遗传和进化机理为基础，从问题域可能潜在解集的一个种群开始，将可能解看做是群体的染色体或个体，并将每一个染色体进行编码操作，模拟生物进化过程对群体反复进行遗传、变异与交叉的操作。根据适应度函数对染色体进行评价，同时以全局并行搜索方式来搜索优化群体中的最优染色体，以求得满足要求的最优解 [ 6 ] 。图2是本研究规划的遗传算法流程图。 3.1. 编码 编码是应用遗传算法时首要解决的问题，针对不同的问题有不同的编码方式。遗传算法的运用过程中，编码的优劣对选择、交叉、变异等遗传操作有直接影响。本研究采序列编码方式：n个工序则编码的长度为n，基因1至n的整数按照顺序排列。如n = 5，则一个合法的染色体可以表示为如图3。 图2. 遗传算法流程图 图3. 一个合法染色体 3.2. 初始种群 根据编码方式，采用随机法产生初始种群。 3.3. 变异 变异是以小概率发生在少数基因位上的变化，它决定了遗传算法的局部搜索性能，避免收敛于某一局部最优解。本研究的变异操作采用两点互易方法，这种变异操作的实现较为简单，能提高运算效率，步骤如下：Step 1，随机产生两个自然数r1、r2；Step 2，交换第r1位和r2位的基因。如r1 = 2，r2 = 4时，染色体变异操作如图4所示。 3.4. 交叉 交叉是按照较大概率从群体中选择两个染色体，交换两个染色体某个或某些部位来获得新染色体，从而推动整个群体进化。本研究采用两点交叉法，步骤如下：step 1，随机选择两个染色体作为父代染色体；step 2，在1到n之间随机产生2个自然数r1和r2；step 3，交换两个父代染色体r1至r2之间的基因片段获到两个子代染色体，为避免子代染色体冲突进行修订处理。如r1 = 2，r2 = 4时，染色体交叉操作如图5所示。 3.5. 解码 本研究解码操作是依据工序约束对编码进行转换，然后依次安装节拍约束将工序填入各个工作站中，将染色体转换成满足适应度的值。问题可以划分两个部分：求工作站数和装配线平衡率，首先在给定生产节拍的条件下，根据工序之间优先关系求得工作站的工序划分，同时得到工作站数；其次，根据工作站上工序的划分求得工作站工作总时间计算装配线平衡率。 图4. 变异操作 图5. 交叉操作 3.6. 适应度函数 适应度函数是用来表明染色体或解的优劣性，是进行自然选择的唯一标准。在遗传算法中针对目标函数制定出适应度函数来衡量群体中各个染色体的优良。适应度较优的染色体遗传到子代的概率就较大；而适应度较差的染色体遗传到子代的概率就相对小。本研究的目的为达到装配线平衡率最优，因此选用上述的装配线平衡率作为适应度函数。 F ( p ) = ∑ i = 1 m T i m × C T × 100 % 3.7. 选择操作 选择操作依据各个染色体的适应度的大小进行选择，本研究的选择策略采用轮盘赌选择法 [ 7 ] ，以各个染色体的选择概率与其适应度成正比为原则，适应度较优的染色体遗传到子代的概率就较大；而适应度较差的染色体遗传到子代的概率就相对小。设种群数量为t，种群中第i个染色体适应度为 f i 。选择操作的步骤如下： Step1，计算种群的总适应度： F = ∑ i = 1 t   f i ； Step2，计算第i个染色体被选择的概率： P i = f i / F ； Step3，计算第i个染色体的概率累积： P s ， i ； Step4，取0至1之间的随机数a，若 a ≤ P s 1 ，则第一个染色体被选择；否则，若 P s ， i − 1 ≤ a ≤ P s ， i ，则第i个染色体被选择； Step5，重复步骤(4)直到选择到第t个染色体结束。
system: [SYSTEM]

user: 4. 某公司装配线平衡现状分析 某生产制造型公司成立于2014年，规模达到中型水平。公司在成立的第二年通过了国际质量体系的相关认证，公司在产品研发上独立自主并拥有该公司自己的设计技术体系，同时具备生产制造能力，公司主要产品是遥控车。该公司通过程序分析将遥控车的装配分为21个步骤，每一个步骤就是一道工序，也是最小的装配作业要素，并规定各工序间的紧前工序，即工序间的优先关系，通过时间测定的方法规定各工序的标准作业时间。遥控车装配步骤与工序时间见表1。 该公司为了完成装配遥控车的一道或者多道工序，在装配线指定一个工作位置，工人在其中进行操作，这个位置就是工作站，每一个工作站由1~3名工人操作。该公司根据遥控车装配线的实际要求，将21道工序分配到8个工作站。该公司的这种分配方式在实际运用中并没有考虑到工站间的生产能力是否均衡，这种规划造成了在生产时间上浪费同时使公司在费用上造成了不必要的开支。所以本研究针对该公司遥控车装配线平衡问题进行分析改善。遥控车装配线工序分配现状与工作站负荷率见表2。 通过表2可知，第1、2、3道工序分配到了1号工作站，负荷率为0.7；第4、5道工序分配到了2号工作站，负荷率为0.823；第6、7、8道工序分配到了3号工作站，负荷率为0.823，同理可知其他工序分配情况与工作站负荷率。 表2中工作站4工作总时间最长，工作总时间为68 s，是遥控车装配线的生产瓶颈同时也是该装配线的生产节拍。生产瓶颈主要指的是在整个生产节拍中耗时最长的工作站，生产瓶颈不仅制约了整个装配线的产出速度，而且会影响道装配线其他工作站装配能力的发挥。 图6中横坐标表示各工作站，纵坐标表示各工作站的负荷率。工作站负荷率越低表示工作站装配能 Table 1 Telecar assembly procedure and process schedul 工序编号 紧前工序 工序描述 工作时间(s) 1 - 装配车轴 25 2 1 装配车圈 12 3 2 装配车胎 12 4 3 组装后发动机箱 36 5 4 装配固定机箱螺丝 20 6 3 装配支架前和发动机 9 7 6 装配箱盖和前轮轴 32 8 7 装配螺丝并将其紧固 15 9 - 装配电池电极片 13 10 9 焊接电极和开关导线 16 11 - 焊接线路板 23 12 11 焊接线路板导线 16 13 10、11 装配开关和螺丝 23 14 12、13 装配线路板 12 15 5、8、14 连接车架箱与前后机 37 16 - 装配两侧挡板 12 17 - 装配后挡板和前后灯 15 18 15、16、17 组合车身 29 19 18 粘图 7 20 19 检测 6 21 20 装配电池壳 3 合计 373 表1. 遥控车装配步骤与工序时间表 Table 2 Status of process distribution and load rate of workstation 工作站编号 总时间/s 工作站内的工序编号 工作站负荷率 1 49 1、2、3 0.7 2 56 4、5 0.823 3 56 6、7、8 0.823 4 68 9、10、11、12 1 5 35 13、14 0.515 6 37 15 0.544 7 27 16、17 0.397 8 45 18、19、20、21 0.332 合计 373 表2. 工序分配现状与工作站负荷率 图6. 工作站负荷率折线图 力发挥空间越大，各工作站之间的负荷率数值越平稳，装配线平衡越理想。通过图6工作站负荷率折线图可看出，各工作站负荷率折线呈现不平稳的状态，说明装配线平衡不理想；第5号、6号、7号、8号工作站负荷率偏低，通过科学的工序与工作站分配，装配能力提升空间很大。 其中各工作站负荷率 = T ( S K ) C T ，( k = 1 ， 2 ， ⋯ ， 8 )。 T ( S K ) 表示第k号工作站总时间，CT表示生产节拍。 依据装配线平衡率公式得到该公司遥控车装配线的装配线平衡率为68.6%，平衡损失率为31.4%。有相关研究 [ 8 ] 表明装配线平衡率和企业利润之间的关系在不同阶段呈现着不同的关系：装配线平衡率在89%之前，利润呈现增长趋势；装配线平衡率在89%至97%之间，利润基本呈现平稳的趋势；装配线平衡率高于97%时，利润会呈现下降趋势。该公司装配线平衡率在68.6%，对企业利润有很大的提升空间。
system: [SYSTEM]

user: 5. 运用Matlab2014a编程求解 本研究基于遗传算法对某公司遥控车生产装配线平衡问题建立模型，在满足工序间优先关系与生产节拍的前提下，对各工序的先后顺序与工作站分配进行调整。并通过Matlab对该模型进行编程求出最优解，在编程中设置生产节拍不大于68，初始种群为50，种群交叉概率为0.7，变异概率为0.1，迭代后结果见图7与图8。 如图7所示，图7的横坐标表示迭代次数，纵坐标表示装配线均衡总目标，即装配线平衡率。本研究运用遗传算法通过设置迭代次数为100，得到装配线平衡率数值为0.933。通过观察图中曲线可知，在迭代次数为18后曲线数值稳定在0.933，达到最优装配线平衡率。 如图8所示，图8的横坐标表示工作站数量与序号，纵坐标表示工作站的工作总时间，直方图中的序号表示21个工序分配到工作站的情况。例如，通过图8可以看出依据遗传算法的执行将21的工序分配到8各工作站中，其中第1号工作站分配有第1、2、9道工序，并且总时间为50 s；第2号工作站分配有第3、4道工序，并且总时间为48s；第3号工作站分配有第6、10、11道工序，并且总时间为48 s。 通过表3分析可知：各工作站得到了重新分配，被安排在各个工作站的工序更加合理。装配线平衡率显著提高，由原来的68.6%上升至93.3%，装配线平衡损失率也由原来的31.4%，下降至6.7%；工 图7. 遗传算法迭代曲线 图8. 工序分配情况 作站总时间最大为50 s，因此改善后的生产节拍由原来的68 s缩短至50 s。根据图9改善后的工作站负荷率折线图可以看出，改善后的各工序负荷率较为平稳。通过对遥控车装配线的优化，各工作站的工作时间与工作站负荷率更加平稳，有效的提高了设备的利用率。
system: [SYSTEM]

user: 6. 结论 本研究基于遗传算法，针对某公司遥控车装配生产线现状建立数学模型，设计了遗传算法流程以及编码、变异、交叉和选择操作，并运用matlab2014a编程求解，使装配线平衡率由原来的68.6%上升到93.3%，装配线损失率由原来的31.4%下降到6.7%，各工作站负荷率更加平稳，优化效果显著；本研究 图9. 改善后的工作站负荷率折线图 Table 3 Improved process allocation and workstation load rat 工作站编号 总时间/s 工作站内的工序编号 工作站负荷率 1 50 1、2、9 1 2 48 3、5 0.96 3 48 6、10、11 0.96 4 48 5、12、16 0.96 5 50 13、14、17 1 6 47 15 0.94 7 37 16、17 0.74 8 45 18、19、20、21 0.9 合计 373 表3. 改善后的工序分配与工作站负荷率 验证了遗传算法结合MATLAB编程在解决装配线平衡问题上的可行性和有效性，同时验证了遗传算法在该领域上具有实际意义；提出的工作站负荷率折线评价指标更好的反映出工作站装配能力的发挥空间和装配线平衡理想程度，为装配线平衡问题的研究提供了参考依据。
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 ] [ 3 ] [ 4 ] 。偏移量是组合系数和为零的同层的相邻几个顶点的线性组合。在已有细分迭代规则的基础上加上偏移量做调整，既可以从逼近型细分构造插值型细分 [ 1 ] [ 2 ] [ 3 ] ，又可以从插值型细分构造逼近型细分 [ 4 ] 。上述工作局限于单个细分格式，没有从理论上对这种方法进行探讨。Luo等 [ 5 ] 指出添加偏移量的操作对应的是在原细分生成函数基础上加上一个Laurent多项式。求和规则是细分的基本性质 [ 6 ] [ 7 ] [ 8 ] ，亓万锋等 [ 8 ] 分析了采用二阶中心差分形式的偏移量所能达到的最高阶求和规则的问题。本文将 [ 8 ] 的结论推广，考虑了高阶中心差分是否能够达到最高阶求和规则的问题。
system: [SYSTEM]

user: 2. 预备知识 细分格式是从给定的初始控制顶点集 P 0 = { p i 0 ∈ ℝ d } i ∈ ℤ 和一个正整数m，利用有限个数 a = { a i } i ∈ ℤ ∈ l 0 ( ℤ ) ， l 0 ( ℤ ) 代表 ℤ 上支撑有限的序列的集合，通过迭代格式 p i k = ∑ j ∈ ℤ a i − m j p j k − 1 ，     i ∈ ℤ ，   k = 1 ， 2 ， ⋯ (1) 迭代k次后得到控制顶点集合 P k = { p i k ∈ ℝ d } i ∈ ℤ 的一种造型方法，其中 m 称为重数， a = { a i } i ∈ ℤ 称为细分掩模，对应的Laurent多项式 a ˜ ( z ) = ∑ i ∈ ℤ a i z i ， z ∈ ( ℂ \ { 0 } ) 称为细分生成函数 [ 8 ] 。 给定m重细分掩模 a = { a i } i ∈ ℤ ∈ l 0 ( ℤ ) ，记 supp     a = { i ∈ ℤ ， a i ≠ 0 } 。在实际应用中，往往使用对称的生成函数即细分掩模满足 a i = a − i ， i ∈ ℤ ，显然生成函数是对称的当且仅当 a ˜ ( z ) = a ˜ ( z − 1 ) 。长度为s的列向量 v = ( v 1 ， v 2 ， ⋯ ， v s ) T 称为中心对称的，如果满足 v i = v s + 1 − i ， i = 1 ， ⋯ ， s 。显然对称的生成函数的系数组成的列向量是中心对称向量。 细分格式可以分为插值型和逼近型两种格式，插值型细分每次迭代时保留上一次的控制顶点，其细分掩模满足 a m i = δ i 。若细分格式不是插值型的，则称为逼近型细分。 称m重细分格式 a ˜ ( z ) 满足k阶求和规则，如果 a ˜ ( 1 ) = m 且 a ˜ ( z ) 含有因式 ( 1 + z + ⋯ + z m − 1 ) k 。易见当 a ˜ ( z ) 满足一阶求和规则时，有 ∑ i a m i = ∑ i a m i + l = 1 ， 0 ≤ l < m 成立。
system: [SYSTEM]

user: 资助信息 国家自然科学基金(No. 61502217)，辽宁省教育厅科研项目(L20168366)。
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 ] 。子宫内膜异位症是育龄期妇女的常见病、多发病，统计资料表明，其发病率近年来呈上升趋势，其病理特点表现为：病变部位多发且广泛、呈局部浸润性，可引起周围组织的粘连，激素依赖且具有正常子宫内膜的分泌功能 [ 3 ] 。异位的子宫内膜最长出现于卵巢，其次多见于子宫直肠陷凹、子宫浆膜层、直肠及乙状结肠、腹膜脏层。目前对于子宫膜异位症的治疗尚无统一的方案且治愈率低，在西医中多采取激素类药物抑制雌激素合成以及手术切除病变部位等方法，而在中医中多采用具有补血调经、活血止痛等功效的方剂以及局部治疗进行治疗 [ 4 ] [ 5 ] 。下面将目前对于子宫内膜异位症的中西医治疗研究现状作一综合性叙述：
system: [SYSTEM]

user: 2. 西医病因病理 关于子宫内膜异位症的西医病因病理尚未明确，目前的研究认为该病主要与免疫因素、血流动力学及血管因素、遗传易感性等有关。临床上较为流行的说法主要包括种植学说、上皮化生学说、诱导学说等 [ 6 ] 。其中种植学说又包括经血逆流学说、血管淋巴管播种以及医源性种植等学说。我国有学者 [ 7 ] 认为子宫内膜的异生主要是由于子宫内膜在机体全身以及局部免疫的影响下，经局部细胞因子、激素、炎症因子的共同作用通过局部黏附的血管迁移而造成种植。除上述学说以外，化学物质如二噁英等化学物质的接触也是子宫内膜异位症的一个可能的病因 [ 8 ] 。
system: [SYSTEM]

user: 3. 中医病因病机 祖国医学的“痛经”、“血瘕”、“不孕”、“月经不调”与子宫内膜异位症的常见临床症状和表现基本相符。早在《素问·平人气象论》中就有关于”寸口脉沉而弱，曰寒热及疝瘕，少腹痛。”的记载了；《金匮要略·妇人杂病脉证并治》将此病记载为：“带下，经水不利，少腹满痛，经一月再见者，土瓜根散主之。”隋代巢元方《诸病源候论》“令人腰痛，不可以俯仰，横骨下有积气，牢如石，小腹里急苦痛，背膂疼，深达腰腹下挛，阴里若生风冷，子门擗，月水不时，乍来乍不来，此病令人无子。”而近代中医理论认为该病的发生发展是多种因素共同引起的，主要包括生活不洁、外感六邪，加之七情所伤、体弱等而患 [ 9 ] [ 10 ] 。 根据子宫内膜异位症的“月经不调”“月经量多”“不孕”等症状，《景岳全书·妇人规》将其将其病机记载为“经行腹痛，证有虚实.......但实中言虚、虚中亦有实”然而“久病必虚”，而子宫内膜异位症患者其病程通常较长，且淤血久聚体内，可见对于子宫内膜异位症的治疗主要以活血化瘀、补血调经等治疗为主 [ 11 ] 。
system: [SYSTEM]

user: 4. 子宫内膜异位症的西医治疗 目前西医中对于子宫内膜异位症的治疗主要以激素药物阻断正常月经周期为主，对于重症患者可采用手术及介入手术治疗。 4.1. 子宫内膜异位症的西医药物治疗 ① 促性腺激素受体激动剂(GnRH-α)是目前最为常用的一种药物，该药是人工合成的一种十肽类化合物，其主要作用为与促性腺激素(GnRH)竞争卵巢上的受体，从而造成GnRH无法产生作用，使卵巢在无GnRH的正常作用下而分泌雌孕激素，从而引起人工闭经，该药的优点在于其对分解该药的酶反应力低下，可长期持续作用，因此该方法又称为药物性卵巢切除。目前国内常用的GnRH-α药物主要包括亮丙瑞林、戈舍瑞林及达菲林，其单独使用即可治疗子宫内膜异位症，但长期应用该药可造成骨钙丢失、阴道干涩、潮热等不良反应，故多应用数月至半年 [ 12 ] 。 ② 米非司酮：米非司酮是一种孕激素受体调节剂，其化学本质为一种甾体，主要作用机制为与孕激素竞争孕激素受体，并且其对下丘脑-垂体-卵巢激素分泌调节轴也具有影响，长期应用该药导致孕激素分泌低下同时造成卵巢萎缩；自该药上个世纪九十年代应用于子宫内膜异位症的治疗以来，由于其不良反应较少、疗效好，被广泛应用于临床，该药无促性腺激素受体激动剂对骨质的影响，常用治疗方案为每日应用米非司酮25~50 mg，连续应用3~6个月，但关于米非司酮的远期疗效以及不良反应和详细机制的研究仍待发掘 [ 13 ] 。 ③ 芳香化酶抑制剂：芳香化酶不仅存在于卵巢，也存在于皮肤、乳腺、脂肪组织中，其作用主要为催化雄烯二酮转化为雌酮，雌酮在相关酶的作用下进一步转化为雌二醇；目前相关研究已表明，在子宫内膜异位症的病人体内，芳香化酶存在着高表达，因此芳香化酶抑制剂对于芳香化酶的抑制可引起患者体内雌激素水平降低，其还能对局部病灶的雌激素生成功能进行一定的抑制，该药多用于绝经后期患者，目前临床上多用的药物包括阿那曲唑和来曲唑 [ 14 ] 。 ④ 其他药物治疗：目前关于其他药物治疗子宫内膜异位症的主要方法还包括抗血管生成药物、抗炎药物、促凋亡和抗增殖药物，还包括一些新型的肝细胞疗法等。 4.2. 子宫内膜异位症的手术治疗 对于子宫内膜异位症的手术治疗，目前有常规手术治疗和介入手术治疗两种方式； ① 常规手术治疗：手术治疗主要适用于重症患者而药物治疗了疗效差或者不能起效者或者是迫切希望保留生育能力者或卵巢囊肿达5 cm以上者。根据手术方式的不同可将手术分为保守性手术、半保留手术、根治性手术。主要术式包括子宫切除术、子宫及附件切除术。前者通过切除子宫但保留了卵巢的内分泌功能，患者仍有生育功能，后者患者无生育功能为根治性手术，根治性手术的优点是疗效好、彻底清除病灶；保留手术虽保留了患者卵巢的排卵及内分泌功能，但术后复发率高，统计资料表明，保守性手术患者术后复发率高达20%~40%。目前腹腔镜对于子宫内膜异位症的手术治疗应用效果较好，具有微创、术后恢复快等优点。值得注意的是仍有部分患者在术后仍需药物治疗 [ 15 ] 。 ② 介入治疗主要针对的是具有囊肿的患者，对患者进行超声引导下抽液硬化治疗，其疗效好、对患者损伤小，但国内有学者 [ 16 ] 报道介入治疗术后仍有20%左右的复发率，但总的说来介入治疗的疗效较好，可值得在临床工作中大力推广应用。 4.3. 子宫内膜异位症的手术联合药物治疗 子宫内膜异位症患者的手术联合药物治疗主要有术前及术后应用药物治疗两种方式，有研究者在对患者术前应用GnRH-α，连用三个月后进行手术，术中患者的出血量均比不用术前GnRH-α治疗的患者更好，术后随访5年内，其复发率也更低 [ 17 ] 。另一方面，由于子宫内膜异位症是一种激素依赖性疾病，故在对患者术后应用抗激素治疗时，可促使残留病灶吸收，降低复发率，因此手术联合药物治疗是一种值得探究及应用的方法。