pajek详细的各功能实现步骤

Benchmark functions.zip 是一些常用的优化算法测试函数，共有17个。

对船舶数学模型的各个推进器参数进行系统辨识，具有实际价值李文华,等:船舶动力定位系统数学模型参数辨识方法研究针对动力定位技术的发展,我国研究人员也进表1离散时间摘要扩展卡尔曼滤波行了积极有益的探索。文献[]用固定增益的卡尔f(k+1)=F((k),(k)+vw(k)曼滤波估计低频运动,而高频运动则用一个参数模系统模型量测|/(k)~N(O,Q(k》)x(k)=H((k)();型来模拟,并用递推增广最小二乘法来估计参数,从u(k)-N(0,Q(k)而估计出船舶的髙频运动。通过控制计算和模拟试初始情况1(0)+5np(0)=验取得了良好的效果。文献[12]提出了水面舰船动力定位控制系统模型参数的离线最速下降寻优的状态估计传递|(k+1)=F((k),()误差协方差传递辨识方法,提高了动力定位系统研制过程的工作效P(h+1)=(h)P()(k)+r(kQ(k )r(k)率。文献[13]在建立船舶三维几何模型基础上,对K()=P(k)HT(k)H()(kH()+R()]-满载船舶从浅水40m到深水500m的水动力系状态估计更新)=()1)((数进行数值计算。利用三维线性势流理论在频域误差协方差更新P)=kk)H(FLK)H()y里研究船舶在浅水中的辐射问题,应用三维源汇分K()R(k)K(k)布法对不同水深下船舶运动的水动力系数,包括附定义φ(k)=0()JH(-)加质量和阻尼系数进行数值计算与分析,得出了有(k)H(k)=0(k)限深水域的附加质量和阻尼系数的渐进特性。文献14]考虑具有修正PM波谱的长峰不规则浪,基于尾部隧道式侧推m,艏部隧道式侧推,艄部方海浪幅值响应算子(RAO研究了船舶在海浪中的六角式推进器。质量阵M可利用文献9]里介绍自由度运动预报模型。为了有效地量化海洋环境对的 Strip Theory计算得到:动力定位船舶的作用,文献[5提出了海洋环境负从/1127400018902-00744载(包括风、海浪和海流)的建模方法,并运用00.07440.1278MATLAB的M文件和SIMUⅠNK分别编制了风干为了得到需要辨识的量需重复进行3项(每项扰力和力矩计算及随机海浪的仿真程序。在三级海2次,共6次)海上试验,以此提高参数估计器的收况下,实现了对海洋环境的仿真,得到了合理的仿真敛性和表现。具体如下结果。文献[1]考虑到船舶的动态特性存在固有的第1项:解耦了的纵荡运动。船舶仅依靠主螺强非线性以及非线性控制改善系统性能和鲁棒性的旋桨山和实现恒速前进,艏向通过艏侧推控制。能力,将非线性控制理论应用到船舶动力定位控制第2项:结耦了的横荡与艏摇运动。通过三个隧道系统的设计中,对某供应船的计算机模型进行仿真,式推进器砌、4、实现两次结耦了的横荡与艏摇运验证了非线性控制系统是有效的。文献[17]提出并动。第3项:在结耦的横荡与艏摇运动中得到方位验证了基于线性核函数在线支持向量回归的模型预角式推进器u的推力系数K6测控制方案。在线支持向量回归算法的引入可以通第1项是为了计算主螺旋桨的推力系数K1和过在线调整,确保预测模型的精确性。Xa,需要的输入量是X本文中X的计算方法是利22船舶数学模型参数辨识用文献[19]里介绍的切片法。第2项是为了计算结文献[18]讨论了使用两个并行测量序列来估计耦了的横荡与艏摇运动的参数数值,可以辨识出的动力定位船舶模型参数的离线并行扩展卡尔曼滤波向量为[ YYNNK3K5]第3项是为了计算全方器算法(O- line parallel extended Kalman filter位推进器的推力系数K6( EKF) Algorithm),见表1。最后采用一项以供给船使用动量方程来代替标准动力学方程,不仅可为对象的全尺度的海上试验来验证提出的参数估计以显著提高状态和参数估计器的性能还具有以下器的收敛性和鲁棒性。优点:实验对象以挪威ABB公司的“ Far Scandia”号供(1)增加数据冗余度;给船为原型。该船总长762m,船宽18:8m,型深(2)降低量测噪声;825m,吃水625m,净吨位4200t,主发动机功率(3)降低环境干扰;3533kW。推进器配置左右舷两个主推进器u1、l2,(4)增加数据记录长度第23卷第3期(总第135期)船羔vd.23N.3012年6月shiP boatJune, 2012(5)以对参数分批进行辨识等手段提高参数辨风。将风速分量定义为识的精确度。L=v,cos(ψ图2显示了实验辨识得到的A和。其中A包(5)W V sin(B-0)含的待求未知量[XyYM而R包含的待求末式中和v分别为风速在X轴和Y轴的分量;v知量是[kk2k3k4k5k6]。和月分别表示风速和风向。如图1所示。假设风速远大于船速,风在纵荡、横荡和艏摇方向的负荷向量可表述为As elements.pAcM(o)V, IV0.5p.A_C-(r )V,V6)0.SpA,Lo C(rm)VV.式中,风的相对角为y=ψp为空气密度,单位e号为kgm3;Lm为船舶总长,单位为m;V为相对风速,103K elements单位为kn;A-和A为正投影面积和侧投影面积,2015单位均为m2;C(y)Cn(y)和C(4)分别为纵荡横荡和艄摇方向的无因次风系数,是通过 Isherwood半经验公式得到的。00003,波浪扰动数学模型波浪干扰力一般分为两种:一种是一阶波浪干图2实验辨识得到的参数曲线扰力,也称高频波浪干扰力。这是在假设波浪为微幅波,未引起船舶大幅摇荡的情况下,认为船舶受到经实验辨识出的动量方程中的量:与波高成线性关系并且与波浪同频率的波浪力。另0.03180种是二阶波浪力,也称波浪漂移力该波浪力与波A000602006l8高平方成比例。0.0075_0.2454这种具有高频率小振幅振荡特性的波浪所产生K=103ding([93,93,20,2.0,28,26]的一阶波浪干扰力最主要是引发船舶的纵摇和垂荡经过计算公式D-M得到运动,对横摇的影响稍次之,而对横荡及艄摇运动的002820影响相对来说就小一些。至于具有慢时变特性的二00.0130475900.081419676阶波浪干扰力,本身同时又是非线性的,它仍然和波写成动力定位模式下的状态空间表达式为:浪的频率有关。波浪的二阶漂移力不但会改变船舶元=AU+Bx(4)疔的航向和航迹,尤其对于在锚泊状态下船舶位置的移动及钻井平台的动力定位系统的工作等均有式中A=MA4M,并且B=MTK。其数值表达式为:重要影响00318000.062800030下面介绍一种估算二阶波浪漂移力方法。19740.0046-0.2428年, Newman提出一种应用频域波浪漂移力系数的0008200082000估算方法。通过把波谱(通常选用PM谱)分为N0∞505-069000108等份,每份有相对应的波浪频率m和波幅A。这样波浪漂移力对横荡、纵荡、艏摇运动的作用力计算公3环境扰动数学模型式为131风扰动数学模型A,(T(W,B=-y)1(W+)风的作用可分为平缓变化的风和快速变化的李文华,等:船舶动力定位系统数学模型参数辨识方法研究[5] Fossen T L. Handbook of Marine Craft Hydrodynamics and式中,T()x0是频域波浪偏移力公式fB是平均Motion Control[M]. Wiley Sons Ltd, 2011: 81-83.波浪方向:是随机的相角。[6] Balchen J G, Jenssen N A, Saelid S Dynamic Positioning可以通过对本估算式进行改变,以避免在数值Using Kalman Filtering and Optimal Control Theory[C]/上产生无物理意义的高频分量。还可对本式进行扩Proceedings of IFAC/IFIP Symposium on Automation in展,用来包括波浪蔓延( wave spreading)。Offshore Oil Field Operation Norway 1976: 183-18633海流扰动数学模型[7]Balchen J G, Jenssen N A Mathisen E, et al. Dynamic作用在海上动力定位船舶上的海流具有方向和Positioning System Based on Kalman Filtering and OptimalControl[J]Modeling, Identification and ControL 1980, 1(3)速度的特征,研究中一般不考虑在大地坐标系下铅135-163垂方向运动。海流分为恒定流和潮汐流。恒定流一般[8] Strand JP, Fossen t inonlinear Passive Observer Design为固定方向和速度的海流,如洋流。潮汐流指海洋for Ships with Adaptive Wave Filtering, In: New Directions因为潮汐运动而引起的海水流动,其典型的表现为in Nonlinear Observer Design(Nijmeijer H, Fossen T L)海流方向的缓慢变化。但对于动力定位来说,海流[M].London: Springer-Verlag London Ld, 1999: 113-134的大小与方向可以认为是确定的,所以海流的模型[9] Guttorm t, Jerome J, Fosset I. Nonlinear Dynamic可以统一按照大小和方向恒定来确立。流的速度分Positioning of Ships with Gain-Scheduled Wave Filtering量表示为5:[C]//The Proceedings of 43rd IEEE Conference orL=V2cos(ψ)Decision and Control, Atlantis, Paradise Island, BahamasDecemher2004:5340-5347ve=y sin(8-n)式中:和v分别为流速在X轴和y轴的分量;V10 i Do K d. Global Robust and Adaptive Output FeedbackDynamic Positioning of Surface Ships[C]/The Proceedings和月分别代表流速和流向。如图1所示。of 2007 IEEE Internati在此没有考虑第摇方向的流速,而海流对水面Automation. Roma, April 2007: 10-14船舶的作用可以通过将各海流速度分量引人到船的1]王晓声船舶动力定位系统设计及试验研究门J国造运动方程中由相对速度向量v=[u-,-a,r丁体现。船,1991(3):12-21[12]边信黔,严渐平,施小成船舶动力定位系统参数辨识4结论方法的研究[J]船舶工程,19994):36-38[13]姜哲,石珦,王磊动力定位船舶水动力参数数值试验本文讨论了船舶及推进器动力学数学模型与船研究[门]实验室研究与搡索,2005(12):14-17.舶外界环境干扰因素数学模型的建模策略。通过对14]李文魁张博田蔚风等.一种波浪中的船舶动力定位已有研究方法的分析研究与总结,有助于建立适用运动建模方法研究[]仪器仪表学报,2007(6):1051于各种海况和操作模式的船舶动力定位系统非线性数学模型。[15]施小成王元慧船舶动力定位海洋环境的建模与仿真J,计算机仿真,2006(11):237-239[16]刘芙蓉陈辉基于非线性控制理论的船舶动力定位控[参考文献制系统的数学模型[〕船海工程,209(5):92-95[1]杜佳璐,张显库汪思源,等船舶动力定位系统的自适[17]邓志良,胡寿松,张军峰船舶动力定位系统的在线模应非线性控制器设计[ C]/proceedings of the2 g chinese型预测控制[门中国造船,2009(6):879Control Conference. Beijing, 2010: 585-589.[2]周利,王磊,陈恒动力定位控制系统研究[船海[18] Fossen T I.Identification of Dynamically Positioned Shipe[].Control Engineering Practice, Volume 4, Issue 3, March程,008,37(2)86-911996:369-376[3]马超庄亚锋陈俊英船舶动力定位系统技术[J中国[19] FaltinsenO M Sea Loads on Ships and Oishore Structures造船,2009,50(增刊):52-57[4]贾欣乐,杨盐生船舶运动数学模型机理建模与数学建[M].Cambridge University Press, 1990:41-45模[M]大连大连海事大学出版社,199:294-356船舶动力定位系统数学模型参数辨识方法研究旧WANFANG DATA文献链接作者:李文华,杜佳璐,张银东,宋健,孙玉清,陈海泉, LI Wen-hua, DU Jia-luZHANG Yin-dong, SONG Jian, SUN Yu-ging, CHEN Hai-quan作者单位李文华,张银东,宋健,孙玉清,陈海泉, LI Wen-hua, ZHANG Yin-dong, SONG Jian, suN Yu-qing, chen Hai-quan(大连海事大学轮机工程学院大连116026),杜佳璐, DU Jia-lu(大连海事大学信息科学技术学院大连116026)刊名:船舶英文刊名:Ship boat年,卷(期):2012,23(3)参考文献(19条1. Balchen J G; Jenssen N A; Mathisen E Dynamic Positioning System Based on Kalmon Filtering andOptimal Control 1980(03)2. Balchen J G; Jenssen N A; Saelid S Dynamic Positioning Using Kalman Filtering and Optimal ControlTheory 19763. Fossen T I Handbook of Marine Craft Hydrodynamics and Motion Control 20114贾欣乐;杨盐生船舶运动数学模型机理建模与数学建模19995.马超;庄亚锋;陈俊英船舶动力定位系统技术2009(增刊)6.周利;王磊;陈恒动力定位控制系统研究[期刊论文]船海工程2008(02)7. Faltinsen 0 M Sea Loads on Ships and Offshore Structures 19908. Fossen t I Identification of Dynamically Positioned Ships 19969.邓志良;胡寿松;张军峰船舶动力定位系统的在线模型预测控制2009(06)10.刘芙蓉;陈辉基于非线性控制理论的船舶动力定位控制系统的数学模型[期刊论文]船海工程2009(05)11.施小成;王元慧船舶动力定位海洋环境的建模与仿真[期刊论文]计算机仿真2006(11)12.李文魁;张博;田蔚风一种波浪中的船舶动力定位运动建模方法硏究[期刊论文]仪器仪表学报2007(06)13.姜哲;石珣;王磊动力定位船舶水动力参数数值试验硏究[期刊论文]实验室硏究与探索2005(12)14.边信黔;严浙平;施小成船舶动力定位系统参数辨识方法的硏究[期刊论文]船舶工程1999(01)15.王晓声船舶动力定位系统设计及试验研究1991(03)Do K d Global robust and Adaptive Output Feedback Dynamic Positioning of Surface Ships 200717. Guttorm T; Jer(o)me J; Fossen T I Nonlinear Dynamic Positioning of Ships with Gain-Scheduled WaveFiltering 200418. Strand J P; Fossen T I Nonlinear Passive Observer Design for Ships with Adaptive Wave Filtering19.杜佳璐;张显库;汪思源船舶动力定位系统的自适应非线性控制器设计2010本文链接http://d.g.wanfangdata.comcn/periodiCalcb201203011.aspx

关于很多高分卫星的数据处理插件，适合ENVI4.8及以上用户，请叫我雷锋哦

三维孔洞储层建模及其地震波场正演模拟，理论讲解很透彻，分析思路清晰1290地球物理学进展26卷预测结果,即从具有确定性资料的控制点(如井点)解释.同时,利用该过程中产生的6冂井的时深关出发,推测出点间(如井间)确定的、唯一的储层参系,用三角剖分网格建立了速度模型并经过井点校数随机模拟是从一个随机函数z(v)中抽取多个可正,实现了对工区构造框架的时深转换通过以波阻能实现,即人工合成反映Z()空间分布的可供选择抗为协变量的孔隙度属性模拟,借助于三维可视化的、等概率的高分辨率实现技术,我们可以大致看到孔洞储集体的形态、分布、对于该工区来说,三维地震资料分辨率较高,对规模及连通性(图1).根据孔隙度发育情况,我们将孔洞储层已经有一定的反映(常表现为低频、不连续强储层分为孔洞欠发育(致密)、孔洞较发育(较致密)振幅反射).通过岩石物理分析又发现孔洞储层低速、和孔洞发育(较疏松)三种类型并分别设计了各自相低密,常规波阻抗反演能够刈其几何形态、空间接触应的弹性参数,同时以模拟出的孔洞形态约束弹性关系定量表征.因此本文将波阻抗数据体作为协变波正演模拟时孔洞体的空间分布量地震属性纵向等值法),采用确定性的协克里金2波动方程正演模拟原理插值算法,对孔洞储层的物性参数进行了三维建模反演所起到的作用,是通过归一化的测井曲线对碳酸盐岩岩溶风化壳孔洞型油气藏属于一种典原始地震数据进行校正,使数据在空间上得到了有型的、复杂的非均质范畴,可以视为由准均匀介质中效的平衡,从而使孔洞反映的更清楚;二是在地质建呈不规则分布的、大小和形状各异的低速体共同组模过程中通过宏观控制,充分利用空间变量的相关成的非层状储集体.在地震剖面上看到的储集体的性,克服低频模型的不足提高属性模拟的分辨率波应是这些低速体的散射(绕射)波.若利用常规波2【V动方程正演模拟方法所使用的均匀介质中的声波方N()程或弹性波方程,难以得到具有复杂非均匀性的孔cline洞型油气藏的地震波场响应2.因此,本文采用非均匀横向各向同性弹性介质中的弹性波波动方程进l()行正演模拟计算,取z轴为垂直对称轴,它可以表示为如下的一阶方程组:a0awta(λ+2naU)+AW(1)7(1(λ1+21)2+λax图1地质模型孔洞储集体俯视图Fig 1 Top view of the cavity reservoiμ(ain geological model其中:(U(x,z,t),W(x,z,t)是速度向量;B(x,z)是密度ρ=g(x,z)的倒数,或者叫疏度;r建模过程中最大的难点是求取准确的速度场,τ(x,z,t),za=x(x,z;t),n=rn(x,z,t)是应木文首先收集整理了工区内6口探井、评价井的钻力张量.A,P/和A1,p分别为水平和垂直方向上井分层数据及多种测井曲线(电阻率、声波时差、密的拉梅系数;为一新的弹性常数可见,在横向各度、自然伽玛等),对其进行了归一化和环境校正,并向同性弹性介质中,独立的弹性常数有五个,它们是制作了合成记录.通过与井旁地震道对比准确标定密度、在垂直方向上的纵、横波速度及纵、横波的各了前中生界侵蚀顶面(15)石炭系双峰灰岩顶面向异性系数,即:(T5。),中下奥陶统顶面(T7)、下奥陶统蓬莱坝组顶面(T7),它们都是区域性的波峰反射在此基础+2uL, Vsi-A pL上采用25m*25m测线密度对该区块6.25km2的3U地震数据体T5°、T5、T74、T78层位进行了精细√λ1+2,CSV闵小刚,等:三维孔洞储层建模其地震波场演模拟1291在具体的有限差分解法上,除了规则网格外,非均匀介质模型的弹性波动方程正演模拟特别是种较为先进的交错网格(图2)最早由 madariaga提当每一个波长中的网格点数多于10个时, Levander出, Virtex在模拟各向同性介质屮SH波和P(1988)2的结果显示,网格色散与网格各向异性均Sⅴ波时也使用了这种差分网格,其精度为o(△2十可忽略不计△x2),在不增加计算工作量和存储容间的前提下,假设U,W分别为介质在x,z两个方向上的速和常规差分网格相比局部精度提高了4倍,且收敛度分量的离散量,R,T,H分别为rxr=和τx的离速度也较快. Levander2又将这种差分网格的精度散量,Lo,M,L1,M1和M2分别为y,kM,A⊥?P1和提高到o(△t2+△x). Crase2则发展了精度可达任g的离散量在(1)式中,各导数项均用中心差分来意阶的高阶交错网格法,但其计算量和内存要求比代替,在如图2所示的一个交错网格中,U,B在节低阶有限差分法大幅度增加.本文使用的是 Virieux点1处计算;W,B在节点2处计算;R,T,M,L,(1986)1的交错网格差分公式,其差分精度为和M1,L1在节点3处计算;H,M2在节点4处计算(△2+△x2), Ikelle l t和 Yung$ K(199)21说这样()式离散为4明该算法可以糈确、稳定地应用于任总复杂变化的=U+B,(R年,-R…)十B,(I1+-H),wH,n-v++B++△(r}一rn)+By2(T+-+),ry=对++(n+2M4)+,△m+-RW(2)+T+,;+(L1+2M1)△tW+U2)+M2△t鲁←z以将震源函数直接赋在rx和n的节点上来模拟震源,即Soure,t)=R(t)t_(source_x, soure_x, t)=R(t.此外,在震源没有激发之前地下介质内部所质点都是静止的,包括质点振动速度为零和所受应力为零.因此,初始条件为图2一个交错树格Fig. 2 A st0,r(x,z,t)=0(t≤0)(3)对于自山表面边界条件,本文采用了模型空间其中,上标k为时间t的离散量,下标i,分别为x的上部加空气的条件,然后再采用吸收边界条件把和z的离散量.△,△x,△z分别为t,x,z的步长空气上边界的弹性波吸收掉,对于空气的下界面,则鉴于 Ricker子波对地震波的分辨率较其它子作正常的分界面来认识,从而获得和实际应用中波函数高,因此,震源选用 Ricker子波,其形式为所采用的地表放炮、地表接收达到一致的效果.R(t)=[1-2Lmf(t-to)] Jexp[-(rf(E to))2]有限差分法在求解波动方程时,会产生不期望式中f表示子波的主频,t为子波持续时间,t为f的数值频散或称网格频散,导致数值模拟结果分辨的函数,在模拟地下激发的地震波时,有限差分交错辛降低2所谓数值频散实质上是一种因离散化求网格中的正应力x和x=是在同一节点上赋值的,解波动方程而产生的伪波动,这种频散既不同于波而vr和vn在节点处的数值并没有参与计算,因此可动方程本身引起的频散,也不同于因波传播的速度1292地球物理学进展26卷频率和角度变化而引起的频散,它是有限差分方法果我们在这里仅分为三种类型:孔涧欠发育(致求解波动方程时所固有的本质特征,无法避免.为了密)、孔洞较发育(较致密)和孔洞发育〔较疏松).消除这种数值频散,前人进行了大量研究,他们的结论是基木一致的,即为了消除数值频散,在使用二阶表1地层框架内各层物性参数有限差分方法时,每个功率对应的波长至少必须使Table I The properties ot each layer用11个网格点,面四阶有限差分则可用二阶差分网in stratigraphie framewor格点数的一半.木文采用的稳定性条件,即计算稳定p(nu/s) v(m/s) (kg/m3)的离散参数区域为151:r4G界面2500三叠系)以Lmd2m≤1(2m-1)fT50界面下石炭系顶)~T46l730≤Ld2m≤T56界面(2m-1)!(巴楚组顶)~T5023102350其中,T74界囿(下奥陶系顶)~T56±8002470T78界面(蓬芠坝组颠)灬T7460002650界面以3702此外,在做波动方程的模型计算时,由于只能在对于试验工区的每条线,其长度均为1625m个有限区域进行,而弹性波在其计算边界上能量为了侏证该区域内均为满叠、孔洞的绕射波收敛以衰诚并不为零,从而产生很强的边界反射,这是模型及边界吸收较为干净,我们在模型的左边延长了计算时所不希望的,需要做人工吸收戌衰减处理,计1200m,右边延长了1575m(延长部分的地层接触算吸收边界的方法有许多种,一般情况下网格周围系并不代表真实情况),即模型总长度为1、4km,的耗散采用质点的速度和应力值乘上一个小于1深度范围为4000~6500m每条线均采用了同样的的因子来平滑的衰减;另种可能性是在网格周围观测系统,具体为:采用零相位对称雷克子波作为震使用低Q值来实现吸收作用,但是后者的吸收效果源(主频40玎z),单边放炮(共20炮,每炮128个检不如前者的吸收效果好,因此本文采用的是第一种波器接收)炮间距160m,检波器间距20m,8次叠方法,具体实珧时釆用了〔 eran等的吸收边界条加,最小偏移距0m,最大偏移距2540m,记录长件实现边界吸收1.6s,Δt=2ms,第一炮的坐标位置为(-1200,0)exp[-a2(I-i)2],1≤i≤1.基于差分稳定条件,取模型中最小介质速度2500m/s其中,I为给定的吸收边界带总节点数;i为吸收边为参考,得到的计算参数为:网格剖分间隔3m界内的节点编号;a为衰减系数,其值的选定与1的3m,时间延拓步长为0.27ms,每个波长(62.5m)大小密切相关,且对吸收效果的影响很大本文中Ⅰ内有20.8个网格.我们一共对33条线进行了正演取为40(即围绕计算区域,再向外设置宽度为40个模拟,图3展示了较为典型的 inline2585线(位于研网格的条形吸收区域)a=0.305/40,i取从0~40究工区的中心部位,地层接触关系以及孔洞体的分节点号.在条形吸收区域中的每个网格结点处,对全布相对比较复杂),从中可以大致看出二维正演模拟部的5个波场量(U,W,R,T,H),在每计算一个时的普遍情况与孔洞体波场响应特征的一般规律问步长后,都做少量的波场减表2展示∫该条线上各孔洞体的几何及物性参3模型计算数,其中④号属于欠发育(充填致密物)类型,①③⑥号属于较发育(充填较致密物)类型,②⑤号属于发在正演之前,我们统计了工区的速度、密度资育(充填较疏松物)类型.此外,建模过程中,我们还料,为了简化模型,并使得孔洞体的地震响应特征更考虑了线与线之间地层起伏渐变、孔洞大小渐变孔具有针对性,我们采用了背景为常速介质、蜜度参数洞物性参数渐变的过程,即所有建模因素都渐变由( arner公式计算的思路(表1).对于孔洞储集的而不是突变,最终保证了三维地震数据体的连体,根据钴井揭示和前面提到的孔隙度属性模拟结续性4期闵小刚,等;三维孔洞储层建模及其地震波场正演模拟1293表2各孔洞储集体的几何及其物性参数最大振幅,且绕射曲率与反射曲率相同,表明二者具Table 2 The geometry and propcrty parameters有不同的传播速度;每个绕射波可分为左右、上下f each cavity reservoir正、反向绕射分支,正向绕射分攴的相位与反射P孔润体尚度宽度vVP波相同,反向绕射分支的相位反转180°,与反射P(m)(m)(m/s)(m/s)(kg/m3)充填物屮心距界面(m)波的相反17396500029002503较致密105弹性波正演模拟生成的炮域合成记录被导人10113480027822470较疏松6FOCUS软件进行常规处理,包括速度谱拾取、动82784500029002500较致密85校、切除、增益、滤波、叠加、偏移和变面积、变密度显①575520030222530致密示等.由于在观测系统中只设计了8次覆盖,为了增⑤18115480027822470较疏松104加速度谱拾取精度,本文采用了由相邻的7个CDP2714850029002500较致密86道集混合构成一个超道集的办法,隔10个CDP拾图4是该模型在590ms时的波场快照,其波场取一个速度文件,并在拾取前先作常规NMO校正清晰,网格频散小,边界吸收较干净这表明,在求解切除,使得原始道集记录能量更强、信噪比史高二维弹性波动方程时,将差分解法和交错网格技术图5、图6分别是TK610井、TK623井所在位置处相结合,通过较好地使用吸收边界和稳定性条件可CDP道集记录及其速度谱,从图中可见各个反射界以显著削弱数值频散,有效地提扃计算精度.同时面的同相轴清晰可辨,对应的能量团集中,而在合成在保证一定的精度前提下,可以采用铰大的空间网记录上T7界面下孔洞所在位置处都有一明显的格间距,提高计算效率.从图巾还可看出,孔洞绕射同相轴,能量团也比较集中,由于TK610升比波和反射波在绕射点处相切,在切点处绕射波具有TK623井孔洞储集体更为发育(尽管二者振幅相1200-80004008001200160020002400280040004505500图3主测线2585地质模型Fig 3 Geclcgical model ol inline258-12004004008001200160045000.10.3图4主测线2585在590ms时波场快照Fig. Snapshot of wave field at 590ms in inlinc25851294地球物理学进展26卷Sg224-230CDP49 SE QNO250030003500400045000.240.60.60.80.8TE1.01. 01.21.2141.4图5TK610井所在位置处CDP道集记录及其速度谱ig. 5 The CDP gather and velocity spectrum at well TK610Sgl58-1640.2ONO250030003500400045000.0.40.60.60.8081.0:1.01212623(2565图6TK623井所在位置处CDP道集记录及共速度谱Fig 6 The CDP gather and velocity spectrum at well TK623当,在地质模型设计时均认为是充填较疏松物,但相消),使得T7界面断断续续,并在该界面下出现TK610井比TK623井在目的层段的厚度要大,横些“短反射”通过仔细分析,我们发现“短反射”中向展布范围也更宽测试产能更高),在合成记录上较强者出现的时间,与孔洞位置相对应.从该模型的孔洞对应的同相轴振幅更强、波形更连续,速度谱上偏移剖面上(图b)可以看到,所有的孔洞体均得到能量团也更强、更集中比铰好的偏移成像,并表现为负正负三个相位的图7是处理完后的叠加和偏移剖面.从叠加剖波形.但鉴于反射波地震勘探的纵向分辨率(大于面上(图a)可以比较清楚的看到孔洞体顶、底的两1/4波长),所有能检测出的孔洞或孔洞组合在叠加组强反射,但是二者之间出现具有绕射特征的弱波剖面上都叠合在T74界面下第一个波峰轴上,在偏代替了成层的背景,这些绕射波的相互下涉(相长、移剖而上都体现在T7界面下第一个黑椭圆体上,4期闵小刚,等:三维孔洞储层建模及其地震波场正演模拟1295601001401802202603003401001401802202603003400.00.00.20.2040.40.60.60.81014露9.926c+089.926e1081022e+091022e+09图?主测线2585对应的叠加剖而(a)和偏移剖面(b)Fig. 7 The stacking section (a) and migration section (b) of inline2585至」其下的“串珠”是孔洞的假象(孔洞组合与围岩(b)之间的多次波及绕射波经偏移归位后形成较强短反L2560L2580L2600射).由于T74界面反射波与沿纵横向有一定分布的孔洞(比较明显的是①、⑤号)的绕射波叠加,使得30001300040孔洞所在位置处T74界面反射波能量变弱,而孔洞底部与下覆围岩之间的正极性反射由于受T7界200600面反射波的负值性续至波叠加,也变得较弱.此外,B40080在构造高点上(④⑥号孔洞体所在位置,④号更为明显),由于孔洞引起的绕射与隆起引起的回转波的相6003600100互丁涉,T56和T7界面不连续,甚至在其间出现空白反射,而实际资料也有这种情况.这说明,对于塔3800800200和油田碳酸盐岩孔洞储集体这类特殊的油气储层来400040001400400说,在解释时遇到层位问断时,不能轻易地开断层,而应该综合考虑构造、孔洞绕射等地震波场特征.这图8联络测线2795实际剖面(a)和正演剖面(b)对比也是塔河油出勘探开发实践中发现“表层弱反射、内Fig 8 Comparison of the actual section (a) and幕强反射”地震特征对应有利储层的一个佐证forward modeling section (b)in crossline2795依据以上思路与工作流程,我们得到了33条沿主测线方向的二维偏移剖面,在并成三维体之前,为而正演模拟釆取的是8次叠加、道间距10m),正演了尽量消除线与线之问因地层起伏造成的不闭合,剖面较好的反映了实际情况.这不仅体现在层位的我们采取先把33条线的速度文件并成三维体,整体形态、分布比较相似(由于速度取了平均,各层的厚平滑后两用每条线对应的、平滑后的速度对其原始度不一致,但不影响我们的主要的,即对孔洞体地共中心点数据进行动校、叠加、偏移的办法,得到33震响应特征的分析),更重要的是,我们所设计的孔条新的二维偏移剖面,再并成一新的三维体,此外,泂体,其位置、形状规模、振幅强弱均与实际地震资由于正演模拟数据体线间距为50m,道间距为料具有相当好的对应关系,这表明我们在止演模拟10m,其空间采样率比实际资料低,本文编制了相和处理时的设计思想和参数选取原则是合理的,这应的算法在频率域对其进行插值,使线间距加密到结果也为我们进行后续工作提供了比较好的数据25m图8是联络测线2795在时间域的实际剖面源由于实质上是2.5维,不是基于面元的真三维,(a)和正演剖面(b)对比,排除二者在采集时的一些所以沿联络测线的剖面上同相轴有抖动现象,这是差异(如实际三维采集资料为24次叠加道间距25m,不可避免的)1296地球物理学进展26卷4结论与建议2]谢桂生,刘洪,赵连功,伪谱法地震波正演模拟的多线程并行计算[冂.地球物理学进展,2005,20(1);17~23.本文从三维角度,建立了与实际资料比较吻合Xie G S, Liu H, Zhao L G. Parallel Algorit hm based on the的孔洞储层模型,并进行了弹性波正演模拟,总结了multithread Technique for pseudospectal modeling of seismic地震响应规律,主要结论如下:wave[J]. Progress in Geophysics(in Chinese), 2005, 20(1)1)结合地震资料建立储层地质模型能够有效[3]刘财,张智,邯志刚,等.线性粘弹体屮地震渡场伪谱法模拟地降低储层模型的不确定性,提高建模精度.同时利技术[门].地球物理学进展。:005,20(3),640~644,用协克里金技术,用波阻抗反演的确定性信息约束Liu C, Zhang Z, Shao Z G, et aL. Pseudo-spectral forward储层的平面非均质性,可以实现孔隙度属性的确定modeling nf seismic wave in linear viscoelasic solid [J]P1性建模),2005,20(3):640~644.4」张智,刘财,邵志刚,伪谱法在常Q粘弹介质地震彼场模拟(2)在求解二维弹性波动方程时,将差分解法和中的应用效果[].地球物理学进展,2005,20(4):945交错网格技术相结合,通过较好地使用吸收边界和949,稳定性条件可以显著削弱数值频散,有效地提高计Zhang Z, Liu C, Shao G. The application of pseudo-spectral算精度.同时,在保证一定的精度前提下,可以采用forward modeling of seismic wave field in constant Q较大的空间网格间距,提高计算效率该方法具有广viscoelastic medium [J]. Progress in Geophysics, 2005,20(4)945~949泛的适用性5]盖良国,马在出,曹景忠,等.一阶弹性波方程交错网格高阶(3)孔洞储集体在偏移剖面上表现为负-正-负差分解法[冂].地球物理学报,200,43(3):411-~419三个相位的波形,但只能确定奥陶系风化面下第Dong LG, Maz T, Cao j Z, et al. A staggered-grid high个负相位是孔洞的发育位置,其下的“串珠”是孔洞order difference method of one-order elastic wave equation]的很象.风化面反射波与沿纵横向有一定分布的孔Chinese J. Geophys. (in Chinese),2000,43(3):411-419洞体的绕射波叠加,使得孔洞所在位置处风化面反[6]董艮国,马在田,曹景忠,一阶弹性波方程交错网格高阶差分解法稳定性研究[门].地球物理学报,200,43(6):856~射能量变弱,而孔洞底部与下覆闱岩之间的正极性反射由于受风化面透射波的负值性续至波叠加,也Dong L G, Ma Z T, Cao J Z. a study on stability of the变得较弱.该结论对于实际地震资料处理、解释以及staggcred-grid high-order difference method of first-order储层预测烃类检测具有普遍的指导意义elastic wave equation. Chinese J. Gcophys. in Chinese)2000。43(6);856~864本文不足之处主要有三点「7]萤良国.复杂地表条件下地震波传播数值模找1.勘探地球(1)在三维孔隙度建模时采用的是常规阻抗信物理进展,2005,28(3);187~194息(约束稀坑脉冲反演),其纵向分辨率不够(只能分Dong L G. Numerical simulation of seismic wave propagation辨1/4波长以上的孔洞储集体),在后续工作中将尝under complex near surface conditions [J]. Progress in试使用地质统计学反演的阻抗体来约束建模以大幅Exploration Geophysics(in Chinese), 2005, 28(3):187--194提高纵向分辨能力[8奚先,姚姚,二维随机介质及波动方程正演模拟[.石油地球物理劫探,2001,36(5):546-552(2)在弹性波正演模拟时采用的是2.5维思想XiX, Yao Y. 2D random media and wavc cquation forward口前正在研制全三维算法有望更逼真的还原孔洞储modeling [J]. Oil Geophysical Prospecting in Chinese集体的真实地下情况001,35(5);546~5523)考虑到缝的各向异性更为复杂,本文尚未涉9]奚先,姚姚,二维粘弹性随机介质中的波场特征分析[刀地及,对于碳酸盐岩中这类油气运移的重要通道,将在球物埋学进展,2004,19(3):608~615今后的工作中进一步研究Xi x, Yao Y. The analysis of the wave field characteristics in2-D viscoelastic random medium LJ. Progress in Geophysics参考文献( References):hinese),2004,19(3):608~[10]奚先,姚姚,二维横各向同性弹性随机介质中的波场特征1]刘文岭.大庆宋芳屯油田芳2区块地震与地质资料综合储层J.地球物理学进展,2004,19(4):924~932地质建模研究(博土论文儿D1.北京:中国地质大学,2002Xi x, Yao Y. The wave field characteristics of 2-DLiu W I. A Study on Reservoir Geological Modeling withclo].ESeismic and Well-log Data in Fang 2 Area of DaqingGeophysics(in Chinese), 2004,19(4):924-932ongfangtun Oil Field (doctor dissertation)(in Chi[111吴永国,贺振华,黄德济.串珠状溶涧模型介质波动方程正Beijing: CUG, 2002.演与偏移[.地球物理学进展,2008,23(2);539~5444期闵小刚,等:三维孔洞储层建模及其地震波场正演模拟1297Wu Y G, He Z H, Huang d J. Wave equation forward[19]肖玉茹,何峰煜,孙义梅,等,古洞穴型碳酸盐岩储层特征modeling and migration for heads-shaped corroded cave model研究一以塔河油田奥陶系古洞穴为例匚门。石油与天然气地EJ]. Progress in Geophysics(in Chinese), 2008, 23(2): 539质,200324(1):71~80.Xiao YR, He f Y, Sun Y M, et al. Reservoir charactetistics12]股文,印兴耀,吴国忧.高特度频率域弹性波方程有限差分of paleocave carhonates-a casc study of Ordovician paleocave方法及波场模拟[」.地球物理学报,2006,49(2):561in tahe oilfield, Tarim basin UJ]. Oil Gas Geology(inChinese),2003:24(1):71-80Yinw, YinXi,WuGC. The method of finite difference of[20]姚蟋,唐文榜.深层碳酸盐岩岩溶风化壳洞缝型油气藏可检high precision elastic wave equations in the frequcncy domain测性的理论研究[门.石油地球物理勘探,2003,38(6):623and wave-field simulation [J. Chinese J, Geophys.629Chinese),2006,49(2):561~568.Yao Y, Tang W B. Theoretical study of detectable cavern[13]马贵,土尚旭,宋建勇.频率域波动方程正演中的多网格Fractured reservoir in weathered Karst of dccp carbonatite迭代箅法[门].石油地球物理勘探,2010,45(1):15[J]. Oil Geophysical Prospecting(in Chinese), 2003,38(6):Ma ZG, Wang S X. Sun J Y. Multigrid iterative algorithm in623~629,domain wave equation forward modeling [J]. Oil [21] Levander A R. Fourth-order finite difference P-SvGeophysical Prospecting(in Chinese ) 2010, 45(1): 1-5seismograms []. Geophysics, 1988, 53(11): 25-36.[14]张金海,王卫民,赵连锋,等.傅里叶有限差分法三维波动[22] Crase e. Iligh- order( space and timc) finite-difference方程正演模拟[.地球物理学报,2007,50(6):1854A, In: 60th SEG Annual1862C].1990:987~991.Zhang j H, Wang W M, Zhao L F, et aL. Modeling 3-D [23] IkelleL T, Yung SK, Daube F. 2-D random media withscalar waves using the Fourier finite-difference method.ellipsoidal autocorrelation function [J]. Geophysics, 199350(6):1854[24]奚先.随机介质模型的构造及其波场模拟(博土论文)[D][15] Qin Z, Lt武汉:中国地质mproved NPML absorbing boundary condition in elastic waveXix. Construction and scismic wave field modeling ofmodeling [J]. Applied Geophysics, 2009. 6(2): 113-121random medium model doctor dissertation ) in Chinese)[16][D].Wwave equation [J. Geophysics, 1986, 51(1): 54-61[25]吴国忱,王华忠.波场模拟中的数值频散分析与校正策略[17] Virieux J. P-Sv wave propagation in heterogeneous mediaLJ.地球物理学进展,2005,20(1):58-65velocity-strcss finite-difference methud LJ]. GeophysicsWu GC, Wang H Z. Analysis of numerical dispersion in1986,V51;889~901.wave-field simulation [J]. Progress in GreaphysiEs ( in18] Igel H, Riollet B. Mora P. Accuracy of staggered 3-D finiteChinese),2005,20(1):58~65difference grids for anisotropie wave propagation [J]. 62th [26] Cerjan C, Kosloff D, Kosloff R, et al. A nonreflectingAnn, Internat, Mtg, Soc. ExpL. Geophys, Expboundary condition for discrete acoustic and elastic- wav1992,1244~1246.equation []. Gcophysics, 1985, 50(4): 705-708.

