建立油气系统模拟用的地质格架时，一般需要一个地表高程图（地貌）。如果探区内探井数量足够多，则可以用探井的海拔数据成图。不过，经常遇到的情况是，探区内的探井数并不多，而且有时用海拔数据成图并不理想。实际上，美国国家地球物理中心（NGDC）已经为我们提供了一套精度为一弧分（1弧度的1/60）的全球地貌数据（ETOP01），用户可以通过http://www.ngdc.noaa.gov/mgg/gdas/gd_designagrid.html，下载相应区域范围的地表高程数据。比如，通过GoogleEarth，我们大概知道松辽盆地的经纬度范围，则可以下载到对应的高程数据，其数据格式为： 经度，纬度，高程。 通过相应的坐标投影转换，则可以得到地质格架所需要的地表高程。比如，这里我们利用Surfer 9.0的坐标变换功能，步骤如下：
在surfer 中打开NGDC上下载来的高程数据，在Data中选择 new projected coordinate，在source columns 中选择 World Geodetic System 1984 (wgs84)投影系统，在target columns 中选择目标区所用的投影参数，比如Beijing 54(北京54)。完成后保存数据，用surfer 成图(网格化时选择 Modified shepard’s method)，结果如下图。
Posted May 10th, 2010 by Trevor ||翻译：Zeng, H.S
冷冰冰的地图如何变得妙趣横生？Zotero的开发人员无时无刻不在寻找对资料库进行创新式的管理与可视化的方法。这样，将资料库反映出的地理信息标在地图上就进入我们的视线。下面，我们将兴奋向Zotero粉丝们宣布Zotero的新插件：地图可视化插件（Zotero maps），这个插件目前可以从这里下载。Zotero地图插件是由Entropy Free LLC 和Zotero项目组共同开发完成的。地图插件可以向用户们展示他们的资料之间的地理信息及其关系，并在地图上标出这些信息和联系。反过来，Zotero也可以用创建的地图作为资料的导航。地图插件的效果可见下图。
下面我将向大家介绍如何在Zotero中使用地图插件：
(more…)
本工具适用于双狐图形数据（2007版Ascii格式）查看，导出；导出的数据可直接用于Trinity,PetroMod等盆地模拟软件中。
软件操作方法：
直接将双狐文件拖到软件界面，df2data将以列表的形式显示双狐图形数据的分层结构，用鼠标右键点击想要导出的数据并选择正确的格式即可。
软件下载：df2data1.4
软件截屏：
官方简介
GCDPlot是一个进行地球化学图解、CIPW计算的 Microsoft Excel VBA 程序，GCDPlot 可以使用 Excel 的数据以标准的 Excel 图表格式绘制各类地球化学图解，另外，GCDPlot还可以对各类图解进行增、删、修改等操作，生成工具所没有附带的图解。功能和内容和 NewPet，MinPet（2个地化数据处理、图解软件）相当。GCDPlot目前为开源软件。
特点
- 纯 Microsoft Excel VBA 加载宏程序
- 直接应用 Excel 数据作图，无须数据转换
- *输出是标准的 Excel 图表，可以在 Excel 中进行修改、修饰、输出、打印等操作
- 程序可以绘制的图解可以任意修改、增删
- 程序提供图解修改工具
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
附件是两个简单的地球化学图版定义文件示例（干酪根显微组分分类三角图.ini和热解S1+S2与有机碳判断源岩质量图.ini），使用时将它们复制到 GCDPlotdig 下面即可，并在GeoChemistryChart.ini中添加相关的图版信息即可，格式可参见文件中的其它定义。自定义一个新的石油地球化学图版可用软件中的config。作图效果见下图：
点击下载GCDPlot修改版包含石油地球化学图版
石油地球化学图版更新状态（2007年11月9日）：
.干酪根显微组分分类三角图
. 甲烷碳同位素-C1/(C1+C2)判断天然气成因（更正：应为 C1/(C2+C3),请自行在Config中修改 2007年12月5日）
.H/C-O/C范氏图,包括三类三分,三类四分,三类五分
.HI-Tmax有机质类型划分图版_三类四分
.热解S1+S2与TOC判断源岩质量图
.有机与无机二氧化碳鉴定图版
.d13CH4(甲烷碳同位素，下同)-CH4(%)识别甲烷成因图版
.d13C1-d13C2-d13C3不同成因有机烷烃气鉴别图版（戴金星，V型鉴别图版）
.d13C2 与d13C2-d13C1天然气成因鉴别图版
(more…)
Simple Ro 是由日本学者Suzuki（1993）[1]提出的计算镜质体反射率的模型。该模型认为在每个温度点下，化 学反应是单一的（相对于平行反应模型），其对应的活化能称为“表观活化能(apparent activation energy，Eapp)”。表观活化能随Ro的增加而增大，其关系可以用Eapp=alpha*ln（Ro）+beta来表示。Simple Ro的优点的在于可以通过拟合不同的alpha和beta，将现有的镜质体反射率模型转换成Simple Ro模型，也可以用实测Ro数据来拟合alpha及Beta得到新的Ro模型。事实上，不同盆地沉积的镜质体由于化学成分和化学结构不 同，其反射率的演化路径也不完全等同于Easy%Ro（目前 广泛应用的Ro校正模型），如富氢镜质体，受到超压作用的镜质体等[2]。Ro模型作为目前热史校正的重要参考，其可靠性直接影响到最终的烃 源岩评价。下面提供的代码是以2C/Ma为加热速率，计算1~300C温度区间Ro的变化，结果输出到文本文件，simpleRo.txt。图1是用OrgPlus（The next release of Z-Kinetics）生成的同等条件下（2C/Ma）Easy%Ro与用Easy%Ro拟合的Simple Ro模型比较，从图上可以看出两者计算结果相差很小。
Ref
[1]Noriyuki Suzuki, Hideki Matsubayashi, et.al., A Simple Kinetic Model of Vitrinite Reflectance, AAPG Bull, 1993, 77:1502-1508.
[2]Buchardt,B., and M.D .Lewan, Reflectance of Vitrinite-like macerals as a thermal maturity index for Cambrian-Ordovician Alum Shale, Southern Scandinavia: AAPG Bull,1990, 74:394-406.
Code here:
#Test on Python 2.6 By zenghuasen
#License： Free
import math
import os
alpha=9.7
beta=54.1
PRO=0.2
ROO=0.2
Ea=0
s=0
F=0.85
FC=0
RO=0.2
oFile=open("simpleRo.txt",'w')
for i in range(1,301):
#Ea=alpha*ln(RO)+Beta
Ea=alpha*math.log(PRO)+beta
Ea= 4186.8*Ea
#duration time (0.5Ma )as second
DT=0.5*1.0e+6*365*24*60*60
s=s+1.0e+13*math.exp(-Ea/(8.314*(i+273.15)))*DT
if(s>500):
FC=F
print(str(i)+ " " + str(FC))
else:
FC=F*(1.0-math.exp(-s))
RO=math.exp(math.log(ROO)+3.7*FC)
oFile.write(str(i)+"\t"+str(RO)+"\n")
PRO=RO
oFile.close()
print("Press any key to quit")
okExit=raw_input()
You can consider OrgPlus as a free alternative for Orgas(BP oil and gas generation analysis kit, now incorporate to MPath) or Kinex(Zetaware Inc.). The tool was compiled based on Net framework(2.x or higher) .
The download link here: get OrgPlus
