A computer-implemented method for determining a lithology of a subterranean formation into which a wellbore has been drilled, comprising: (a) receiving a set of measurement logs comprising a plurality of measurement logs, each including a measured characteristic of the wellbore plotted according to depth, the set further including a cuttings percentage log, a gamma ray log a formation strength log and total gas in mud log over a region of the wellbore, (b) processing the cuttings percentage log to segment the region into a plurality of zones, with zone boundaries being established based on appearances or disappearances of rock types in the cuttings percentage log, (c) processing the cuttings percentage log to compute a total percentage of each rock type in each of the plurality of zones and further computing an aggregated length of layers of each rock type in each of the plurality of zones, wherein said length is defined as a change in said depth; (d) processing the gamma ray log, the formation strength log, and the total gas log to compute locations of clay and shale, limestone, and sandstone in each of the plurality of zones, wherein said processing the gamma ray log, the formation strength log, and the total gas in mud log includes: (i) processing the gamma ray log to identify a local gamma ray maximum and identifying a depth interval of clay and shale at the local gamma ray maximum; (ii) processing the formation strength log to identify a local formation strength maximum and identifying a depth interval of limestone at the local formation strength maximum; and (iii) processing the total gas in mud log to identify a local total gas maximum and identifying a depth interval of sandstone at the local total gas maximum. (a) receiving a set of measurement logs comprising a plurality of measurement logs, each including a measured characteristic of the wellbore plotted according to depth, the set further including a cuttings percentage log, a gamma ray log a formation strength log and total gas in mud log over a region of the wellbore, (b) processing the cuttings percentage log to segment the region into a plurality of zones, with zone boundaries being established based on appearances or disappearances of rock types in the cuttings percentage log, (c) processing the cuttings percentage log to compute a total percentage of each rock type in each of the plurality of zones and further computing an aggregated length of layers of each rock type in each of the plurality of zones, wherein said length is defined as a change in said depth; (d) processing the gamma ray log, the formation strength log, and the total gas log to compute locations of clay and shale, limestone, and sandstone in each of the plurality of zones, wherein said processing the gamma ray log, the formation strength log, and the total gas in mud log includes: (i) processing the gamma ray log to identify a local gamma ray maximum and identifying a depth interval of clay and shale at the local gamma ray maximum; (ii) processing the formation strength log to identify a local formation strength maximum and identifying a depth interval of limestone at the local formation strength maximum; and (iii) processing the total gas in mud log to identify a local total gas maximum and identifying a depth interval of sandstone at the local total gas maximum. (i) processing the gamma ray log to identify a local gamma ray maximum and identifying a depth interval of clay and shale at the local gamma ray maximum; (ii) processing the formation strength log to identify a local formation strength maximum and identifying a depth interval of limestone at the local formation strength maximum; and (iii) processing the total gas in mud log to identify a local total gas maximum and identifying a depth interval of sandstone at the local total gas maximum.