In this study, we present a numerical simulation approach for designing Fresnel zone plate (FZP) patterns. By optimizing surface phase parameters using desired merit functions in ray-tracing software, the obtained surface phase was converted into an FZP pattern through a 5-step procedure. A comparison between our numerical simulation approach and the traditional analytical method showed a negligible zone size difference of 0.606 nm and nearly absolute agreement of 17.549 μm in focal spot size. The FZP pattern was experimentally verified by an expected focal spot size of 18.55 μm. Our approach demonstrated design flexibility and has potential applications in simulating various functionalities in FZP patterns and refractive-diffractive hybrid lenses to address specific optical challenges. The surface phase can be freely modified based on optimization objectives that cannot be achieved using the analytical approach, ensuring high-precision design for accurate extraction.