工具栏图标、文件图标素材打包下载。程序设计用图标

xilinx的fpga+pcie数据采集卡，包括linux驱动以及测试程序

UDS_PCAN_APIA应用程序，整车网络诊断应用程序，超值！（PEAK CAN UDS Application Programming InterfaceUser Manual.pdf）PCAN-UDS APi- User ManualContents1 PCAN-UDS API Documentation2 Introduction2.1 Understanding PCAN-UDS2.2 Using PCAN-UDS2.3 Features7888992.4 System Requi rements2.5 Scope of supply3 DLL API Reference3.1 Namespaces103.1.1 Peak Can uds3.2 Units3.21 PuDs Unit3.3 Classes3.3.1 UDSApi3.3.2 TUDSApi3. 4 structures1022334553.4.1 TPUDSMsg3.4.2 TPUDSSessionInfo3.43 TPUDSNetAddrinfo3.5 Types213.5.1 TPUDSCANHand]e223.5.2 TPUDSstatus233.5.3 TPUDSBaudrate253.5.4 TPUDSHWType283.5.5 TPUDSResult303.5.6 TPUDSParameter313.5.7 TPUDSService393.5.8 TPUDSAddress423.5.9 TPUDSCanId443.5.10 TPUDSProtoco l463.5.11 TPUDSAddressingType483.5.12 TPUDSMessageType493.5.13 TPUDSSVCParamDSC503.5.14 TPUDSSVCParamER513.5.15 TPUDSSVCParamcc533.5.16 TPUDSSVCParamTP543.5.17 TPUDSSVCParamcdTCS543.5.18 TPUDSSvCParamROE553.5.19 TPUDSSvCParamROERe commendedserviceID573.5.20 TPUDSSVCParamLC583.5.21 TPUDSSvcParamLCBaudrateidentifier593.5.22 TPUDSSVCParamDI603.5.23 TPUDSSVCParamRDBPI643.5.24 TPUDSSVCParamDDDI653,525 TPUDSSyCParamRDTCI66PCAN-UDS APi- User Manual3.5.26 TPUDSSVCParamRDTCI DTCSVM6935.27 TPUDSSYCParamIOCBI703.5.28 TPUDSSvCParamRC3.5.29 TPUDSSVCParaMRC RID723.6 Methods733.6.1 Initialize753.6.2 Initialize(TpudsCanhandle, tpudsbaudrate)3.6.3 Initialize(TPUdsCANhandle, TPUdSBaudrate, TPudSHWType, UInt32,UInt16)83.6.4 Uninitialize813.6.5 Setvalue843.6.6 Setvalue (TPUdsCanhandle, tpudsparameter, UInt32, uint32)843.6Setvalue (TPUdSCaNHandle, TPUDSParameter, stringBufferUint32)873.6.8 Setvalue (TPUDSANHandle, TPUDSParameter, Byte[], Uint32)883.6.9 Setvalue(Tpudscanhand le, tpudsparameter, IntPtr, UInt32)3.6.10 Getvalue933.6.11 Getvalue (TPUDSCANHandle, TPUDSParameter, StringBufferUint32)933.6. 12 Getvalue (TPUDSCANHandle, tpudsparameter, uint32, Uint32)963.6.13 Getvalue (TPUDsCaNHandle, TPUDSParameter, Byte l], UInt32)993.6. 14 Getvalue (TPUdSCAnhandle, tpudSParameter, Intptr, UInt32)1013.6.15 Getstatus1043.6.16Read1073.6.17 Write3.6.18 Reset1143.6.19 WaitForsing lemessage1163. 6.20 WaitFormultiplemessage1203.6.21 Waitforseryice1263.6.22 WaitForservicefunctional1303.6.23 ProcessResponse1333.6. 24 SvCDiagnosticsessioncontro l1383.6.25 SVCECUReset1413.6.26 SvcSecuri tyAccess1453.6.27 SvCCommunicationControl1483.6.28 SvcTesterpresent1523.6.29 SvcsecuredDataTransmission1553.6.30 SvcControlDTCSetting1583.6.31 SvcResponseonEvent1623.6.32 SVCLinkcontrol1663.6.33 SVCReaddatabyidentifier1703.6. 34 SvcReadMemory ByAddress1733.6.35 SvcReadscal ingdatabyidentifier1773.6. 36 SvcReadDataByperiodicIdentifier1803.6.37 SvcDynamicallydefinedataIdentifierDBID1843.6.38 SvcDynamicall ydefineDataIdentifierDBMA1883.6. 39 SvcDynamical lyDefineDataIdentifierCDDDI1933.6.40 SvcWri teDataByidentifier1973.6. 41 Svcwri teMemory byaddress2003.6.42 SvcClearDi agnosticInformation2053. 6. 43 SVCReadDTCInformation2083.6.44 SvCReadDTCInformationRDTCSSBDTC2113. 6. 45 SvCReaddTCInformationRDTCSSBRN215PCAN-UDS APi- User Manual3. 6.46 SVcReadDTCInformationReportExtended2183.6. 47 SvcReadDTCInformationReportseverity2213,648 SvcReaddTCInformationrsIodtc2253. 6.49 SvCReadDTCInformationNoParam2283.6.50 SvcInputout put contro byidentifier2323. 6.51 SyCRoutineControl2363.6.52 SvCReques tOwn load2393.6.53 SvcRequestUp load2433. 6.54 SVCTransferData2483.6.55 SvCRequestTransferExit2513.7 Functions2563.7.1 UDS Initialize2583.7.2 UDs Uninitialize2593.7.3 UDs Setvalue2603.7.4 UDs Getvalue2613.7.5 UDS Getstatus2623.7.6 UDS Read2643.7.7 UDs Write2653.7.8 UDs Reset2663.7.9 UDS_WaitForsinglemessage2673.7.10 UDS_waitForMultipleMessage2693.7.11 UDs Wai ce2723.7.12 UDS WaitForserviceFunctional2733.7.13 UDS_ Processresponse2753.7.14 UDS_SvcDiagnosticSessionControl2773.7.15 UDS SVCECUReset2783.7.16 DS_SVCSecuri tyAccess2803.7.17 UDS SVCCommunicationcontrol2813.7.18 UDs SvCTesterpresent2833719 UDS SvCSecuredDatatransmission2843.7.20 UDS_SvCControlDTCSetting2863.7.21 UDS_SVCResponseonEvent2873,7.22 UDs SVCLinkcontrol2893.7.23 UDS_SvcReaddatabyidentifier2913.7.24 UDS_SvcReadMemory byAddress2923.7.25 uDs_ SvcReadscalingdatabyidentifier2943.7.26 UDS_SvCReadDataBy Periodi iDentifier2953.7. 27 UDS_SVcDynamical l yDefineDataIdentifierDBID2973.7.28 UDS_SvcDynami call ydefinedataIdentifierDBMa2993.7.29 UDS_SvcDynami cal l yDefineDataIdentifierCDDDI3013. 7.30 UDS_SvcWriteDataByIdentifier3023,7.31 UDs SvcWri teMemorybyaddress3033.7. UDS_SvcClearDiagnosticInformation3053.7.33 UDS SVCReadDTCInformation3073.7. UDs SyCReadDTCInformationRdtCSSBDTC3093.7.35 uDs SvCReadDTCInformationRdtcssbrn3103.7.36 UDS_ SvCReadDTCInformationReportExtended3113.7.37 UDS_SvcReadDTCInformationReportseverity3133.7.38 UDS SVCReadDTCInformationRSIODTC3153,739 UDS SVCReadDTCInformationNoParam3163. 7.40 UDS_SvcInputoutput contro l byIdentifier3,7. 41 UDs SyCRoutinecontrol319PCAN-UDS APi- User Manual3.7.42 UDS_SvcRequestDown load3213.7.43 UDS_ SVCRequestupload32337.44 UDS SyCTransferData3253.7.45 UDS_SVCRequestTransferExit3263.8 Definitions3293.8.1 PCAN-UDS Handle Definitions3293.8.2 Parameter value defintions3313.8.3 TPUDSMsg Member value Definitions3323.8.4 PCAN-UDs Service parameter Definitions3334 Additional Information3354.1 PCAn Fundamentals33542 PCAN-Basic3364.3 UDS and ISO-TP Network Addressing Information3384.3.1 ISO-TP network addressing format3384.4 USing Events3405 License Information3426PCAN-UDS APi- User Manual1 PCAn-UDS APi DocumentationWelcome to the documentation of PCan-UD APl, a PEAK CAN API that implements ISo 15765-3, UDS in CANan international standard that allows a diagnostic tester(client) to control diagnostic functions in an on-vehicleElectronic Control Unit(ECU or serveIn the following chapters you will find all the information needed to take advantage of this aPlIntroduction on page 8DLL API Reference on page 10Additional Information on page 335PCAN-UDS APi- User Manual2 IntroductionPCAN-UDS is a simple programming interface intended to support windows automotive applications that usePEAK-Hardware to communicate with Electronic Control Units(ECU) connected to the bus systems of a car, formaintenance purpose2.1 Understanding PCAN-UDSUDS stands for Unified Diagnostic Services and is a communication protocol of the automotive industry. thisprotocol is described in the norm iSo 14229-1The UDS protocol is the result of 3 other standardized diagnostic communication protocolsIS0 14230-3, as known as Keyword 2000 Protocol(KWP2000L IS0 15765-3, as known as diagnostic on CANISo 15765-2, as known as ISo-TPThe idea of this protocol is to contact all electronic data units installed andCAN OBDninterconnected in a car, in order to provide maintenance, as checking for errors,actualizing of firmware, etcUDS is a Client/Server oriented protocol. In a UDS session(diagnostic session ),aprogram application on a computer constitutes the client(within UDS, it is calledPCAN-UDSTester), the server is the ecu being tested and the diagnostic requests from client toserver are called services. The client always starts with a request and this ends with apositive or negative response from the server(ECuSince the transport protocol of UDS is done using ISo-TP, an international standardPCAN ISOTPfor sending data packets over a CAN Bus, the maximum data length that can betransmitted in a single data-block is 4095 bytes.PCAN-UDS API is an implementation of the Uds on CAN standard the physicalcommunication is carried out by PCAN-Hardware (PCAN-USB, PCAN-PCI etc )throughPCAN-Basithe pCAN-ISo-TP and PCAN-Basic API (free CAN APls from PEAK-System). Because ofthis it is necessary to have also the pCAN-1S0-tP and PCAN-Basic APls(PCAN-ISO-TP. dll and PCAN Basic. dll) present on the working computer where UdS is intended tobe used. PCAN-UDS, PCAN-ISO-TP and PCan-Basic apis are free and available for allFigure 1: Relationship of thepeople that acquire a pCAn-hardware2.2 Using PCAN-UDSSince PCAN-UDS API is built on top of the PCAN-1So-TP API and PCAN-Basic APls, it shares similar functions. Itoffers the possibility to use several PCAn-UDS (PUds) channels within the same application in an easy way. Thecommunication process is divided in 3 phases: initialization interaction and finalization of a puds-channelInitialization In order to do UDS on CAN communication using a channel, it is necessary to initialize it first. Thisis done by making a call to the function UDS_ Initialize (class- method: InitializePCAN-UDS APi- User ManualInteraction: After a successful initialization a channel is ready to communicate with the connected can bus.Further configuration is not needed the 24 functions starting with UDS Svc(class-methods: starting with Svccan be used to transmit UdS requests and the utility functions starting with Uds WaitFor(class- methodsstarting with WaitFor) are used to retrieve the results of a previous request. the Uds read and UDS Write(class-methods: Read and Write are lower level functions to read and write UDs messages from scratch. Ifdesired, extra configuration can be made to improve a communication session, like service request timeouts orISo-TP parametersFinalization: When the communication is finished, the function UDS_ Uninitialize(class-method: Uninitializeshould be called in order to release the puds-channel and the resources allocated for it. In this way thechannel is marked as free"and can be used from other applications23 FeaturesI mplementation of the UDS protocol(iSo 14229-1)for the communication with control unitsWindows DLLs for the development of 32-bit and 64-bit applicationsPhysical communication via Can using a Can interface of the pcan seriesUses the pcan-Basic programming interface to access the can hardware in the computerUses the pCAn-ISo-TP programming interface(iso 15765-2)for the transfer of data packages up to 4095bytes via the can bus2.4 System Requi rementsL- Windows 10, 8.1, 7(32/64-bitAt least 512 Mb ram and 1 GHz CPUPC CAN interface from peak-SystemPCAN-Basic APlL PCAN-SO-TP API2.5 Scope of supplyInterface DLL, examples, and header files for all common programming languagesDocumentation in pdf formatDocumentation in HTML Help formatPCAN-UDS APi- User Manual3 DLL API ReferenceThis section contains information about the data types (classes, structures, types, defines enumerations)andAPI functions which are contained in the pcan-uds api3.1 NamespacesPEAK offers the implementation of some specific programming interfaces as namespaces for the. NEtFramework programming environment. The following namespaces are available:NamespacesNameDescription}PeakContains all namespaces that are part of the managed programming environment fromPEAK-SystemPeak CanContains types and classes for using the PCan aPi from PEAK-SystemPeak Can. LightContains types and classes for using the PCAn-Light API from PEAK-SystemPeak Can basicContains types and classes for using the pcan-Basic APl from PEAK-SystemPeak Can CcpContains types and classes for using the CCP API implementation from PEAK-SystemPeak Can XcpContains types and classes for using the XcP aPi implementation from PEAK-SystemPeak Can. Iso TpContains types and classes for using the pCAN-IS0-TP aPl implementation from PEAKSystelPeak Can, UdsContains types and classes for using the PCan-UDS API implementation from PEAK-SystemPeakCan.Obdll Contains types and classes for using the PCAN-OBDIll API implementation from PEAKSystemt}Peak. LinContains types and classes used to handle with lin devices from PEAK-Systemt}Peak. RP1210AContains types and classes used to handle with can devices from PEak-System through theTMC Recommended Practices 1210, version A, as known as RP1210(A3.1.1 Peak Can UdsThe peak Can. Uds namespace contains types and classes to use the pcan-UdS aPi within the. NET Frameworkprogramming environment and handle pcan devices from peak-SystemRemarks: Under the delphi environment, these elements are enclosed in the puds-Unit. the functionality of allelements included here is just the same. the difference between this namespace and the delphi unit consists inthe fact that delphi accesses the Windows api directly it is not managed code)AliasesAliasDescriptionTPUDSCANHandle Represents a pCAn-UDS channel handleClassesClassDescription像曰UDSApiDefines a class which represents the PCAN-UDS API10

本程序用于求解多车型多目标下的车辆路线问题，程序中考虑了两种车型，建立的目标函数是车辆总运营成本最小，考虑的约束有容量约束、最大行驶距离约束和时间窗约束，采用的优化算法是遗传算法，程序内部有详细的注释，方便修改。

本文档是根据网上零散的FreeRtos嵌入式系统资料，整理成一个完整的FreeRtos系统嵌入式开发的文档，总共含附件分38章节，上传分享给大家

里面包含随机森林的matlab实现代码，并且有简单的训练数据集和测试数据集