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

user: 6. M角数的研究现状和展望 比较文献 [ 2 ] [ 4 ] [ 5 ] [ 6 ] [ 7 ] 和表2，容易发现，当 取适当的有理数时，就是文献 [ 2 ] [ 4 ] [ 5 ] [ 6 ] [ 7 ] 研究的内容，它们都属于求解二元二次Diophtantus方程。其表达式和递推公式只是表2中数组 的两个特殊性质。 另一方面，由表2的规律性数字和上面给出的拓扑空间可以直观的想象出，满足条件为 的数列在表2的“反面”，属于复数域，故不存在满足条件的项为自然数的递推公式。因此，满足条件为 的M角数关系式只能在不同 的角数中寻找。例如：我们计算很少几个有理数组 ，就可以得到3角数和4角数的关系式(Plutarch定理 [ 1 ] ) 它是一个二次的代数关系式，但它不是满足条件 的递推公式。 众所周知，有一些数学问题用代数的方法解决很困难，另有一些数学问题用几何的方法解决很困难，还有一些数学问题仅用代数方法或仅用几何方法解决都是困难的。但很少有文献阐述出现这种现象的根源是什么？ 本文的逆定理明确指出了产生这种困难的原因是高次方程多根的选择问题，由上面的讨论我们又可以看出，高次方程多根的选择问题可以扩充公式(4)中自然数 为有理数，并用拓扑的性质来替换。也就是说，拓扑几何把方程的求根和根的选择这两个步骤合并在了一起。 因此，涉及无穷概念(这里指偏序集)的连续问题，用几何的方法更直接；涉及无穷概念的离散问题，用代数的方法更直接。交替涉及连续和离散的问题，仅用代数方法或仅用几何方法解决都是困难的。 另一方面，高次方程多根的选择和公式(4)中自然数 扩充为有理数是等价的，都属于偏序集研究的范畴。由此可知，Diophtantus方程中未解决的离散问题(其中绝大部分都可以归类为分析学的问题)的解决方案，应该主要在偏序集理论中去寻找，直观的拓扑几何理论和许多代数理论很难解决这类问题，即使使用上奇异微积分理论也是如此。遗憾的是，偏序集理论远未达到完善的程度。 继续讨论上述的解决方案，不可避免的会涉及多个数学分支。不再这里赘述。 实际上，尽管本文给出了M角数的离散代数性质和拓扑几何性质，并对亏格为1的Diophtantus方程建立了代数和几何的联系。但M角数更多重要性质还远远没有被挖掘出来。
system: [SYSTEM]

user: 摘 要 社会信用体系是金融基础设施建设的重要组成部分，对实现普惠金融具有重要意义。为推进农村普惠金融发展，提高农村金融服务的覆盖率和可得性，贵州省兴仁县农村商业银行结合当地的实际情况积极开展农村信用体系建设，探索出有特色的农村信用体系建设模式。本文通过总结兴仁县农村商业银行开展信用体系建设的具体实践，分析其在信用体系建设过程中取得的成效，进而对在全国范围内推进社会信用体系建设提出政策建议。 关键词 :社会信用体系建设，普惠金融，兴仁县农村商业银行 Copyright © 2017 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.1. 普惠金融的现实困境 根据世界银行的定义，普惠金融是指“立足商业可持续的原则，保证社会各阶层，特别是弱势群体以及低收入人群的金融需求能够有尊严地以合理的成本得到满足”，党的十八届三中全会在《关于全面深化改革若干重大问题的决定》中把推进“普惠金融”作为重要的工作任务；但截至2014年末，全国仍然存在1570个金融服务空白乡镇 [ 1 ] ，农村地区金融排斥现象依然严重。所以，普惠金融的难点在农村，农村普惠金融的推进直接影响整个社会普惠金融目标的实现。农户贷款可得性是农村普惠金融的一项重要内容。截至2014年末，金融机构农户贷款余额为5.4万亿元，仅占全国各项贷款比重的6.4% [ 1 ] ，且绝大部分贷款为担保贷款，极少对农户发放信用贷款，而农户恰恰是最缺乏担保物或担保人的群体。因此，如何开发农户信用资产，使金融机构敢于发放农户信用贷款，是解决农户贷款难问题的重要途径。但是，目前社会信用体系建设的严重滞后导致农户违约成本低，农户信用不值钱，金融机构对信用贷款望而却步，因此，完善社会信用体系对于推进农村普惠金融事业具有重要意义。 1.2. 社会信用体系建设是实现普惠金融的有效途径 国务院于2014年6月发布的《社会信用体系建设纲要(2014~2020)》指出，社会信用体系是一种社会机制，它的核心作用在于记录社会主体的信用状况，揭示社会主体的信用优劣，并整合全社会力量褒扬诚信，惩戒失信。建设社会信用体系的主要目的是要提高全社会的诚信意识和信用水平，为发展普惠金融事业奠定良好的社会基础。农村社会信用体系建设是我国社会信用体系建设工作的重点和难点，也是整个社会普惠金融目标实现的关键环节。建立和完善农村社会信用体系，有利于增强农民的信用意识，提高农村金融服务的覆盖率和可得性，推动农村普惠金融的实现。在社会信用体系缺位的情况下，失信惩戒机制不能有效建立，农户失信行为难以受到惩罚，农户的守信行为也得不到更好的激励，进而导致农村金融生态恶化。通过采集农户的信用信息，建立所有农户电子化信用信息档案，是建立农村社会信用体系工作的重要组成部分。在此基础上，金融机构可以根据农户的信用等级和信用记录为贷款农户提供不同的利率优惠政策，让农户逐步意识到诚实守信的价值，进而增强农民的信用意识；同时，电子化的信用信息也可以帮助金融机构了解农户的基本情况，有效判断农户的信用水平，降低金融机构和农户之间的信息不对称，简化贷前调查的程序，节约放贷成本，提高金融机构为农村地区服务的积极性。
system: [SYSTEM]

