約翰留著長長的鬍子 作品

第三十八章,照著答案挑解題過程

    1963年,碩士畢業的約翰·斯旺森進入西屋公司宇航實驗室工作。他發現西屋公司的工程師雖然身邊就有計算機,但在進行有限元分析的時候都是用手工計算,而沒有想到利用計算機輔助完成。斯旺森自己琢磨,打算開發一種通用工程軟件來進行自動化的有限元分析,但西屋公司拒絕了約翰·斯旺森的想法。



    覺得呆在西屋公司沒搞頭,斯旺森在1969年離職,自己創辦了一個公司:alysis Sys系統公司(SASI)。



    1970年,約翰·斯旺森推出了初代有限元分析軟件。那時候SASI公司窮得還買不起大型計算機,有限元分析軟件開發成功之後,試運行是找IBM按小時租用的計算機來運行的。不過SASI馬上就成功地把這個軟件賣回給了西屋公司,賺到了老東家的錢。



    SASI也就是後來的ANSYS(安思)。



    有限元分析軟件和集成電路設計軟件是1953年交給上海計算機所軟件中心的兩個重點項目。本來集成電路設計軟件是頭號重點,但回國路上看一列列的油罐車備受刺激,唐華必須得提前用上有限元分析。



    唐華拿有限元分析軟件用來……學術造假。



    當然學術造假這說法略有些誇張,唐華只是拿到了標準答案,然後照著答案挑選解題過程而已。



    地質學這個學科直到21世紀還是大片的未知之地。學地質的,在做碩士畢業設計時就會遇到導師不知道、學校也不知道,上網一查全世界都沒人研究過的領域。



    在50年代,關於中國大地構造和成油理論有兩個學說分支,到21世紀,分支成功地增加到了兩位數……



    所以現在唐華整理整理手頭的書籍論文,就能發掘出好幾種模型構造和推演方法。甚至只按照李四光30年代到50年代的著作,就能總結出4種分析法。



    唐華一種方法一種方法地試,用超算先跑個低精度的運算,等三四個小時,結果出來之後看中國地圖上的柱狀圖。當跑到第四種方法時,柱狀圖的分佈、柱子長短和他記憶中幾大油田的分佈非常吻合。唐華當場就確定了第四種方法“很有希望”。



    這就是照著正確答案挑選解題過程(也就是分析方法)——應該說挑出來的分析方法有較大概率是正確方法。



    唐華現在給李四光展示的,是選定了第四種方法之後,用高精度數值重新運算一次得出的結果,整整用了27個小時。如果華夏構造體系的模型更精細一點,運算時間可能還要翻幾倍。



    “有柱子說明這個地區有油田?松遼盆地這個柱子挺高的啊。怎麼山東這兒也有一個很高的柱子。”李四光看屏幕上的運算結果,若有所思地說,“我看看你們運算的時候用的是什麼公式。”



    唐華坐到李四光旁邊,把源文件調了出來,一行行地解釋這些運算都是用的什麼公式、怎麼計算的。



    “停一下,這裡不對。”李四光說。



    唐華:“這一行?盆地演化數值分析?”



    “對。盆地模擬的主要目的是計算資源量,你這裡是做到了演算每個地層的壓實過程,計算出烴源層的排烴量,但是,不能由此計算出油氣的聚集數量。因為石油到底是怎麼排烴而生成的,這其實不僅和壓實過程有關,還和浮力有關,就是所謂的浮力排烴猜想。我直覺覺得這裡要改。”



    李四光找了支筆,在紙上寫下一個新的公式,“這一部分的運算,你換這個公式試試看。”



    唐華笑容略顯僵硬……還是照著李四光的公式把這一段源文件給改了。



    “這個地方,你用的參數是過時的,好像是我1934年出的那本書裡的吧。1947年美國地質學會刊的論文裡提出了一種新的算法,參數應該是這樣的。我猜你是隻收集和石油有關的論文,沒有注意到會刊裡的這篇。這種算法非常重要,提出這種思想的人很厲害。”



    唐華笑容略顯苦澀……但李四光在一旁看著,只好乖乖地換上新參數。



    ……



    李四光用了一個小時,把源文件重審了一遍,改動了三四個地方。或者是換用另一個公式進行運算,或者是用新的參數。



    “李部長,那我們再運行一次,用高精度數值進行運算?”



    李四光:“對,儘量算得準確一點。”



    唐華檢視源文件,這樣改了幾個地方後,運算量似乎反而有所下降。具體降低多少目測不出來,上一次運算是27小時,這次是24小時左右?



    “李部長,開始運算了,按運算能力預估,大概要24小時得出結果。機房保持開機運行,維護性工作有計算中心的管理員來幫忙,明天的這個時候我們再回來就可以了。”