近地层风沙颗粒跃移运动学：机理、特征与应用的深度剖析

近地层风沙颗粒跃移运动学：机理、特征与应用的深度剖析一、引言1.1研究背景与意义风沙问题是一个全球性的环境难题，对生态系统、人类生活和经济发展都产生了深远的影响。据统计，全球约有40%的陆地面积受到风沙活动的威胁，而中国也是受风沙危害较为严重的国家之一，约1/3的国土面积受到风沙的影响。风沙灾害不仅导致土地沙漠化、土壤侵蚀和植被退化，还会引发沙尘暴等极端天气事件，对农业、畜牧业、交通运输、能源等行业造成巨大损失。例如，2000年春季，中国北方地区连续发生了12次沙尘暴，波及范围包括西北、华北和东北等地，造成了严重的经济损失和生态破坏，据估算，直接经济损失超过540亿元。此外，风沙活动还会对人体健康造成危害，沙尘暴中的沙尘颗粒可携带大量的细菌、病毒和有害物质，被人体吸入后可能引发呼吸道疾病、心血管疾病等。近地层风沙颗粒跃移运动是风沙活动的基本过程，也是风沙物理学研究的核心内容之一。在风力作用下，沙粒从地表起动并以跳跃的方式在近地层空气中运动，这一过程涉及到复杂的空气动力学、颗粒力学和多相流相互作用等问题。了解近地层风沙颗粒跃移运动的规律，对于深入理解风沙现象的形成机制、预测风沙灾害的发生发展以及制定有效的防治措施具有重要的科学意义和实际应用价值。从科学研究的角度来看，近地层风沙颗粒跃移运动学的研究有助于揭示风沙物理过程的本质，丰富和完善风沙物理学的理论体系。通过对风沙颗粒跃移运动的观测和实验研究，可以获取沙粒的起动条件、运动轨迹、速度、加速度等关键参数，为建立准确的风沙运动模型提供数据支持。同时，研究风沙颗粒与气流之间的相互作用机制，以及颗粒间的碰撞、摩擦等过程，有助于深入理解风沙流的动力学特性和能量传输规律。此外，近地层风沙颗粒跃移运动还与土壤风蚀、风沙地貌的形成演化等密切相关，对这些领域的研究也具有重要的推动作用。在实际应用方面，研究近地层风沙颗粒跃移运动学对于风沙灾害的防治具有重要的指导意义。风沙灾害的防治需要深入了解风沙运动的规律，以便采取针对性的措施。例如，通过研究风沙颗粒的跃移运动，可以确定风沙流的输沙能力和危害范围，为风沙防护工程的设计和建设提供科学依据。在沙漠地区的铁路、公路建设中，需要考虑风沙对路基和轨道的侵蚀危害，通过对风沙颗粒跃移运动的研究，可以优化防护措施，减少风沙对交通设施的破坏。此外，研究近地层风沙颗粒跃移运动还可以为农业生产、生态环境建设等提供参考。在干旱半干旱地区，风沙活动会导致土壤肥力下降、农作物减产，通过了解风沙颗粒的运动规律，可以采取合理的农业措施，如种植防风固沙植物、设置沙障等，减少风沙对农田的危害，保护生态环境。综上所述，近地层风沙颗粒跃移运动学的研究具有重要的理论和实际意义。通过深入研究这一领域，可以为风沙灾害的防治、生态环境的保护和可持续发展提供科学依据和技术支持，对于改善人类生存环境、促进经济社会的发展具有重要的作用。1.2国内外研究现状近地层风沙颗粒跃移运动学的研究历史悠久，国内外学者在这一领域开展了大量的研究工作，取得了丰硕的成果。国外方面，早期的研究主要集中在风沙颗粒的起动条件和运动轨迹的初步探索。Bagnold在1941年发表的《ThePhysicsofBlownSandandDesertDunes》一书中，对风沙运动进行了系统的研究，提出了许多经典的理论和观点，如沙粒的起动风速公式、跃移运动的基本模型等，为后续的研究奠定了坚实的基础。此后，众多学者在此基础上不断深入研究。White等研究了马格努斯力对沙粒跃移轨迹的影响，发现马格努斯力会使沙粒的运动轨迹发生偏移。Ungar等模拟了沙粒与空气之间的耦合作用，得出了跃移达到稳定之后的风速廓线和输沙通量分布，进一步揭示了风沙流内部的动力学机制。随着科技的不断进步，实验技术和数值模拟方法在风沙颗粒跃移运动研究中得到了广泛应用。风洞实验可以精确控制实验条件，获取详细的风沙运动数据。高速摄像机、激光粒度仪等先进设备的应用，使得对沙粒运动参数的测量更加准确和细致。在数值模拟方面，计算流体力学（CFD）方法以及离散单元法（DEM）等被用于模拟风沙颗粒的运动过程，能够对复杂的多相流相互作用进行深入分析。国内在近地层风沙颗粒跃移运动学的研究起步相对较晚，但发展迅速。朱震达等早在20世纪60年代就对塔克拉玛干沙漠沙漠南缘策勒的沙丘变化进行了研究，并于80年代出版了《塔克拉玛干风沙地貌研究》一书，为我国风沙地貌研究提供了重要的基础资料。此后，我国学者在风沙颗粒起动、跃移运动特征、风沙流结构等方面开展了大量的研究工作。董治宝等对风沙颗粒的起动条件、受力分析及其影响因子进行了深入研究，提出了多种起动学说，如风压起动说、升力起动说、冲击起动说等。邹学勇等通过风洞试验与野外观测试验，研究了跃移输沙量随高度的分布规律，发现其通常遵循指数函数或幂函数分布，但受下垫面特点及气象条件的影响，也会出现偏离。李建国等研究了不同角度和不同高度防沙堤周围的流场分布，为风沙防护工程的设计提供了理论依据。近年来，我国学者还在风沙颗粒运动的数值模拟方面取得了重要进展，结合我国的实际情况，建立了一系列适合我国风沙环境的数值模型。尽管国内外在近地层风沙颗粒跃移运动学方面取得了众多成果，但仍存在一些不足之处。一方面，现有的研究大多基于理想条件，如均匀的沙粒粒径、平坦的下垫面等，而实际风沙环境复杂多变，沙粒粒径分布广泛，下垫面起伏不平且存在植被覆盖等情况，这些复杂因素对风沙颗粒跃移运动的影响尚未得到充分研究。另一方面，风沙颗粒与气流之间的相互作用机制，以及颗粒间的碰撞、摩擦等过程的研究还不够深入，目前的理论模型和数值模拟方法在描述这些复杂过程时仍存在一定的局限性。此外，不同研究方法和实验条件下得到的结果存在一定的差异，缺乏统一的标准和对比分析，这也给研究成果的应用和推广带来了困难。针对现有研究的不足，本文将综合运用实验观测、理论分析和数值模拟等方法，深入研究近地层风沙颗粒跃移运动学。在实验方面，将设计更加贴近实际风沙环境的实验方案，考虑沙粒粒径分布、下垫面特征和植被覆盖等因素的影响，获取更加准确和全面的实验数据。在理论分析方面，将进一步完善风沙颗粒运动的理论模型，深入研究风沙颗粒与气流之间的相互作用机制，以及颗粒间的碰撞、摩擦等过程，提高理论模型的准确性和适用性。在数值模拟方面，将开发更加先进的数值模拟方法，结合实验数据进行验证和优化，实现对近地层风沙颗粒跃移运动的高精度模拟。通过以上研究，旨在揭示近地层风沙颗粒跃移运动的规律，为风沙灾害的防治和生态环境的保护提供更加科学的依据和技术支持。1.3研究内容与方法1.3.1研究内容风沙颗粒跃移的基本原理：深入剖析风沙颗粒跃移的物理过程，对沙粒的起动条件展开研究，明确在不同风力、地表状况及沙粒特性下，沙粒从静止状态转变为运动状态所需的临界风速等条件。同时，对作用于风沙颗粒的各种力进行详细分析，如重力、风力、气动阻力、升力、马格努斯力、萨夫曼力等，研究这些力在沙粒运动过程中的相互作用及对沙粒运动状态的影响机制。此外，对风沙颗粒的运动形式进行深入研究，除了经典的跃移运动，还包括蠕移和悬移等运动形式，探讨它们在不同条件下的发生机制、相互转化关系以及各自对风沙输移的贡献。风沙颗粒跃移的运动特征：精确测量风沙颗粒跃移的运动轨迹，运用高速摄像机、粒子图像测速技术（PIV）等先进设备，获取不同条件下沙粒在空间中的运动路径信息，分析运动轨迹随时间的变化规律，以及与风速、风向、沙粒粒径等因素的关联。深入研究风沙颗粒的运动速度，包括起跳速度、上升速度、下降速度以及水平方向的运动速度等，分析速度在跃移过程中的变化情况，探究影响速度大小和方向的各种因素，如风力大小、沙粒质量、空气阻力等。此外，研究风沙颗粒的运动加速度，分析加速度在不同阶段的变化规律，以及对沙粒运动轨迹和速度的影响。风沙颗粒跃移的影响因素：研究风速和风向对风沙颗粒跃移的影响，通过实验和数值模拟，分析不同风速和风向条件下沙粒的起动、运动轨迹、速度和输沙率的变化规律，探讨风速和风向的波动对风沙颗粒跃移的影响机制。分析沙粒粒径和形状对其跃移运动的影响，不同粒径和形状的沙粒在风力作用下的受力情况和运动特性存在差异，研究这些差异对沙粒跃移的起动、轨迹、速度和输沙率的影响，为风沙运动的理论研究和数值模拟提供基础数据。此外，研究地表粗糙度和植被覆盖对风沙颗粒跃移的影响，地表粗糙度和植被覆盖会改变近地层气流的结构和特性，进而影响沙粒的起动和运动，分析它们对风沙颗粒跃移的抑制或促进作用，以及在不同地表条件下风沙颗粒跃移的变化规律。1.3.2研究方法实验观测：在实验室风洞中开展风沙颗粒跃移实验，风洞实验能够精确控制风速、风向、温度、湿度等环境参数，为研究风沙颗粒跃移提供了理想的实验条件。通过在风洞中铺设不同类型的沙床，模拟真实的风沙环境，利用高速摄像机、激光粒度仪、热膜风速仪等先进设备，对风沙颗粒的运动轨迹、速度、加速度、粒径分布等参数进行精确测量。同时，在野外风沙活动频繁的地区进行实地观测，获取实际风沙环境下的风沙颗粒跃移数据。在沙漠地区设置观测站点，安装风速仪、风向仪、集沙仪、颗粒物监测仪等设备，长期监测风沙颗粒的运动情况、输沙率、风沙流结构等参数，分析风沙颗粒跃移在自然条件下的变化规律，以及与气象条件、地形地貌、植被覆盖等因素的关系。数值模拟：运用计算流体力学（CFD）方法，模拟近地层气流的运动特性，建立三维的气流运动模型，考虑空气的粘性、湍流等因素，求解Navier-Stokes方程，得到近地层气流的速度场、压力场等信息。将风沙颗粒视为离散相，与连续相的气流进行耦合计算，考虑风沙颗粒与气流之间的相互作用力，如气动阻力、升力等，以及颗粒间的碰撞、摩擦等过程，采用离散单元法（DEM）或颗粒轨道模型（PTM）等方法，模拟风沙颗粒在气流中的运动轨迹、速度、加速度等参数。此外，结合实验数据对数值模拟结果进行验证和优化，通过对比实验测量值和数值模拟值，调整模型参数和计算方法，提高数值模拟的准确性和可靠性。理论分析：基于经典力学和空气动力学原理，建立风沙颗粒跃移的理论模型，考虑风沙颗粒的受力情况、运动方程以及与气流的相互作用，推导风沙颗粒的起动条件、运动轨迹方程、速度和加速度的计算公式等。运用数学方法对理论模型进行求解和分析，采用解析法或数值解法，得到风沙颗粒跃移的关键参数随时间和空间的变化规律，分析各种因素对风沙颗粒跃移的影响机制，为实验观测和数值模拟提供理论指导。此外，对不同的理论模型进行对比和分析，探讨它们的优缺点和适用范围，结合实际情况选择合适的理论模型，或者对现有模型进行改进和完善。二、近地层风沙颗粒跃移运动的基本原理2.1风沙颗粒的起动条件2.1.1起动学说概述风沙颗粒的起动是风沙运动的起始环节，多年来，众多学者针对这一现象提出了多种起动学说，每种学说都从特定角度对风沙颗粒的起动机制进行了阐释。风压起动说由拜格诺（Bagnold）和伊万诺夫（Ivanov）提出。拜格诺通过风洞实验认为，表面沙颗粒是在风压力直接作用下发生起动的，当风速达到一定程度，风对沙粒表面产生的压力足以克服沙粒与地表之间的摩擦力时，沙粒便开始运动。伊万诺夫则认为沙颗粒运动时，颗粒与接触表面的摩擦力和风压力不平衡时产生反力矩，从而导致沙粒起动。然而，该学说存在一定局限性，它难以解释一些粒径较小的沙粒在相对较低风速下的起动现象，且对于复杂地表条件下的沙粒起动情况，风压起动说的解释力相对薄弱。升力起动说指出，由于地表上沙土颗粒上下存在风速差，根据伯努利定律，颗粒上部风速大、压力小，颗粒下部速度小、压力大，从而使颗粒受到一个向上的压差作用力，即升力，当升力大于沙粒重力时，沙粒起动上升而离开地面。这一学说在解释沙粒起动时考虑了气流速度分布对沙粒的作用，但实际风沙环境中，沙粒形状不规则，气流的紊流特性也较为复杂，使得升力的计算和作用机制变得难以准确把握，这限制了该学说在实际应用中的准确性。冲击起动说认为，冲击力是沙粒起动的一个主要起动力。拜格诺根据风洞顶上供应时，可以降低沙面的起动流速实验得出此结论，伊万诺夫也认为沙粒脱离地表的主要升力是冲击力。但董治宝对该学说提出质疑，实际风沙运动中，冲击起动说难以全面解释不同粒径沙粒在各种环境下的起动情况，且对于冲击力的量化以及与其他作用力的协同关系研究尚不完善。除上述学说外，还有压差起动说、振动起动说、猝发起动说、滑移起动说、湍流起动说、斜面飞升说、负压起动说和涡旋起动说等。这些学说从不同方面，如气流的湍流特性、沙粒间的相互作用、地表的振动等，对风沙颗粒的起动进行了探讨，但它们也都存在各自的局限性，无法完全统一地解释风沙颗粒的起动现象。由于风沙起动过程同时受多种力的作用，且不同性质（如不同密度、粒径和形状等）沙粒的起动过程中各种力所起的作用也可能不同，因此，试图仅仅靠某一种学说或某一种力来解释风沙的起动是不够完善的。2.1.2受力分析风沙颗粒的起动和运动受到多种力的综合作用，这些力相互影响，共同决定了沙粒的运动状态。重力是作用在风沙颗粒上的基本力之一，其大小与沙粒质量和重力加速度有关，方向竖直向下。重力对沙粒的运动起着阻碍作用，尤其是在沙粒的上升阶段，重力会使沙粒的上升速度逐渐减小。在相同风力条件下，质量较大的沙粒受到的重力影响更为显著，其起动相对困难，运动轨迹也会因重力作用而更加贴近地面。浮力是由于空气对沙粒的作用而产生的，方向与重力相反，竖直向上。在一般情况下，浮力相对较小，对沙粒运动的影响不如重力明显，但在研究非常细小的沙粒或考虑高精度的风沙运动模型时，浮力的作用不可忽视。对于粒径较小的沙尘颗粒，浮力与重力的比值相对较大，可能会对其悬浮和长距离输送产生一定影响。气动阻力是风沙颗粒在气流中运动时，气流对沙粒表面产生的摩擦力和压力差所引起的阻力。气动阻力的方向与沙粒运动方向相反，其大小与沙粒的形状、粒径、运动速度以及空气的密度和粘性等因素密切相关。根据流体力学原理，当沙粒在气流中运动时，其周围的气流会发生绕流，形成复杂的流场结构，从而产生气动阻力。在风沙运动中，随着沙粒运动速度的增加，气动阻力也会迅速增大，对沙粒的运动起到显著的阻碍作用。压力梯度力是由于气流中存在压力梯度而作用在沙粒上的力。在近地层气流中，由于地形、地表粗糙度等因素的影响，气流的压力分布不均匀，从而产生压力梯度。当沙粒处于这样的压力梯度场中时，会受到压力梯度力的作用，其方向指向压力降低的方向。压力梯度力在风沙颗粒起动和运动过程中可能会起到重要作用，特别是在复杂地形和近地表气流变化剧烈的区域，它可能会改变沙粒的受力平衡，影响沙粒的起动和运动轨迹。此外，风沙颗粒还可能受到马格努斯力、萨夫曼力、静电力等其他力的作用。马格努斯力是由于沙粒的旋转而产生的，当沙粒在气流中运动时，如果发生旋转，就会受到马格努斯力的作用，其方向垂直于沙粒的运动方向和旋转轴。萨夫曼力则是由于沙粒在速度梯度场中运动而产生的，方向与速度梯度方向和沙粒运动方向垂直。静电力是由于沙粒表面带电而产生的，在风沙流中，沙粒之间的摩擦、碰撞等过程可能会使沙粒带电，从而受到静电力的作用。这些力在特定条件下可能对风沙颗粒的运动产生重要影响，但目前对它们的研究还相对较少，其作用机制和影响程度尚有待进一步深入探讨。2.1.3影响因子风沙颗粒的起动受到多种因素的综合影响，这些因素相互关联，共同决定了风沙运动的起始和发展。风速是影响风沙颗粒起动的关键因素之一，它直接决定了作用在沙粒上的风力大小。一般来说，风速越大，作用在沙粒上的力就越大，当风力超过沙粒与地表之间的摩擦力等阻力时，沙粒就会起动。研究表明，存在一个临界起动风速，只有当实际风速达到或超过这个临界值时，风沙颗粒才会开始运动。不同粒径的沙粒，其临界起动风速也不同，粒径越大，所需的临界起动风速越高。在沙漠地区，当风速达到8-10m/s时，粒径为0.1-0.25mm的沙粒就容易被起动，而对于粒径更大的沙粒，则需要更高的风速。此外，风速的变化特征，如阵风、风速的脉动等，也会对风沙颗粒的起动产生影响。阵风的突然出现和风速的脉动会使沙粒受到的风力瞬间增大，更容易满足起动条件，从而增加风沙起动的可能性。沙粒粒径和形状对其起动也有着重要影响。粒径不同的沙粒，其重力、惯性以及与气流的相互作用特性都存在差异。较小粒径的沙粒，由于其质量较轻，惯性小，在相对较低的风速下就能够被起动。但同时，小粒径沙粒更容易受到气流中紊流和其他微小作用力的影响，其运动轨迹可能更加复杂。而较大粒径的沙粒，虽然需要更高的风速才能起动，但一旦起动，由于其较大的惯性，运动相对较为稳定。沙粒的形状也会影响其起动，不规则形状的沙粒与地表的接触面积和摩擦力分布不均匀，起动时所需克服的阻力也更为复杂。扁平状的沙粒在相同风速下，可能比球状沙粒更难起动，因为其与地表的摩擦力更大，且气流对其作用力的方向和大小也会因形状的不规则而发生变化。地表粗糙度是影响风沙颗粒起动的重要下垫面因素。地表粗糙度的存在会改变近地层气流的结构和速度分布，进而影响沙粒的受力情况。粗糙的地表会使气流产生更多的紊流和漩涡，增加气流对沙粒的作用力。同时，地表粗糙度还会影响沙粒与地表的接触状态，增加沙粒起动的难度。在戈壁地区，地表粗糙度较大，风沙颗粒的起动风速相对较高；而在平坦的沙漠地区，地表粗糙度较小，沙粒更容易被起动。此外，地表的微地形起伏、植被覆盖等也会对地表粗糙度产生影响，从而间接影响风沙颗粒的起动。沙丘的迎风坡和背风坡，由于地形的变化，地表粗糙度不同，风沙颗粒的起动情况也会有明显差异。土壤含水率对风沙颗粒的起动也有显著影响。当土壤含水率较高时，水分会在沙粒之间形成水膜，增加沙粒之间的粘结力。这使得沙粒需要更大的风力才能克服这种粘结力而起动，从而提高了风沙颗粒的起动风速。研究表明，土壤含水率每增加1%，风沙颗粒的起动风速可能会增加0.1-0.3m/s。在湿润的沙质土壤地区，由于土壤含水率较高，风沙活动相对较少；而在干旱地区，土壤含水率低，沙粒容易起动，风沙活动频繁。此外，土壤中的盐分含量、有机质含量等也会影响土壤的物理性质，进而对风沙颗粒的起动产生一定的影响。2.2跃移运动的基本模型2.2.1模型假设在风沙颗粒跃移运动的研究中，为了建立合理的理论模型，通常会基于一些假设条件。这些假设在一定程度上简化了复杂的实际情况，使得理论分析和数学求解成为可能。在大多数经典的跃移运动模型中，常假设沙粒为球形。这一假设主要是为了简化对沙粒受力和运动的分析。球形沙粒的几何形状规则，其在气流中的受力特性相对易于描述。在计算气动阻力时，基于球形假设可以运用较为成熟的流体力学公式，如斯托克斯定律适用于小雷诺数下球形颗粒的阻力计算。然而，实际的沙粒形状是极其复杂多样的，并非理想的球形。不规则形状的沙粒在气流中会受到更复杂的气动力作用，其阻力系数、升力系数等与球形沙粒存在显著差异。扁平状或长条状的沙粒，其在气流中的取向不同，所受到的气动力大小和方向也会发生变化，这使得基于球形假设的模型在描述实际沙粒运动时存在一定的局限性。假设风速均匀也是常见的模型假设之一。均匀风速的假设便于建立简单的数学模型，分析风沙颗粒在稳定气流中的运动规律。在一些早期的研究中，通过假设近地层风速为均匀分布，推导出了沙粒的起动风速公式和跃移轨迹方程。但在实际的风沙环境中，风速并非均匀不变。近地层气流受到地表粗糙度、地形起伏、植被覆盖等多种因素的影响，会产生复杂的湍流结构。在沙漠地区，沙丘的存在会导致气流在沙丘迎风坡加速、背风坡减速并形成漩涡，使得风速在空间上呈现出不均匀的分布。此外，大气边界层中的风速还存在着明显的垂直梯度，随着高度的增加，风速逐渐增大。因此，均匀风速假设在实际应用中需要谨慎考虑，尤其是对于复杂地形和地表条件下的风沙运动研究，该假设可能会导致较大的误差。沙粒与空气的相互作用遵循牛顿定律也是跃移运动模型的重要假设。根据牛顿第二定律，可以建立沙粒的运动方程，分析沙粒在各种力作用下的加速度、速度和位移。在考虑气动阻力时，认为阻力与沙粒相对空气的速度成正比，方向相反，这是基于牛顿力学中力与运动的关系。然而，当涉及到微小颗粒或高速气流时，一些微观效应和非牛顿流体行为可能会对沙粒与空气的相互作用产生影响。在纳米级别的沙尘颗粒运动中，布朗运动等微观效应不能被忽略，此时牛顿定律可能无法准确描述沙粒的运动。此外，在高雷诺数的气流中，空气的粘性和湍流特性会变得更加复杂，可能会出现一些不符合牛顿定律简单假设的现象。2.2.2模型参数跃移运动模型中涉及多个关键参数，这些参数对于准确描述风沙颗粒的运动状态和规律起着至关重要的作用。沙粒直径是一个基本的参数，它直接影响沙粒的重力、惯性以及与气流的相互作用。较小直径的沙粒，其重力相对较小，在气流中更容易被起动和携带。粒径为0.05-0.1mm的粉沙颗粒，在较低的风速下就能够被扬起并参与跃移运动。而较大直径的沙粒，如粒径大于0.5mm的粗沙，由于其重力较大，需要更高的风速才能使其起动。沙粒直径还会影响沙粒在气流中的运动轨迹和速度。较小粒径的沙粒，受到气流中微小扰动的影响更大，其运动轨迹可能更加曲折，速度变化也更为复杂。沙粒密度也是一个重要参数，它决定了沙粒的质量和惯性。不同类型的沙粒，其密度存在差异，常见的石英沙密度约为2650kg/m³。密度较大的沙粒，在相同风力条件下，其运动速度相对较慢，起动也更为困难。在相同风速下，密度为3000kg/m³的沙粒比密度为2650kg/m³的沙粒更难被起动，且起动后的运动速度也会更低。沙粒密度还会影响沙粒与其他颗粒或物体碰撞时的能量传递和反弹特性。密度大的沙粒在碰撞时具有更大的动量，可能会对被碰撞物体造成更大的冲击。风速是决定风沙颗粒跃移运动的关键因素之一。风速的大小直接影响作用在沙粒上的风力大小。随着风速的增加，作用在沙粒上的气动阻力和升力增大，当风力超过沙粒与地表之间的摩擦力等阻力时，沙粒就会起动并开始跃移运动。在沙漠地区，当风速达到6-8m/s时，粒径为0.1-0.25mm的沙粒就容易被起动。风速的方向也会影响沙粒的运动轨迹。在风向变化频繁的区域，沙粒的运动轨迹会更加复杂，可能会出现曲折、迂回的路径。空气密度也是跃移运动模型中不可忽视的参数。空气密度与气压、温度等因素有关，在标准大气压和常温下，空气密度约为1.29kg/m³。空气密度的变化会影响气动阻力和升力的大小。在高海拔地区，由于气压较低，空气密度减小，相同风速下作用在沙粒上的气动阻力和升力也会相应减小，这使得沙粒的起动和运动特性发生变化。在海拔4000m的高原地区，空气密度约为0.81kg/m³，相比平原地区，沙粒在该地区的起动风速可能会更高，运动轨迹也会有所不同。起跳角和冲击角也是描述沙粒跃移运动的重要参数。起跳角是沙粒起跳时速度方向与水平面的夹角，它影响沙粒的初始运动方向和上升高度。较小的起跳角通常会使沙粒的运动轨迹更加贴近地面，而较大的起跳角则可能使沙粒上升到较高的高度。冲击角是沙粒撞击地面时速度方向与地面的夹角，它对沙粒撞击地面后的反弹和溅起其他沙粒的情况有重要影响。当冲击角较小时，沙粒可能会在地面上滚动一段距离后再次起跳；而当冲击角较大时，沙粒可能会强烈反弹并溅起更多的沙粒。2.2.3模型求解方法为了求解跃移运动模型方程，通常采用数值方法来获得近似解。不同的数值方法具有各自的特点和适用范围，在风沙颗粒跃移运动研究中发挥着重要作用。有限差分法是一种常用的数值方法。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。在求解风沙颗粒跃移运动方程时，通过将方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。对于描述沙粒运动轨迹的微分方程，利用泰勒级数展开将其转化为差分形式，然后通过迭代求解代数方程组得到沙粒在各个网格节点上的位置、速度等参数。有限差分法的优点是数学概念直观，表达简单，易于编程实现。它在处理规则区域的问题时具有较高的计算效率。在简单的平板边界层内的风沙颗粒跃移模拟中，有限差分法能够快速准确地得到沙粒的运动轨迹和速度分布。但该方法也存在一些缺点，如对网格划分的依赖性较强。若网格划分不合理，可能导致数值解的精度下降，甚至出现错误。当网格间距过大时，可能会丢失一些沙粒运动的细节信息；而网格间距过小时，又会增加计算量和计算时间。有限差分法在处理复杂边界条件和不规则区域时相对困难，需要采用特殊的数值技巧来确保边界条件的正确实现。有限元法也是求解跃移运动模型的重要方法之一。其基本思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。在风沙颗粒跃移运动研究中，有限元法可以用于处理复杂的几何形状和边界条件。对于具有不规则地形的风沙区域，通过将其划分为三角形或四边形等单元，可以准确地模拟沙粒在复杂地形下的运动。有限元法的优点是对计算域（几何区域）复杂度适应性好，能够处理各种复杂的边界条件和不规则形状的计算区域。它在求解结构强度、结构动响应等问题中得到了广泛应用，在风沙运动研究中也逐渐受到重视。但有限元法的计算过程相对复杂，需要进行单元剖分、插值函数选择、总体合成等多个步骤，计算量较大。在大规模的风沙流模拟中，有限元法的计算时间可能较长，对计算机的硬件要求也较高。边界元法是一种基于边界积分方程的数值方法。该方法将求解域内的问题转化为边界上的积分方程，通过对边界进行离散化，将积分方程转化为代数方程组进行求解。在风沙颗粒跃移运动模型求解中，边界元法可以利用边界条件简化计算，减少计算量。在模拟风沙颗粒在固定边界（如防风墙、建筑物等）附近的运动时，边界元法可以只对边界进行离散，而不需要对整个求解域进行离散，从而提高计算效率。边界元法的优点是可以降低问题的维数，减少计算量，尤其适用于处理具有无限域或半无限域的问题。在研究风沙颗粒在大气边界层中的长距离传输时，边界元法可以有效地处理无限远边界条件。但边界元法也存在一些局限性，如需要求解奇异积分，对边界的光滑性要求较高。当边界条件复杂或边界存在不连续时，边界元法的计算精度和效率可能会受到影响。三、近地层风沙颗粒跃移运动的特征3.1运动轨迹特征3.1.1轨迹形状风沙颗粒在近地层的跃移轨迹形状是研究其运动学的关键内容之一，它受到多种因素的综合影响，呈现出复杂的变化规律。通过实验观测和数值模拟等方法，众多研究对风沙颗粒跃移轨迹形状进行了深入探究。在早期的研究中，Bagnold通过风洞实验观察到风沙颗粒的跃移轨迹近似为抛物线。他认为沙粒在气流中主要受到重力和风力的作用，在这两种力的共同影响下，沙粒从地面起跳后，先做上升运动，随着高度增加，重力作用逐渐显现，沙粒速度减小，达到一定高度后开始下降，最终落回地面，其运动轨迹类似于平抛运动的抛物线。这一观点为后续研究奠定了基础，但实际风沙环境远比理想实验条件复杂，抛物线模型并不能完全准确地描述风沙颗粒的跃移轨迹。随着研究的深入，学者们发现风沙颗粒的跃移轨迹并非严格的抛物线，而是存在一定的偏差和不规则性。White等通过实验研究发现，马格努斯力会对沙粒跃移轨迹产生显著影响。当沙粒在气流中运动时，如果发生旋转，就会受到马格努斯力的作用，该力垂直于沙粒的运动方向和旋转轴，使得沙粒的运动轨迹发生偏移。对于表面粗糙或形状不规则的沙粒，在气流中更容易发生旋转，马格努斯力对其轨迹的影响更为明显。一些研究表明，风速的脉动和紊流特性也是导致沙粒跃移轨迹不规则的重要因素。近地层气流并非均匀稳定，存在着各种尺度的涡旋和速度脉动，这些都会使作用在沙粒上的风力发生变化，从而导致沙粒的运动轨迹出现波动和曲折。在强风条件下，气流的紊流强度增大，沙粒跃移轨迹的不规则性更加显著。沙粒粒径和形状对跃移轨迹形状也有着重要影响。不同粒径的沙粒，其重力、惯性以及与气流的相互作用特性存在差异，导致运动轨迹不同。粒径较小的沙粒，质量轻，惯性小，更容易受到气流中微小扰动的影响，其运动轨迹相对更加曲折。粒径为0.05-0.1mm的粉沙颗粒，在气流中可能会出现频繁的上下波动和横向偏移，轨迹较为复杂。而粒径较大的沙粒，如粒径大于0.5mm的粗沙，由于其重力较大，惯性也大，运动相对较为稳定，轨迹相对较为平滑。沙粒的形状同样会影响其运动轨迹。不规则形状的沙粒在气流中的受力情况复杂，其与气流的接触面积和受力方向不断变化，导致运动轨迹不规则。扁平状的沙粒在运动过程中，可能会因为与气流的不同夹角而受到不同方向和大小的作用力，从而使轨迹发生扭曲。3.1.2运动高度与距离风沙颗粒跃移的高度和距离是衡量风沙运动强度和范围的重要指标，它们受到多种因素的共同影响，对风沙地貌的形成和演变以及风沙灾害的发生发展具有重要意义。研究表明，风沙颗粒跃移高度与风速密切相关。一般来说，风速越大，作用在沙粒上的风力就越大，沙粒获得的动能也越大，从而能够跃移到更高的高度。在风洞实验中，当风速从5m/s增加到10m/s时，沙粒的跃移高度明显增加。通过对大量实验数据的分析，发现风沙颗粒跃移高度与风速之间存在幂函数关系，即跃移高度随着风速的增加而呈幂次增长。不同粒径的沙粒，其跃移高度随风速的变化规律也有所不同。粒径较小的沙粒，在较低风速下就能够跃移到一定高度，且随着风速增加，跃移高度增长较为迅速；而粒径较大的沙粒，需要更高的风速才能达到相同的跃移高度，且增长速度相对较慢。沙粒粒径对跃移高度也有着显著影响。粒径较小的沙粒，质量轻，惯性小，在相同风力条件下更容易被扬起并跃移到较高的高度。粒径为0.1-0.25mm的中沙和细沙，在一般起沙风条件下，跃移高度通常在0.5-1.5m之间；而粒径大于0.5mm的粗沙，由于重力较大，跃移高度一般较低，大多在0.5m以下。研究还发现，沙粒的形状也会对跃移高度产生一定影响。不规则形状的沙粒，其在气流中的受力情况复杂，可能会导致跃移高度的变化。扁平状或长条状的沙粒，相较于球形沙粒，在相同条件下跃移高度可能会有所不同。地表粗糙度是影响风沙颗粒跃移高度的重要下垫面因素。粗糙的地表会使气流产生更多的紊流和漩涡，增加气流对沙粒的作用力，但同时也会使沙粒更容易与地表发生碰撞，从而限制沙粒的跃移高度。在戈壁地区，地表粗糙度较大，风沙颗粒的跃移高度相对较低；而在平坦的沙漠地区，地表粗糙度较小，沙粒更容易跃移到较高的高度。植被覆盖也会对风沙颗粒跃移高度产生影响。植被可以增加地表粗糙度，降低风速，同时还能阻挡沙粒的运动，减少沙粒的跃移高度。在有植被覆盖的区域，风沙颗粒的跃移高度通常比无植被覆盖区域低。风沙颗粒跃移的距离同样受到多种因素的影响。风速是决定跃移距离的关键因素之一，风速越大，沙粒在水平方向上的运动速度越快，跃移距离也就越远。在野外风沙观测中发现，当风速达到15m/s以上时，沙粒的跃移距离明显增加，一些沙粒甚至可以跃移数十米远。沙粒粒径也会影响跃移距离，粒径较大的沙粒，由于其惯性大，在水平方向上的运动相对稳定，跃移距离相对较远；而粒径较小的沙粒，容易受到气流扰动的影响，运动轨迹不稳定，跃移距离相对较短。地表粗糙度和植被覆盖等下垫面因素也会对跃移距离产生影响。粗糙的地表和植被覆盖会增加沙粒与地表的摩擦力和碰撞概率，从而缩短沙粒的跃移距离。在沙丘迎风坡，由于地表坡度和粗糙度的变化，风沙颗粒的跃移距离会受到影响，可能会出现沙粒在较短距离内就发生沉积的情况。3.2运动速度与加速度特征3.2.1速度变化规律风沙颗粒在跃移过程中的速度变化规律是其运动学研究的关键内容，深入了解这一规律对于揭示风沙运动的本质具有重要意义。众多学者通过实验观测和理论分析等方法，对风沙颗粒跃移速度的变化进行了广泛研究。在风沙颗粒的起跳阶段，起跳速度是一个重要参数，它受到多种因素的影响。风速是决定起跳速度的关键因素之一，风速越大，作用在沙粒上的风力就越大，沙粒获得的初始动能也越大，起跳速度也就越高。通过风洞实验发现，当风速从5m/s增加到10m/s时，沙粒的起跳速度明显增大。沙粒粒径也会对起跳速度产生影响。粒径较小的沙粒，质量轻，惯性小，在相同风力条件下更容易被加速，起跳速度相对较高。粒径为0.1-0.25mm的中沙和细沙，起跳速度一般在1-3m/s之间；而粒径大于0.5mm的粗沙，由于重力较大，起跳速度相对较低，大多在1m/s以下。此外，沙粒的形状和地表状况也会影响起跳速度。不规则形状的沙粒与地表的摩擦力和作用力更为复杂，可能导致起跳速度的变化。粗糙的地表会使沙粒起跳时受到更大的阻力，从而降低起跳速度。在跃移的上升阶段，风沙颗粒的速度会发生变化。随着沙粒上升，重力作用逐渐显现，沙粒的上升速度逐渐减小。同时，气动阻力也会对沙粒的运动产生阻碍作用，进一步降低沙粒的速度。研究表明，上升速度与起跳速度、风速以及沙粒的受力情况密切相关。当起跳速度较大时，沙粒在上升阶段能够保持较高的速度，但随着高度增加，速度下降也会更快。风速的大小和方向也会影响上升速度，在逆风条件下，沙粒的上升速度会受到更大的抑制。在下降阶段，风沙颗粒在重力作用下加速下降，速度逐渐增大。但同时，气动阻力也会随着速度的增加而增大，当气动阻力与重力达到平衡时，沙粒会达到一个稳定的下降速度。这个稳定下降速度与沙粒的粒径、密度以及空气的密度和粘性等因素有关。粒径较大、密度较大的沙粒，其稳定下降速度相对较高。研究还发现，沙粒在下降过程中，可能会受到气流中涡旋和其他沙粒的影响，导致速度出现波动。在强风条件下，气流中的涡旋较多，沙粒下降速度的波动会更加明显。水平方向上，风沙颗粒的运动速度也会受到多种因素的影响。风速是决定水平速度的主要因素，风速越大，沙粒在水平方向上的运动速度就越快。在野外风沙观测中发现，当风速达到15m/s以上时，沙粒在水平方向上的运动速度可达到5-8m/s。风向的变化也会导致沙粒水平运动方向和速度的改变。此外，沙粒之间的相互碰撞以及与地表的碰撞也会影响水平速度。沙粒在碰撞过程中会发生动量交换，导致速度大小和方向的变化。在风沙流中，大量沙粒相互碰撞，使得水平速度呈现出复杂的分布特征。3.2.2加速度变化规律风沙颗粒在跃移过程中的加速度变化规律是其运动学研究的重要方面，它反映了沙粒在各种力作用下运动状态的改变。通过理论分析和实验研究，学者们对风沙颗粒跃移加速度的变化进行了深入探讨。在起跳瞬间，风沙颗粒受到风力、重力、气动阻力等多种力的作用，产生一个较大的加速度。此时，风力是使沙粒获得加速度的主要动力，而重力和气动阻力则起到阻碍作用。风速越大，作用在沙粒上的风力就越大，起跳加速度也就越大。研究表明，起跳加速度与风速的平方成正比。沙粒粒径也会影响起跳加速度，粒径较小的沙粒，质量轻，在相同风力条件下更容易被加速，起跳加速度相对较大。在风洞实验中，粒径为0.1-0.25mm的沙粒，起跳加速度可达10-20m/s²；而粒径大于0.5mm的沙粒，起跳加速度一般在5-10m/s²之间。此外，沙粒的形状和地表状况也会对起跳加速度产生影响。不规则形状的沙粒与地表的摩擦力和作用力更为复杂，可能导致起跳加速度的变化。粗糙的地表会使沙粒起跳时受到更大的阻力，从而降低起跳加速度。在上升阶段，风沙颗粒受到重力和气动阻力的作用，加速度方向与运动方向相反，沙粒做减速运动，加速度逐渐减小。随着沙粒上升高度的增加，重力作用逐渐增强，气动阻力也会随着速度的减小而减小，但重力的影响更为显著。研究发现，上升阶段的加速度与沙粒的运动速度、高度以及空气的密度和粘性等因素密切相关。当沙粒运动速度较大时，气动阻力较大，加速度也较大；随着高度增加，重力作用增强，加速度逐渐减小。在风速为10m/s的条件下，粒径为0.2mm的沙粒在上升初期的加速度约为-15m/s²，随着上升高度的增加，加速度逐渐减小，当达到最大上升高度时，加速度趋近于0。在下降阶段，风沙颗粒受到重力和气动阻力的作用，加速度方向与运动方向相同，沙粒做加速运动。重力是使沙粒加速的主要动力，而气动阻力则起到阻碍作用。随着沙粒下降速度的增加，气动阻力逐渐增大，当气动阻力与重力达到平衡时，加速度减小为0，沙粒达到稳定下降速度。稳定下降速度与沙粒的粒径、密度以及空气的密度和粘性等因素有关。粒径较大、密度较大的沙粒，其稳定下降速度相对较高，达到稳定下降速度所需的时间也相对较短。研究还发现，沙粒在下降过程中，可能会受到气流中涡旋和其他沙粒的影响，导致加速度出现波动。在强风条件下，气流中的涡旋较多，沙粒下降加速度的波动会更加明显。此外，风沙颗粒在跃移过程中，还可能受到其他力的作用，如马格努斯力、萨夫曼力等，这些力也会对加速度产生影响。马格努斯力是由于沙粒的旋转而产生的，当沙粒在气流中运动时，如果发生旋转，就会受到马格努斯力的作用，其方向垂直于沙粒的运动方向和旋转轴。马格努斯力会使沙粒的运动轨迹发生偏移，同时也会改变沙粒的加速度。萨夫曼力则是由于沙粒在速度梯度场中运动而产生的，方向与速度梯度方向和沙粒运动方向垂直。萨夫曼力对沙粒加速度的影响相对较小，但在某些情况下也不可忽视。在研究风沙颗粒在复杂气流环境中的运动时，需要综合考虑这些力对加速度的影响，以更准确地描述沙粒的运动状态。3.3跃移频率与输沙量特征3.3.1跃移频率风沙颗粒的跃移频率是衡量风沙运动活跃程度的重要指标，它反映了沙粒在单位时间内的跃移次数，对于深入理解风沙运动过程和输沙机制具有关键作用。众多研究表明，风沙颗粒跃移频率与风速之间存在密切的正相关关系。随着风速的增加，作用在沙粒上的风力增大，更多的沙粒获得足够的能量起动并参与跃移运动，从而导致跃移频率显著增加。在风洞实验中，当风速从5m/s逐渐增加到10m/s时，风沙颗粒的跃移频率呈现出明显的上升趋势。通过对大量实验数据的统计分析发现，跃移频率与风速之间近似满足幂函数关系，即跃移频率随风速的增大而呈幂次增长。这是因为风速的增加不仅提高了沙粒的起动概率，还加快了沙粒的运动速度，使得沙粒在单位时间内能够完成更多次的跃移。沙粒粒径对跃移频率也有着重要影响。粒径较小的沙粒，质量轻，惯性小，在相同风力条件下更容易被起动和加速，因此其跃移频率相对较高。粒径为0.1-0.25mm的中沙和细沙，在一般起沙风条件下，跃移频率通常在10-20次/秒之间；而粒径大于0.5mm的粗沙，由于重力较大，起动和加速相对困难，跃移频率一般较低，大多在5-10次/秒之间。研究还发现，沙粒形状的不规则性也会对跃移频率产生一定影响。不规则形状的沙粒与气流的相互作用更为复杂，其在气流中的受力情况不断变化，可能导致跃移频率的波动。扁平状或长条状的沙粒，相较于球形沙粒，在相同条件下跃移频率可能会有所不同。地表粗糙度是影响风沙颗粒跃移频率的重要下垫面因素。粗糙的地表会使气流产生更多的紊流和漩涡，增加气流对沙粒的作用力，但同时也会使沙粒更容易与地表发生碰撞，从而改变沙粒的运动状态。在戈壁地区，地表粗糙度较大，风沙颗粒在运动过程中与地表的碰撞概率增加，跃移频率可能会相对降低；而在平坦的沙漠地区，地表粗糙度较小，沙粒受到的阻碍较小，跃移频率相对较高。此外，植被覆盖也会对跃移频率产生影响。植被可以增加地表粗糙度，降低风速，同时还能阻挡沙粒的运动，减少沙粒的跃移频率。在有植被覆盖的区域，风沙颗粒的跃移频率通常比无植被覆盖区域低。3.3.2输沙量近地层风沙颗粒跃移的输沙量是衡量风沙活动强度和风沙灾害程度的关键指标，它反映了风沙流中携带的沙粒总量，对于研究风沙地貌的形成、土地沙漠化的发展以及风沙灾害的防治具有重要意义。大量研究表明，输沙量与风速之间存在显著的正相关关系。风速是决定输沙量的关键因素之一，随着风速的增大，作用在沙粒上的风力增强，更多的沙粒被起动并参与跃移运动，输沙量也随之急剧增加。在野外风沙观测中发现，当风速从8m/s增加到12m/s时，输沙量可增加数倍甚至数十倍。通过对大量实验数据和野外观测数据的分析，发现输沙量与风速之间通常满足幂函数关系，即输沙量随风速的增大而呈幂次增长。这是因为风速的增加不仅提高了沙粒的起动概率和跃移频率，还增大了沙粒的运动速度和跃移距离，使得风沙流能够携带更多的沙粒。沙粒粒径对输沙量也有着重要影响。不同粒径的沙粒在风沙流中的运动特性和输移能力存在差异，从而导致输沙量的变化。粒径较大的沙粒，虽然起动相对困难，但一旦起动，由于其较大的惯性和质量，在跃移过程中能够携带更多的能量，对输沙量的贡献也较大。粒径大于0.5mm的粗沙，虽然在风沙流中的含量相对较少，但由于其单个沙粒的质量较大，在输沙量中所占的比例可能并不低。而粒径较小的沙粒，虽然起动容易且跃移频率较高，但由于其质量轻，单个沙粒对输沙量的贡献相对较小。粒径为0.05-0.1mm的粉沙颗粒，在风沙流中含量较多，但由于其质量小，总体输沙量可能不如粒径较大的沙粒。研究还发现，沙粒的粒径分布对输沙量也有重要影响。当风沙流中存在多种粒径的沙粒时，不同粒径沙粒之间的相互作用和协同效应会影响输沙量的大小。如果风沙流中粗细沙粒搭配合理，可能会增加沙粒之间的摩擦力和碰撞概率，从而提高输沙量；反之，如果风沙流中沙粒粒径过于单一，输沙量可能会受到一定限制。地表粗糙度和植被覆盖是影响输沙量的重要下垫面因素。地表粗糙度的增加会使气流产生更多的紊流和漩涡，增加气流对沙粒的作用力，但同时也会使沙粒更容易与地表发生碰撞，从而降低沙粒的运动速度和跃移距离，减少输沙量。在戈壁地区，地表粗糙度较大，风沙颗粒在运动过程中与地表的碰撞频繁，输沙量相对较低；而在平坦的沙漠地区，地表粗糙度较小，沙粒受到的阻碍较小，输沙量相对较高。植被覆盖可以增加地表粗糙度，降低风速，同时还能阻挡沙粒的运动，减少沙粒的起动和跃移，从而显著降低输沙量。在有植被覆盖的区域，风沙颗粒的输沙量通常比无植被覆盖区域低很多。研究表明，植被覆盖度每增加10%，输沙量可能会降低20%-30%。此外，地表的微地形起伏、土壤含水率等因素也会对输沙量产生影响。沙丘的迎风坡和背风坡，由于地形的变化，输沙量会有明显差异；土壤含水率较高时，沙粒之间的粘结力增大，起动风速提高，输沙量也会相应减少。四、近地层风沙颗粒跃移运动的影响因素4.1气象因素4.1.1风速与风向风速与风向是影响近地层风沙颗粒跃移运动的关键气象因素，它们对风沙颗粒的起动、运动轨迹和输沙量有着显著的影响。风速是决定风沙颗粒能否起动以及起动后运动状态的首要因素。当风速达到一定阈值时，作用在沙粒上的风力足以克服沙粒与地表之间的摩擦力、粘结力等阻力，沙粒开始起动并进入跃移运动。研究表明，临界起动风速与沙粒粒径密切相关，粒径越大，所需的临界起动风速越高。粒径为0.1-0.25mm的沙粒，其临界起动风速一般在4-6m/s左右；而粒径大于0.5mm的沙粒，临界起动风速则可能超过8m/s。随着风速的进一步增大，更多的沙粒被起动，跃移运动变得更加活跃，输沙量也随之增加。风速与输沙量之间存在幂函数关系，输沙量随风速的增大而迅速增加。当风速从8m/s增加到12m/s时，输沙量可能会增加数倍甚至数十倍。这是因为风速的增大不仅提高了沙粒的起动概率，还增大了沙粒的运动速度和跃移距离，使得风沙流能够携带更多的沙粒。风速的变化特征，如阵风、风速的脉动等，也会对风沙颗粒跃移运动产生重要影响。阵风的突然出现会使沙粒受到瞬间增大的风力作用，更容易满足起动条件，从而增加风沙起动的可能性。风速的脉动会导致沙粒受到的风力不稳定，运动轨迹出现波动和曲折。在强风条件下，风速的脉动更为明显，沙粒的运动轨迹更加复杂，输沙量的变化也更加剧烈。风向的变化会直接影响风沙颗粒的运动方向。在单一风向的条件下，风沙颗粒的运动方向相对较为稳定，呈现出与风向一致的趋势。但在实际的风沙环境中，风向往往是多变的。当风向发生改变时，风沙颗粒的运动方向也会随之改变，这可能导致沙粒的运动轨迹出现曲折、迂回的情况。在沙漠地区，由于地形和大气环流的影响，风向经常发生变化，使得风沙颗粒的运动轨迹变得复杂多样。风向的变化还会影响风沙颗粒的堆积和分布。在不同风向的交替作用下，沙粒可能会在不同的区域堆积，形成复杂的风沙地貌。此外，风速和风向的协同作用也会对风沙颗粒跃移运动产生影响。在风向变化频繁的情况下，不同方向的风力对沙粒的作用相互叠加，可能会使沙粒的运动状态更加复杂。在一个区域内，同时存在来自不同方向的风，且风速大小不同，沙粒可能会受到多个方向风力的拉扯，其运动轨迹和速度会发生频繁的变化。这种复杂的风速和风向条件会导致输沙量的分布更加不均匀，对风沙地貌的形成和演变产生重要影响。4.1.2温度与湿度温度和湿度是影响近地层风沙颗粒跃移运动的重要气象因素，它们通过改变沙粒的物理性质和运动状态，对风沙运动产生显著影响。温度的变化会影响沙粒的物理性质，进而影响其跃移运动。一方面，温度升高会导致沙粒的热膨胀，使沙粒之间的接触状态发生变化，从而影响沙粒与地表之间的摩擦力和粘结力。当温度升高时，沙粒的体积膨胀，与地表的接触面积减小，摩擦力和粘结力降低，沙粒更容易起动。在炎热的夏季，沙漠地区的地表温度可达50℃以上，此时沙粒更容易被风吹起，风沙活动相对频繁。另一方面，温度还会影响空气的密度和粘性，进而影响风沙颗粒与气流之间的相互作用。随着温度升高，空气密度减小，粘性降低，作用在沙粒上的气动阻力和升力也会发生变化。在高温条件下，空气密度减小，相同风速下作用在沙粒上的气动阻力减小，沙粒的运动速度可能会增加。湿度对风沙颗粒跃移运动的影响主要体现在沙粒的吸湿和粘结作用上。当空气湿度较大时，沙粒表面会吸附水分，形成水膜，增加沙粒之间的粘结力。这使得沙粒需要更大的风力才能克服粘结力而起动，从而提高了风沙颗粒的起动风速。研究表明，土壤含水率每增加1%，风沙颗粒的起动风速可能会增加0.1-0.3m/s。在湿润的沙质土壤地区，由于土壤含水率较高，风沙活动相对较少。湿度还会影响沙粒的表面性质和静电特性，进而影响沙粒之间的相互作用。在高湿度环境下，沙粒表面的水膜可能会改变其表面电荷分布，导致沙粒之间的静电力发生变化，影响沙粒的团聚和分散状态。温度和湿度的变化还会对风沙流的结构和输沙量产生影响。在高温低湿的条件下，风沙流中的沙粒更容易起动和运动，输沙量相对较大。而在低温高湿的条件下，沙粒的起动受到抑制，输沙量相对较小。温度和湿度的日变化和季节变化也会导致风沙颗粒跃移运动的周期性变化。在一天中，白天温度较高，湿度较低，风沙活动相对较强；夜晚温度降低，湿度增加，风沙活动减弱。在季节上，春季和冬季温度较低，湿度较小，风沙活动较为频繁；夏季和秋季温度较高，湿度较大，风沙活动相对较少。4.2地表因素4.2.1地表粗糙度地表粗糙度是影响近地层风沙颗粒跃移运动的重要地表因素之一，它对风沙颗粒的起动、运动轨迹和输沙量等方面都有着显著的影响。地表粗糙度的增加会改变近地层气流的结构。当气流流经粗糙的地表时，会在地表附近产生更多的紊流和漩涡。这些紊流和漩涡会使气流的速度分布变得不均匀，增加了气流的能量耗散。在戈壁地区，地表布满了大小不一的砾石，地表粗糙度较大，气流流经时会产生强烈的紊流。通过实验观测和数值模拟发现，在这种粗糙地表条件下，近地层气流的平均速度会降低，且速度的脉动增强。这种气流结构的改变会对风沙颗粒的运动产生重要影响。在风沙颗粒的起动阶段，地表粗糙度的增加会使沙粒的起动风速增大。粗糙的地表会增加沙粒与地表之间的摩擦力和碰撞概率，使得沙粒需要更大的风力才能克服这些阻力而起动。研究表明，当地表粗糙度增大时，临界起动风速会相应提高。在地表粗糙度较大的区域，沙粒的起动风速可能比平坦地表高出1-2m/s。这是因为粗糙地表上的沙粒受到的气流作用力更加分散，且更容易与地表的障碍物发生碰撞，从而增加了起动的难度。对于风沙颗粒的运动轨迹，地表粗糙度的变化也会产生影响。粗糙的地表会使风沙颗粒在运动过程中更容易与地表发生碰撞，导致运动轨迹的改变。在风洞实验中，当在沙床上设置不同高度的障碍物以模拟不同的地表粗糙度时，发现随着地表粗糙度的增加，风沙颗粒的运动轨迹变得更加曲折，跃移高度和距离也会减小。这是因为沙粒在与地表障碍物碰撞后，会损失部分能量，改变运动方向，从而影响其后续的运动轨迹。地表粗糙度还会对风沙颗粒的输沙量产生影响。由于地表粗糙度的增加会使气流的能量耗散增大，风速降低，从而减少了风沙流中能够携带的沙粒数量。在地表粗糙度较大的区域，输沙量通常比平坦地表低。研究表明，地表粗糙度每增加一定程度，输沙量可能会降低20%-40%。这是因为风速的降低使得沙粒的起动和跃移受到抑制，同时也减少了沙粒在气流中的悬浮时间和运动距离，进而降低了输沙量。4.2.2植被覆盖植被覆盖是影响近地层风沙颗粒跃移运动的另一个重要地表因素，它通过多种方式对风沙运动起到抑制作用。植被可以显著降低风速，从而减少风沙颗粒的起动和跃移。植被的枝叶会对气流产生阻挡作用，使气流在穿过植被时发生分流和减速。在有植被覆盖的区域，风速会明显降低。通过实地观测发现，在植被覆盖率为30%-50%的草原地区，近地层风速比无植被覆盖的裸地降低了30%-50%。风速的降低使得作用在沙粒上的风力减小，当风力不足以克服沙粒与地表之间的摩擦力等阻力时，沙粒就难以起动，从而抑制了风沙运动。植被覆盖还能增加地表粗糙度。植被的根系、茎干和枝叶等都会增加地表的粗糙度，改变近地层气流的结构。植被的根系可以固定土壤，防止土壤颗粒被风吹起。植物的茎干和枝叶在地面上形成了复杂的障碍物，使气流在流经时产生更多的紊流和漩涡。在沙漠边缘的绿洲地区，种植的防护林带使得地表粗糙度大幅增加，近地层气流的速度分布变得更加不均匀。这种增加的地表粗糙度会使风沙颗粒的起动风速增大，运动轨迹变得更加复杂，从而减少了风沙颗粒的跃移运动。植被还可以直接阻挡风沙颗粒的运动。植物的枝叶能够拦截风沙流中的沙粒，使其沉积在植被表面或附近。在风沙流经过植被时，部分沙粒会被枝叶捕获，从而减少了风沙流中的沙粒含量。研究表明，植被的覆盖度越高，对风沙颗粒的拦截效果越好。在植被覆盖率达到70%以上的区域，风沙流中的沙粒含量可减少50%以上。此外，植被还可以通过改变地表的微地形，如形成植被沙丘等，进一步阻挡风沙颗粒的运动，降低风沙灾害的影响。4.3沙粒自身因素4.3.1粒径与形状沙粒的粒径和形状是影响近地层风沙颗粒跃移运动的重要自身因素，它们对沙粒的起动、运动轨迹和输沙量等方面都有着显著的影响。不同粒径的沙粒在起动过程中表现出明显的差异。粒径较小的沙粒，由于其质量较轻，惯性小，与地表之间的摩擦力和粘结力相对较小，因此在相对较低的风速下就能够被起动。粒径为0.05-0.1mm的粉沙颗粒，在风速达到3-4m/s时就可能开始起动。而粒径较大的沙粒，如粒径大于0.5mm的粗沙，其质量较大，惯性也大，与地表的摩擦力和粘结力较强，需要更高的风速才能克服这些阻力而起动。研究表明，沙粒的起动风速与粒径之间存在一定的正相关关系，粒径每增大一个量级，起动风速可能会增加1-2m/s。在跃移运动过程中，沙粒粒径对运动轨迹和速度也有着重要影响。粒径较小的沙粒，受到气流中微小扰动的影响更大，其运动轨迹相对更加曲折。这些小粒径沙粒在气流中容易受到紊流和漩涡的作用，导致其运动方向和速度频繁变化。粒径为0.1-0.25mm的中沙和细沙，在跃移过程中可能会出现多次起伏和横向偏移，运动轨迹较为复杂。而粒径较大的沙粒，由于其惯性大，运动相对较为稳定，轨迹相对较为平滑。大粒径沙粒在运动过程中能够保持相对稳定的速度和方向，不易受到气流中微小扰动的影响。沙粒的形状同样会对跃移运动产生重要影响。实际的沙粒形状并非理想的球形，而是多种多样，如扁平状、长条状、不规则多面体等。不规则形状的沙粒在气流中的受力情况复杂，其与气流的接触面积和受力方向不断变化，导致运动轨迹不规则。扁平状的沙粒在运动过程中，由于其与气流的夹角不同，受到的气动力大小和方向也会发生变化，从而使轨迹发生扭曲。长条状的沙粒在气流中可能会发生旋转，进一步增加了运动的复杂性。沙粒形状还会影响其与其他沙粒或物体碰撞时的反弹和溅射情况。不规则形状的沙粒在碰撞时，可能会产生不同的反弹角度和溅射效果，从而影响风沙流中沙粒的分布和运动。4.3.2矿物成分与密度沙粒的矿物成分和密度是影响其跃移运动的重要自身因素，它们通过改变沙粒的物理性质和受力情况，对风沙颗粒的跃移运动产生显著影响。不同的矿物成分决定了沙粒的密度、硬度、表面性质等物理特性。石英是沙粒中常见的矿物成分，其密度相对较大，约为2650kg/m³，硬度较高。富含石英的沙粒在跃移运动中，由于其较大的密度和硬度，具有较大的惯性，运动相对较为稳定。而一些含有云母等轻质矿物的沙粒，其密度相对较小，在相同风力条件下，更容易被起动和加速，运动轨迹也可能更加复杂。云母的密度约为2700-3100kg/m³，但由于其片状结构，在气流中更容易受到气动力的作用而发生旋转和摆动。沙粒的密度直接影响其在跃移运动中的受力和运动状态。密度较大的沙粒，在相同风力条件下，受到的重力作用更大，起动相对困难。在风洞实验中，当风速为8m/s时，密度为2650kg/m³的沙粒可能刚刚达到起动条件，而密度为2000kg/m³的沙粒可能已经开始活跃的跃移运动。密度还会影响沙粒在气流中的运动速度和轨迹。密度大的沙粒，由于其惯性大，在运动过程中速度变化相对较小，轨迹相对较为平滑；而密度小的沙粒，容易受到气流扰动的影响，速度变化较大，轨迹也更加曲折。矿物成分和密度还会影响沙粒之间的相互作用。不同矿物成分的沙粒，其表面电荷分布和化学性质可能存在差异，从而影响沙粒之间的静电作用力和粘结力。在风沙流中，沙粒之间的相互碰撞和摩擦会导致电荷的转移和积累，不同矿物成分的沙粒可能会因为电荷差异而发生团聚或分散。密度不同的沙粒在碰撞时，会发生动量的交换，影响彼此的运动方向和速度。密度大的沙粒在与密度小的沙粒碰撞时，可能会使密度小的沙粒获得更大的动量，从而改变其运动轨迹。五、近地层风沙颗粒跃移运动的研究方法5.1实验观测方法5.1.1风洞实验风洞实验是研究近地层风沙颗粒跃移运动的重要手段之一，它能够在可控的实验条件下，精确模拟风沙运动过程，为深入了解风沙颗粒的运动特性提供关键数据。风洞实验装置主要由风洞本体、供沙系统、测量系统等部分组成。风洞本体是产生稳定气流的核心部件，通常采用直流式或回流式结构。直流式风洞结构简单，气流通过风机直接进入实验段，但其能耗较大，且气流稳定性相对较差；回流式风洞则通过回流管道使气流循环流动，能耗较低，气流稳定性好，在风沙颗粒跃移运动研究中应用更为广泛。实验段是进行风沙实验的区域，其尺寸和形状根据研究需求而定，一般要求实验段的气流均匀性好，湍流度低。供沙系统用于向实验段提供沙粒，常见的供沙方式有重力供沙、机械供沙和气动供沙等。重力供沙是利用沙粒的重力作用，使其自然下落进入实验段；机械供沙则通过机械装置，如螺旋输送机、振动筛等，将沙粒输送到实验段；气动供沙是利用气流将沙粒吹入实验段，能够实现连续稳定的供沙。测量系统用于测量风沙颗粒的运动参数，如风速、沙粒速度、运动轨迹、输沙率等。常用的测量仪器包括风速仪、激光粒度仪、高速摄像机、集沙仪等。风速仪用于测量实验段内的风速，常见的风速仪有热线风速仪、超声风速仪等；激光粒度仪可以测量沙粒的粒径分布；高速摄像机能够拍摄沙粒的运动轨迹，通过图像分析获取沙粒的速度和加速度等参数；集沙仪则用于收集风沙流中的沙粒，测量输沙率。在进行风洞实验时，首先需要根据研究目的和需求，确定实验方案，包括实验段的风速、沙粒的粒径和供给量、实验时间等参数。然后，安装和调试实验装置，确保风洞本体、供沙系统和测量系统正常工作。在实验过程中，启动风机，使实验段内产生稳定的气流，同时通过供沙系统向实验段提供沙粒，模拟风沙运动过程。利用测量系统实时测量风沙颗粒的运动参数，并记录数据。在研究不同风速对风沙颗粒跃移运动的影响时，可以设置多个不同的风速值，依次进行实验，测量每个风速下沙粒的运动轨迹、速度、输沙率等参数。实验结束后，对采集到的数据进行整理和分析，通过数据处理和统计分析方法，提取有价值的信息，揭示风沙颗粒跃移运动的规律。数据采集是风洞实验的关键环节，直接影响研究结果的准确性和可靠性。为了确保数据的质量，需要合理选择测量仪器，并进行校准和标定。风速仪在使用前需要进行校准，以保证测量的风速准确可靠。高速摄像机的帧率和分辨率也需要根据实验要求进行设置，以确保能够清晰拍摄沙粒的运动轨迹。在数据采集过程中，要注意避免外界干扰，确保实验环境的稳定性。同时，为了提高数据的代表性，需要进行多次重复实验，对实验数据进行统计分析，减小实验误差。通过对多次实验数据的平均值和标准差进行计算，可以评估实验数据的可靠性和重复性。5.1.2野外观测野外观测是研究近地层风沙颗粒跃移运动的另一种重要方法，它能够获取实际风沙环境下的第一手数据，弥补风洞实验在真实场景模拟方面的不足。野外观测的方法主要包括定点观测和移动观测。定点观测是在选定的观测站点上，设置各种观测仪器，长期监测风沙颗粒的运动情况。观测站点通常选择在风沙活动频繁、具有代表性的地区，如沙漠边缘、戈壁滩等。在观测站点上，安装风速仪、风向仪、集沙仪、颗粒物监测仪等设备，实时测量风速、风向、输沙率、风沙颗粒的浓度和粒径分布等参数。移动观测则是利用移动观测平台，如车辆、无人机等，在不同的区域进行风沙观测。通过移动观测，可以获取不同地形、地貌和下垫面条件下的风沙运动数据，了解风沙运动的空间变化规律。利用搭载了风速仪、集沙仪和高清摄像机的无人机，在沙漠地区进行飞行观测，获取不同区域的风沙运动信息。常用的野外观测设备有多种。集沙仪是测量输沙率的重要设备，它通过收集风沙流中的沙粒，来计算单位时间内通过单位面积的沙量。常见的集沙仪有垂直集沙仪、水平集沙仪和旋转集沙仪等。垂直集沙仪的集沙口垂直于地面，能够收集垂直方向上的沙粒；水平集沙仪的集沙口平行于地面，适用于测量水平方向上的输沙率；旋转集沙仪则通过旋转的集沙杯来收集沙粒，能够在不同风向条件下进行测量。风速仪和风向仪用于测量近地层的风速和风向，常见的风速仪有三杯式风速仪、螺旋桨式风速仪和超声风速仪等，风向仪则有风向标式和电子罗盘式等。传感器在野外观测中也起着重要作用，如颗粒物监测仪可以实时监测风沙颗粒的浓度和粒径分布，土壤湿度传感器可以测量地表土壤的含水率，这些数据对于分析风沙颗粒跃移运动与环境因素的关系具有重要意义。野外观测对于研究近地层风沙颗粒跃移运动具有重要意义。它能够真实反映实际风沙环境下的风沙运动情况，为理论研究和数值模拟提供可靠的数据支持。通过野外观测，可以获取不同季节、不同天气条件下的风沙运动数据，了解风沙运动的时间变化规律。在春季和冬季，风沙活动较为频繁，通过野外观测可以分析不同季节风沙颗粒跃移运动的差异。野外观测还可以研究不同地形、地貌和下垫面条件对风沙颗粒跃移运动的影响。在沙丘地区，通过观测不同部位的风沙运动情况，可以了解沙丘的形态变化和风沙输移规律。此外，野外观测还能够为风沙灾害的监测和预警提供依据，通过实时监测风沙颗粒的运动参数，及时发现风沙灾害的迹象，采取相应的防护措施。5.2数值模拟方法5.2.1计算流体力学（CFD）方法计算流体力学（CFD）方法在模拟近地层风沙颗粒跃移运动中具有重要作用，其原理基于流体力学的基本守恒方程，通过数值计算的方式求解这些方程，从而得到流场的相关信息。CFD方法的基本原理是将连续的流体空间离散为有限个计算单元，利用数值方法求解Navier-Stokes方程等控制方程。Navier-Stokes方程描述了流体的质量守恒、动量守恒和能量守恒，在近地层风沙颗粒跃移运动模拟中，通过对该方程的求解，可以得到气流的速度场、压力场等信息。对于二维不可压缩粘性流体，其Navier-Stokes方程的分量形式为：\begin{cases}\frac{\partialu}{\partialx}+\frac{\partialv}{\partialy}=0\\\rho\left(\frac{\partialu}{\partialt}+u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}\right)=-\frac{\partialp}{\partialx}+\mu\left(\frac{\partial^{2}u}{\partialx^{2}}+\frac{\partial^{2}u}{\partialy^{2}}\right)\\\rho\left(\frac{\partialv}{\partialt}+u\frac{\partialv}{\partialx}+v\frac{\partialv}{\partialy}\right)=-\frac{\partialp}{\partialy}+\mu\left(\frac{\partial^{2}v}{\partialx^{2}}+\frac{\partial^{2}v}{\partialy^{2}}\right)\end{cases}其中，u和v分别是x和y方向的速度分量，\rho是流体密度，p是压力，\mu是动力粘性系数，t是时间。在实际计算中，需要对这些方程进行离散化处理，常用的离散方法有有限差分法、有限元法和有限体积法等。有限差分法是将方程中的导数用差商代替，通过在网格节点上建立代数方程组来求解；有限元法是将求解区域划分为有限个单元，通过插值函数将方程离散为单元节点上的代数方程组；有限体积法是将控制方程在有限大小的控制体积上进行积分，从而得到离散方程。在模拟近地层风沙颗粒跃移运动时，CFD方法通常将气流视为连续相，而将风沙颗粒视为离散相。通过耦合连续相和离散相的相互作用，来模拟风沙颗粒在气流中的运动。在考虑风沙颗粒与气流之间的相互作用力时，常用的模型有单向耦合模型和双向耦合模型。单向耦合模型只考虑气流对风沙颗粒的作用力，而不考虑风沙颗粒对气流的反作用；双向耦合模型则同时考虑了两者之间的相互作用。在双向耦合模型中，风沙颗粒对气流的反作用会影响气流的速度场和压力场，从而进一步影响风沙颗粒的运动。例如，风沙颗粒在气流中运动会产生尾流，尾流会改变周围气流的速度分布，这种变化又会反过来影响风沙颗粒的运动轨迹。CFD方法在模拟近地层风沙颗粒跃移运动方面具有诸多优点。它能够处理复杂的几何形状和边界条件，对于实际的风沙环境，如沙丘、建筑物等复杂地形条件下的风沙颗粒跃移运动，CFD方法可以通过建立相应的几何模型，准确地模拟气流在这些复杂地形中的流动特性，进而分析风沙颗粒的运动情况。在研究沙漠地区的风沙流经过沙丘时，CFD方法可以精确地模拟气流在沙丘迎风坡的加速和背风坡的减速过程，以及沙丘周围的气流漩涡结构，从而深入了解风沙颗粒在这些区域的运动规律。CFD方法可以提供详细的流场信息，包括气流速度、压力、温度等参数的分布，以及风沙颗粒的运动轨迹、速度、加速度等信息。这些信息对于深入研究风沙颗粒跃移运动的机制和规律具有重要意义。通过CFD模拟，可以得到不同高度和位置处的气流速度和压力分布，分析这些参数对风沙颗粒起动和运动的影响。此外，CFD方法还具有高效性和灵活性，能够快速地进行大量的数值实验，节省时间和成本。在研究不同风速、风向、沙粒粒径等因素对风沙颗粒跃移运动的影响时，可以通过改变CFD模型的输入参数，快速地得到不同工况下的模拟结果，从而进行系统的分析和研究。然而，CFD方法也存在一些缺点。