user: 2. 兴仁县农村商业银行信用体系建设实践 兴仁县农村人口占总人口的94%，是典型的山区农业县。该县于2011年正式成为贵州省农村社会信用体系建设试点县，并于2014年被中国人民银行总行列为贵州省唯一全国信用体系建设试验区。兴仁县农村商业银行(以下简称“农商行”)是兴仁县农村社会信用体系建设的开拓者与实践者。农商行在兴仁县委、县政府的领导下以及县金融办、人民银行等的支持下，积极开展社会农村信用体系建设工作，已经探索出有特色的社会信用体系建设模式。 2.1. 根据评定标准，稳步推进信用评级工作 兴仁县农商行根据政府出台的《兴仁县农村信用等级评价暂行办法》稳步推进信用评级工作。该工作主要分为两部分：一是开展农户的信用评级工作。农商行负责人逐户进入农户家中进行调查，根据农户的家庭经济水平、偿债能力、信用状况以及村内其他人员对农户的评价等多个维度进行综合评价，将农户划分为“特优”、“优秀”、“较好”、“一般”以及等外级五个等级。评级完成后，农商行对评级结果进行公示。二是开展信用组、信用村、信用乡镇评定工作。农商行在农户信用评级的基础上建立了“申请、评审、授牌、年检”环环相扣的信用组、信用村、信用乡镇评定制度。农户的信用等级每一季度更新一次，信用组、信用村、信用乡镇每年考核评定一次。农商行针对不同的信用等级制定相应的利率优惠政策，让农户切身感受良好信用所带来的更加实惠的金融服务。信用评级工作的开展提高了农户的违约成本，降低了农户的不良贷款率。农村地区是一个熟人社会，农户的生产与生活活动具有一定的社区性，农商行公示农户信用等级的做法会提高农户的信息透明度。信用状况较差的农户在重复交易的过程中会受到农商行以及其他交易对手的排斥，从而提高其违约成本。因此，社会信用体系建设会促进农户自我规范、自我约束，降低道德风险的发生概率，进而降低农户贷款的不良率。 2.2. 创新金融产品和金融服务，培育农户的信用意识 金融知识贫乏、诚信意识差是偏远地区农户面临的普遍问题。为普及金融知识、培养诚信意识，兴仁县农商行在县委领导以及金融办等多方支持下于2015年4月在村级行政单位设立金融服务站，并于2015年6月逐步开展金融夜校工作。金融服务站的设立架起了农商行和农户之间金融服务的桥梁，有利于推进农村信用体系建设的步伐。一方面，金融服务站在为农户提高服务的过程中，可以采集到农户更多的信息，提高信用评级的准确性和有效性。另一方面，服务站通过举办“金融知识大讲堂”等开展金融知识宣传，不仅有针对性的为农户提供知识、技术、信息等方面的服务，同时也可以提高农户对农商行服务以及信用体系的认识。金融夜校开展的工作始终结合各个农村地区的实际情况，以更好地实现对农户宣传教育的目的。夜校主要利用农闲时间向农户宣传金融知识和诚实守信的重要性，有时也会利用地方党政部门、村委会召开的各类宣传会、动员会、推进会等对农户进行宣传或培训。金融夜校的宣传和培训都始终围绕农户关心的金融问题展开阐述，采用通俗易懂的语言表达，通过列举典型案例、现场演示以及有奖问答等方式让农户真正坐得住、听得懂、记得牢、用得好。夜校工作的开展帮助农户了解了金融政策，熟悉了农商行的各项服务，增长了农户的金融知识和信用意识。截至2016年6月末，兴仁县农商行共开展金融夜校291次，培训人次16，506人，覆盖了全县162个行政村。金融服务站和金融夜校对普及金融知识、培育农村信用主体的信用意识具有积极作用。农商行面对面的宣传培育方式，不仅可以拉近自身与农户之间的联系，强化农民对农商行金融产品及服务的了解，同时也可以提升农民防范和识别金融风险的能力，培育农户的信用意识。 2.3. 建立合理的奖惩机制，改善农村信用环境 农商行为营造“守信激励，失信惩戒”的社会氛围，建立了合理的贷款制度和奖惩机制，让农户切身感受诚实守信的重要性，进而达到改善农村信用环境的目的。农商行在为农户评级授信后，会与农户签订《农户小额信用贷款最高限额循环借款合同》。合同中规定，在农户没有逾期未还的贷款余额或利息时，农商行会按照“余额控制、周转使用、随用随贷、动态管理”的原则为农户在最高授信额度内提供贷款。该原则简化了农户贷款的流程，有利于降低农户贷款成本。但是，当农户的贷款余额出现逾期时，农商行也会严格按照合同中约定的逾期还款的惩罚利率收取农户的罚息。当然，在农户还清贷款以及交清罚息前，农商行不会为农户提供贷款。但是，只要农户可以还清之前的欠账，农商行依然会一如既往的为农户提供所需贷款，不会因之前的不良记录对农户产生偏见，从而促进了农商行与农户之间的良性互动。
system: [SYSTEM]

user: 3. 兴仁农商行推进农村社会信用体系建设取得的成效 兴仁县农商行在推进社会信用体系建设的基础上，通过创新金融产品和金融服务，建立了有效的信用激励机制，明显改善了农村信用环境。同时，社会信用体系建设工作的开展也缓解了农户的融资瓶颈，提高了便民利民的服务水平，显著推进了普惠金融的发展。 3.1. 农户贷款难题明显改善 兴仁农商行在社会信用体系建设过程中，逐渐提高农户和农商行之间的信息透明度，改变了传统的农户贷款以担保贷款为主的方式，有效缓解了农户贷款难的问题。截至2016年6月末，兴仁县共创建信用乡镇10个，信用村126个，信用组1492个，为92139户农户进行了建档评级授信，授信额度为336498万元，户均授信额度为36，520元。同一时间，农商行累计为农户提供18464笔贷款，贷款农户的平均贷款额度为53，848元。其次，农户贷款覆盖面有明显的提高，由2010年的46.89%上升到2016的68.54%。最后，根据在兴仁县2个乡镇4个行政村的调研发现，约有85%的农户表示在2015年到2016年6月底期间有借款行为发生，其中71.5%的贷款都是从兴仁农商行获得的。 3.2. 农村信用环境明显改善，农民信用意识显著提高 农商行的宣传教育工作提高了农民的信用意识，使整个农村的信用环境也得到了极大的改善。自兴仁县农商行开展农村信用体系建设以来，信用乡镇、信用村以及信用组的个数得到了成倍的增长，信用乡镇个数由3个上升到10个，信用村的个数由13个上升到126个、信用组的个数由497个上升到1492个。农户信用证持证比例由原来70.23%上升到89.45%，持证比例逐年上升。而且，建档农户的不良贷款率逐渐降低，截止到2016年6月末，不良贷款率由最初的2.38%降到1.71%。在创建“信用村”、“信用组”活动中，村干部积极协助金融机构清收不良贷款，参与金融机构的贷后管理；一些有逾期贷款的农户或农户子女主动归还欠款，积极参与信用等级评定，逐渐形成了诚实守信的良好风气。 3.3. 显著提高了金融服务的覆盖率和可得性 农村社会信用体系的建立，缓解了农村地区农商行与农户之间信息不对称的“瓶颈”，进一步激发农商行创新金融产品和金融服务的积极性。农商行先后在农村地区设立了便民服务点、金融服务站以及惠农“村村通”，显著提高了农村地区金融服务的覆盖率和可得性。金融产品和金融服务的创新有效解决了广大贫困地区农户、返乡农民工等的信贷需求，打通了金融服务的“最后一公里”。截止2016年6月，金融服务站登记外出务工人员861人，存款合计达到442万元。另外，兴仁县162个行政村共布设了186台“村村通”助农服务终端，实现了助农服务的全覆盖，有效的解决了外出农民工往返携带现金的不便，实现了农民工在外打工挣钱，在家门口取款的愿望。惠农“村村通”等的设立不仅为农户提供了便利，而且也降低了农户取款与转账的成本。
system: [SYSTEM]

user: 4. 政策建议 兴仁县农商行在积极推进农村社会信用体系建设实践中显著推进了当地农村普惠金融的发展，其开展信用体系建设工作的具体实践，对在全国范围内推进社会信用体系建设具有重要的借鉴意义。 4.1. 加大政府支持力度，构建多方配合的工作体系 农村社会信用体系建设的涉及面非常广，是一个系统的工程 [ 2 ] 。地方政府的大力支持以及各级乡镇政府、村委会的积极配合是保证信用体系建设工作有序推进的关键因素。一方面，地方政府不仅要提供适当的经费支持，以确保信用体系工作的顺利开展；也要为信用体系建设出台一系列措施办法，为信用体系建设工作的顺利开展提供有力的保障。另一方面，当地的金融机构作为信用体系建设工作的具体实践者，要与其他金融相关部门以及各级乡镇政府、村委会建立密切联系。同时，金融机构作为政府与农户之间的联结者也需要根据政策目标细化政府的各项政策措施，以确保信用体系建设工作的稳步推进。 4.2. 逐步建立信息系统，保持信息系统的动态更新 农村社会信用体系的建设不是一次性工作，它需要长期持久的维护，并要定期更新信息系统，只有这样才可以保证信息系统的有效性。信用体系建设初期需要做到应评尽评，通过进村入户、扫街扫巷的形式，对辖内农户逐户开展建档、评级、授信工作。同时，要通过多方位多角度对农户的信息进行评价以保证信息真实可靠。信息系统的后期维护也是信用体系建设不可或缺的部分，要做到定期检查，及时更新。信息系统的更新需要结合线上线下两种审核方式。线上审核是指要根据农户的贷款信息、还款情况，定期审核农户的信用状况；线下审核可以采取兴仁农商行责任负责制的做法，即1名信贷员负责1个村庄，信贷员每天主要的工作任务就是了解村内农户的基本情况，与农户建立密切的联系。线下审核方式是线上审核的一种重要补充，可以更深入了解农户的信息。因此，线上和线下两种审核方式的结合可以保证信息系统动态更新的准确性。 4.3. 开展宣传教育工作，培育农户信用意识 提高农户的信用水平是社会信用体系建设的重要工作之一。因此，应大力开展信用知识以及金融知识等的宣传教育工作，增强农户的信用意识，提高农户的信用水平。在农村地区开展金融知识的宣传教育工作时，要结合农民自身特点，针对不同的受众群体采取不同方式开展。一方面，可以通过广播、横幅和宣传标语等形式进行集中宣传 [ 3 ] ；另一方面，要结合当地的特色在每天的农闲时间以不同的方式吸引农户参加培训活动。总之，要通过开展宣传教育工作努力营造“守信受益、失信惩戒”的社会氛围，让农户深刻认识到诚实守信的价值，进而增强农户的信用意识，改善农村的信用环境。 4.4. 充分利用信用信息，创新金融产品和服务 在农村社会信用体系建设过程中要充分利用农户的信用信息创新金融产品和金融服务，以促进农村信用信息系统的进一步完善。针对不同信用等级的农户给予有差别的优惠政策与服务，对信用好的农户在利率、服务等方面给予更优惠的政策，让农户深刻认识到诚实守信的价值。另外，金融机构应加强对信用信息产品和服务的开发，将信用信息产品和服务应用到金融服务空白的乡镇和村庄，促使金融机构为更多的农户服务，从而进一步推进农村信用体系的建设。
system: [SYSTEM]

user: 摘 要 Lattice Boltzmann Method (LBM)是一种近年来发展的基于分子动理论的数值方法，它具有计算简单、边界易于处理的特点。本文基于采用D2Q9模型的LBM方法对NACA0012翼型进行了数值模拟与分析，数值结果表明：对于简单翼型的低雷诺数、非定常湍流流场，LBM方法的数值结果与CFL3D参考值可以很好的吻合。因而，该方法对于低雷诺数下简单翼型流场计算是一种有效的计算手段。 关键词 :LBM方法，D2Q9模型，翼型，低雷诺数 Copyright © 2017 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. 引言 Lattice Boltzmann Method (LBM)是在过去的20多年里逐渐发展起来的一种新的不同于传统计算流体力学(CFD)的计算方法 [ 1 ] 。相比于采用有限差分、有限体积等传统的计算方法有很大的不同，LBM基于分子动理论，具有比较清楚的物理背景。传统的计算方法是在求解离散方程时直接对宏观的控制方程进行离散化处理，然后使用某种数值方法进行求解，最后得到宏观的各个物理量。而LBM从微观动力学的角度出发，把宏观物理量当作微观量统计平均的结果 [ 2 ] 。 LBM方法是由格子气自动机(Lattice Gas Automata， LGA)发展而来 [ 3 ] 。在LBM方法中，流体在介观层次上被流散成流体粒子的同时物理区域也被离散为一系列的粒子(微团)，时间被离散成一系列的时间步长。其中，描述流体粒子运动的方程成为Boltzmann方程或者相应的离散形式 [ 4 ] 。在LBM方法中，速度、密度等宏观物理量不作为变量，基本变量是细观上的粒子分布函数。LBM方法通过粒子分布函数的演变过程，从不同的角度建立起宏观与微观、连续与离散之间的联系，按照简化以后的动力学模型流体粒子被约束在有限的格线上运动，需要对这些相关粒子的相关特征值做平均才可以获得宏观层次的密度、速度等参数 [ 5 ] [ 6 ] 。 本文针对低雷诺数下的NACA0012翼型进行了数值模拟，通过数值模拟结果与相关文献参考数据的对比，分析并验证LBM方法在数值模拟低雷诺数下翼型绕流时是一种有效、精准的计算方法。
system: [SYSTEM]

user: 2. 数值方法和几何模型 2.1. 数值方法 LBM法将粒子速度离散为有限速度空间 ，N为离散后速度的种类。速度离散后的玻尔兹曼输运方程为 (1) 式中， 为t时刻r点在α方向上的粒子分布函数， 为其离散速度， 为碰撞算子， 为外力在α方向上的分力。 宏观物理量由下式确定： (2) (3) 为求解式(1)，需对其时间及空间进行离散，得到完全离散的玻尔兹曼方程，并对其进行滤波： (4) 求解(1)式中必须同时离散速度空间和物理空间，最大的困难在于速度 的处理，由于宏观量是微观量的积分形式，He [ 7 ] 用高斯积分，对速度空间离散，得到如下的方程(2D)： (5) 其中： (6) (7) (8) 因为是在2维空间且具有9个离散速度，因此将这种速度离散记为D2Q9。 其实，LBM最初起源于LGA (Lattice gas automata)。经过Higuera和Jimenez提出的一种将碰撞算子线性化的简化模型 [ 8 ] 以及Higuera等进一步提出的增加模型的数值稳定性的强化碰撞算子方法 [ 9 ] ，LBM去除了统计噪音并去处理碰撞算子指数复杂性的缺点。由于改进后的LBM方法使用的仍然是Fermi-Dirac平衡态分布函数，LGA的其他缺点仍然存在。为了进一步简化碰撞算子，Chen、Qian提出了用SRT (single relaxation time)过程代替碰撞项的格子BGK (LBGK)模型 [ 10 ] [ 11 ] [ 12 ] ，相比之前提出的模型，碰撞算子不会因粒子种类的增加而变得复杂，主要是由于LBGK模型中碰撞算子的本身会发生改变。 通过对速度空间进行离散，宏观量为： (9) 通过多尺度展开，假定密度变化很小，可由(5)和(9)得到Navier-Stokes方程： (10) (11) 这里粘性系数 ，压力 。 经过上述处理后，LBM方法完全克服了格子气自动机的一系列缺点。 边界处理是LBM方法进行运用中一项非常重要的内容，它对数值计算的精度、计算的稳定性及其计算的效率有很大的影响。通过LBM方法进行数模拟计算时，由于边界节点上的部分分布函数是未知的，需要根据已知的宏观边界条件确定边界节点的分布函数，因此处理不同边边界条件时需要不同的边界处理格式。目前，将现有的边界处理格式主要分为启发式格式、动力学格式、外推格式以及复杂边界处理格式等四类。相对于其他三种格式，启发式格式不需要较复杂的数学推导的公式求解，一般只要根据边界条件的物理本质通过微观粒子的运动规则直接确定边界节点上的未知分布函数；动力学格式主要是通过边界宏观物理量或者宏观物理量梯度值的定义，直接求解边界节点上未知分布函数的方程组： ， (12) 以获得边界节点上待定的分布函数。虽然启发式和动力学格式可以处理一些特定的边界条件，但是存在很大的局限性，并且在过程中需要对边界条件的一些物理性质做出假设，很难推广到更一般的边界上。因此，借鉴传统计算流体力学(CFD)方法中的边界处理方法来构造LBM方法边界条件处理格式提出了外推式格式。复杂格式综合了其他三种边界处理格式，主要用于处理曲线边界和运动边界。四类处理格式具有不同的使用范围，其对模拟结果的精度以及稳定影响也差别很大，因此在模拟计算时需要根据具体的物理问题选择合适的边界处理格式。随着LBM方法的不断发展以及解决问题的需要，最近几年多位学者提出了用以处理复杂边界条件的不同格式，如连续反弹格式(continuous bounce-back scheme) [ 13 ] 、体积法(volumetric method) [ 14 ] [ 15 ] 以及浸润边界(immersed boundary) [ 16 ] [ 17 ] [ 18 ] 。 网格划分方面，非均匀网格、多重网格、无网格技术等网格处理技术的出现延伸了LBM方法的应用范围，使得LBM方法逐渐在更多的领域得到实际的应用 [ 4 ] 。 2.2. 几何模型 本文选取了NACA0012翼型作为研究对象，它作为一款常用翼型研究数据较多，因此可以方便的进行对比研究。本文算例计算区域大小为60 c × 40 c，计算区域外围为远场边界条件。为了提高数值计算的精度及其可靠性，本文LBM方法所采用的网格生成方法为无网格生成技术。格点在物理量变化较为剧烈的部分分布较密，而变化平缓的部分分布比较稀疏。远场区域解析度为2.56 m，翼型及其尾迹解析度为0.005 m，如图1所示，仅在交界处使用插值的方法，不仅节省计算资源并且保证了计算的精度 [ 19 ] 。
system: [SYSTEM]

user: 3. 数值模拟 本文用LBM方法对雷诺数为 下NACA0012翼型进行数值模拟，将结果与参考文献 [ 20 ] [ 21 ] 所提供参考值进行对比分析，说明LBM方法在计算低雷诺数下翼型绕流的有效性。为了达到基于弦长等于500的雷诺数，根据表1设定了模拟参数。由于来流为层流，在翼型边界层不应用壁面函数模型。 图2为NACA0012翼型在一流场下不同解析度对于精度的影响的仿真结果。通过不同解析度下LBM法计算得到的阻力系数与CFL3D参考值的对比，结果表明当解析度为0.005 m时，计算结果与文献参考值可以很好的吻合。在翼型及尾迹解析度为0.005 m时，其升力系数、阻力系数如表2所示。 将LBM方法数值模拟流场中五个垂直于翼型截面取无量纲化的X和Y方向速度分量以及压力系数与CFL3D数据进行对比，如图3所示在弦长( )的五个不同位置。 图1. NACA0012翼型格点分布图 图2. 解析度对阻力系数的影响 图3. 数据点线 图4为LBM法和CFL3D在X与Y方向无量纲化速度分布对比。通过分析我们可以发现：LBM计算所得X方向上的速度分量与CFL3D参考值十分吻合。但LBM法在X方向上的分量与CFL3D参考值在靠近翼型壁面时存在一些偏差，特别是在弦长位置 。这可能是由于计算过程中靠近壁面的粒子层大小的设置不够精准。 压力系数 被定义为 ，其中 为静压。图5为LBM法和CFL3D在翼型弦长五个不 同位置压力系数分布对比。通过对比分析可以得知，在计算压力系数时LBM法计算所得结果相对于CFL3D的参考值总体趋势比较吻合，但是在弦长 与 位置靠近壁面处的存在一定的差异。造成壁面处速度的Y方向差异的原因与速度X方向存在差异的原因一样，主要是靠近壁面附近的粒子层精度的设置问题。 图4. 翼型不同位置X与Y方向的速度分量。(a) 处X方向速度分量；(b) 处Y方向速度分量；(c) 处X方向速度分量；(d) 处Y方向速度分量；(e) 处X方向速度分量；(f) 处Y方向速度分量；(g) 处X方向速度分量；(h) 处Y方向速度分量；(i) 处X方向速度分量；(j) 处Y方向速度分量 Table 1 Simulation paramete 自由来流速度 50 m∙s −1 密度 ρ 1 kg∙m − 3 动态粘度 μ 0.1 pa∙s 弦长 c 1 m 雷诺数 Re 500 攻角 α 0度 表1. 模拟参数 Table 2 NACA0012 lift coefficient and drag coefficien 方法 升力系数 阻力系数 LBM 0.1734 CFL3D 0.1741 理想结果 0 0.1741 误差 - 0.402% 表2. NACA0012升力系数、阻力系数 综合上述低雷诺数下的NACA0012翼型绕流数值模拟结果的对比分析，可以发现LBM方在对简单对称翼型进行数值模拟具有很高的精准度，模拟结果与文献参考值可以很好地吻合。同时，这也验证了LBM方法在计算低雷诺数下简单对称翼型是一种有效的计算方法。在计算的过程中，对同一流场采用分区域的划分，可以提高区域格点的分布质量，从而在一定程度上能够提高计算的精度。但是，LBM方法在处理翼型边界时采用阶梯近似的方法会对翼型表面物理量的准确性计算产生一定的误差。 图5. 翼型不同位置的压力系数。(a) ; (b) ; (c) ; (d) ; (e)
system: [SYSTEM]

user: 4. 结论 Lattice Boltzmann Method (LBM)是最近发展起来的一种新的计算方法，其所具物理图像清晰、边界条件容易处理以并行性能好的特点，使LBM 在一些传统方法研究受到一定限制的领域取得了开拓性的进展。本文采用LBM方法对低雷诺数下的NACA0012翼型进行了模拟，结果表明LBM方法在计算低雷诺数下的简单对称翼型是一种有效的计算方法。未来的工作将以本文为基础，对高雷诺数、复杂边界条件下的翼型进行数值研究，说明LBM方法在更复杂条件下的应用。
system: [SYSTEM]

user: 摘 要 新疆油田冲积扇砂砾岩储层岩石类型多样，颗粒大小混杂，分选极差，储层孔隙结构具有典型的多种模态特征，即单模态、双模态与复模态等。利用铸体薄片、压汞、恒速压汞资料定量表征了七东1区砂砾岩储层复杂模态孔隙特征以及连通特征，利用多元回归分析法建立复杂模态孔隙结构储层的识别模型。以此为基础，分析了冲积扇砂砾岩储层孔隙结构复杂模态成因机制。研究表明：单模态储层主要以Ⅰ类储层为主，双模态储层以Ⅱ、Ⅲ类储层为主，复模态储层以Ⅲ类储层为主。从单模态类型到复模态孔喉组合类型，孔隙发育由好变差，由连通较好的网状到连通较差的星点状，面孔率降低，平均孔径减小，孔喉配位数减小。利用已有147块样品统计发现，七东1区下克拉玛依组储层中单模态储层占19%、双模态储层占48%、复模态储层占33%。由此提出不同模态判断标准，建立了储层孔隙结构预测模型，揭示了储层构型、成岩后作用及注水(聚)开发中后期的系列变化等多因素是导致形成复杂模态的原因。 关键词 :冲积扇，砂砾岩，复杂模态，孔隙结构，定量表征 Copyright © 2017 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 ] [ 6 ] [ 7 ] 、碳酸岩岩储层、低渗致密储层、火山岩储层、页岩及煤层等类型 [ 8 ] [ 9 ] [ 10 ] ，其复杂性主要体现在沉积机制的多样而引起的多种岩石相类型，至少包括砾岩相、砂砾岩相、小砾岩相、粗砂岩相、中–细砂岩相、粉砂岩相、泥质砂岩相、泥岩相等8类 [ 11 ] ，不同的岩石相类型其颗粒分选、颗粒含量、排列方式、杂基含量、胶结物胶结方式及含量具有明显差异性，因此各类岩石相还可以进一步细分亚类，学者们 [ 12 ] [ 13 ] [ 14 ] 认识到这些是导致孔隙结构模态(包括孔隙、喉道类型及其连通关系等方面)差异的关键因素。 学者早期对储层的研究发现，宏观尺度的沉积相分布和后期埋深、压力、温度等因素控制了成岩作用，进而控制了储层质量的差异 [ 15 ] ；而微观尺度的颗粒分选、颗粒含量、排列方式，杂基含量及胶结方式等影响了孔隙结构 [ 16 ] 。对于碎屑岩中分选差、粒度粗、粒度分布范围大的砂砾岩储层来说，由于砾岩储层孔隙空间容易被细的砂岩填充，基于此学者 [ 12 ] 提出了砂砾岩双模态特征。之后，有学者 [ 13 ] [ 14 ] 研究发现在砂岩、砾岩、泥岩3类岩性同时存在的砂砾岩储层中，受沉积相类型及成岩作用方式的控制，其储层孔隙结构具有复模态特征。然而，另有学者 [ 17 ] [ 18 ] 认为，储层复模态的形成不仅与岩石粒度分选有关，同时还与粒径比和颗粒相对含量相关，并建立了碎屑岩储层孔隙结构模态的定量数学模型；同时，提出了砂岩储层也可以形成复模态的观点，但并未进一步深入揭示砂岩储层孔隙复模态的成因。 新疆油田砂砾岩储层属于典型冲积扇沉积，岩石颗粒快速沉积，颗粒分选差，大小混杂，岩石颗粒结构具有典型的复杂模态特征 [ 19 ] 。几十年的勘探开发及试验研究表明，不同的储层孔隙结构模态对开发方式的适应性具有显著差异，是制约砾岩油气藏高效开发的关键因素。因此，研究储层孔隙结构复杂模态对提高砾岩油藏采收率具有重要的实践意义，同时，对深化砂砾岩储层地质学具有重要的理论意义。以新疆油田七东1区砂砾岩为对象，利用铸体薄片、压汞与恒速压汞等方法，对砂砾岩孔隙结构进行定量表征，拟建立砂砾岩模态判断方程和物性预测模型，并深入探讨砂砾岩复杂模态孔隙结构的成因。
system: [SYSTEM]

user: 2. 研究区地质概况 2.1. 研究区位置及油藏特征 七东1区位于新疆克拉玛依市白碱滩地区，距克拉玛依市以东约30 km，地表平坦，交通便利。三叠系下克拉玛依组(简称克下组)油藏含油面积7.0 km 2 ，地质储量1010 × 10 4 t，油藏埋深1100~1400 m。构造上处于准噶尔盆地西北缘克–乌逆掩断裂带白碱滩段的下盘，是一个四周被断裂切割成似菱形的封闭断块油藏。 2.2. 地层与沉积储层特征 研究区七东1区钻遇了第四系、新近系、古近系、白垩系、侏罗系、三叠系、石炭系。目的层克下组不整合沉积在石炭系之上，克下组划分为2个砂层组，细分为12个单层(小层) (图1)。 研究区砾岩类占岩石总类的64.82%，砂岩类占35.18%。其中，含砾粗砂岩占32.21%，砂砾岩占23.72%，小砾岩占21.15%，不等粒砾岩占18.58%，粗砂岩占2.96%，细砾岩占1.38%。储层岩石中砾石成分复杂，以花岗岩、凝灰岩为主，砂岩颗粒主要为石英、长石和岩屑；碎屑颗粒分选中等～差，磨圆度差，以次棱角状～次圆状为主；岩石碎屑为颗粒支撑，骨架颗粒以点接触为主。储层矿物以石英、钾长石为主，其次为斜长石，少量碳酸盐矿物。黏土矿物(粒径 < 8 µm)绝对体积分数为0.86%~9.23% (平均3.62%)。岩石黏土含量低，与油田长期注水开发，部分黏土微粒被运移有关。主要黏土矿物为高岭石，平均体积分数为67.09%，其次为伊利石、伊–蒙间层等。整体上成分成熟度和结构成熟度低，反映了近源短距离搬运和快速堆积的沉积特征。泥岩以深灰色和棕红色为主，反映了干旱氧化的沉积环境。储层非均质性极其严重，其储层内部孔隙结构模态多样。 图1. 研究区目的层柱状图(ES7008井)
system: [SYSTEM]

user: 摘 要 本文提出分布意义下的股票价格动量并研究其在证券投资中的应用。首先假设每支股票为一个多粒子系统，引入股票交易筹码分布，在实证检验中利用正态分布对其进行估计；其次，本文给出一个迭代算法用来估计股票的成本分布；最后，基于成本分布定义新股票价格动量，并利用动量/反转策略对其进行实证分析，结果表明：新动量在很多方面优于其它动量。 关键词 :股票交易筹码分布，成本分布，价格动量，动量/反转策略 Copyright © 2017 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. 引言 近年来，经济物理吸引了大量专家学者的关注。Jegadeesh和Titman [ 1 ] 于1993年研究股票市场的动量效应，他们通过使用基于收益率的动量/反转策略构造投资组合，最终得到了显著异常收益。此后，国内外很多学者将该策略应用到不同股票市场以检验是否存在套利机会。王永宏和赵学军 [ 2 ] 对深沪两市1993年以前上市的所有股票进行研究，实证结果表明深沪股票市场存在明显的反转效应；朱战宇、吴冲锋和王承炜 [ 3 ] 使用1995年至2001年交易数据(该时间段经历了牛市和熊市)，利用重叠抽样价格动量检验方法研究中国股市，结果表明：在以月为周期的检验中不存在显著动量效应，而使用以周为周期的动量策略可获得超额收益；杨德勇和王家庆 [ 4 ] 通过研究发现：在我国A股市场中存在显著的短期反转效应，而对于中长期，A股股票具有显著的动量效应；冯科和郑琛 [ 5 ] 证实中小板市场存在短期反转效应，而且动态投资策略能够实现更为可观的收益，但收益波动也更为剧烈；许强和顾海峰 [ 6 ] 发现我国创业板市场存在短期反转效应，中期动量效应以及长期反转效应。 Choi [ 7 ] 从经典物理的角度出发，定义了10种股票动量。他选择对数价格作为股票的位置，持有期收益率 和对数收益率 作为速度，交易量，每天交易金额和波动率倒数作为质量。他将经典物理动量嵌入到动量/反转策略 [ 1 ] 中(代替原有的持有期收益)来研究美国和韩国股市的均衡性，最终得到如下结果：大部分经典动量策略是有效的，且优于传统动量策略。但动量 在美国股市中表现最好，在韩国股市中表现最差。这个差异表明也许存在更科学更符合实际的股票动量定义。如果我们把每支股票的每一股份看作一个粒子，那么每支股票可被视为热力学或统计物理中的一个多粒子系统 [ 8 ] ，而估计所有投资者的成本就成为一个重要任务。本文提出股票交易筹码分布和成本分布，并基于这两个分布定义新股票价格动量，然后利用反转策略将新定义和其他定义作对比，发现新股票价格动量对应的投资组合可以带来更高的收益回报率和夏普比率，因而本文构建出更好的股票动量。
system: [SYSTEM]

user: 2. 股票交易筹码分布和成本分布 正如引言所述，每支股票可以被当做一个系统，它有许多“粒子”(股份)组成。在实际交易中，每一笔交易仅涉及部分“粒子”。很明显，交易的部分股份不能代表该股票所有股份，但是如果想了解股票系统的准确状态，我们需要所有股份的信息。因此，本文提出股票交易筹码分布和成本分布来推断股票的状态。本文把每支股票的所有股份当做一个集合，将其看作一个样本空间 ，每一股份看作空间中的一个样本点 。给定交易时间区间 ，交易股份作为 的一个子集，用 表示。 注：由于股票价格服从对数正态分布，所以对数价格服从正态分布 [ 7 ] 。下面所提及的价格都是指对数价格。 在时间段 中的交易价格 是空间 上的随机变量，称它的分布为股票交易筹码分布，该分布表示在时间段 内以价格 交易的股份比例。而不同的持有者在时刻 有不同的成本，因为他们有不同的购买价格和持有时间。相似地，如果把成本 作为空间 上的随机变量，则把它的分布称为成本分布。该分布表示在某一时刻 ，成本是 的股份比例。具体表示如下： 令 表示在时间段 内股票的交易价格， 表示在此期间与价格 对应的交易量， 表示此期间内总的交易量，即 。因此，在 期间的股票交易筹码分布为 。 本文以深圳万科(000002)2011年1月6号的数据为例来计算它的交易筹码分布。当天总交易量 是77189390股，价格 ，其对应的交易量 为2052074股。因此，股票交易筹码分布在该价格上的值为 。 与股票交易筹码分布类似，令 表示股票的成本，令 表示在时刻 ，成本为 的股份总数。 假设 为流通股总数，则该股票的成本分布在时刻 ，价格 处的值为 。
system: [SYSTEM]

user: 3. 股票价格动量 令 表示某支股票在时刻 的成本分布，令 表示在时间段 内股票的交易状态，它被定义为 (1) 因此， ， ，其中 表示流通股总数， 表示在该时间段内的交易量。 在中国股市中，交易规则是 规则，即当天买入的股票当天不能卖出。 如果令 表示 时的成本，则 (2) 其中 表示交易价格，即股票交易后成本 为 。 股票速度为 (3) 由于投资者不知道买家和卖家，即 和 的关系未知，为了利用 和 推导出 ，本文假设随机变量 和 独立，即每一股份等可能被交易。 令 且 ，则可以得到 (4) 股票 的动量 定义为 (5) 其中质量 为总流通股数 ， ， ，则每股平均动量可表示为 (6) 现在本文已经给每支股票的每一股份定义了成本，速度以及动量，它们都是样本空间 上的随机变量，并且都有各自的概率分布，因此新的动量具有分布意义。
system: [SYSTEM]

user: 4. 迭代算法 通常来说，投资者只能得到每天的交易数据，如开盘价，收盘价，最高价，最低价和交易量，现在的问题就是如何根据以上几种数据得到股票交易筹码分布和成本分布。 正态分布是概率论中最重要的分布之一，也是自然和社会科学中常被用来估计实值随机变量分布最常用分布之一。如果已知一随机变量的均值和方差，现要推测其分布，根据最大熵原理 [ 9 ] ，可知其服从正态分布。 因此，正态分布是理想选择。本文利用正态分布来推断股票交易筹码分布和成本分布，下面给出计算动量的迭代算法： (1) 估计股票交易筹码分布 令 ， ， 和 分别表示一交易日价格的最大值，最小值，开盘价和收盘价。平均价格由下式估计： (7) 对于正态分布，它有一重要性质： 原则。该性质表明正态分布随机变量 的取值落在均值左右3倍标准差范围内的概率接近于1，即 ，因此本文依此给 一个估计： (8) 估计出 和 ，就可以通过正态分布 得到股票交易筹码分布，如下图1。 在图1中，柱状图表示由深圳万科2011年1月6日实际交易数据得到的股票交易筹码分布，曲线表示其概率密度函数 ，其参数由式(7)和(8)估计。 (2) 估计成本分布 令 表示时刻 对应的股票交易筹码分布， 代表时刻 对应的成本分布，则时刻 对应的成本分布由下式估计： (9) 其中 代表时刻 对应的股票换手率。 图1. 深圳万科(000002)的股票交易筹码分布 在图2中，柱状图表示由深圳万科实际交易数据得到的2011年3月29日成本分布，曲线表示其成本分布的概率密度函数。 (3) 计算动量 利用式(6)计算股票动量 。 使用上述迭代算法可根据实际数据估计股票交易筹码分布和成本分布，计算股票动量。 至此，基于以上讨论，我们可以看出新动量定义不同于之前的，因为之前动量定义是依据经典物理，而新动量具有分布意义，所以新定义可以被认为是对先前的扩展。
system: [SYSTEM]

user: 5. 实证检验 5.1. 动量/反转策略 为了研究哪种动量定义在股票市场中更有效，Choi [ 7 ] 给出了基于经典物理的股票动量定义。他以两种收益率 和 定义股票速度 ，又分别以换手率 ，每日交易金额与市场资本比值 和波动倒 数 作为股票质量 。他还定义三种动量计算方式，即 ， ， ，并利用文献 [ 1 ] 中的动量策略比较11种动量定义，最终得到如下结果：和传统动量(收益率) 相比，几乎所有动量定义都是有效的，并且在韩国股市中，动量 的效果优于其他动量，而在美国股市中，动量 表现最好。动量策略描述如下：首先计算所有股票在 至 期间的动量，然后将动量从小到大排序，并将其分10组(分组数目也可根据具体情况稍作改变，如股票总数较少，可选择将 图2. 深圳万科(000002)的成本分布 其分为5组 [ 10 ] )，称第一组为loser，称最后一组为winner。对于winner和loser中的股票，按照每支股票购买相同数量的原则，买入winner，卖出loser来构建投资组合。不妨设每支股票购买数量为1股，并持有该投资组合 天。在持有到期时卖出winner，买入loser，反转策略与之相反。 本文利用动量/反转策略对12种动量进行实证检验，以比较新旧定义的优劣。如果新动量策略对应的收益优于包括传统策略在内的其他策略，则很显然，基于股票交易筹码分布和成本分布的新动量策略在寻找套利机会方面有很大优势。 5.2. 实验步骤 本文使用MATLAB软件编写程序进行实证检验，具体步骤如下： 第一步：下载股票交易数据并进行简单整理； 第二步：计算股票交易筹码分布和成本分布； 第三步：计算股票动量并排序，构建投资组合； 第四步：计算持有期收益率及风险指标。 5.3. 数据集 本文所用数据均来自国泰安和锐思数据库。选取上海证券交易所上市的203支股票，时间段选为2001年1月至2010年12月。由于起始成本分布未知，首先利用2001年数据得到股票交易筹码分布，然后使用迭代算法得到初始成本分布，并将2002年1月4日作为起始日。本文使用反转策略，取 ， ，即选择6周作为股票回看期，1周作为股票持有期。由于数据区间较长，故采用无重叠投资策略。 5.4. 实证结果 附录结果显示：所有基于物理动量的反转策略都优于传统反转策略，其中 是最好的，它的平均收益为0.93%，波动是3.17%，而传统动量策略的平均收益为0.13%，标准差为2.01%。新策略 的收益是 的7倍。策略 的收益为0.72%，优于 和 的，略差于 。 以上所有物理动量投资组合的表现都是源于较强的反转现象，即反转策略对应的投资组合收益均为正值。首先，所有物理动量策略中loser的收益均高于传统策略的。例如策略 中loser的平均收益为0.19%，在所有投资组合中是最高的；winner的收益为-0.73%。此外，策略 中loser的波动小于传统动量策略。而所有物理动量策略对应的winner收益较接近，均差于 的。 在所有动量策略中，策略 的夏普比率是最大的，几乎是传统策略 的6倍，对于风险指标 、 和 ， 的表现比 好，略差于策略 ， 和 。 的夏普比率高于 和 ，但是低于 。对于指标 、 和 ， 和 ， 及 相比，均表现出较低风险。 每个组(winner和loser)对应的风险和分组的目的是一致的。策略 的loser有最高的夏普比率和较低的风险指标，而winner有更高的风险。策略 中loser对应的夏普比率大于 ， 和 的。对于指标 和 ， 对应的loser (winner)小于(大于) ， 和 的。 各策略投资组合的累计收益率在图3中，其中 (绿色)， (蓝色)， (红色)， (紫色)， (浅蓝色)。
system: [SYSTEM]

user: 6. 结论 本文首先提出股票交易筹码分布和成本分布，其中股票交易筹码分布用来表示既定时间区间内股票的交易状态，成本分布用来表示某一时刻股票的状态，然后又提出一个估计成本分布的迭代算法。基于 图3. 反转策略对应的累计收益率 上述两个分布，得到分布意义下的股票价格动量。经过在中国股票市场的实证检验， 结果表明中国股市存在强的反转效应，并且新动量/反转策略的效果明显优于其他动量策略。
system: [SYSTEM]

user: 稳定碳、氧同位素地球化学特征与 成岩环境分析
system: [SYSTEM]

user: 摘 要 四川盆地东北部上二叠统长兴组白云岩为重要的优质储层，近年来，长兴组白云岩成岩环境或成因问题成为研究者讨论的热点话题。在岩石学研究的基础上，运用稳定碳、氧同位素，对研究区长兴组白云岩进行综合研究表明，白云石化具有多期次的特点，按照成岩环境划分为蒸发环境下的准同生白云岩和埋藏条件下的早、中、晚期成岩白云岩。并划分为早、中、晚成岩阶段，结合碳、氧同位素分析结果进一步表明晚成岩阶段白云岩是构造热液白云石化的结果。 关键词 :碳、氧同位素，白云岩，上二叠统长兴组，白云石化 Copyright © 2017 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 ] ，其余各类白云岩，如颗粒白云岩、晶粒白云岩和生物礁白云岩皆为多期次埋藏白云石化作用的产物。以前人研究为基础，笔者拟从地质综合研究角度出发，以岩石学为基础，结合研究区构造和沉积相特征，主要是利用稳定碳、氧同位素等地球化学测试数据，试对研究区各种类型白云岩进行全面的成岩环境分析，增加对川东北长兴组白云岩形成环境的认识，希望能够为盆地白云岩储层预测提供地球化学上的依据。
system: [SYSTEM]

user: 2. 区域地质背景 四川盆地已经成为我国重要的天然气富集盆地之一，盆地东北部上广泛发育的碳酸盐岩地层具有巨大的勘探潜力。四川盆地经过复杂的构造演化，不同阶段、不同性质的构造叠加，形成了四川盆地东部具有构造叠加性质的川东弧形褶皱带(图1)。中二叠世发生的“峨眉地裂运动”，使川东地区形成了张性断块，其原有的古构造基底格局被彻底改变。随后，中二叠世茅口期和龙潭期之间的东吴运动，使四川盆地整体抬升，在川东北地区形成了一个北东向的古斜坡。构造格局的改变，使西部的峨嵋山和康滇古陆一带为陆源剥蚀区。自龙潭期，沉积古地理格局自西向东依次为海陆交互相、浅海相和深海相。至长兴期，川东地区便以广阔的碳酸盐台地沉积为主，并伴有深海相“台盆”，即“开江–梁平”海槽和“城口–鄂西”海槽 [ 4 ] 。研究区即位于开江–梁平海槽东段(图2)，在区域构造上，开江–梁平海槽位于川东弧形褶皱带的东北缘，川东弧形褶皱带的空间展布与形成演化对研究区的构造演化具有极其重要影响，进而控制了该区的沉积 [ 5 ] 。研究区上二叠统长兴组主要为一套海相碳酸盐岩沉积，其中白云岩构成主要的储集岩类。 图1. 研究区构造位置图 图2. 开江–梁平海槽东段长兴组构造–古地理格局示意图
system: [SYSTEM]

user: 3. 白云岩岩石学特征 四川盆地东北部上二叠统长兴组白云岩分布较为广泛，岩石类型复杂多样。笔者在研究过程中，将原岩沉积特征和白云岩结构特征以及所具有的特殊构造为主要依据，将研究区长兴组白云岩分为泥–微晶白云岩、粉晶白云岩、细晶白云岩、中-粗晶白云岩和构造碎裂化白云岩等5种类型，各类型白云岩岩石学特征如下。 3.1. 泥–微晶白云岩 泥-微晶白云岩以具有泥-微晶结构为主要特征，晶体形态以他形为主，晶体之间接触紧密呈镶嵌状(图3(a))。常具有泥质条带、鸟眼和石膏假晶等潮汐和暴露成因标志。该类白云岩一般认为与准同生白云石化作用有关，即蒸发条件下的潟湖–潮坪环境中，高Mg/Ca比的卤水交代灰泥形成。 3.2. 粉晶白云岩 该类白云岩的原始沉积构造保存较为完整，有两种主要岩石类型：灰质粉晶颗粒或礁白云岩(图3(b))和粉-细晶颗粒或礁白云岩(图3(c))；颗粒组分以生物碎屑为主，前者晶体形态以半自形为主，少数他形，粉晶结构，白云石的晶面较脏，颗粒和礁结构皆保存较为完好；后者在晶体大小和形态上与前者较为一致，颗粒或生物礁结构少量破坏，且几乎不含残余灰质组分。该类白云岩应与生物碎屑或造礁生物被优先白云石化有关。根据其结构特征，表明该类白云岩为埋藏阶段早期白云石化作用的产物。 3.3. 细晶白云岩 该类白云岩晶体大小在0.05~0.25 mm之间，部分为中晶级(0.25~0.5 mm)，白云石晶面较脏，晶形较好，且多见雾心亮边结构，白云石化不均匀并伴有强烈的重结晶作用。在该类白云岩的某些岩石类型中，如具不明显残余结构的粉–细晶颗粒白云岩(图3(d))或礁白云岩或部分中–细晶白云岩(图3(e))，礁白云岩中的晶间孔和晶间溶孔十分发育，在较大的溶孔内往往有白云石晶簇分布，除局部充填有少量沥青外，一般不含其他的外来物质，此为深部溶解作用的典型标志，并确定为该类白云岩为埋藏阶段中期白云石化作用形成。 3.4. 中–粗晶白云岩 该类白云岩中晶体大小以中晶级(0.25~0.5 mm)为主，部分为粗晶级(0.5~1.0 mm)，晶形较好，部分白云石具有雾心亮边结构(图3(f)、图3(g))。值得注意的是，在粗晶白云石中常见晶体呈马鞍状结构特征。与之相应的是研究区长兴组分布在台内海槽的东段，表明构造活动相伴生的热液流体对长兴组白云岩具有一定影响，该类具马鞍状结构的白云石属于埋藏阶段晚期构造热液活动成因 [ 6 ] [ 7 ] [ 8 ] ，下文会利用地球化学数据进一步分析论证。 3.5. 构造碎裂化白云岩 该类白云岩以发生构造碎裂化作用为特征(图3(h)、图3(i))，伴随有强烈的溶蚀现象，发育呈网状分布的裂缝，导致岩石的原始结构完全消失，溶蚀作用沿裂缝表现较为强烈，孔、洞、缝连通性好并在其内充填有炭质沥青。该类白云岩一般认为是成岩阶段晚期白云岩发生碎裂化作用所导致的。 图3. 研究区长兴组常见的白云岩类型图谱